NFTest: automated testing of Nextflow pipelines.

MOTIVATION: The ongoing expansion in the volume of biomedical data has contributed to a growing complexity in the tools and technologies used in research with an increased reliance on complex workflows written in orchestration languages such as Nextflow to integrate algorithms into processing pipelines. The growing use of workflows involving various tools and algorithms has led to increased scrutiny of software development practices to avoid errors in individual tools and in the connections between them. RESULTS: To facilitate test-driven development of Nextflow pipelines, we created NFTest, a framework for automated pipeline testing and validation with customizability options for Nextflow features. It is open-source, easy to initialize and use, and customizable to allow for testing of complex workflows with test success configurable through a broad range of assertions. NFTest simplifies the testing burden on developers by automating tests once defined and providing a flexible interface for running tests to validate workflows. This reduces the barrier to rigorous biomedical workflow testing and paves the way toward reducing computational errors in biomedicine. AVAILABILITY AND IMPLEMENTATION: NFTest is an open-source Python framework under the GPLv2 license and is freely available at https://github.com/uclahs-cds/tool-NFTest. The call-sSNV Nextflow pipeline is available at: https://github.com/uclahs-cds/pipeline-call-sSNV.

The genome assembly of Island Oak (Quercus tomentella), a relictual island tree species.

Island oak (Quercus tomentella) is a rare relictual island tree species that exists only on six islands off the coast of California and Mexico, but was once widespread throughout mainland California. Currently, this species is endangered by threats such as non-native plants, grazing animals, and human removal. Efforts for conservation and restoration of island oak currently underway could benefit from information about its range-wide genetic structure and evolutionary history. Here we present a high-quality genome assembly for Q. tomentella, assembled using PacBio HiFi and Omni-C sequencing, developed as part of the California Conservation Genomics Project (CCGP). The resulting assembly has a length of 781 Mb, with a contig N50 of 22.0 Mb and a scaffold N50 of 63.4 Mb. This genome assembly will provide a resource for genomics-informed conservation of this rare oak species. Additionally, this reference genome will be the first one available for a species in Quercus section Protobalanus, a unique oak clade present only in western North America.

Human T cell generation is restored in CD3δ severe combined immunodeficiency through adenine base editing.

CD3δ SCID is a devastating inborn error of immunity caused by mutations in CD3D, encoding the invariant CD3δ chain of the CD3/TCR complex necessary for normal thymopoiesis. We demonstrate an adenine base editing (ABE) strategy to restore CD3δ in autologous hematopoietic stem and progenitor cells (HSPCs). Delivery of mRNA encoding a laboratory-evolved ABE and guide RNA into a CD3δ SCID patient's HSPCs resulted in a 71.2% ± 7.85% (n = 3) correction of the pathogenic mutation. Edited HSPCs differentiated in artificial thymic organoids produced mature T cells exhibiting diverse TCR repertoires and TCR-dependent functions. Edited human HSPCs transplanted into immunodeficient mice showed 88% reversion of the CD3D defect in human CD34+ cells isolated from mouse bone marrow after 16 weeks, indicating correction of long-term repopulating HSCs. These findings demonstrate the preclinical efficacy of ABE in HSPCs for the treatment of CD3δ SCID, providing a foundation for the development of a one-time treatment for CD3δ SCID patients.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks.

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

Ancient Introgression Between Distantly Related White Oaks (Quercus sect. Quercus) Shows Evidence of Climate-Associated Asymmetric Gene Exchange.

Ancient introgression can be an important source of genetic variation that shapes the evolution and diversification of many taxa. Here, we estimate the timing, direction, and extent of gene flow between two distantly related oak species in the same section (Quercus sect. Quercus). We estimated these demographic events using genotyping by sequencing data, which generated 25 702 single nucleotide polymorphisms for 24 individuals of California scrub oak (Quercus berberidifolia) and 23 individuals of Engelmann oak (Quercus engelmannii). We tested several scenarios involving gene flow between these species using the diffusion approximation-based population genetic inference framework and model-testing approach of the Python package DaDi. We found that the most likely demographic scenario includes a bottleneck in Q. engelmannii that coincides with asymmetric gene flow from Q. berberidifolia into Q. engelmannii. Given that the timing of this gene flow coincides with the advent of a Mediterranean-type climate in the California Floristic Province, we propose that changing precipitation patterns and seasonality may have favored the introgression of climate-associated genes from the endemic into the non-endemic California oak.

Landscape genomics of Quercus lobata reveals genes involved in local climate adaptation at multiple spatial scales.

Understanding how the environment shapes genetic variation provides critical insight about the evolution of local adaptation in natural populations. At multiple spatial scales and multiple geographic contexts within a single species, such information could address a number of fundamental questions about the scale of local adaptation and whether or not the same loci are involved at different spatial scales or geographic contexts. We used landscape genomic approaches from three local elevational transects and rangewide sampling to (a) identify genetic variation underlying local adaptation to environmental gradients in the California endemic oak, Quercus lobata; (b) examine whether putatively adaptive SNPs show signatures of selection at multiple spatial scales; and (c) map putatively adaptive variation to assess the scale and pattern of local adaptation. Of over 10 k single-nucleotide polymorphisms (SNPs) generated with genotyping-by-sequencing, we found signatures of natural selection by climate or local environment at over 600 SNPs (536 loci), some at multiple spatial scales across multiple analyses. Candidate SNPs identified with gene-environment tests (LFMM) at the rangewide scale also showed elevated associations with climate variables compared to the background at both rangewide and elevational transect scales with gradient forest analysis. Some loci overlap with those detected in other oak species, raising the question of whether the same loci might be involved in local climate adaptation in different congeneric species that inhabit different geographic contexts. Mapping landscape patterns of adaptive versus background genetic variation identified regions of marked local adaptation and suggests nonlinear association of candidate SNPs and environmental variables. Taken together, our results offer robust evidence for novel candidate genes for local climate adaptation at multiple spatial scales.

An epigenetic gene silencing pathway selectively acting on transgenic DNA in the green alga Chlamydomonas.

Silencing of exogenous DNA can make transgene expression very inefficient. Genetic screens in the model alga Chlamydomonas have demonstrated that transgene silencing can be overcome by mutations in unknown gene(s), thus producing algal strains that stably express foreign genes to high levels. Here, we show that the silencing mechanism specifically acts on transgenic DNA. Once a permissive chromatin structure has assembled, transgene expression can persist even in the absence of mutations disrupting the silencing pathway. We have identified the gene conferring the silencing and show it to encode a sirtuin-type histone deacetylase. Loss of gene function does not appreciably affect endogenous gene expression. Our data suggest that transgenic DNA is recognized and then quickly inactivated by the assembly of a repressive chromatin structure composed of deacetylated histones. We propose that this mechanism may have evolved to provide protection from potentially harmful types of environmental DNA.

An in situ high-throughput screen identifies inhibitors of intracellular Burkholderia pseudomallei with therapeutic efficacy.

Burkholderia pseudomallei (Bp) and Burkholderia mallei (Bm) are Tier-1 Select Agents that cause melioidosis and glanders, respectively. These are highly lethal human infections with limited therapeutic options. Intercellular spread is a hallmark of Burkholderia pathogenesis, and its prominent ties to virulence make it an attractive therapeutic target. We developed a high-throughput cell-based phenotypic assay and screened ∼220,000 small molecules for their ability to disrupt intercellular spread by Burkholderia thailandensis, a closely related BSL-2 surrogate. We identified 268 hits, and cross-species validation found 32 hits that also disrupt intercellular spread by Bp and/or Bm Among these were a fluoroquinolone analog, which we named burkfloxacin (BFX), which potently inhibits growth of intracellular Burkholderia, and flucytosine (5-FC), an FDA-approved antifungal drug. We found that 5-FC blocks the intracellular life cycle at the point of type VI secretion system 5 (T6SS-5)-mediated cell-cell spread. Bacterial conversion of 5-FC to 5-fluorouracil and subsequently to fluorouridine monophosphate is required for potent and selective activity against intracellular Burkholderia In a murine model of fulminant respiratory melioidosis, treatment with BFX or 5-FC was significantly more effective than ceftazidime, the current antibiotic of choice, for improving survival and decreasing bacterial counts in major organs. Our results demonstrate the utility of cell-based phenotypic screening for Select Agent drug discovery and warrant the advancement of BFX and 5-FC as candidate therapeutics for melioidosis in humans.

Landscape genomics provides evidence of climate-associated genetic variation in Mexican populations of Quercus rugosa.

Local adaptation is a critical evolutionary process that allows plants to grow better in their local compared to non-native habitat and results in species-wide geographic patterns of adaptive genetic variation. For forest tree species with a long generation time, this spatial genetic heterogeneity can shape the ability of trees to respond to rapid climate change. Here, we identify genomic variation that may confer local environmental adaptations and then predict the extent of adaptive mismatch under future climate as a tool for forest restoration or management of the widely distributed high-elevation oak species Quercus rugosa in Mexico. Using genotyping by sequencing, we identified 5,354 single nucleotide polymorphisms (SNPs) genotyped from 103 individuals across 17 sites in the Trans-Mexican Volcanic Belt, and, after controlling for neutral genetic structure, we detected 74 F ST outlier SNPs and 97 SNPs associated with climate variation. Then, we deployed a nonlinear multivariate model, Gradient Forests, to map turnover in allele frequencies along environmental gradients and predict areas most sensitive to climate change. We found that spatial patterns of genetic variation were most strongly associated with precipitation seasonality and geographic distance. We identified regions of contemporary genetic and climatic similarities and predicted regions where future populations of Q. rugosa might be at risk due to high expected rate of climate change. Our findings provide preliminary details for future management strategies of Q. rugosa in Mexico and also illustrate how a landscape genomic approach can provide a useful tool for conservation and resource management strategies.

Complete genomic sequences of Propionibacterium freudenreichii phages from Swiss cheese reveal greater diversity than Cutibacterium (formerly Propionibacterium) acnes phages.

BACKGROUND: A remarkable exception to the large genetic diversity often observed for bacteriophages infecting a specific bacterial host was found for the Cutibacterium acnes (formerly Propionibacterium acnes) phages, which are highly homogeneous. Phages infecting the related species, which is also a member of the Propionibacteriaceae family, Propionibacterium freudenreichii, a bacterium used in production of Swiss-type cheeses, have also been described and are common contaminants of the cheese manufacturing process. However, little is known about their genetic composition and diversity. RESULTS: We obtained seven independently isolated bacteriophages that infect P. freudenreichii from Swiss-type cheese samples, and determined their complete genome sequences. These data revealed that all seven phage isolates are of similar genomic length and GC% content, but their genomes are highly diverse, including genes encoding the capsid, tape measure, and tail proteins. In contrast to C. acnes phages, all P. freudenreichii phage genomes encode a putative integrase protein, suggesting they are capable of lysogenic growth. This is supported by the finding of related prophages in some P. freudenreichii strains. The seven phages could further be distinguished as belonging to two distinct genomic types, or 'clusters', based on nucleotide sequences, and host range analyses conducted on a collection of P. freudenreichii strains show a higher degree of host specificity than is observed for the C. acnes phages. CONCLUSIONS: Overall, our data demonstrate P. freudenreichii bacteriophages are distinct from C. acnes phages, as evidenced by their higher genetic diversity, potential for lysogenic growth, and more restricted host ranges. This suggests substantial differences in the evolution of these related species from the Propionibacteriaceae family and their phages, which is potentially related to their distinct environmental niches.

High-throughput sequencing of the chloroplast and mitochondrion of Chlamydomonas reinhardtii to generate improved de novo assemblies, analyze expression patterns and transcript speciation, and evaluate diversity among laboratory strains and wild isolates.

Chlamydomonas reinhardtii is a unicellular chlorophyte alga that is widely studied as a reference organism for understanding photosynthesis, sensory and motile cilia, and for development of an algal-based platform for producing biofuels and bio-products. Its highly repetitive, ~205-kbp circular chloroplast genome and ~15.8-kbp linear mitochondrial genome were sequenced prior to the advent of high-throughput sequencing technologies. Here, high coverage shotgun sequencing was used to assemble both organellar genomes de novo. These new genomes correct dozens of errors in the prior genome sequences and annotations. Genome sequencing coverage indicates that each cell contains on average 83 copies of the chloroplast genome and 130 copies of the mitochondrial genome. Using protocols and analyses optimized for organellar transcripts, RNA-Seq was used to quantify their relative abundances across 12 different growth conditions. Forty-six percent of total cellular mRNA is attributable to high expression from a few dozen chloroplast genes. RNA-Seq data were used to guide gene annotation, to demonstrate polycistronic gene expression, and to quantify splicing of psaA and psbA introns. In contrast to a conclusion from a recent study, we found that chloroplast transcripts are not edited. Unexpectedly, cytosine-rich polynucleotide tails were observed at the 3'-end of all mitochondrial transcripts. A comparative genomics analysis of eight laboratory strains and 11 wild isolates of C. reinhardtii identified 2658 variants in the organellar genomes, which is 1/10th as much genetic diversity as is found in the nucleus.

First Draft Assembly and Annotation of the Genome of a California Endemic Oak Quercus lobata Née (Fagaceae).

Oak represents a valuable natural resource across Northern Hemisphere ecosystems, attracting a large research community studying its genetics, ecology, conservation, and management. Here we introduce a draft genome assembly of valley oak (Quercus lobata) using Illumina sequencing of adult leaf tissue of a tree found in an accessible, well-studied, natural southern California population. Our assembly includes a nuclear genome and a complete chloroplast genome, along with annotation of encoded genes. The assembly contains 94,394 scaffolds, totaling 1.17 Gb with 18,512 scaffolds of length 2 kb or longer, with a total length of 1.15 Gb, and a N50 scaffold size of 278,077 kb. The k-mer histograms indicate an diploid genome size of ∼720-730 Mb, which is smaller than the total length due to high heterozygosity, estimated at 1.25%. A comparison with a recently published European oak (Q. robur) nuclear sequence indicates 93% similarity. The Q. lobata chloroplast genome has 99% identity with another North American oak, Q. rubra Preliminary annotation yielded an estimate of 61,773 predicted protein-coding genes, of which 71% had similarity to known protein domains. We searched 956 Benchmarking Universal Single-Copy Orthologs, and found 863 complete orthologs, of which 450 were present in > 1 copy. We also examined an earlier version (v0.5) where duplicate haplotypes were removed to discover variants. These additional sources indicate that the predicted gene count in Version 1.0 is overestimated by 37-52%. Nonetheless, this first draft valley oak genome assembly represents a high-quality, well-annotated genome that provides a tool for forest restoration and management practices.

Ni induces the CRR1-dependent regulon revealing overlap and distinction between hypoxia and Cu deficiency responses in Chlamydomonas reinhardtii.

The selectivity of metal sensors for a single metal ion is critical for cellular metal homeostasis. A suite of metal-responsive regulators is required to maintain a prescribed balance of metal ions ensuring that each apo-protein binds the correct metal. However, there are cases when non-essential metals ions disrupt proper metal sensing. An analysis of the Ni-responsive transcriptome of the green alga Chlamydomonas reinhardtii reveals that Ni artificially turns on the CRR1-dependent Cu-response regulon. Since this regulon also responds to hypoxia, a combinatorial transcriptome analysis was leveraged to gain insight into the mechanisms by which Ni interferes with the homeostatic regulation of Cu and oxygen status. Based on parallels with the effect of Ni on the hypoxic response in animals, we propose that a possible link between Cu, oxygen and Ni sensing is an as yet uncharacterized prolyl hydroxylase that regulates a co-activator of CRR1. This analysis also identified transcriptional responses to the pharmacological activation of the Cu-deficiency regulon. Although the Ni-responsive CRR1 regulon is composed of 56 genes (defined as the primary response), 259 transcripts responded to Ni treatment only when a copy of the wild-type CRR1 gene was present. The genome-wide impact of CRR1 target genes on the transcriptome was also evident from the 210 transcripts that were at least 2-fold higher in the crr1 strain, where the abundance of many CRR1 targets was suppressed. Additionally, we identified 120 transcripts that responded to Ni independent of CRR1 function. The putative functions of the proteins encoded by these transcripts suggest that high Ni results in protein damage.

An Indexed, Mapped Mutant Library Enables Reverse Genetics Studies of Biological Processes in Chlamydomonas reinhardtii.

The green alga Chlamydomonas reinhardtii is a leading unicellular model for dissecting biological processes in photosynthetic eukaryotes. However, its usefulness has been limited by difficulties in obtaining mutants in specific genes of interest. To allow generation of large numbers of mapped mutants, we developed high-throughput methods that (1) enable easy maintenance of tens of thousands of Chlamydomonas strains by propagation on agar media and by cryogenic storage, (2) identify mutagenic insertion sites and physical coordinates in these collections, and (3) validate the insertion sites in pools of mutants by obtaining >500 bp of flanking genomic sequences. We used these approaches to construct a stably maintained library of 1935 mapped mutants, representing disruptions in 1562 genes. We further characterized randomly selected mutants and found that 33 out of 44 insertion sites (75%) could be confirmed by PCR, and 17 out of 23 mutants (74%) contained a single insertion. To demonstrate the power of this library for elucidating biological processes, we analyzed the lipid content of mutants disrupted in genes encoding proteins of the algal lipid droplet proteome. This study revealed a central role of the long-chain acyl-CoA synthetase LCS2 in the production of triacylglycerol from de novo-synthesized fatty acids.

Chlamydomonas Genome Resource for Laboratory Strains Reveals a Mosaic of Sequence Variation, Identifies True Strain Histories, and Enables Strain-Specific Studies.

Chlamydomonas reinhardtii is a widely used reference organism in studies of photosynthesis, cilia, and biofuels. Most research in this field uses a few dozen standard laboratory strains that are reported to share a common ancestry, but exhibit substantial phenotypic differences. In order to facilitate ongoing Chlamydomonas research and explain the phenotypic variation, we mapped the genetic diversity within these strains using whole-genome resequencing. We identified 524,640 single nucleotide variants and 4812 structural variants among 39 commonly used laboratory strains. Nearly all (98.2%) of the total observed genetic diversity was attributable to the presence of two, previously unrecognized, alternate haplotypes that are distributed in a mosaic pattern among the extant laboratory strains. We propose that these two haplotypes are the remnants of an ancestral cross between two strains with ∼2% relative divergence. These haplotype patterns create a fingerprint for each strain that facilitates the positive identification of that strain and reveals its relatedness to other strains. The presence of these alternate haplotype regions affects phenotype scoring and gene expression measurements. Here, we present a rich set of genetic differences as a community resource to allow researchers to more accurately conduct and interpret their experiments with Chlamydomonas.

Genome-wide gene order distances support clustering the gram-positive bacteria.

Initially using 143 genomes, we developed a method for calculating the pair-wise distance between prokaryotic genomes using a Monte Carlo method to estimate the conservation of gene order. The method was based on repeatedly selecting five or six non-adjacent random orthologs from each of two genomes and determining if the chosen orthologs were in the same order. The raw distances were then corrected for gene order convergence using an adaptation of the Jukes-Cantor model, as well as using the common distance correction D' = -ln(1-D). First, we compared the distances found via the order of six orthologs to distances found based on ortholog gene content and small subunit rRNA sequences. The Jukes-Cantor gene order distances are reasonably well correlated with the divergence of rRNA (R (2) = 0.24), especially at rRNA Jukes-Cantor distances of less than 0.2 (R (2) = 0.52). Gene content is only weakly correlated with rRNA divergence (R (2) = 0.04) over all distances, however, it is especially strongly correlated at rRNA Jukes-Cantor distances of less than 0.1 (R (2) = 0.67). This initial work suggests that gene order may be useful in conjunction with other methods to help understand the relatedness of genomes. Using the gene order distances in 143 genomes, the relations of prokaryotes were studied using neighbor joining and agreement subtrees. We then repeated our study of the relations of prokaryotes using gene order in 172 complete genomes better representing a wider-diversity of prokaryotes. Consistently, our trees show the Actinobacteria as a sister group to the bulk of the Firmicutes. In fact, the robustness of gene order support was found to be considerably greater for uniting these two phyla than for uniting any of the proteobacterial classes together. The results are supportive of the idea that Actinobacteria and Firmicutes are closely related, which in turn implies a single origin for the gram-positive cell.

Systems-level analysis of nitrogen starvation-induced modifications of carbon metabolism in a Chlamydomonas reinhardtii starchless mutant.

To understand the molecular basis underlying increased triacylglycerol (TAG) accumulation in starchless (sta) Chlamydomonas reinhardtii mutants, we undertook comparative time-course transcriptomics of strains CC-4348 (sta6 mutant), CC-4349, a cell wall-deficient (cw) strain purported to represent the parental STA6 strain, and three independent STA6 strains generated by complementation of sta6 (CC-4565/STA6-C2, CC-4566/STA6-C4, and CC-4567/STA6-C6) in the context of N deprivation. Despite N starvation-induced dramatic remodeling of the transcriptome, there were relatively few differences (5 × 10(2)) observed between sta6 and STA6, the most dramatic of which were increased abundance of transcripts encoding key regulated or rate-limiting steps in central carbon metabolism, specifically isocitrate lyase, malate synthase, transaldolase, fructose bisphosphatase and phosphoenolpyruvate carboxykinase (encoded by ICL1, MAS1, TAL1, FBP1, and PCK1 respectively), suggestive of increased carbon movement toward hexose-phosphate in sta6 by upregulation of the glyoxylate pathway and gluconeogenesis. Enzyme assays validated the increase in isocitrate lyase and malate synthase activities. Targeted metabolite analysis indicated increased succinate, malate, and Glc-6-P and decreased Fru-1,6-bisphosphate, illustrating the effect of these changes. Comparisons of independent data sets in multiple strains allowed the delineation of a sequence of events in the global N starvation response in C. reinhardtii, starting within minutes with the upregulation of alternative N assimilation routes and carbohydrate synthesis and subsequently a more gradual upregulation of genes encoding enzymes of TAG synthesis. Finally, genome resequencing analysis indicated that (1) the deletion in sta6 extends into the neighboring gene encoding respiratory burst oxidase, and (2) a commonly used STA6 strain (CC-4349) as well as the sequenced reference (CC-503) are not congenic with respect to sta6 (CC-4348), underscoring the importance of using complemented strains for more rigorous assignment of phenotype to genotype.

Amino acid signatures of salinity on an environmental scale with a focus on the Dead Sea.

The increase of the acidic nature of proteins as an adaptation to hypersalinity has been well documented within halophile isolates. Here we explore the effect of salinity on amino acid preference on an environmental scale. Via pyrosequencing, we have obtained two distinct metagenomic data sets from the Dead Sea, one from a 1992 archaeal bloom and one from the modern Dead Sea. Our data, along with metagenomes from environments representing a range of salinities, show a strong linear correlation (R(2) = 0.97) between the salinity of an environment and the ratio of acidic to basic amino acids encoded by its inhabitants. Using the amino acid composition of putative protein-encoding reads and the results of 16S rRNA amplicon sequencing, we differentiate recovered sequences representing microorganisms indigenous to the Dead Sea from lateral gene transfer events and foreign DNA. Our methods demonstrate lateral gene transfer events between a halophilic archaeon and relatives of the thermophilic bacterial genus Thermotoga and suggest the presence of indigenous Dead Sea representatives from 10 traditionally non-hyperhalophilic bacterial lineages. The work suggests the possibility that amino acid bias of hypersaline environments might be preservable in fossil DNA or fossil amino acids, serving as a proxy for the salinity of an ancient environment. Finally, both the amino acid profile of the 2007 Dead Sea metagenome and the V9 amplicon library support the conclusion that the dominant microorganism inhabiting the Dead Sea is most closely related to a thus far uncultured relative of an alkaliphilic haloarchaeon.

Terminal addition, the Cambrian radiation and the Phanerozoic evolution of bilaterian form.

We examine terminal addition, the process of addition of serial elements in a posterior subterminal growth zone during animal development, across modern taxa and fossil material. We argue that terminal addition was the basal condition in Bilateria, and that modification of terminal addition was an important component of the rapid Cambrian evolution of novel bilaterian morphology. We categorize the often-convergent modifications of terminal addition from the presumed ancestral condition. Our focus on terminal addition and its modification highlights trends in the history of animal evolution evident in the fossil record. These trends appear to be the product of departure from the initial terminal addition state, as is evident in evolutionary patterns within-fossil groups such as trilobites, but is also more generally related to shifts in types of morphologic change through the early Phanerozoic. Our argument is contingent on dates of metazoan divergence that are roughly convergent with the first appearance of metazoan fossils in the latest Proterozoic and Cambrian, as well as on an inference of homology of terminal addition across bilaterian Metazoa.