The first high-quality genome assembly and annotation of Anthocidaris crassispina.

Anthocidaris crassispina is a very popular edible sea urchin distributed along the coast of the South China Sea. In this study, we performed whole-genome sequencing and generated a chromosome-level assembly of this species. The total length of the genomic contig sequence was 891.02 Mb, and contig N50 was 808.15 kb when Hifiasm was used for assembly. The Hi-C library was constructed and sequenced, yielding approximately 68.61 Gb of data. After Hi-C assembly, approximately 886.72 Mb of sequence was able to be mapped onto 21 chromosomes, accounting for 99.52% of the total genome length. Among the sequences located on the chromosomes, those for which the order and direction could be determined accounted for approximately 826.82 Mb, or 93.24% of the total length. These results provide valuable resources for further study of A. crassispina at the genetic level.

Combinatorial regulatory states define cell fate diversity during embryogenesis.

Cell fate specification occurs along invariant species-specific trajectories that define the animal body plan. This process is controlled by gene regulatory networks that regulate the expression of the limited set of transcription factors encoded in animal genomes. Here we globally assess the spatial expression of ~90% of expressed transcription factors during sea urchin development from embryo to larva to determine the activity of gene regulatory networks and their regulatory states during cell fate specification. We show that >200 embryonically expressed transcription factors together define >70 cell fates that recapitulate the morphological and functional organization of this organism. Most cell fate-specific regulatory states consist of ~15-40 transcription factors with similarity particularly among functionally related cell types regardless of developmental origin. Temporally, regulatory states change continuously during development, indicating that progressive changes in regulatory circuit activity determine cell fate specification. We conclude that the combinatorial expression of transcription factors provides molecular definitions that suffice for the unique specification of cell states in time and space during embryogenesis.

Morphological and Molecular Evidence of an Intergeneric Host-Range in Clavisodalis sentifer (Crustacea: Copepoda: Taeniacanthidae) Associated with Diadematid Sea Urchins from the Western Pacific.

Sea urchins have a wide variety of symbionts on their body surfaces and inside their bodies. Copepods of the genus Clavisodalis (Taeniacanthidae) collected from the esophagus of sea urchins of the genera Diadema and Echinothrix in southern Japan were identified based on their morphological characteristics, and molecular analysis was conducted to determine whether genetic variation occurs in copepods from different localities and hosts. Morphological observations identified individuals from southern Japan as Clavisodalis sentifer Dojiri and Humes, 1982, making this the first record of this species in the northern hemisphere and the first record of its genus in Japan. Morphological and molecular analysis suggested that the copepod specimens collected from multiple hosts across two genera would be the same species. Considering the typically observed high level of host specificity among taeniacanthid copepods, the utilization of hosts from two genera by C. sentifer is noteworthy.

Stable germline transgenesis using the Minos Tc1/mariner element in the sea urchin Lytechinus pictus.

Stable transgenesis is a transformative tool in model organism biology. Although the sea urchin is one of the oldest animal models in cell and developmental biology, studies in this animal have largely relied on transient manipulation of wild animals, without a strategy for stable transgenesis. Here, we build on recent progress to develop a more genetically tractable sea urchin species, Lytechinus pictus, and establish a robust transgene integration method. Three commonly used transposons (Minos, Tol2 and piggyBac) were tested for non-autonomous transposition, using plasmids containing a polyubiquitin promoter upstream of a H2B-mCerulean nuclear marker. Minos was the only transposable element that resulted in significant expression beyond metamorphosis. F0 animals were raised to sexual maturity, and spawned to determine germline integration and transgene inheritance frequency, and to characterize expression patterns of the transgene in F1 progeny. The results demonstrate transgene transmission through the germline, the first example of a germline transgenic sea urchin and, indeed, of any echinoderm. This milestone paves the way for the generation of diverse transgenic resources that will dramatically enhance the utility, reproducibility and efficiency of sea urchin research.

Chromosome-level genome assembly and annotation of the black sea urchin Arbacia lixula (Linnaeus, 1758).

The black sea urchin (Arbacia lixula) is a keystone species inhabiting the coastal shallow waters of the Mediterranean Sea, which is a key driver of littoral communities' structure. Here, we present the first genome assembly and annotation of this species, standing as the first Arbacioida genome, including both nuclear and mitochondrial genomes. To obtain a chromosome-level assembly, we used a combination of PacBio high fidelity (HiFi) reads and chromatin capture reads (Omni-C). In addition, we generated a high-quality nuclear annotation of both coding and non-coding genes, by using published RNA-Seq data from several individuals of A. lixula and gene models from closely related species. The nuclear genome assembly has a total span of 607.91 Mb, being consistent with its experimentally estimated genome size. The assembly contains 22 chromosome-scale scaffolds (96.52% of the total length), which coincides with its known karyotype. A total of 72,767 transcripts were predicted from the nuclear genome, 24,171 coding, and 48,596 non-coding that included lncRNA, snoRNA, and tRNAs. The circularized mitochondrial genome had 15,740 bp comprising 13 protein-coding genes, 2 rRNA, and 22 tRNA. This reference genome will enhance ongoing A. lixula studies and benefit the wider sea urchin scientific community.

Transgenerational acclimation to acidified seawater and gene expression patterns in a sea urchin.

Transgenerational responses of susceptible calcifying organisms to progressive ocean acidification are an important issue in reducing uncertainty of future predictions. In this study, a two-generation rearing experiment was conducted using mature Mesocentrotus nudus, a major edible sea urchin that occurs along the coasts of northern Japan. Morphological observations and comprehensive gene expression analysis (RNA-seq) of resulting larvae were performed to examine transgenerational acclimation to acidified seawater. Two generations of rearing experiments showed that larvae derived from parents acclimated to acidified seawater tended to have higher survival and show less reduction in body size when exposed to acidified seawater of the same pH, suggesting that a positive carry-over effect occurred. RNA-seq analysis showed that gene expression patterns of larvae originated from both acclimated and non-acclimated parents to acidified seawater tended to be different than control condition, and the gene expression pattern of larvae originated from acclimated parents was substantially different than that of larvae of non-acclimated and control parents.

Genomic signatures of exceptional longevity and negligible aging in the long-lived red sea urchin.

The red sea urchin (Mesocentrotus franciscanus) is one of the Earth's longest-living animals, reported to live more than 100 years with indeterminate growth, life-long reproduction, and no increase in mortality rate with age. To understand the genetic underpinnings of longevity and negligible aging, we constructed a chromosome-level assembly of the red sea urchin genome and compared it to that of short-lived sea urchin species. Genome-wide syntenic alignments identified chromosome rearrangements that distinguish short- and long-lived species. Expanded gene families in long-lived species play a role in innate immunity, sensory nervous system, and genome stability. An integrated network of genes under positive selection in the red sea urchin was involved in genomic regulation, mRNA fidelity, protein homeostasis, and mitochondrial function. Our results implicated known longevity genes in sea urchin longevity but also revealed distinct molecular signatures that may promote long-term maintenance of tissue homeostasis, disease resistance, and negligible aging.

CRISPR-Cas9-Mediated Gene Knockout in a Non-Model Sea Urchin, Heliocidaris crassispina.

Sea urchins have been used as model organisms in developmental biology research and the genomes of several sea urchin species have been sequenced. Recently, genome editing technologies have become available for sea urchins, and methods for gene knockout using the CRISPRCas9 system have been established. Heliocidaris crassispina is an important marine fishery resource with edible gonads. Although H. crassispina has been used as a biological research material, its genome has not yet been published, and it is a non-model sea urchin for molecular biology research. However, as recent advances in genome editing technology have facilitated genome modification in non-model organisms, we applied genome editing using the CRISPR-Cas9 system to H. crassispina. In this study, we targeted genes encoding ETS transcription factor (HcEts) and pigmentation-related polyketide synthase (HcPks1). Gene fragments were isolated using primers designed by inter-specific sequence comparisons within Echinoidea. When Ets gene was targeted using two sgRNAs, one successfully introduced mutations and impaired skeletogenesis. In the Pks1 gene knockout, when two sgRNAs targeting the close vicinity of the site corresponding to the target site that showed 100% mutagenesis efficiency of the Pks1 gene in Hemicentrotus pulcherrimus, mutagenesis was not observed. However, two other sgRNAs targeting distant sites efficiently introduced mutations. In addition, Pks1 knockout H. crassispina exhibited an albino phenotype in the pluteus larvae and adult sea urchins after metamorphosis. This indicates that the CRISPRCas9 system can be used to modify the genome of the non-model sea urchin H. crassispina.

Transglobal spread of an ecologically relevant sea urchin parasite.

Mass mortality of the dominant coral reef herbivore Diadema antillarum in the Caribbean in the early 1980s contributed to a persistent phase shift from coral- to algal-dominated reefs. In 2022, a scuticociliate most closely related to Philaster apodigitiformis caused further mass mortality of D. antillarum across the Caribbean, leading to >95% mortality at affected sites. Mortality was also reported in the related species Diadema setosum in the Mediterranean in 2022, though the causative agent of the Mediterranean outbreak has not yet been determined. In April 2023, mass mortality of Diadema setosum occurred along the Sultanate of Oman's coastline. Urchins displayed signs compatible with scuticociliatosis including abnormal behavior, drooping and loss of spines, followed by tissue necrosis and death. Here we report the detection of an 18S rRNA gene sequence in abnormal urchins from Muscat, Oman, that is identical to the Philaster strain responsible for D. antillarum mass mortality in the Caribbean. We also show that scuticociliatosis signs can be elicited in Diadema setosum by experimental challenge with the cultivated Philaster strain associated with Caribbean scuticociliatosis. These results demonstrate the Philaster sp. associated with D. antillarum mass mortality has rapidly spread to geographically distant coral reefs, compelling global-scale awareness and monitoring for this devastating condition through field surveys, microscopy, and molecular microbiological approaches, and prompting investigation of long-range transmission mechanisms.

Local Genomic Instability of the SpTransformer Gene Family in the Purple Sea Urchin Inferred from BAC Insert Deletions.

The SpTransformer (SpTrf) gene family in the purple sea urchin, Strongylocentrotus purpuratus, encodes immune response proteins. The genes are clustered, surrounded by short tandem repeats, and some are present in genomic segmental duplications. The genes share regions of sequence and include repeats in the coding exon. This complex structure is consistent with putative local genomic instability. Instability of the SpTrf gene cluster was tested by 10 days of growth of Escherichia coli harboring bacterial artificial chromosome (BAC) clones of sea urchin genomic DNA with inserts containing SpTrf genes. After the growth period, the BAC DNA inserts were analyzed for size and SpTrf gene content. Clones with multiple SpTrf genes showed a variety of deletions, including loss of one, most, or all genes from the cluster. Alternatively, a BAC insert with a single SpTrf gene was stable. BAC insert instability is consistent with variations in the gene family composition among sea urchins, the types of SpTrf genes in the family, and a reduction in the gene copy number in single coelomocytes. Based on the sequence variability among SpTrf genes within and among sea urchins, local genomic instability of the family may be important for driving sequence diversity in this gene family that would be of benefit to sea urchins in their arms race with marine microbes.

Genome-Wide Comparative Analysis of SRCR Gene Superfamily in Invertebrates Reveals Massive and Independent Gene Expansions in the Sponge and Sea Urchin.

Without general adaptative immunity, invertebrates evolved a vast number of heterogeneous non-self recognition strategies. One of those well-known adaptations is the expansion of the immune receptor gene superfamily coding for scavenger receptor cysteine-rich domain containing proteins (SRCR) in a few invertebrates. Here, we investigated the evolutionary history of the SRCR gene superfamily (SRCR-SF) across 29 metazoan species with an emphasis on invertebrates. We analyzed their domain architectures, genome locations and phylogenetic distribution. Our analysis shows extensive genome-wide duplications of the SRCR-SFs in Amphimedon queenslandica and Strongylocentrotus purpuratus. Further molecular evolution study reveals various patterns of conserved cysteines in the sponge and sea urchin SRCR-SFs, indicating independent and convergent evolution of SRCR-SF expansion during invertebrate evolution. In the case of the sponge SRCR-SFs, a novel motif with seven conserved cysteines was identified. Exon-intron structure analysis suggests the rapid evolution of SRCR-SFs during gene duplications in both the sponge and the sea urchin. Our findings across nine representative metazoans also underscore a heightened expression of SRCR-SFs in immune-related tissues, notably the digestive glands. This observation indicates the potential role of SRCR-SFs in reinforcing distinct immune functions in these invertebrates. Collectively, our results reveal that gene duplication, motif structure variation, and exon-intron divergence might lead to the convergent evolution of SRCR-SF expansions in the genomes of the sponge and sea urchin. Our study also suggests that the utilization of SRCR-SF receptor duplication may be a general and basal strategy to increase immune diversity and tissue specificity for the invertebrates.

Echinobase: a resource to support the echinoderm research community.

Echinobase (www.echinobase.org) is a model organism knowledgebase serving as a resource for the community that studies echinoderms, a phylum of marine invertebrates that includes sea urchins and sea stars. Echinoderms have been important experimental models for over 100 years and continue to make important contributions to environmental, evolutionary, and developmental studies, including research on developmental gene regulatory networks. As a centralized resource, Echinobase hosts genomes and collects functional genomic data, reagents, literature, and other information for the community. This third-generation site is based on the Xenbase knowledgebase design and utilizes gene-centric pages to minimize the time and effort required to access genomic information. Summary gene pages display gene symbols and names, functional data, links to the JBrowse genome browser, and orthology to other organisms and reagents, and tabs from the Summary gene page contain more detailed information concerning mRNAs, proteins, diseases, and protein-protein interactions. The gene pages also display 1:1 orthologs between the fully supported species Strongylocentrotus purpuratus (purple sea urchin), Lytechinus variegatus (green sea urchin), Patiria miniata (bat star), and Acanthaster planci (crown-of-thorns sea star). JBrowse tracks are available for visualization of functional genomic data from both fully supported species and the partially supported species Anneissia japonica (feather star), Asterias rubens (sugar star), and L. pictus (painted sea urchin). Echinobase serves a vital role by providing researchers with annotated genomes including orthology, functional genomic data aligned to the genomes, and curated reagents and data. The Echinoderm Anatomical Ontology provides a framework for standardizing developmental data across the phylum, and knowledgebase content is formatted to be findable, accessible, interoperable, and reusable by the research community.

Spotting disease disrupts the microbiome of infected purple sea urchins, Strongylocentrotus purpuratus.

BACKGROUND: Spotting disease infects a variety of sea urchin species across many different marine locations. The disease is characterized by discrete lesions on the body surface composed of discolored necrotic tissue that cause the loss of all surface appendages within the lesioned area. A similar, but separate disease of sea urchins called bald sea urchin disease (BSUD) has overlapping symptoms with spotting disease, resulting in confusions in distinguishing the two diseases. Previous studies have focus on identifying the underlying causative agent of spotting disease, which has resulted in the identification of a wide array of pathogenic bacteria that vary based on location and sea urchin species. Our aim was to investigate the spotting disease infection by characterizing the microbiomes of the animal surface and various tissues. RESULTS: We collected samples of the global body surface, the lesion surface, lesioned and non-lesioned body wall, and coelomic fluid, in addition to samples from healthy sea urchins. 16S rRNA gene was amplified and sequenced from the genomic DNA. Results show that the lesions are composed mainly of Cyclobacteriaceae, Cryomorphaceae, and a few other taxa, and that the microbial composition of lesions is the same for all infected sea urchins. Spotting disease also alters the microbial composition of the non-lesioned body wall and coelomic fluid of infected sea urchins. In our closed aquarium systems, sea urchins contracted spotting disease and BSUD separately and therefore direct comparisons could be made between the microbiomes from diseased and healthy sea urchins. CONCLUSION: Results show that spotting disease and BSUD are separate diseases with distinct symptoms and distinct microbial compositions.

microRNA-1 regulates sea urchin skeletogenesis by directly targeting skeletogenic genes and modulating components of signaling pathways.

microRNAs are evolutionarily conserved non-coding RNAs that direct post-transcriptional regulation of target transcripts. In vertebrates, microRNA-1 (miR-1) is expressed in muscle and has been found to play critical regulatory roles in vertebrate angiogenesis, a process that has been proposed to be analogous to sea urchin skeletogenesis. Results indicate that both miR-1 inhibitor and miR-1 mimic-injected larvae have significantly less F-actin enriched circumpharyngeal muscle fibers and fewer gut contractions. In addition, miR-1 regulates the positioning of skeletogenic primary mesenchyme cells (PMCs) and skeletogenesis of the sea urchin embryo. Interestingly, the gain-of-function of miR-1 leads to more severe PMC patterning and skeletal branching defects than its loss-of-function. The results suggest that miR-1 directly suppresses Ets1/2, Tbr, and VegfR7 of the skeletogenic gene regulatory network, and Nodal, and Wnt1 signaling components. This study identifies potential targets of miR-1 that impacts skeletogenesis and muscle formation and contributes to a deeper understanding of miR-1's function during development.

Characterization of p53/p63/p73 and Myc expressions during embryogenesis of the sea urchin.

BACKGROUND: Some marine invertebrate organisms are considered not to develop tumors due to unknown mechanisms. To gain an initial insight into how tumor-related genes may be expressed and function during marine invertebrate development, we here leverage sea urchin embryos as a model system and characterize the expressions of Myc and p53/p63/p73 which are reported to function synergistically in mammalian models as an oncogene and tumor suppressor, respectively. RESULTS: During sea urchin embryogenesis, a combo gene of p53/p63/p73 is found to be maternally loaded and decrease after fertilization both in transcript and protein, while Myc transcript and protein are zygotically expressed. p53/p63/p73 and Myc proteins are observed in the cytoplasm and nucleus of every blastomere, respectively, throughout embryogenesis. Both p53/p63/p73 and Myc overexpression results in compromised development with increased DNA damage after the blastula stage. p53/p63/p73 increases the expression of parp1, a DNA repair/cell death marker gene, and suppresses endomesoderm gene expressions. In contrast, Myc does not alter the expression of specification genes or oncogenes yet induces disorganized morphology. CONCLUSIONS: p53/p63/p73 appears to be important for controlling cell differentiation, while Myc induces disorganized morphology yet not through conventional oncogene regulations or apoptotic pathways during embryogenesis of the sea urchin.

Parentage influence on gene expression under acidification revealed through single-embryo sequencing.

The dissolution of anthropogenic carbon dioxide (CO2 ) in seawater has altered its carbonate chemistry in the process of ocean acidification (OA). OA affects the viability of marine species. In particular, calcifying organisms and their early planktonic larval stages are considered vulnerable. These organisms often utilize energy reserves for metabolism rather than growth and calcification as supported by bulk RNA-sequencing (RNA-seq) experiments. Yet, transcriptomic profiling of a bulk sample reflects the average gene expression of the population, neglecting the variations between individuals, which forms the basis for natural selection. Here, we used single-embryo RNA-seq on larval sea urchin Heliocidaris crassispina, which is a commercially and ecologically valuable species in East Asia, to document gene expression changes to OA at an individual and family level. Three paternal half-sibs groups were fertilized and exposed to 3 pH conditions (ambient pH 8.0, 7.7 and 7.4) for 12 h prior to sequencing and oxygen consumption assay. The resulting transcriptomic profile of all embryos can be distinguished into four clusters, with differences in gene expressions that govern biomineralization, cell differentiation and patterning, as well as metabolism. While these responses were influenced by pH conditions, the male identities also had an effect. Specifically, a regression model and goodness of fit tests indicated a significant interaction between sire and pH on the probability of embryo membership in different clusters of gene expression. The single-embryo RNA-seq approach is promising in climate stressor research because not only does it highlight potential impacts before phenotypic changes were observed, but it also highlights variations between individuals and lineages, thus enabling a better determination of evolutionary potential.

Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins.

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species.

Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.

Single-Cell Transcriptomic Analysis Reveals the Molecular Profile of Go-Opsin Photoreceptor Cells in Sea Urchin Larvae.

The ability to perceive and respond to light stimuli is fundamental not only for spatial vision but also to many other light-mediated interactions with the environment. In animals, light perception is performed by specific cells known as photoreceptors and, at molecular level, by a group of GPCRs known as opsins. Sea urchin larvae possess a group of photoreceptor cells (PRCs) deploying a Go-Opsin (Opsin3.2) which have been shown to share transcription factors and morphology with PRCs of the ciliary type, raising new questions related to how this sea urchin larva PRC is specified and whether it shares a common ancestor with ciliary PRCs or it if evolved independently through convergent evolution. To answer these questions, we combined immunohistochemistry and fluorescent in situ hybridization to investigate how the Opsin3.2 PRCs develop in the sea urchin Strongylocentrotus purpuratus larva. Subsequently, we applied single-cell transcriptomics to investigate the molecular signature of the Sp-Opsin3.2-expressing cells and show that they deploy an ancient regulatory program responsible for photoreceptors specification. Finally, we also discuss the possible functions of the Opsin3.2-positive cells based on their molecular fingerprint, and we suggest that they are involved in a variety of signaling pathways, including those entailing the thyrotropin-releasing hormone.

Widespread priming of transcriptional regulatory elements by incipient accessibility or RNA polymerase II pause in early embryos of the sea urchin Strongylocentrotus purpuratus.

Transcriptional regulatory elements (TREs) are the primary nodes that control developmental gene regulatory networks. In embryo stages, larvae, and adult differentiated red spherule cells of the sea urchin Strongylocentrotus purpuratus, transcriptionally engaged TREs are detected by Precision Run-On Sequencing (PRO-seq), which maps genome-wide at base pair resolution the location of paused or elongating RNA polymerase II (Pol II). In parallel, TRE accessibility is estimated by the Assay for Transposase-Accessible Chromatin using Sequencing (ATAC-seq). Our analysis identifies surprisingly early and widespread TRE accessibility in 4-cell cleavage embryos that is not necessarily followed by concurrent or subsequent transcription. TRE transcriptional differences identified by PRO-seq provide more contrast among embryonic stages than ATAC-seq accessibility differences, in agreement with the apparent excess of accessible but inactive TREs during embryogenesis. Global TRE accessibility reaches a maximum around the 20-hour late blastula stage, which coincides with the consolidation of major embryo regionalizations and peak histone variant H2A.Z expression. A transcriptional potency model based on labile nucleosome TRE occupancy driven by DNA sequences and the prevalence of histone variants is proposed in order to explain the basal accessibility of transcriptionally inactive TREs during embryogenesis. However, our results would not reconcile well with labile nucleosome models based on simple A/T sequence enrichment. In addition, a large number of distal TREs become transcriptionally disengaged during developmental progression, in support of an early Pol II paused model for developmental gene regulation that eventually resolves in transcriptional activation or silencing. Thus, developmental potency in early embryos may be facilitated by incipient accessibility and transcriptional pause at TREs.