Connexins evolved after early chordates lost innexin diversity

Gap junction channels are formed by two unrelated protein families. Non-chordates use the primordial innexins, while chordates use connexins that superseded the gap junction function of innexins. Chordates retained innexin-homologs, but N-glycosylation prevents them from forming gap junctions. It is puzzling why chordates seem to exclusively use the new gap junction protein and why no chordates should exist that use non-glycosylated innexins to form gap junctions. Here, we identified glycosylation sites of 2270 innexins from 152 non-chordate and 274 chordate species. Among all chordates, we found not a single innexin without glycosylation sites. Surprisingly, the glycosylation motif is also widespread among non-chordate innexins indicating that glycosylated innexins are not a novelty of chordates. In addition, we discovered a loss of innexin diversity during the early chordate evolution. Most importantly, the most basal living chordates, which lack connexins, exclusively possess innexins with glycosylation sites. A bottleneck effect might thus explain why connexins have become the only protein used to form chordate gap junctions.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

Cardiopharyngeal deconstruction and ancestral tunicate sessility

A key problem in understanding chordate evolution has been the origin of sessility of ascidians, and whether the appendicularian free-living style represents a primitive or derived condition of tunicates. To address this problem, we performed comprehensive developmental and genomic comparative analyses of the cardiopharyngeal gene regulatory network (GRN) between appendicularians and ascidians. Our results reveal that the cardiopharyngeal GRN has suffered a process of evolutionary deconstruction with massive ancestral losses of genes (Mesp, Ets1/2, Gata4/5/6, Mek1/2, Tbx1/10, and RA- and FGF-signaling related genes) and subfunctions (e.g. FoxF, Islet, Ebf, Mrf, Dach and Bmp signaling). These losses have led to the deconstruction of two modules of the cardiopharyngeal GRN that in ascidians are related to early and late multipotent state cells involved in lineage fate determination towards first and secondary heart fields, and siphon muscle. Our results allow us to propose an evolutionary scenario, in which the evolutionary deconstruction of the cardiopharyngeal GRN has had an adaptive impact on the acceleration of the developmental cardiac program, the redesign of the cardiac architecture into an open-wide laminar structure, and the loss of pharyngeal muscle. Our findings, therefore, provide evidence supporting that the ancestral tunicate had a sessile ascidian-like lifestyle, and points to the deconstruction of the cardiopharyngeal GRN in appendicularians as a key event that facilitated the evolution of their pelagic free-living style connected to the innovation of the house.

Functional requirement of terminal inverted repeats for efficient ProtoRAG activity reveals the early evolution of V(D)J recombination

The discovery of ProtoRAG in amphioxus indicated that vertebrate RAG recombinases originated from an ancient transposon. However, the sequences of ProtoRAG terminal inverted repeats (TIRs) were obviously dissimilar to the consensus sequence of mouse 12/23RSS and recombination mediated by ProtoRAG or RAG made them incompatible with each other. Thus, it is difficult to determine whether or how 12/23RSS persisted in the vertebrate RAG system that evolved from the TIRs of ancient RAG transposons. Here, we found that the activity of ProtoRAG is highly dependent on its asymmetric 5′TIR and 3′TIR, which are composed of conserved TR1 and TR5 elements and a partially conserved TRsp element of 27/31 bp to separate them. Similar to the requirements for the recombination signal sequences (RSSs) of RAG recombinase, the first CAC in TR1, the three dinucleotides in TR5 and the specific length of the partially conserved TRsp are important for the efficient recombination activity of ProtoRAG. In addition, the homologous sequences flanking the signal sequences facilitate ProtoRAG- but not RAG-mediated recombination. In addition to the diverged TIRs, two differentiated functional domains in BbRAG1L were defined to coordinate with the divergence between TIRs and RSSs. One of these is the CTT* domain, which facilitates the specific TIR recognition of the BbRAGL complex, and the other is NBD*, which is responsible for DNA binding and the protein stabilization of the BbRAGL complex. Thus, our findings reveal that the functional requirement for ProtoRAG TIRs is similar to that for RSS in RAG-mediated recombination, which not only supports the common origin of ProtoRAG TIRs and RSSs from the asymmetric TIRs of ancient RAG transposons, but also reveals the development of RAG and RAG-like machineries during chordate evolution.

The evolution of the 9aaTAD domain in Sp2 proteins: inactivation with valines and intron reservoirs

The Sp1 transcription factor has been defined as glutamine-rich activator. The Nine amino acid TransActivation Domains (9aaTAD) have been identified in numerous transcription activators. Here, we identified the conserved 9aaTAD motif in the Sp1 and in all nine members of SP family with broad natural 9aaTAD variations. We showed by the amino acid substitutions that the glutamine residues are completely dispensable for 9aaTADs function. We described the 9aaTAD domains’ origin and evolutionary history. The ancestral Sp2 gene with inactive 9aaTAD has duplicated in early chordates and created new paralogs Sp1, Sp3 and Sp4. We discovered that the accumulation of valines in the 9aaTADs correlated with the domain inactivation. The Sp2 activation domain, whose dormancy have lasted over 100 million years during chordate evolution, enabled later diversification in the Sp1-4 clade, including both repressors and activators. The new paralogs Sp1 and Sp3 activation domains have regained their original activator function by loss of valines in their 9aaTADs.