Parallel evolution of amphioxus and vertebrate small-scale gene duplications

Characterization of Myxovirus resistance (Mx) gene from Chinese seabass Lateolabrax maculatus: Insights into the evolution and function of Mx genes

Chinese seabass (Lateolabrax maculatus) stands out as one of the most sought-after and economically significant species in aquaculture within China. Diseases of L. maculatus occur frequently due to the degradation of the germplasm, the aggravation of environmental pollution of water, and the reproduction of pathogenic microorganisms, inflicting considerable economic losses on the Chinese seabass industry. The Myxovirus resistance (Mx) gene plays pivotal roles in the antiviral immune response ranging from mammals to fish. However, the function of the Mx gene in L. maculatus is still unknown. Firstly, the origin and evolutionary history of Mx proteins was elucidated in this study. Subsequently, an Mx gene from L. maculatus (designed as LmMxA gene) was identified, and its functions in combating antiviral and antibacterial threats were investigated. Remarkably, our findings suggested that while Mx group genes were present in chordates, DYN group genes were present in everything from single-celled animals to humans. Furthermore, our investigation revealed that the LmMxA mRNA level increased in the kidney, spleen and liver subsequent to Vibrio anguillarum and poly(I:C) challenged. Immunofluorescence analysis indicated that LmMxA is predominantly localization in the nucleus and the cytoplasm. Notably, the expression of MAVS, IFN1 and Mx1 increased when LmMxA was overexpression within the EPC cells. Moreover, through assessment via cytopathic effect (CPE), virus titer, and antibacterial activity, it becomes evident that LmMxA exerts a dual role in bolstering both antiviral and antibacterial immune responses. These compelling findings laid the foundation for further exploring the mechanism of LmMxA in response to innate immunity of L. maculatus.

A novel C-type lectin protein (BjCTL5) interacts with apoptosis stimulating proteins of p53 (ASPP) to activate NF-κB signaling pathway in primitive chordate

C-type lectin proteins (CTLs), a group of pattern recognition receptors (PRRs), play pivotal roles in immune responses. However, the signal transduction and regulation of CTLs in cephalochordates have yet to be explored. In this study, we examined the composition of CTLs in Branchiostoma japonicum, identifying a total of 272 CTLs. These CTLs underwent further analysis concerning domain arrangement, tandem and segmental duplication events. A multidomain C-type lectin gene, designated as BjCTL5, encompassing CLECT, KR, CUB, MAM, and SR domains, was the focal point of our investigation. BjCTL5 exhibits ubiquitous expression across all detected tissues and is responsive to stimulation by LPS, mannose, and poly (I:C). The recombinant protein of BjCTL5 can bind to Escherichia coli and Staphylococcus aureus, inducing their agglutination and inhibiting the proliferation of S. aureus. Yeast two-hybrid, CoIP, and confocal immunofluorescence experiments revealed the interaction between BjCTL5 and apoptosis-stimulating proteins of p53, BjASPP. Intriguingly, BjCTL5 was observed to induce the luciferase activity of the NF-κB promoter in HEK293T cells. These results suggested a potential interaction between BjCTL5 and BjASPP, implicating that they involve in the activation of the NF-κB signaling pathway, which provides an evolutionary viewpoint on NF-κB signaling pathway in primitive chordate.

Conservation of cis-Regulatory Syntax Underlying Deuterostome Gastrulation

Throughout embryonic development, the shaping of the functional and morphological characteristics of embryos is orchestrated by an intricate interaction between transcription factors and cis-regulatory elements. In this study, we conducted a comprehensive analysis of deuterostome cis-regulatory landscapes during gastrulation, focusing on four paradigmatic species: the echinoderm Strongylocentrotus purpuratus, the cephalochordate Branchiostoma lanceolatum, the urochordate Ciona intestinalis, and the vertebrate Danio rerio. Our approach involved comparative computational analysis of ATAC-seq datasets to explore the genome-wide blueprint of conserved transcription factor binding motifs underlying gastrulation. We identified a core set of conserved DNA binding motifs associated with 62 known transcription factors, indicating the remarkable conservation of the gastrulation regulatory landscape across deuterostomes. Our findings offer valuable insights into the evolutionary molecular dynamics of embryonic development, shedding light on conserved regulatory subprograms and providing a comprehensive perspective on the conservation and divergence of gene regulation underlying the gastrulation process.

Genomic Resources and Annotations for a Colonial Ascidian, the Light-Bulb Sea Squirt Clavelina lepadiformis

Ascidian embryos have been studied since the birth of experimental embryology at the end of the 19th century. They represent textbook examples of mosaic development characterized by a fast development with very few cells and invariant cleavage patterns and lineages. Ascidians belong to tunicates, the vertebrate sister group, and their study is essential to shed light on the emergence of vertebrates. Importantly, deciphering developmental gene regulatory networks has been carried out mostly in two of the three ascidian orders, Phlebobranchia and Stolidobranchia. To infer ancestral developmental programs in ascidians, it is thus essential to carry out molecular embryology in the third ascidian order, the Aplousobranchia. Here, we present genomic resources for the colonial aplousobranch Clavelina lepadiformis: a transcriptome produced from various embryonic stages, and an annotated genome. The assembly consists of 184 contigs making a total of 233.6 Mb with a N50 of 8.5 Mb and a L50 of 11. The 32,318 predicted genes capture 96.3% of BUSCO orthologs. We further show that these resources are suitable to study developmental gene expression and regulation in a comparative framework within ascidians. Additionally, they will prove valuable for evolutionary and ecological studies.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Insights into RAG Evolution from the Identification of "Missing Link" Family A RAGL Transposons

A series of "molecular domestication" events are thought to have converted an invertebrate RAG-like (RAGL) transposase into the RAG1-RAG2 (RAG) recombinase, a critical enzyme for adaptive immunity in jawed vertebrates. The timing and order of these events are not well understood, in part because of a dearth of information regarding the invertebrate RAGL-A transposon family. In contrast to the abundant and divergent RAGL-B transposon family, RAGL-A most closely resembles RAG and is represented by a single orphan RAG1-like (RAG1L) gene in the genome of the hemichordate Ptychodera flava (PflRAG1L-A). Here, we provide evidence for the existence of complete RAGL-A transposons in the genomes of P. flava and several echinoderms. The predicted RAG1L-A and RAG2L-A proteins encoded by these transposons intermingle sequence features of jawed vertebrate RAG and RAGL-B transposases, leading to a prediction of DNA binding, catalytic, and transposition activities that are a hybrid of RAG and RAGL-B. Similarly, the terminal inverted repeats (TIRs) of the RAGL-A transposons combine features of both RAGL-B transposon TIRs and RAG recombination signal sequences. Unlike all previously described RAG2L proteins, RAG2L-A proteins contain an acidic hinge region, which we demonstrate is capable of efficiently inhibiting RAG-mediated transposition. Our findings provide evidence for a critical intermediate in RAG evolution and argue that certain adaptations thought to be specific to jawed vertebrates (e.g. the RAG2 acidic hinge) actually arose in invertebrates, thereby focusing attention on other adaptations as the pivotal steps in the completion of RAG domestication in jawed vertebrates.

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

Amphioxus as a model to study the evolution of development in chordates

Cephalochordates and tunicates represent the only two groups of invertebrate chordates, and extant cephalochordates - commonly known as amphioxus or lancelets - are considered the best proxy for the chordate ancestor, from which they split around 520 million years ago. Amphioxus has been an important organism in the fields of zoology and embryology since the 18th century, and the morphological and genomic simplicity of cephalochordates (compared to vertebrates) makes amphioxus an attractive model for studying chordate biology at the cellular and molecular levels. Here we describe the life cycle of amphioxus, and discuss the natural histories and habitats of the different species of amphioxus. We also describe their use as laboratory animal models, and discuss the techniques that have been developed to study different aspects of amphioxus.

Genome-wide transcriptomics and microRNAomics analyses uncover multi-faceted mechanisms to cope with copper stress in ancient macrobenthos amphioxus

The mechanisms underlying the toxicity of environmental stress are unclear for marine macrobenthos. Copper/Cu has posed the most serious threats to amphioxus, an ancient and model benthic cephalochordate. Herein, a dynamic change in the physiological parameters (GR, SOD, ATP, and MDA) was detected with ROS accumulation in Branchiostoma belcheri exposed to 0.3 mg·L-1 Cu. Transcriptomes and microRNAomes of B. belcheri were generated to investigate the molecular mechanisms by which this amphioxus copes with Cu exposure. Time-specific genes identified at different time points after exposure were involved in the stimulus and immune response, detoxification and ionic homeostasis, aging and the nervous system, sequentially, with prolongation of exposure time, forming a dynamic process of molecular response to Cu stress. In total, 57 differentially expressed miRNAs were identified under Cu stress. Transcriptomics-miRNAomics analyses indicate that these miRNAs targeted genes associated with many key biological processes such as xenobiotics degradation, oxidative stress, and energy metabolism. The constructed miRNA-mRNA-pathway network uncovered a broad post-transcriptional regulatory mechanism in B. belcheri to cope with Cu stress. Overall, this integrated analyses show that enhanced defense response, accelerated ROS elimination, and repressed ATP production constitute a comprehensive strategy to cope with Cu toxicity in the ancient macrobenthos.

Insights into RAG evolution from the identification of "missing link" family A RAGL transposons

A series of "molecular domestication" events are thought to have converted an invertebrate RAG-like (RAGL) transposase into the RAG1-RAG2 (RAG) recombinase, a critical enzyme for adaptive immunity in jawed vertebrates. The timing and order of these events is not well understood, in part because of a dearth of information regarding the invertebrate RAGL-A transposon family. In contrast to the abundant and divergent RAGL-B transposon family, RAGL-A most closely resembles RAG and is represented by a single orphan RAG1-like (RAG1L) gene in the genome of the hemichordate Ptychodera flava (PflRAG1L-A). Here, we provide evidence for the existence of complete RAGL-A transposons in the genomes of P. flava and several echinoderms. The predicted RAG1L-A and RAG2L-A proteins encoded by these transposons intermingle sequence features of jawed vertebrate RAG and RAGL-B transposases, leading to a prediction of DNA binding, catalytic, and transposition activities that are a hybrid of RAG and RAGL-B. Similarly, the terminal inverted repeats (TIRs) of the RAGL-A transposons combine features of both RAGL-B transposon TIRs and RAG recombination signal sequences. Unlike all previously described RAG2L proteins, PflRAG2L-A and echinoderm RAG2L-A contain an acidic hinge region, which we demonstrate is capable of efficiently inhibiting RAG-mediated transposition. Our findings provide evidence for a critical intermediate in RAG evolution and argue that certain adaptations thought to be specific to jawed vertebrates (e.g., the RAG2 acidic hinge) actually arose in invertebrates, thereby focusing attention on other adaptations as the pivotal steps in the completion of RAG domestication in jawed vertebrates.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Retinoic Acid and POU Genes in Developing Amphioxus: A Focus on Neural Development

POU genes are a family of evolutionarily conserved transcription factors with key functions in cell type specification and neurogenesis. In vitro experiments have indicated that the expression of some POU genes is controlled by the intercellular signaling molecule retinoic acid (RA). In this work, we aimed to characterize the roles of RA signaling in the regulation of POU genes in vivo. To do so, we studied POU genes during the development of the cephalochordate amphioxus, an animal model crucial for understanding the evolutionary origins of vertebrates. The expression patterns of amphioxus POU genes were assessed at different developmental stages by chromogenic in situ hybridization and hybridization chain reaction. Expression was further assessed in embryos subjected to pharmacological manipulation of endogenous RA signaling activity. In addition to a detailed description of the effects of these treatments on amphioxus POU gene expression, our survey included the first description of Pou2 and Pou6 expression in amphioxus embryos. We found that Pit-1, Pou2, Pou3l, and Pou6 expression are not affected by alterations of endogenous RA signaling levels. In contrast, our experiments indicated that Brn1/2/4 and Pou4 expression are regulated by RA signaling in the endoderm and the nerve cord, respectively. The effects of the treatments on Pou4 expression in the nerve cord revealed that, in developing amphioxus, RA signaling plays a dual role by (1) providing anteroposterior patterning information to neural cells and (2) specifying neural cell types. This finding is coherent with a terminal selector function of Pou4 for GABAergic neurons in amphioxus and represents the first description of RA-induced changes in POU gene expression in vivo.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications

BACKGROUND

Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates.

RESULTS

We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements.

CONCLUSIONS

In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.