Genomic profiles and clinical presentation of chordoma.

Chordoma is a rare bone cancer with variable clinical outcomes. Here, we recruited 184 sporadic chordoma patients from the US and Canada and collected their clinical and treatment data. The average age at diagnosis was 45.5 years (Range 5-78) and the chordoma site distribution was 49.2% clivus, 26.2% spinal, and 24.0% sacral. Most patients (97.5%) received surgery as the primary treatment, among whom 85.3% also received additional treatment. Except for the most prevalent cancers like prostate, lung, breast, and skin cancer, there was no discernible enrichment for any specific cancer type among patients or their family members. Among a subset of patients (N = 70) with tumor materials, we conducted omics analyses and obtained targeted panel sequencing and SNP array genotyping data for 51 and 49 patients, respectively. The most recurrent somatic driver mutations included PIK3CA (12%), followed by chromatin remodeling genes PBRM1 and SETD2. Amplification of the 6q27 region, containing the chordoma susceptibility gene TBXT, was detected in eight patients (16.3%). Clival patients appeared to be less likely to carry driver gene mutations, chromosome arm level deletion events (e.g., 5p, 5p, and 9p), or 6q27 amplification compared to sacral patients. After adjusting for age, sex, tumor site, and additional treatment, patients with somatic deletions of 14q (OR = 13.73, 95% CI 1.96-96.02, P = 0.008) and 18p (OR = 13.68, 95% CI 1.77-105.89, P = 0.012) were more likely to have persistent chordoma. The study highlights genomic heterogeneity in chordoma, potentially linked to location and clinical progression.

ALCAM-mediated cDC1 CD8 T cells interactions are suppressed in advanced lung tumors.

Conventional type 1 dendritic cells (cDC1) are critical regulators of anti-tumoral T-cell responses. The structure and abundance of intercellular contacts between cDC1 and CD8 T cells in cancer tissues is important to determine the outcome of the T-cell response. However, the molecular determinants controlling the stability of cDC1-CD8 interactions during cancer progression remain poorly investigated. Here, we generated a genetic model of non-small cell lung cancer crossed to a fluorescent cDC1 reporter (KP-XCR1venus) to allow the detection of cDC1-CD8T cell clusters in tumor tissues across tumor stages. We found that cDC1-CD8 clusters are abundant and productive at the early stages of tumor development but progressively diminish in advanced tumors. Transcriptional profiling and flow cytometry identified the adhesion molecule ALCAM/CD166 (Activated Leukocyte Cell Adhesion Molecule, ligand of CD6) as highly expressed by lung cDC1 and significantly downregulated in advanced tumors. Analysis of human datasets indicated that ALCAM is downregulated in non-small cell lung cancer and its expression correlates to better prognosis. Mechanistically, triggering ALCAM on lung cDC1 induces cytoskeletal remodeling and contact formation whereas its blockade prevents T-cell activation. Together, our results indicate that ALCAM is important to stabilize cDC1-CD8 interactions at early tumor stages, while its loss in advanced tumors contributes to immune evasion.

Longitudinal assessment of urinary ALCAM, HPX, and PRDX6 in Korean patients with systemic lupus erythematosus: implications for disease activity monitoring and treatment response.

Introduction: This study aimed to demonstrate the potential of activated leukocyte cell adhesion molecule (ALCAM), hemopexin (HPX), and peroxiredoxin 6 (PRDX6) as urine biomarkers for systemic lupus erythematosus (SLE). Methods: Urine samples were collected from 138 Korean patients with SLE from the Ajou Lupus Cohort and 39 healthy controls (HC). The concentrations of urine biomarkers were analyzed using enzyme-linked immunosorbent assay kits specific for ALCAM, HPX, and PRDX6, respectively. Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic utility, and Pearson's correlation analysis was conducted to assess the relationships between the disease activity and urine biomarkers. Results: Patients with SLE and patients with lupus nephritis (LN) showed significantly elevated ALCAM, HPX, and PRDX6 levels compared with HCs. ALCAM, HPX, and PRDX6 showed significant diagnostic values, especially for lupus nephritis (LN), with areas under the receiver operating characteristic curve for LN was 0.850 for ALCAM (95% CI, 0.778-0.921), 0.781 for HPX (95% CI, 0.695-0.867), and 0.714 for PRDX6 (95% CI, 0.617-0.812). Correlation analysis revealed that all proteins were significantly associated with anti-double stranded DNA antibody (ALCAM, r = 0.350, p < 0.001; HPX, r = 0.346, p < 0.001; PRDX6, r = 0.191, p = 0.026) and SLEDAI (ALCAM, r = 0.526, p < 0.001; HPX, r = 0.479, p < 0.001; PRDX6, r = 0.262, p = 0.002). Results from the follow-up of the three biomarker levels in these patients revealed a significant decrease, showing a positive correlation with changes in SLEDAI-2k scores (ALCAM, r = 0.502, p < 0.001; HPX, r = 0.475, p < 0.001; PRDX6, r = 0.245, p = 0.026), indicating their potential as indicators for tracking disease activity. Discussions: Urinary ALCAM, HPX, and PRDX6 levels have diagnostic value and reflect disease activity in Korean patients with SLE, emphasizing their potential for non-invasive monitoring and treatment response evaluation.

Chordoma: Genetics and Contemporary Management.

Chordomas, arising from notochord remnants, are rare neoplasms with aggressive growth patterns despite their histologically low-grade nature. This review explores their embryological origins, molecular markers like brachyury, and genetic alterations driving pathogenesis. Diagnosis relies on advanced imaging and biopsy confirmation due to overlapping features with chondrosarcoma. The WHO classification distinguishes conventional, dedifferentiated, and poorly differentiated chordomas, each with distinct prognostic implications. Recent genomic analyses uncovered recurrent mutations in PI3K signaling pathways and chromatin remodeling genes, informing prognostic models. Surgery remains the cornerstone of treatment, though adjuvant radiation complements surgical resection. Although chordomas are generally considered refractory to medical therapy, emerging targeted molecular strategies show potential promise in ongoing trials. This review aims to provide a concise yet comprehensive overview of chordomas, guiding clinicians in diagnosis, treatment, and prognostication for improved patient outcomes.

Unraveling molecular advancements in chordoma tumorigenesis and treatment response: a review of scientific discoveries and clinical implications.

Chordomas are tumors thought to originate from notochordal remnants that occur in midline structures from the cloves of the skull base to the sacrum. In adults, the most common location is the sacrum, followed by the clivus and then mobile spine, while in children a clival origin is most common. Most chordomas are slow growing. Clinical presentation of chordomas tend to occur late, with local invasion and large size often complicating surgical intervention. Radiation therapy with protons has been proven to be an effective adjuvant therapy. Unfortunately, few adjuvant systemic treatments have demonstrated significant effectiveness, and chordomas tend to recur despite intensive multimodal care. However, insight into the molecular underpinnings of chordomas may guide novel therapeutic approaches including selection for immune and molecular therapies, individualized prognostication of outcomes, and real-time noninvasive assessment of disease burden and evolution. At the genomic level, elevated levels of brachyury stemming from duplications and mutations resulting in altered transcriptional regulation may introduce druggable targets for new surgical adjuncts. Transcriptome and epigenome profiling have revealed promoter- and enhancer-dependent mechanisms of protein regulation, which may influence therapeutic response and long-term disease history. Continued scientific and clinical advancements may offer further opportunities for treatment of chordomas. Single-cell transcriptome profiling has further provided insight into the heterogeneous molecular pathways contributing to chordoma propagation. New technologies such as spatial transcriptomics and emerging biochemical analytes such as cell-free DNA have further augmented the surgeon-clinician's armamentarium by facilitating detailed characterization of intra- and intertumoral biology while also demonstrating promise for point-of-care tumor quantitation and assessment. Recent and ongoing clinical trials highlight accelerating interest to translate laboratory breakthroughs in chordoma biology and immunology into clinical care. In this review, the authors dissect the landmark studies exploring the molecular pathogenesis of chordoma. Incorporating this into an outline of ongoing clinical trials and discussion of emerging technologies, the authors aimed to summarize recent advancements in understanding chordoma pathogenesis and how neurosurgical care of chordomas may be augmented by improvements in adjunctive treatments.

Brachyury promotes extracellular matrix synthesis through transcriptional regulation of Smad3 in nucleus pulposus.

Metabolic dysfunction of the extracellular matrix (ECM) is one of the primary causes of intervertebral disc degeneration (IVDD). Previous studies have demonstrated that the transcription factor Brachyury (Bry) has the potential to promote the synthesis of collagen II and aggrecan, while the specific mechanism is still unknown. In this study, we used a lipopolysaccharide (LPS)-induced model of nucleus pulposus cell (NPC) degeneration and a rat acupuncture IVDD model to elucidate the precise mechanism through which Bry affects collagen II and aggrecan synthesis in vitro and in vivo. First, we confirmed Bry expression decreased in degenerated human nucleus pulposus (NP) cells (NPCs). Knockdown of Bry exacerbated the decrease in collagen II and aggrecan expression in the lipopolysaccharide (LPS)-induced NPCs degeneration in vitro model. Bioinformatic analysis indicated that Smad3 may participate in the regulatory pathway of ECM synthesis regulated by Bry. Chromatin immunoprecipitation followed by quantitative polymerase chain reaction (ChIP-qPCR) and luciferase reporter gene assays demonstrated that Bry enhances the transcription of Smad3 by interacting with a specific motif on the promoter region. In addition, Western blot and reverse transcription-qPCR assays demonstrated that Smad3 positively regulates the expression of aggrecan and collagen II in NPCs. The following rescue experiments revealed that Bry-mediated regulation of ECM synthesis is partially dependent on Smad3 phosphorylation. Finally, the findings from the in vivo rat acupuncture-induced IVDD model were consistent with those obtained from in vitro assays. In conclusion, this study reveals that Bry positively regulates the synthesis of collagen II and aggrecan in NP through transcriptional activation of Smad3.NEW & NOTEWORTHY Mechanically, in the nucleus, Bry enhances the transcription of Smad3, leading to increased expression of Smad3 protein levels; in the cytoplasm, elevated substrate levels further lead to an increase in the phosphorylation of Smad3, thereby regulating collagen II and aggrecan expression. Further in vivo experiments provide additional evidence that Bry can alleviate IVDD through this mechanism.

Germline findings in cancer predisposing genes from a small cohort of chordoma patients.

INTRODUCTION: Chordoma is a rare slow-growing tumor that occurs along the length of the spinal axis and arises from primitive notochordal remnants (Stepanek et al., Am J Med Genet 75:335-336, 1998). Most chordomas are sporadic, but a small percentage of cases are due to hereditary cancer syndromes (HCS) such as tuberous sclerosis 1 and 2 (TSC1/2), or constitutional variants in the gene encoding brachyury T (TBXT) (Pillay et al., Nat Genet 44:1185-1187, 2012; Yang et al., Nat Genet 41:1176-1178, 2009). PURPOSE: The genetic susceptibility of these tumors is not well understood; there are only a small number of studies that have performed germline genetic testing in this population. METHODS: We performed germline genetic in chordoma patients using genomic DNA extracted by blood or saliva. CONCLUSION: We report here a chordoma cohort of 24 families with newly found germline genetic mutations in cancer predisposing genes. We discuss implications for genetic counseling, clinical management, and universal germline genetic testing for cancer patients with solid tumors.

Basal delamination during mouse gastrulation primes pluripotent cells for differentiation.

The blueprint of the mammalian body plan is laid out during gastrulation, when a trilaminar embryo is formed. This process entails a burst of proliferation, the ingression of embryonic epiblast cells at the primitive streak, and their priming toward primitive streak fates. How these different events are coordinated remains unknown. Here, we developed and characterized a 3D culture of self-renewing mouse embryonic cells that captures the main transcriptional and architectural features of the early gastrulating mouse epiblast. Using this system in combination with microfabrication and in vivo experiments, we found that proliferation-induced crowding triggers delamination of cells that express high levels of the apical polarity protein aPKC. Upon delamination, cells become more sensitive to Wnt signaling and upregulate the expression of primitive streak markers such as Brachyury. This mechanistic coupling between ingression and differentiation ensures that the right cell types become specified at the right place during embryonic development.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona.

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

Snail Transcriptionally Represses Brachyury to Promote the Mesenchymal-Epithelial Transition in Ascidian Notochord Cells.

Mesenchymal-epithelial transition (MET) is a widely spread and evolutionarily conserved process across species during development. In Ciona embryogenesis, the notochord cells undergo the transition from the non-polarized mesenchymal state into the polarized endothelial-like state to initiate the lumen formation between adjacent cells. Based on previously screened MET-related transcription factors by ATAC-seq and Smart-Seq of notochord cells, Ciona robusta Snail (Ci-Snail) was selected for its high-level expression during this period. Our current knockout results demonstrated that Ci-Snail was required for notochord cell MET. Importantly, overexpression of the transcription factor Brachyury in notochord cells resulted in a similar phenotype with failure of lumen formation and MET. More interestingly, expression of Ci-Snail in the notochord cells at the late tailbud stage could partially rescue the MET defect caused by Brachyury-overexpression. These results indicated an inverse relationship between Ci-Snail and Brachyury during notochord cell MET, which was verified by RT-qPCR analysis. Moreover, the overexpression of Ci-Snail could significantly inhibit the transcription of Brachyury, and the CUT&Tag-qPCR analysis demonstrated that Ci-Snail is directly bound to the upstream region of Brachyury. In summary, we revealed that Ci-Snail promoted the notochord cell MET and was essential for lumen formation via transcriptionally repressing Brachyury.

Development of a Small Molecule Downmodulator for the Transcription Factor Brachyury.

Brachyury is an oncogenic transcription factor whose overexpression drives chordoma growth. The downmodulation of brachyury in chordoma cells has demonstrated therapeutic potential, however, as a transcription factor it is classically deemed "undruggable". Given that direct pharmacological intervention against brachyury has proven difficult, attempts at intervention have instead targeted upstream kinases. Recently, afatinib, an FDA-approved kinase inhibitor, has been shown to modulate brachyury levels in multiple chordoma cell lines. Herein, we use afatinib as a lead to undertake a structure-based drug design approach, aided by mass-spectrometry and X-ray crystallography, to develop DHC-156, a small molecule that more selectively binds brachyury and downmodulates it as potently as afatinib. We eliminated kinase-inhibition from this novel scaffold while demonstrating that DHC-156 induces the post-translational downmodulation of brachyury that results in an irreversible impairment of chordoma tumor cell growth. In doing so, we demonstrate the feasibility of direct brachyury modulation, which may further be developed into more potent tool compounds and therapies.

Integrating Clinical Cancer and PTM Proteomics Data Identifies a Mechanism of ACK1 Kinase Activation.

Beyond the most common oncogenes activated by mutation (mut-drivers), there likely exists a variety of low-frequency mut-drivers, each of which is a possible frontier for targeted therapy. To identify new and understudied mut-drivers, we developed a machine learning (ML) model that integrates curated clinical cancer data and posttranslational modification (PTM) proteomics databases. We applied the approach to 62,746 patient cancers spanning 84 cancer types and predicted 3,964 oncogenic mutations across 1,148 genes, many of which disrupt PTMs of known and unknown function. The list of putative mut-drivers includes established drivers and others with poorly understood roles in cancer. This ML model is available as a web application. As a case study, we focused the approach on nonreceptor tyrosine kinases (NRTK) and found a recurrent mutation in activated CDC42 kinase-1 (ACK1) that disrupts the Mig6 homology region (MHR) and ubiquitin-association (UBA) domains on the ACK1 C-terminus. By studying these domains in cultured cells, we found that disruption of the MHR domain helps activate the kinase while disruption of the UBA increases kinase stability by blocking its lysosomal degradation. This ACK1 mutation is analogous to lymphoma-associated mutations in its sister kinase, TNK1, which also disrupt a C-terminal inhibitory motif and UBA domain. This study establishes a mut-driver discovery tool for the research community and identifies a mechanism of ACK1 hyperactivation shared among ACK family kinases. IMPLICATIONS: This research identifies a potentially targetable activating mutation in ACK1 and other possible oncogenic mutations, including PTM-disrupting mutations, for further study.

Methylated ctDNA Quantification: Noninvasive Approach to Monitoring Hepatocellular Carcinoma Burden.

BACKGROUND: Noninvasive, precision monitoring of hepatocellular carcinoma (HCC) treatment efficacy would greatly facilitate personalized therapy and improve patient outcomes. We hypothesize that quantifying methylated circulating tumor DNA (ctDNA) can be used to effectively monitor HCC burden without the need for biopsy. STUDY DESIGN: Blood samples were collected from 25 patients, 21 with HCC and 4 with benign liver masses, at various timepoints throughout the course of treatment at a high-volume academic medical center. Quantification of methylated ctDNA molecules assessed CpG sites on more than 550 preselected cancer-specific amplicons. The tumor methylation score (TMS) was calculated by measuring the difference between the amount of methylation in the plasma and buffy coat with a normal cutoff value of 120 or less. RESULTS: Among 10 patients with surgical HCC (5 surgical resections and 5 liver transplants), TMS revealed a statistically significant, rapid postoperative decline in 9. One patient who had a persistently elevated TMS on postoperative day 1 was subsequently found to have had metastatic disease. Patients in the negative control cohort all had normal-range pre- and postoperative TMS. Preoperative TMS correlated moderately with tumor burden on pathology (Spearman r = 0.54) of surgical specimens. From 11 subjects undergoing systemic therapy or Y90 radioembolization, analysis of 16 time periods demonstrated that the change in TMS (ΔTMS) was better associated with tumor progression than the change in Δalpha-fetoprotein (area under the curve 0.800 and 0.783, respectively). A composite score combining ΔTMS and Δalpha-fetoprotein further improved performance for detecting tumor progression with an area under the curve of 0.892. CONCLUSIONS: These findings indicate that ctDNA methylation scores can effectively evaluate changes in tumor burden without the need for tumor biopsy.

Conserved enhancers control notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

Genome-wide identification of notochord enhancers comprising the regulatory landscape of the brachyury locus in mouse.

The node and notochord are important signaling centers organizing the dorso-ventral patterning of cells arising from neuro-mesodermal progenitors forming the embryonic body anlage. Owing to the scarcity of notochord progenitors and notochord cells, a comprehensive identification of regulatory elements driving notochord-specific gene expression has been lacking. Here, we have used ATAC-seq analysis of FACS-purified notochord cells from Theiler stage 12-13 mouse embryos to identify 8921 putative notochord enhancers. In addition, we established a new model for generating notochord-like cells in culture, and found 3728 of these enhancers occupied by the essential notochord control factors brachyury (T) and/or Foxa2. We describe the regulatory landscape of the T locus, comprising ten putative enhancers occupied by these factors, and confirmed the regulatory activity of three of these elements. Moreover, we characterized seven new elements by knockout analysis in embryos and identified one new notochord enhancer, termed TNE2. TNE2 cooperates with TNE in the trunk notochord, and is essential for notochord differentiation in the tail. Our data reveal an essential role of Foxa2 in directing T-expressing cells towards the notochord lineage.

Fusion crystallization reveals the behavior of both the 1TEL crystallization chaperone and the TNK1 UBA domain.

Human thirty-eight-negative kinase-1 (TNK1) is implicated in cancer progression. The TNK1 ubiquitin-associated (UBA) domain binds polyubiquitin and plays a regulatory role in TNK1 activity and stability. No experimentally determined molecular structure of this unusual UBA domain is available. We fused the UBA domain to the 1TEL variant of the translocation ETS leukemia protein sterile alpha motif (TELSAM) crystallization chaperone and obtained crystals diffracting as far as 1.53 Å. GG and GSGG linkers allowed the UBA to reproducibly find a productive binding mode against its host 1TEL polymer and crystallize at protein concentrations as low as 0.2 mg/mL. Our studies support a mechanism of 1TEL fusion crystallization and show that 1TEL fusion crystals require fewer crystal contacts than traditional protein crystals. Modeling and experimental validation suggest the UBA domain may be selective for both the length and linkages of polyubiquitin chains.

Clinical impact of myosteatosis measured by magnetic resonance imaging on long-term outcomes of hepatocellular carcinoma after radical hepatectomy.

AIMS: A variety of factors have been reported to affect long-term outcomes after radical resection of hepatocellular carcinoma (HCC). However, the indicators remain controversial. The purpose of this study was to evaluate the relationship between myosteatosis of the multifidus muscle and long-term outcomes after radical surgery for HCC. METHODS: We retrospectively analyzed clinicopathological data for 187 patients with HCC who underwent radical surgery at Tokushima University between January 2009 and December 2020 and measured the density of fat in the multifidus muscle at L3 on their preoperative magnetic resonance images (MRI). Associations of myosteatosis and clinicopathological factors with long-term outcomes were evaluated. RESULTS: The patients were divided into a myosteatosis-negative group (n = 122) and a myosteatosis-positive group (n = 65). The cancer-specific survival rate after hepatectomy was significantly worse in the myosteatosis-positive group than in the myosteatosis-negative group (p = 0.03). Univariate analysis identified multiple tumors, stage III/IV disease, an alfa-fetoprotein level ≥ 10 ng/ml, PIVKA-II ≥ 400 AU/ml, vp(+) status, and myosteatosis to be prognostic factors for cancer-specific survival. Multivariate analysis revealed multiple tumors, an alfa-fetoprotein level ≥ 10 ng/ml, and myosteatosis to be independent prognostic factors. CONCLUSIONS: Myosteatosis measured by MRI is a simple and useful predictor of the long-term outcome after radical surgery for HCC.

Eomes restricts Brachyury functions at the onset of mouse gastrulation.

Mammalian specification of mesoderm and definitive endoderm (DE) is instructed by the two related Tbx transcription factors (TFs) Eomesodermin (Eomes) and Brachyury sharing partially redundant functions. Gross differences in mutant embryonic phenotypes suggest specific functions of each TF. To date, the molecular details of separated lineage-specific gene regulation by Eomes and Brachyury remain poorly understood. Here, we combine mouse embryonic and stem-cell-based analyses to delineate the non-overlapping, lineage-specific transcriptional activities. On a genome-wide scale, binding of both TFs overlaps at promoters of target genes but shows specificity for distal enhancer regions that is conferred by differences in Tbx DNA-binding motifs. The unique binding to enhancer sites instructs the specification of anterior mesoderm (AM) and DE by Eomes and caudal mesoderm by Brachyury. Remarkably, EOMES antagonizes BRACHYURY gene regulatory functions in coexpressing cells during early gastrulation to ensure the proper sequence of early AM and DE lineage specification followed by posterior mesoderm derivatives.

Transcriptomic evidence for Brachyury expression in the caudal tip region of adult Ptychodera flava (Hemichordata).

Most metazoans have a single copy of the T-box transcription factor gene Brachyury. This gene is expressed in cells of the blastopore of late blastulae and the archenteron invagination region of gastrulae. It appears to be crucial for gastrulation and mesoderm differentiation of embryos. Although this expression pattern is shared by most deuterostomes, Brachyury expression has not been reported in adult stages. Here we show that Brachyury of an indirect developer, the hemichordate acorn worm Ptychodera flava, is expressed not only in embryonic cells, but also in cells of the caudal tip (anus) region of adults. This spatially restricted expression, shown by whole-mount in situ hybridization, was confirmed by Iso-Seq RNA sequencing and single-cell RNA-seq (scRNA-seq) analysis. Iso-Seq analysis showed that gene expression occurs only in the caudal region of adults, but not in anterior regions, including the stomochord. scRNA-seq analysis showed a cluster that contained Brachyury-expressing cells comprising epidermis- and mesoderm-related cells, but which is unlikely to be associated with the nervous system or muscle. Although further investigation is required to examine the roles of Brachyury in adults, this study provides important clues for extending studies on Brachyury expression involved in development of the most posterior region of deuterostomes.

The evolutionary history of Brachyury genes in Hydrozoa involves duplications, divergence, and neofunctionalization.

Brachyury, a member of T-box gene family, is widely known for its major role in mesoderm specification in bilaterians. It is also present in non-bilaterian metazoans, such as cnidarians, where it acts as a component of an axial patterning system. In this study, we present a phylogenetic analysis of Brachyury genes within phylum Cnidaria, investigate differential expression and address a functional framework of Brachyury paralogs in hydrozoan Dynamena pumila. Our analysis indicates two duplication events of Brachyury within the cnidarian lineage. The first duplication likely appeared in the medusozoan ancestor, resulting in two copies in medusozoans, while the second duplication arose in the hydrozoan ancestor, resulting in three copies in hydrozoans. Brachyury1 and 2 display a conservative expression pattern marking the oral pole of the body axis in D. pumila. On the contrary, Brachyury3 expression was detected in scattered presumably nerve cells of the D. pumila larva. Pharmacological modulations indicated that Brachyury3 is not under regulation of cWnt signaling in contrast to the other two Brachyury genes. Divergence in expression patterns and regulation suggest neofunctionalization of Brachyury3 in hydrozoans.