Detection and phylogenetic assessment of conserved synteny derived from whole genome duplications

A Revised Perspective on the Evolution of Troponin I and Troponin T Gene Families in Vertebrates

The troponin (Tn) complex, responsible for the Ca2+ activation of striated muscle, is composed of three interacting protein subunits: TnC, TnI, and TnT, encoded by TNNC, TNNI, and TNNT genes. TNNI and TNNT are sister gene families, and in mammals the three TNNI paralogs (TNNI1, TNNI2, TNNI3), which encode proteins with tissue-specific expression, are each in close genomic proximity with one of the three TNNT paralogs (TNNT2, TNNT3, TNNT1, respectively). It has been widely presumed that all vertebrates broadly possess genes of these same three classes, although earlier work has overlooked jawless fishes (cyclostomes) and cartilaginous fishes (chimeras, rays, and sharks), which are distantly related to other jawed vertebrates. With a new phylogenetic and synteny analysis of a diverse array of vertebrates including these taxonomic groups, we define five distinct TNNI classes (TNNI1-5), with TNNI4 and TNNI5 being only present in non-amniote vertebrates and typically found in tandem, and four classes of TNNT (TNNT1-4). These genes are located in four genomic loci that were generated by the 2R whole-genome duplications. TNNI3, encoding "cardiac TnI" in tetrapods, was independently lost in cartilaginous and ray-finned fishes. Instead, ray-finned fishes predominantly express TNNI1 in the heart. TNNI5 is highly expressed in shark hearts and contains a N-terminal extension similar to that of TNNI3 found in tetrapod hearts. Given that TNNI3 and TNNI5 are distantly related, this supports the hypothesis that the N-terminal extension may be an ancestral feature of vertebrate TNNI and not an innovation unique to TNNI3, as has been commonly believed.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions

Background

Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown.

Results

Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels.

Conclusions

Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08304-2.

Whole-Genome Duplications and the Diversification of the Globin-X Genes of Vertebrates

Globin-X (GbX) is an enigmatic member of the vertebrate globin gene family with a wide phyletic distribution that spans protostomes and deuterostomes. Unlike canonical globins such as hemoglobins and myoglobins, functional data suggest that GbX does not have a primary respiratory function. Instead, evidence suggests that the monomeric, membrane-bound GbX may play a role in cellular signaling or protection against the oxidation of membrane lipids. Recently released genomes from key vertebrates provide an excellent opportunity to address questions about the early stages of the evolution of GbX in vertebrates. We integrate bioinformatics, synteny, and phylogenetic analyses to characterize the diversity of GbX genes in nonteleost ray-finned fishes, resolve relationships between the GbX genes of cartilaginous fish and bony vertebrates, and demonstrate that the GbX genes of cyclostomes and gnathostomes derive from independent duplications. Our study highlights the role that whole-genome duplications (WGDs) have played in expanding the repertoire of genes in vertebrate genomes. Our results indicate that GbX paralogs have a remarkably high rate of retention following WGDs relative to other globin genes and provide an evolutionary framework for interpreting results of experiments that examine functional properties of GbX and patterns of tissue-specific expression. By identifying GbX paralogs that are products of different WGDs, our results can guide the design of experimental work to explore whether gene duplicates that originate via WGDs have evolved novel functional properties or expression profiles relative to singleton or tandemly duplicated copies of GbX.

Early vertebrate origin of CTCFL, a CTCF paralog, revealed by proximity-guided shark genome scaffolding

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in diverse animals. The gene encoding this protein has a paralog which was first identified to be expressed exclusively in the testis in mammals and designated as CTCFL (also called BORIS). CTCFL orthologs were reported only among amniotes, and thus CTCFL was once thought to have arisen in the amniote lineage. In this study, we identified elasmobranch CTCFL orthologs, and investigated its origin with the aid of a shark genome assembly improved by proximity-guided scaffolding. Our analysis employing evolutionary interpretation of syntenic gene location suggested an earlier timing of the gene duplication between CTCF and CTCFL than previously thought, that is, around the common ancestor of extant vertebrates. Also, our transcriptomic sequencing revealed a biased expression of the catshark CTCFL in the testis, suggesting the origin of the tissue-specific localization in mammals more than 400 million years ago. To understand the historical process of the functional consolidation of the long-standing chromatin regulator CTCF, its additional paralogs remaining in some of the descendant lineages for spatially restricted transcript distribution should be taken into consideration.

The EJC component Magoh in non-vertebrate chordates

Earliest craniates possess a newly enlarged, elaborated forebrain with new cell types and neuronal networks. A key question in vertebrate evolution is when and how this cerebral expansion took place. The exon-junction complex (EJC) plays an essential role in mRNA processing of all Eukarya. Recently, it has been proposed that the EJC represses recursive RNA splicing in Deuterostomes, with implication in human brain diseases like microcephaly and depression. However, the EJC or EJC subunit contribution to brain development in non-vertebrate Deuterostomes remained unknown. Being interested in the evolution of chordate characters, we focused on the model species, Branchiostoma lanceolatum (Cephalochordata) and Ciona robusta (Tunicata), with the aim to investigate the ancestral and the derived expression state of Magoh orthologous genes. This study identifies that Magoh is part of a conserved syntenic group exclusively in vertebrates and suggests that Magoh has experienced duplication and loss events in mammals. During early development in amphioxus and ascidian, maternal contribution and zygotic expression of Magoh genes in various types of progenitor cells and tissues are consistent with the condition observed in other Bilateria. Later in development, we also show expression of Magoh in the brain of cephalochordate and ascidian larvae. Collectively, these results provide a basis to further define what functional role(s) Magoh exerted during nervous system development and evolution.

Inference of the ancestral vertebrate phenotype through vestiges of the whole-genome duplications

Inferring the phenotype of the last common ancestor of living vertebrates is a challenging problem because of several unresolvable factors. They include the lack of reliable out-groups of living vertebrates, poor information about less fossilizable organs and specialized traits of phylogenetically important species, such as lampreys and hagfishes (e.g. secondary loss of vertebrae in adult hagfishes). These factors undermine the reliability of ancestral reconstruction by traditional character mapping approaches based on maximum parsimony. In this article, we formulate an approach to hypothesizing ancestral vertebrate phenotypes using information from the phylogenetic and functional properties of genes duplicated by genome expansions in early vertebrate evolution. We named the conjecture as ‘chronological reconstruction of ohnolog functions (CHROF)’. This CHROF conjecture raises the possibility that the last common ancestor of living vertebrates may have had more complex traits than currently thought.

Asymmetric paralog evolution between the "cryptic" gene Bmp16 and its well-studied sister genes Bmp2 and Bmp4

The vertebrate gene repertoire is characterized by “cryptic” genes whose identification has been hampered by their absence from the genomes of well-studied species. One example is the Bmp16 gene, a paralog of the developmental key genes Bmp2 and -4. We focus on the Bmp2/4/16 group of genes to study the evolutionary dynamics following gen(om)e duplications with special emphasis on the poorly studied Bmp16 gene. We reveal the presence of Bmp16 in chondrichthyans in addition to previously reported teleost fishes and reptiles. Using comprehensive, vertebrate-wide gene sampling, our phylogenetic analysis complemented with synteny analyses suggests that Bmp2, -4 and -16 are remnants of a gene quartet that originated during the two rounds of whole-genome duplication (2R-WGD) early in vertebrate evolution. We confirm that Bmp16 genes were lost independently in at least three lineages (mammals, archelosaurs and amphibians) and report that they have elevated rates of sequence evolution. This finding agrees with their more “flexible” deployment during development; while Bmp16 has limited embryonic expression domains in the cloudy catshark, it is broadly expressed in the green anole lizard. Our study illustrates the dynamics of gene family evolution by integrating insights from sequence diversification, gene repertoire changes, and shuffling of expression domains.

Duplication and diversification of insulin genes in ray-finned fish

Insulin is a key hormone for the regulation of metabolism in vertebrates. Insulin is produced by pancreatic islet cells in response to elevated glucose levels and leads to the uptake of glucose by tissues such as liver and adipose tissue to store energy. Insulin also has additional functions in regulating development. Previous work has shown that the proglucagon gene, which encodes hormones counter regulating insulin, is duplicated in teleost fish, and that the peptide hormones encoded by these genes have diversified in function. I sought to determine whether similar processes have occurred to insulin genes in these species. Searches of fish genomes revealed an unexpected diversity of insulin genes. A triplication of the insulin gene occurred at the origin of teleost fish, however one of these three genes, insc, has been lost in most teleost fish lineages. The two other insulin genes, insa and insb, have been retained but show differing levels of selective constraint suggesting that they might have diversified in function. Intriguingly, a duplicate copy of the insa gene, which I named insab, is found in many fish. The coding sequence encoded by insab genes is under weak selective constraint, with its predicted protein sequences losing their potential to be processed into a two-peptide hormone. However, these sequences have retained perfectly conserved cystine residues, suggesting that they maintain insulin's three-dimensional structure and therefore might modulate the processing and secretion of insulin produced by the other genes.

Diversification of the functions of proglucagon and glucagon receptor genes in fish

The teleost fish-specific genome duplication gave rise to a great number of species inhabiting diverse environments with different access to nutrients and life histories. This event produced duplicated gcg genes, gcga and gcgb, for proglucagon-derived peptides, glucagon and GLP-1 and duplicated gcgr receptor genes, gcgra and gcgrb, which play key roles connecting the consumption of nutrients with glucose metabolism. We conducted a systematic survey of the genomes from 28 species of fish (24 bony (Superclass Osteichthyes), 1 lobe-finned (Class Sarcoperygii), 1 cartilaginous (Superclass Chondrichthyes), and 2 jawless (Superclass Agnatha)) and find that almost all surveyed ray-finned fish contain gcga and gcgb genes with different coding potential and duplicated gcgr genes, gcgra and gcgrb that form two separate clades in the phylogenetic tree consistent with the accepted species phylogeny. All gcgb genes encoded only glucagon and GLP-1 and gcga genes encoded glucagon, GLP-1, and GLP-2, indicating that gcga was subfunctionalized to produce GLP-2. We find a single glp2r, but no glp1r suggesting that duplicated gcgrb was neofunctionalized to bind GLP-1, as demonstrated for the zebrafish gcgrb (Oren et al., 2016). In functional experiments with zebrafish gcgrb and GLP-1 from diverse fish we find that anglerfish GLP-1a, encoded by gcga, is less biologically active than the gcgb anglerfish GLP-1b paralog. But some other fish (zebrafish, salmon, and catfish) gcga GLP-1a display similar biological activities, indicating that the regulation of glucose metabolism by GLP-1 in ray-finned fish is species-specific. Searches of genomes in cartilaginous fish identified a proglucagon gene that encodes a novel GLP-3 peptide in addition to glucagon, GLP-1, and GLP-2, as well as a single gcgr, glp2r, and a new glucagon receptor-like receptor whose identity still needs to be confirmed. The sequence of the shark GLP-1 contained an N-terminal mammalian-like extension that in mammals undergoes a proteolytic cleavage to release biologically active GLP-1. Our results indicate that early in vertebrate evolution diverse regulatory mechanisms emerged for the control of glucose metabolism by proglucagon-derived peptides and their receptors and that in ray-finned fish they included subfunctionalization and neofunctionalization of these genes.

Evolution of the vertebrate insulin receptor substrate (Irs) gene family

BACKGROUND

Insulin receptor substrate (Irs) proteins are essential for insulin signaling as they allow downstream effectors to dock with, and be activated by, the insulin receptor. A family of four Irs proteins have been identified in mice, however the gene for one of these, IRS3, has been pseudogenized in humans. While it is known that the Irs gene family originated in vertebrates, it is not known when it originated and which members are most closely related to each other. A better understanding of the evolution of Irs genes and proteins should provide insight into the regulation of metabolism by insulin.

RESULTS

Multiple genes for Irs proteins were identified in a wide variety of vertebrate species. Phylogenetic and genomic neighborhood analyses indicate that this gene family originated very early in vertebrae evolution. Most Irs genes were duplicated and retained in fish after the fish-specific genome duplication. Irs genes have been lost of various lineages, including Irs3 in primates and birds and Irs1 in most fish. Irs3 and Irs4 experienced an episode of more rapid protein sequence evolution on the ancestral mammalian lineage. Comparisons of the conservation of the proteins sequences among Irs paralogs show that domains involved in binding to the plasma membrane and insulin receptors are most strongly conserved, while divergence has occurred in sequences involved in interacting with downstream effector proteins.

CONCLUSIONS

The Irs gene family originated very early in vertebrate evolution, likely through genome duplications, and in parallel with duplications of other components of the insulin signaling pathway, including insulin and the insulin receptor. While the N-terminal sequences of these proteins are conserved among the paralogs, changes in the C-terminal sequences likely allowed changes in biological function.

Towards the bridging of molecular genetics data across Xenopus species

BACKGROUND

The clawed African frog Xenopus laevis has been one of the main vertebrate models for studies in developmental biology. However, for genetic studies, Xenopus tropicalis has been the experimental model of choice because it shorter life cycle and due to a more tractable genome that does not result from genome duplication as in the case of X. laevis. Today, although still organized in a large number of scaffolds, nearly 85% of X. tropicalis and 89% of X. laevis genomes have been sequenced. There is expectation for a comparative physical map that can be used as a Rosetta Stone between X. laevis genetic studies and X. tropicalis genomic research.

RESULTS

In this work, we have mapped using coarse-grained alignment the 18 chromosomes of X. laevis, release 9.1, on the 10 reference scaffolds representing the haploid genome of X. tropicalis, release 9.0. After validating the mapping with theoretical data, and estimating reference averages of genome sequence identity, 37 to 44% between the two species, we have carried out a synteny analysis for 2,112 orthologous genes. We found that 99.6% of genes are in the same organization.

CONCLUSIONS

Taken together, our results make possible to establish the correspondence between 62 and 65.5% of both genomes, percentage of identity, synteny and automatic annotation of transcripts of both species, providing a new and more comprehensive tool for comparative analysis of these two species, by allowing to bridge molecular genetics data among them.

Incorporating tree-thinking and evolutionary time scale into developmental biology

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

CRISPR multitargeter: a web tool to find common and unique CRISPR single guide RNA targets in a set of similar sequences

Genome engineering has been revolutionized by the discovery of clustered regularly interspaced palindromic repeats (CRISPR) and CRISPR-associated system genes (Cas) in bacteria. The type IIB Streptococcus pyogenes CRISPR/Cas9 system functions in many species and additional types of CRISPR/Cas systems are under development. In the type II system, expression of CRISPR single guide RNA (sgRNA) targeting a defined sequence and Cas9 generates a sequence-specific nuclease inducing small deletions or insertions. Moreover, knock-in of large DNA inserts has been shown at the sites targeted by sgRNAs and Cas9. Several tools are available for designing sgRNAs that target unique locations in the genome. However, the ability to find sgRNA targets common to several similar sequences or, by contrast, unique to each of these sequences, would also be advantageous. To provide such a tool for several types of CRISPR/Cas system and many species, we developed the CRISPR MultiTargeter software. Similar DNA sequences in question are duplicated genes and sets of exons of different transcripts of a gene. Thus, we implemented a basic sgRNA target search of input sequences for single-sgRNA and two-sgRNA/Cas9 nickase targeting, as well as common and unique sgRNA target searches in 1) a set of input sequences; 2) a set of similar genes or transcripts; or 3) transcripts a single gene. We demonstrate potential uses of the program by identifying unique isoform-specific sgRNA sites in 71% of zebrafish alternative transcripts and common sgRNA target sites in approximately 40% of zebrafish duplicated gene pairs. The design of unique targets in alternative exons is helpful because it will facilitate functional genomic studies of transcript isoforms. Similarly, its application to duplicated genes may simplify multi-gene mutational targeting experiments. Overall, this program provides a unique interface that will enhance use of CRISPR/Cas technology.

Evolution of the vertebrate Pax4/6 class of genes with focus on its novel member, the Pax10 gene

The members of the paired box (Pax) family regulate key developmental pathways in many metazoans as tissue-specific transcription factors. Vertebrate genomes typically possess nine Pax genes (Pax1-9), which are derived from four proto-Pax genes in the vertebrate ancestor that were later expanded through the so-called two-round (2R) whole-genome duplication. A recent study proposed that pax6a genes of a subset of teleost fishes (namely, acanthopterygians) are remnants of a paralog generated in the 2R genome duplication, to be renamed pax6.3, and reported one more group of vertebrate Pax genes (Pax6.2), most closely related to the Pax4/6 class. We propose to designate this new member Pax10 instead and reconstruct the evolutionary history of the Pax4/6/10 class with solid phylogenetic evidence. Our synteny analysis showed that Pax4, -6, and -10 originated in the 2R genome duplications early in vertebrate evolution. The phylogenetic analyses of relationships between teleost pax6a and other Pax4, -6, and -10 genes, however, do not support the proposed hypothesis of an ancient origin of the acanthopterygian pax6a genes in the 2R genome duplication. Instead, we confirmed the traditional scenario that the acanthopterygian pax6a is derived from the more recent teleost-specific genome duplication. Notably, Pax6 is present in all vertebrates surveyed to date, whereas Pax4 and -10 were lost multiple times in independent vertebrate lineages, likely because of their restricted expression patterns: Among Pax6-positive domains, Pax10 has retained expression in the adult retina alone, which we documented through in situ hybridization and quantitative reverse transcription polymerase chain reaction experiments on zebrafish, Xenopus, and anole lizard.

The Dlx genes as clues to vertebrate genomics and craniofacial evolution

The group of Dlx genes belongs to the homeobox-containing superfamily, and its members are involved in various morphogenetic processes. In vertebrate genomes, Dlx genes exist as multiple paralogues generated by tandem duplication followed by whole genome duplications. In this review, we provide an overview of the Dlx gene phylogeny with an emphasis on the chordate lineage. Referring to the Dlx gene repertoire, we discuss the establishment and conservation of the nested expression patterns of the Dlx genes in craniofacial development. Despite the accumulating genomic sequence resources in diverse vertebrates, embryological analyses of Dlx gene expression and function remain limited in terms of species diversity. By supplementing our original analysis of shark embryos with previous data from other osteichthyans, such as mice and zebrafish, we support the previous speculation that the nested Dlx expression in the pharyngeal arch is likely a shared feature among all the extant jawed vertebrates. Here, we highlight several hitherto unaddressed issues regarding the evolution and function of Dlx genes, with special reference to the craniofacial development of vertebrates.

Saltatory evolution of the ectodermal neural cortex gene family at the vertebrate origin

The ectodermal neural cortex (ENC) gene family, whose members are implicated in neurogenesis, is part of the kelch repeat superfamily. To date, ENC genes have been identified only in osteichthyans, although other kelch repeat-containing genes are prevalent throughout bilaterians. The lack of elaborate molecular phylogenetic analysis with exhaustive taxon sampling has obscured the possible link of the establishment of this gene family with vertebrate novelties. In this study, we identified ENC homologs in diverse vertebrates by means of database mining and polymerase chain reaction screens. Our analysis revealed that the ENC3 ortholog was lost in the basal eutherian lineage through single-gene deletion and that the triplication between ENC1, -2, and -3 occurred early in vertebrate evolution. Including our original data on the catshark and the zebrafish, our comparison revealed high conservation of the pleiotropic expression pattern of ENC1 and shuffling of expression domains between ENC1, -2, and -3. Compared with many other gene families including developmental key regulators, the ENC gene family is unique in that conventional molecular phylogenetic inference could identify no obvious invertebrate ortholog. This suggests a composite nature of the vertebrate-specific gene repertoire, consisting not only of de novo genes introduced at the vertebrate origin but also of long-standing genes with no apparent invertebrate orthologs. Some of the latter, including the ENC gene family, may be too rapidly evolving to provide sufficient phylogenetic signals marking orthology to their invertebrate counterparts. Such gene families that experienced saltatory evolution likely remain to be explored and might also have contributed to phenotypic evolution of vertebrates.