Comparative Developmental Transcriptomics Reveals Rewiring of a Highly Conserved Gene Regulatory Network during a Major Life History Switch in the Sea Urchin Genus Heliocidaris.

The ecologically significant shift in developmental strategy from planktotrophic (feeding) to lecithotrophic (nonfeeding) development in the sea urchin genus Heliocidaris is one of the most comprehensively studied life history transitions in any animal. Although the evolution of lecithotrophy involved substantial changes to larval development and morphology, it is not known to what extent changes in gene expression underlie the developmental differences between species, nor do we understand how these changes evolved within the context of the well-defined gene regulatory network (GRN) underlying sea urchin development. To address these questions, we used RNA-seq to measure expression dynamics across development in three species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph H. tuberculata, and an outgroup planktotroph Lytechinus variegatus. Using well-established statistical methods, we developed a novel framework for identifying, quantifying, and polarizing evolutionary changes in gene expression profiles across the transcriptome and within the GRN. We found that major changes in gene expression profiles were more numerous during the evolution of lecithotrophy than during the persistence of planktotrophy, and that genes with derived expression profiles in the lecithotroph displayed specific characteristics as a group that are consistent with the dramatically altered developmental program in this species. Compared to the transcriptome, changes in gene expression profiles within the GRN were even more pronounced in the lecithotroph. We found evidence for conservation and likely divergence of particular GRN regulatory interactions in the lecithotroph, as well as significant changes in the expression of genes with known roles in larval skeletogenesis. We further use coexpression analysis to identify genes of unknown function that may contribute to both conserved and derived developmental traits between species. Collectively, our results indicate that distinct evolutionary processes operate on gene expression during periods of life history conservation and periods of life history divergence, and that this contrast is even more pronounced within the GRN than across the transcriptome as a whole.

Contingent interactions among biofilm-forming bacteria determine preservation or decay in the first steps toward fossilization of marine embryos.

Fossils of soft tissues provide important records of early animals and embryos, and there is substantial evidence for a role for microbes in soft tissue fossilization. We are investigating the initial events in interactions of bacteria with freshly dead tissue, using marine embryos as a model system. We previously found that microbial invasion can stabilize embryo tissue that would otherwise disintegrate in hours or days by generating a bacterial pseudomorph, a three dimensional biofilm that both replaces the tissue and replicates its morphology. In this study, we sampled seawater at different times and places near Sydney, Australia, and determined the range and frequency of different taphonomic outcomes. Although destruction was most common, bacteria in 35% of seawater samples yielded morphology­preserving biofilms. We could replicate the taphonomic pathways seen with seawater bacterial communities using single cultured strains of marine gammaproteobacteria. Each given species reproducibly generated a consistent taphonomic outcome and we identified species that yielded each of the distinct pathways produced by seawater bacterial communities. Once formed,bacterial pseudomorphs are stable for over a year and resist attack by other bacteria and destruction by proteases and other lytic enzymes. Competition studies showed that the initial action of a pseudomorphing strain can be blocked by a strain that destroys tissues. Thus embryo preservation in nature may depend on contingent interactions among bacterial species that determine if pseudomorphing occurs.We used Artemia nauplius larvae to show that bacterial biofilm replacement of tissue is not restricted to embryos, but is relevant for preservation of small multicellular organisms. We present a model for bacterial self­assembly of large­scale three­dimensional tissue pseudomorphs, based on smallscaleinteractions among individual bacterial cells to form local biofilms at structural boundaries within the tissue. Localbiofilms then conjoin to generate the pseudomorph.

Axoneme beta-tubulin sequence determines attachment of outer dynein arms.

Axonemes of motile eukaryotic cilia and flagella have a conserved structure of nine doublet microtubules surrounding a central pair of microtubules. Outer and inner dynein arms on the doublets mediate axoneme motility [1]. Outer dynein arms (ODAs) attach to the doublets at specific interfaces [2-5]. However, the molecular contacts of ODA-associated proteins with tubulins of the doublet microtubules are not known. We report here that attachment of ODAs requires glycine 56 in the beta-tubulin internal variable region (IVR). We show that in Drosophila spermatogenesis, a single amino acid change at this position results in sperm axonemes markedly deficient in ODAs. Moreover, we found that axonemal beta-tubulins throughout the phylogeny have invariant glycine 56 and a strongly conserved IVR, whereas nonaxonemal beta-tubulins vary widely in IVR sequences. Our data reveal a deeply conserved physical requirement for assembly of the macromolecular architecture of the motile axoneme. Amino acid 56 projects into the microtubule lumen [6]. Imaging studies of axonemes indicate that several proteins may interact with the doublet-microtubule lumen [3, 4, 7, 8]. This region of beta-tubulin may determine the conformation necessary for correct attachment of ODAs, or there may be sequence-specific interaction between beta-tubulin and a protein involved in ODA attachment or stabilization.

Embryo fossilization is a biological process mediated by microbial biofilms.

Fossilized embryos with extraordinary cellular preservation appear in the Late Neoproterozoic and Cambrian, coincident with the appearance of animal body fossils. It has been hypothesized that microbial processes are responsible for preservation and mineralization of organic tissues. However, the actions of microbes in preservation of embryos have not been demonstrated experimentally. Here, we show that bacterial biofilms assemble rapidly in dead marine embryos and form remarkable pseudomorphs in which the bacterial biofilm replaces and exquisitely models details of cellular organization and structure. The experimental model was the decay of cleavage stage embryos similar in size and morphology to fossil embryos. The data show that embryo preservation takes place in 3 distinct steps: (i) blockage of autolysis by reducing or anaerobic conditions, (ii) rapid formation of microbial biofilms that consume the embryo but form a replica that retains cell organization and morphology, and (iii) bacterially catalyzed mineralization. Major bacterial taxa in embryo decay biofilms were identified by using 16S rDNA sequencing. Decay processes were similar in different taphonomic conditions, but the composition of bacterial populations depended on specific conditions. Experimental taphonomy generates preservation states similar to those in fossil embryos. The data show how fossilization of soft tissues in sediments can be mediated by bacterial replacement and mineralization, providing a foundation for experimentally creating biofilms from defined microbial species to model fossilization as a biological process.

Genetic basis for divergence in developmental gene expression in two closely related sea urchins.

The genetic basis for divergence in developmental gene expression among species is poorly understood, despite growing evidence that such changes underlie many interesting traits. Here we quantify transcription in hybrids of Heliocidaris tuberculata and Heliocidaris erythrogramma, two closely related sea urchins with highly divergent developmental gene expression and life histories. We find that most expression differences between species result from genetic influences that affect one stage of development, indicating limited pleiotropic consequences for most mutations that contribute to divergence in gene expression. Activation of zygotic transcription is broadly delayed in H. erythrogramma, the species with the derived life history, despite its overall faster premetamorphic development. Altered expression of several terminal differentiation genes associated with the derived larval morphology of H. erythrogramma is based largely on differences in the expression or function of their upstream regulators, providing insights into the genetic basis for the evolution of key life history traits.

Cooperativity between the beta-tubulin carboxy tail and the body of the molecule is required for microtubule function.

Using Drosophila spermatogenesis as a model, we show that function of the beta-tubulin C-terminal tail (CTT) is not independent of the body of the molecule. For optimal microtubule function, the beta-tubulin CTT and body must match. beta2 is the only beta-tubulin used in meiosis and spermatid differentiation. beta1-tubulin is used in basal bodies, but beta1 cannot replace beta2. However, when beta1 is co-expressed with beta2, both beta-tubulins are equally incorporated into all microtubules, and males exhibit near wild type fertility. In contrast, co-expression of beta2beta1C and beta1beta2C, two reciprocal chimeric molecules with bodies and tails swapped, results in defects in meiosis, cytoskeletal microtubules, and axonemes; males produce few functional sperm and few or no progeny. In these experiments, all the same beta-tubulin parts are present, but unlike the co-assembled native beta-tubulins, the "trans" configuration of the co-assembled chimeras is poorly functional. Our data thus reveal essential intra-molecular interactions between the CTT and other parts of the beta-tubulin molecule, even though the CTT is a flexible surface feature of tubulin heterodimers and microtubules. In addition, we show that Drosophila sperm tail length depends on the total tubulin pool available for axoneme assembly and spermatid elongation. D. melanogaster and other Drosophila species have extraordinarily long sperm tails, the length of which is remarkably constant in wild type flies. We show that in males of experimental genotypes that express wild type tubulins but have half the amount of the normal tubulin pool size, sperm tails are substantially shorter than wild type.

Deciphering the fossil record of early bilaterian embryonic development in light of experimental taphonomy.

Experimental analyses of decay in a tunicate deuterostome and three lophotrochozoans indicate that the controls on decay and preservation of embryos, identified previously based on echinoids, are more generally applicable. Four stages of decay are identified regardless of the environment of death and decay. Embryos decay rapidly in oxic and anoxic conditions, although the gross morphology of embryos is maintained for longer under anoxic conditions. Under anoxic reducing conditions, the gross morphology of the embryos is maintained for the longest period of time, compatible with the timescale required for bacterially mediated mineralization of soft tissues. All four stages of decay were encountered under all environmental conditions, matching the spectrum of preservational qualities encountered in all fossil embryo assemblages. The preservation potential of embryos of deuterostomes and lophotrochozoans is at odds with the lack of such embryos in the fossil record. Rather, the fossil record of embryos, as sparse as it is, is dominated by forms interpreted as ecdysozoans, cnidarians, and stem-metazoans. The dearth of deuterostome and lophotrochozoan embryos may be explained by the fact that ecdysozoans, at least, tend to deposit their eggs in the sediment rather than through broadcast spawning. However, fossil embryos remain very rare and the main controlling factor on their fossilization may be the unique conspiracy of environmental conditions at a couple of sites. The preponderance of fossilized embryos of direct developers should not be used in evidence against the existence of indirect development at this time in animal evolutionary history.

Drosophila KAP interacts with the kinesin II motor subunit KLP64D to assemble chordotonal sensory cilia, but not sperm tails.

BACKGROUND: Kinesin II-mediated anterograde intraflagellar transport (IFT) is essential for the assembly and maintenance of flagella and cilia in various cell types. Kinesin associated protein (KAP) is identified as the non-motor accessory subunit of Kinesin II, but its role in the corresponding motor function is not understood. RESULTS: We show that mutations in the Drosophila KAP (DmKap) gene could eliminate the sensory cilia as well as the sound-evoked potentials of Johnston's organ (JO) neurons. Ultrastructure analysis of these mutants revealed that the ciliary axonemes are absent. Mutations in Klp64D, which codes for a Kinesin II motor subunit in Drosophila, show similar ciliary defects. All these defects are rescued by exclusive expression of DmKAP and KLP64D/KIF3A in the JO neurons of respective mutants. Furthermore, reduced copy number of the DmKap gene was found to enhance the defects of hypomorphic Klp64D alleles. Unexpectedly, however, both the DmKap and the Klp64D mutant adults produce vigorously motile sperm with normal axonemes. CONCLUSIONS: KAP plays an essential role in Kinesin II function, which is required for the axoneme growth and maintenance of the cilia in Drosophila type I sensory neurons. However, the flagellar assembly in Drosophila spermatids does not require Kinesin II and is independent of IFT.

Tinkering: new embryos from old--rapidly and cheaply.

Marine embryos and larvae reflect distinct life histories and body plans from their adults, and are relatively simple in morphology and genetic regulation. Evolution of development to produce highly modified larval forms can be rapid among closely related species. We have studied the mechanisms by which the non-feeding direct-developing larva of the sea urchin Heliocidaris erythrogramma evolved from an indirect-developing feeding larva, the pluteus. H. erythrogramma diverged within 4 million years from its sister species, H. tuberculata, which has a typical pluteus larva. Radical evolution of H. erythrogramma early development allows it to reach metamorphosis in three to four days versus the several weeks required for the pluteus. We have used embryology, cross species hybrids, and manipulation of gene expression in embryos to dissect developmental changes and the genic controls that underlie these changes. Evolution of a new larval form resulted largely from several heterochronies in which conserved regulatory pathways are shifted in timing, producing new temporal relationships to other developmental events. Other changes in gene regulation also have contributed to rapid evolution of larval features, including the origin of unexpected and novel tissue identities that transcend changes within homologous features.

Adaptive evolution of bindin in the genus Heliocidaris is correlated with the shift to direct development.

Sea urchins are widely used to study both fertilization and development. In this study we combine the two fields to examine the evolution of reproductive isolation in the genus Heliocidaris. Heliocidaris tuberculata develops indirectly via a feeding larva, whereas the only other species in the genus, H. erythrogramma, has evolved direct development through a nonfeeding larva. We estimated the time of divergence between H. erythrogramma and H. tuberculata from mitochondrial DNA divergence, quantified levels of gametic compatibility between the two species in cross-fertilization assays, and examined the mode of evolution of the sperm protein bindin by sequencing multiple alleles of the two species. Bindin is the major component of the sea urchin sperm acrosomal vesicle, and is involved in sperm-egg attachment and fusion. Based on our analyses, we conclude that: the two species of Heliocidaris diverged less than five million years ago, indicating that direct development can evolve rapidly in sea urchins; since their divergence, the two species have become gametically incompatible; Heliocidaris bindin has evolved under positive selection; and this positive selection is concentrated on the branch leading to H. erythrogramma. Three hypotheses can explain the observed pattern of selection on bindin: (1) it is a correlated response to the evolution of direct development in H. erythrogramma; (2) it is the result of an intraspecific process acting in H. erythrogramma but not in H. tuberculata; or (3) it is the product of reinforcement on the species that invests more energy into each egg to avoid hybridization.

Natural hybridization in the sea urchin genus Pseudoboletia between species without apparent barriers to gamete recognition.

Marine species with high dispersal potential often have huge ranges and minimal population structure. Combined with the paucity of geographic barriers in the oceans, this pattern raises the question as to how speciation occurs in the sea. Over the past 20 years, evidence has accumulated that marine speciation is often linked to the evolution of gamete recognition proteins. Rapid evolution of gamete recognition proteins in gastropods, bivalves, and sea urchins is correlated with gamete incompatibility and contributes to the maintenance of species boundaries between sympatric congeners. Here, we present a counterexample to this general pattern. The sea urchins Pseudoboletia indiana and P. maculata have broad ranges that overlap in the Indian and Pacific oceans. Cytochrome oxidase I sequences indicated that these species are distinct, and their 7.3% divergence suggests that they diverged at least 2 mya. Despite this, we suspected hybridization between them based on the presence of morphologically intermediate individuals in sympatric populations at Sydney, Australia. We assessed the opportunity for hybridization between the two species and found that (1) individuals of the two species occur within a meter of each other in nature, (2) they have overlapping annual reproductive cycles, and (3) their gametes cross-fertilize readily in the laboratory and in the field. We genotyped individuals with intermediate morphology and confirmed that many were hybrids. Hybrids were fertile, and some female hybrids had egg sizes intermediate between the two parental species. Consistent with their high level of gamete compatibility, there is minimal divergence between P. indiana and P. maculata in the gamete recognition protein bindin, with a single fixed amino acid difference between the two species. Pseudoboletia thus provides a well-characterized exception to the idea that broadcast spawning marine species living in sympatry develop and maintain species boundaries through the divergence of gamete recognition proteins and the associated evolution of gamete incompatibility.

Axoneme-dependent tubulin modifications in singlet microtubules of the Drosophila sperm tail.

Drosophila melanogaster sperm tubulins are posttranslationally glutamylated and glycylated. We show here that axonemes are the substrate for these tubulin C-terminal modifications. Axoneme architecture is required, but full length, motile axonemes are not necessary. Tubulin glutamylation occurs during or shortly after assembly into the axoneme; only glutamylated tubulins are glycylated. Tubulins in other testis microtubules are not modified. Only a small subset of total Drosophila sperm axoneme tubulins have these modifications. Biochemical fractionation of Drosophila sperm showed that central pair and accessory microtubules have the majority of poly-modified tubulins, whereas doublet microtubules have only small amounts of mono- and oligo-modified tubulins. Glutamylation patterns for different beta-tubulins experimentally assembled into axonemes were consistent with utilization of modification sites corresponding to those identified in other organisms, but surrounding sequence context was also important. We compared tubulin modifications in the 9 + 9 + 2 insect sperm tail axonemes of Drosophila with the canonical 9 + 2 axonemes of sperm of the sea urchin Lytichinus pictus and the 9 + 0 motile sperm axonemes of the eel Anguilla japonica. In contrast to Drosophila sperm, L. pictus sperm have equivalent levels of modified tubulins in both doublet and central pair microtubule fractions, whereas the doublets of A. japonica sperm exhibit little glutamylation but extensive glycylation. Tubulin C-terminal modifications are a prevalent feature of motile axonemes, but there is no conserved pattern for placement or amount of these

Axoneme specialization embedded in a "generalist" beta-tubulin.

The relationship between the primary structure of the beta-tubulin C-terminal tail (CTT) and axoneme structure and function is explored using the spermatogenesis-specific beta2-tubulin of Drosophila. We previously showed that all beta-tubulins used for motile 9 + 2 axonemes contain a conserved sequence motif in the proximal part of the CTT, the beta-tubulin axoneme motif. The differential ability of tubulin isoforms and abilities of beta2-tubulin C-terminal truncations to form axonemes led us to hypothesize that the axoneme motif is essential for axoneme formation and the distal half of the CTT was less important. The studies we report here indicate that it is not that simple. Unexpectedly, some changes in the core sequence of the axoneme motif did not disrupt formation of motile axonemes. And, while deletion of the distal CTT did not disrupt the ability to produce functional sperm [Popodi et al., Cell Motil Cytoskeleton 2005;62:48-64], changing the amino acid sequence in this region can. Thus both regions are important. The deep conservation of the axoneme motif in all eukaryotic groups implies that the presence of the sequence motif confers a functional advantage. The central pair is the axoneme structure most sensitive to perturbations in tubulin molecules; we hypothesize central pair assembly is facilitated by the presence of this motif. Our data reveal that beta2-tubulin has robust properties for axoneme assembly, and that axonemal specializations are embedded in both the CTT and the body of the beta2 molecule.

Experimental taphonomy shows the feasibility of fossil embryos.

The recent discovery of apparent fossils of embryos contemporaneous with the earliest animal remains may provide vital insights into the metazoan radiation. However, although the putative fossil remains are similar to modern marine animal embryos or larvae, their simple geometric forms also resemble other organic and inorganic structures. The potential for fossilization of animals at such developmental stages and the taphonomic processes that might affect preservation before mineralization have not been examined. Here, we report experimental taphonomy of marine embryos and larvae similar in size and inferred cleavage mode to presumptive fossil embryos. Under conditions that prevent autolysis, embryos within the fertilization envelope can be preserved with good morphology for sufficiently long periods for mineralization to occur. The reported fossil record exhibits size bias, but we show that embryo size is unlikely to be a major factor in preservation. Under some conditions of death, fossilized remains will not accurately reflect the cell structure of the living organism. Although embryos within the fertilization envelope have high preservation potential, primary larvae have negligible preservation potential. Thus the paleo-embryological record may have strong biases on developmental stages preserved. Our data provide a predictive basis for interpreting the fossil record to unravel the evolution of ontogeny in the origin of metazoans.

Cellular and subcellular structure of neoproterozoic animal embryos.

Stereoblastic embryos from the Doushantuo Formation of China exhibit occasional asynchronous cell division, with diminishing blastomere volume as cleavage proceeded. Asynchronous cell division is common in modern embryos, implying that sophisticated mechanisms for differential cell division timing and embryonic cell lineage differentiation evolved before 551 million years ago. Subcellular structures akin to organelles, coated yolk granules, or lipid vesicles occur in these embryos. Paired reniform structures within embryo cells may represent fossil evidence of cells about to undergo division. Embryos exhibit no evidence of epithelial organization, even in embryos composed of approximately 1000 cells. Many of these features are compatible with metazoans, but the absence of epithelialization is consistent only with a stem-metazoan affinity for Doushantuo embryos.

The proximal region of the beta-tubulin C-terminal tail is sufficient for axoneme assembly.

We have used Drosophila testis-specific beta2-tubulin to determine sequence requirements for different microtubules. The beta2-tubulin C-terminal tail has unique sperm-specific functions [Dev Biol 158:267-286 (2003)] and is also important for forming stable heterodimers with alpha-tubulin, a general function common to all microtubules [Mol Biol Cell 12(7):2185-2194 (2001)]. beta-tubulins utilized in motile 9 + 2 axonemes contain a C-terminal sequence "axoneme motif" [Science 275 (1997) 70-73]. C-terminal truncated beta2-tubulin cannot form the sperm tail axoneme. Here we show that a partially truncated beta2-tubulin (beta2Delta7) containing only the proximal portion of the C-terminal tail, including the axoneme motif, can support production of functional motile sperm. We conclude that these proximal eight amino acids specify the binding site for protein(s) essential to support assembly of the motile axoneme. Males that express beta2Delta7, although they are fertile, produce fewer sperm than wild type males. Beta2Delta7 causes a slightly increased error rate in spermatogenesis attributable to loss of stabilizing properties intrinsic to the full-length C-terminal tail. Therefore, beta2Delta7 males would be at a selective disadvantage and it is likely that the full-length C-terminus would be essential in the wild and in evolution.

The best of all worlds or the best possible world? Developmental constraint in the evolution of beta-tubulin and the sperm tail axoneme.

Through evolutionary history, some features of the phenotype show little variation. Stabilizing selection could produce this result, but the possibility also exists that a feature is conserved because it is developmentally constrained--only one or a few developmental mechanisms can produce that feature. We present experimental data documenting developmental constraint in the assembly of the motile sperm tail axoneme. The 9+2 microtubule architecture of the eukaryotic axoneme has been deeply conserved. We argue that the quality of motility supported by axonemes with this morphology explains their long conservation, rather than a developmental necessity for the 9+2 architecture. However, our functional tests in Drosophila spermatogenesis reveal considerable constraint in the coevolution of testis-specific beta-tubulin and the sperm tail axoneme. The evolution of testis beta-tubulins used in insect sperm tail axonemes is highly punctuated, indicating some pressure acting on their evolution. We provide a mechanistic explanation for their punctuated evolution by testing structure-function relationships between testis beta-tubulin and the motile axoneme in D. melanogaster. We discovered that a highly conserved sequence feature of beta-tubulins used in motile axonemes is needed to specify central pair formation. Second, our data suggest that cooperativity in the function of internal beta-tubulin amino acids is needed to support the long axonemes characteristic of Drosophila sperm tails. Thus, central pair formation constrains the evolution of the axoneme motif, and intramolecular cooperativity makes the evolution of the internal residues path dependent, which slows their evolution. Our results explain why a highly specialized beta-tubulin is needed to construct the Drosophila sperm tail axoneme. We conclude that these constraints have fixed testis-specific beta-tubulin identity in Drosophila.

Regulatory punctuated equilibrium and convergence in the evolution of developmental pathways in direct-developing sea urchins.

We made hybrid crosses between closely and distantly related sea urchin species to test two hypotheses about the evolution of gene regulatory systems in the evolution of ontogenetic pathways and larval form. The first hypothesis is that gene regulatory systems governing development evolve in a punctuational manner during periods of rapid morphological evolution but are relatively stable over long periods of slow morphological evolution. We compared hybrids between direct and indirect developers from closely and distantly related families. Hybrids between eggs of the direct developer Heliocidaris erythrogramma and sperm of the 4-million year distant species H. tuberculata, an indirect developer, restored feeding larval structures and paternal gene expression that were lost in the evolution of the direct-developing maternal parent. Hybrids resulting from the cross between eggs of H. erythrogramma and sperm of the 40-million year distant indirect-developer Pseudoboletia maculata are strikingly similar to hybrids between the congeneric hybrids. The marked similarities in ontogenetic trajectory and morphological outcome in crosses of involving either closely or distantly related indirect developing species indicates that their regulatory mechanisms interact with those of H. erythrogramma in the same way, supporting remarkable conservation of molecular control pathways among indirect developers. Second, we tested the hypothesis that convergent developmental pathways in independently evolved direct developers reflect convergence of the underlying regulatory systems. Crosses between two independently evolved direct-developing species from two 70-million year distant families, H. erythrogramma and Holopneustes purpurescens, produced harmoniously developing hybrid larvae that maintained the direct mode of development and did not exhibit any obvious restoration of indirect-developing features. These results are consistent with parallel evolution of direct-developing features in these two lineages.