The zebrafish genome in context: ohnologs gone missing

Some zebrafish genes appear to lack an ortholog in the human genome and researchers often call them "novel" genes. The origin of many so-called "novel" genes becomes apparent when considered in the context of genome duplication events that occurred during evolution of the phylum Chordata, including two rounds at about the origin of the subphylum Vertebrata (R1 and R2) and one round before the teleost radiation (R3). Ohnologs are paralogs stemming from such genome duplication events, and some zebrafish genes said to be "novel" are more appropriately interpreted as "ohnologs gone missing", cases in which ohnologs are preserved differentially in different evolutionary lineages. Here we consider ohnologs present in the zebrafish genome but absent from the human genome. Reasonable hypotheses are that lineage-specific loss of ohnologs can play a role in establishing lineage divergence and in the origin of developmental innovations. How does the evolution of ohnologs differ from the evolution of gene duplicates arising from other mechanisms, such as tandem duplication or retrotransposition? To what extent do different major vertebrate lineages or different teleost lineages differ in ohnolog content? What roles do differences in ohnolog content play in the origin of developmental mechanisms that differ among lineages? This review explores these questions.     

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.     

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal–gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.     

Multiple Kisspeptin Receptors in Early Osteichthyans Provide New Insights into the Evolution of This Receptor Family

Deorphanization of GPR54 receptor a decade ago led to the characterization of the kisspeptin receptor (Kissr) in mammals and the discovery of its major role in the brain control of reproduction. While a single gene encodes for Kissr in eutherian mammals including human, other vertebrates present a variable number of Kissr genes, from none in birds, one or two in teleosts, to three in an amphibian, xenopus. In order to get more insight into the evolution of Kissr gene family, we investigated the presence of Kissr in osteichthyans of key-phylogenetical positions: the coelacanth, a representative of early sarcopterygians, the spotted gar, a non-teleost actinopterygian, and the European eel, a member of an early group of teleosts (elopomorphs). We report the occurrence of three Kissr for the first time in a teleost, the eel. As measured by quantitative RT-PCR, the three eel Kissr were differentially expressed in the brain-pituitary-gonadal axis, and differentially regulated in experimentally matured eels, as compared to prepubertal controls. Subfunctionalisation, as shown by these differences in tissue distribution and regulation, may have represented significant evolutionary constraints for the conservation of multiple Kissr paralogs in this species. Furthermore, we identified four Kissr in both coelacanth and spotted gar genomes, providing the first evidence for the presence of four Kissr in vertebrates. Phylogenetic and syntenic analyses supported the existence of four Kissr paralogs in osteichthyans and allowed to propose a clarified nomenclature of Kissr (Kissr-1 to -4) based on these paralogs. Syntenic analysis suggested that the four Kissr paralogs arose through the two rounds of whole genome duplication (1R and 2R) in early vertebrates, followed by multiple gene loss events in the actinopterygian and sarcopterygian lineages. Due to gene loss there was no impact of the teleost-specific whole genome duplication (3R) on the number of Kissr paralogs in current teleosts.     

The natural history of class I primate alcohol dehydrogenases includes gene duplication, gene loss, and gene conversion

Background

Gene duplication is a source of molecular innovation throughout evolution. However, even with massive amounts of genome sequence data, correlating gene duplication with speciation and other events in natural history can be difficult. This is especially true in its most interesting cases, where rapid and multiple duplications are likely to reflect adaptation to rapidly changing environments and life styles. This may be so for Class I of alcohol dehydrogenases (ADH1s), where multiple duplications occurred in primate lineages in Old and New World monkeys (OWMs and NWMs) and hominoids.

Methodology/Principal Findings

To build a preferred model for the natural history of ADH1s, we determined the sequences of nine new ADH1 genes, finding for the first time multiple paralogs in various prosimians (lemurs, strepsirhines). Database mining then identified novel ADH1 paralogs in both macaque (an OWM) and marmoset (a NWM). These were used with the previously identified human paralogs to resolve controversies relating to dates of duplication and gene conversion in the ADH1 family. Central to these controversies are differences in the topologies of trees generated from exonic (coding) sequences and intronic sequences.

Conclusions/Significance

We provide evidence that gene conversions are the primary source of difference, using molecular clock dating of duplications and analyses of microinsertions and deletions (micro-indels). The tree topology inferred from intron sequences appear to more correctly represent the natural history of ADH1s, with the ADH1 paralogs in platyrrhines (NWMs) and catarrhines (OWMs and hominoids) having arisen by duplications shortly predating the divergence of OWMs and NWMs. We also conclude that paralogs in lemurs arose independently. Finally, we identify errors in database interpretation as the source of controversies concerning gene conversion. These analyses provide a model for the natural history of ADH1s that posits four ADH1 paralogs in the ancestor of Catarrhine and Platyrrhine primates, followed by the loss of an ADH1 paralog in the human lineage.     

Cis-regulatory characterization of sequence conservation surrounding the Hox4 genes

Hox genes are key regulators of anterior-posterior axis patterning and have a major role in hindbrain development. The zebrafish Hox4 paralogs have strong overlapping activities in hindbrain rhombomeres 7 and 8, in the spinal cord and in the pharyngeal arches. With the aim to predict enhancers that act on the hoxa4a, hoxb4a, hoxc4a and hoxd4a genes, we used sequence conservation around the Hox4 genes to analyze all fish:human conserved non-coding sequences by reporter assays in stable zebrafish transgenesis. Thirty-four elements were functionally tested in GFP reporter gene constructs and more than 100 F1 lines were analyzed to establish a correlation between sequence conservation and cis-regulatory function, constituting a catalog of Hox4 CNEs. Sixteen tissue-specific enhancers could be identified. Multiple alignments of the CNEs revealed paralogous cis-regulatory sequences, however, the CNE sequence similarities were found not to correlate with tissue specificity. To identify ancestral enhancers that direct Hox4 gene activity, genome sequence alignments of mammals, teleosts, horn shark and the cephalochordate amphioxus, which is the most basal extant chordate possessing a single prototypical Hox cluster, were performed. Three elements were identified and two of them exhibited regulatory activity in transgenic zebrafish, however revealing no specificity. Our data show that the approach to identify cis-regulatory sequences by genome sequence alignments and subsequent testing in zebrafish transgenesis can be used to define enhancers within the Hox clusters and that these have significantly diverged in their function during evolution.     

Evolution of developmental roles of Pax2/5/8 paralogs after independent duplication in urochordate and vertebrate lineages

Background

Gene duplication provides opportunities for lineage diversification and evolution of developmental novelties. Duplicated genes generally either disappear by accumulation of mutations (nonfunctionalization), or are preserved either by the origin of positively selected functions in one or both duplicates (neofunctionalization), or by the partitioning of original gene subfunctions between the duplicates (subfunctionalization). The Pax2/5/8 family of important developmental regulators has undergone parallel expansion among chordate groups. After the divergence of urochordate and vertebrate lineages, two rounds of independent gene duplications resulted in the Pax2, Pax5, and Pax8 genes of most vertebrates (the sister group of the urochordates), and an additional duplication provided the pax2a and pax2b duplicates in teleost fish. Separate from the vertebrate genome expansions, a duplication also created two Pax2/5/8 genes in the common ancestor of ascidian and larvacean urochordates.

Results

To better understand mechanisms underlying the evolution of duplicated genes, we investigated, in the larvacean urochordate Oikopleura dioica, the embryonic gene expression patterns of Pax2/5/8 paralogs. We compared the larvacean and ascidian expression patterns to infer modular subfunctions present in the single pre-duplication Pax2/5/8 gene of stem urochordates, and we compared vertebrate and urochordate expression to infer the suite of Pax2/5/8 gene subfunctions in the common ancestor of olfactores (vertebrates + urochordates). Expression pattern differences of larvacean and ascidian Pax2/5/8 orthologs in the endostyle, pharynx and hindgut suggest that some ancestral gene functions have been partitioned differently to the duplicates in the two urochordate lineages. Novel expression in the larvacean heart may have resulted from the neofunctionalization of a Pax2/5/8 gene in the urochordates. Expression of larvacean Pax2/5/8 in the endostyle, in sites of epithelial remodeling, and in sensory tissues evokes like functions of Pax2, Pax5 and Pax8 in vertebrate embryos, and may indicate ancient origins for these functions in the chordate common ancestor.

Conclusion

Comparative analysis of expression patterns of chordate Pax2/5/8 duplicates, rooted on the single-copy Pax2/5/8 gene of amphioxus, whose lineage diverged basally among chordates, provides new insights into the evolution and development of the heart, thyroid, pharynx, stomodeum and placodes in chordates; supports the controversial conclusion that the atrial siphon of ascidians and the otic placode in vertebrates are homologous; and backs the notion that Pax2/5/8 functioned in ancestral chordates to engineer epithelial fusions and perforations, including gill slit openings.     

Genome size reduction in the chicken has involved massive loss of ancestral protein-coding genes

Both mean genomes size and the variance in genome size among species are smaller on average in birds (class Aves) than in the other tetrapod classes. In order to test whether loss of protein-coding genes has contributed to genome size reduction in birds, we compared the chicken genome and five mammalian genomes. Numbers of members (paralogs) were significantly lower in the chicken gene families than in the corresponding mammalian families. Phylogenetic analyses of chicken, mammal, and fish paralogs supported the hypothesis that chicken-specific loss of paralogs occurred much more frequently than mammal-specific gene duplications. Moreover, the phylogenetic analyses supported the hypothesis that a substantial majority of the paralogs lost in chicken originated from duplications prior to the most recent common ancestor of tetrapods and bony fishes. In addition to loss of paralogs, numerous gene families present in the mammalian genomes were missing in the chicken genome; over 1,000 of these families were found in bony fishes, implying presence of the family in the tetrapod ancestor. In the set of families with more members on average in the mammals than in the chicken, immune system function was associated with a greater degree of gene family size reduction in the chicken, consistent with other evidence that immune system gene families have become particularly compact in birds.     

Genetic redundancy in soybean photoresponses associated with duplication of the phytochrome A gene

Gene and genome duplications underlie the origins of evolutionary novelty in plants. Soybean, Glycine max, is considered to be a paleopolyploid species with a complex genome. We found multiple homologs of the phytochrome A gene (phyA) in the soybean genome and determined the DNA sequences of two paralogs designated GmphyA1 and GmphyA2. Analysis of the GmphyA2 gene from the lines carrying a recessive allele at a photoperiod insensitivity locus, E4, revealed that a Ty1/copia-like retrotransposon was inserted in exon 1 of the gene, which resulted in dysfunction of the gene. Mapping studies suggested that GmphyA2 is encoded by E4. The GmphyA1 gene was mapped to a region of linkage group O, which is homeologous to the region harboring E4 in linkage group I. Plants homozygous for the e4 allele were etiolated under continuous far red light, but the de-etiolation occurred partially, indicating that the mutation alone did not cause a complete loss of phyA function. The genetic redundancy suggests that the presence of duplicated copies of phyA genes accounts for the generation of photoperiod insensitivity, while protecting against the deleterious effects of mutation. Thus, this phenomenon provides a link between gene duplication and establishment of an adaptive response of plants to environments.     

Molecular Cloning,Functional Characterization,and Evolutionary Analysis of Vitamin D Receptors Isolated from Basal Vertebrates

The vertebrate genome is a result of two rapid and successive rounds of whole genome duplication, referred to as 1R and 2R. Furthermore, teleost fish have undergone a third whole genome duplication (3R) specific to their lineage, resulting in the retention of multiple gene paralogs. The more recent 3R event in teleosts provides a unique opportunity to gain insight into how genes evolve through specific evolutionary processes. In this study we compare molecular activities of vitamin D receptors (VDR) from basal species that diverged at key points in vertebrate evolution in order to infer derived and ancestral VDR functions of teleost paralogs. Species include the sea lamprey (Petromyzon marinus), a 1R jawless fish; the little skate (Leucoraja erinacea), a cartilaginous fish that diverged after the 2R event; and the Senegal bichir (Polypterus senegalus), a primitive 2R ray-finned fish. Saturation binding assays and gel mobility shift assays demonstrate high affinity ligand binding and classic DNA binding characteristics of VDR has been conserved across vertebrate evolution. Concentration response curves in transient transfection assays reveal EC₅₀ values in the low nanomolar range, however maximum transactivational efficacy varies significantly between receptor orthologs. Protein-protein interactions were investigated using co-transfection, mammalian 2-hybrid assays, and mutations of coregulator activation domains. We then combined these results with our previous study of VDR paralogs from 3R teleosts into a bioinformatics analysis. Our results suggest that 1, 25D₃ acts as a partial agonist in basal species. Furthermore, our bioinformatics analysis suggests that functional differences between VDR orthologs and paralogs are influenced by differential protein interactions with essential coregulator proteins. We speculate that we may be observing a change in the pharmacodynamics relationship between VDR and 1, 25D₃ throughout vertebrate evolution that may have been driven by changes in protein-protein interactions between VDR and essential coregulators.     

Comprehensive Evolutionary and Expression Analysis of FCS-Like Zinc finger Gene Family Yields Insights into Their Origin,Expansion and Divergence

Plant evolution is characterized by frequent genome duplication events. Expansion of habitat resulted in the origin of many novel genes and genome duplication events which in turn resulted in the expansion of many regulatory gene families. The plant-specific FCS-Like Zinc finger (FLZ) gene family is characterized by the presence of a FCS-Like Zinc finger (FLZ) domain which mediates the protein-protein interaction. In this study, we identified that the expansion of FLZ gene family size in different species is correlated with ancestral and lineage-specific whole genome duplication events. The subsequent gene loss found to have a greater role in determining the size of this gene family in many species. However, genomic block duplications played the significant role in the expansion of FLZ gene family in some species. Comparison of Arabidopsis thaliana and Oryza sativa FLZ gene family revealed monocot and dicot specific evolutionary trends. The FLZ genes were found to be under high purifying selection. The spatiotemporal expression analyses of Arabidopsis thaliana FLZ gene family revealed that majority of the members are highly expressed in reproductive organs. FLZ genes were also found to be highly expressed during vegetative-to-reproductive phase transition which is correlated with the proposed role of this gene family in sugar signaling. The comparison of sequence, structural and expression features of duplicated genes identified lineage-specific redundancy and divergence. This extensive evolutionary analysis and expression analysis of Arabidopsis thaliana FLZ genes will pave the way for further functional analysis of FLZ genes.     

Comparative genomics using <Emphasis Type="Italic">Fugu</Emphasis> reveals insights into regulatory subfunctionalization

Background  

A major mechanism for the preservation of gene duplicates in the genome is thought to be mediated via loss or modification of cis-regulatory subfunctions between paralogs following duplication (a process known as regulatory subfunctionalization). Despite a number of gene expression studies that support this mechanism, no comprehensive analysis of regulatory subfunctionalization has been undertaken at the level of the distal cis-regulatory modules involved. We have exploited fish-mammal genomic alignments to identify and compare more than 800 conserved non-coding elements (CNEs) that associate with genes that have undergone fish-specific duplication and retention.     

The Origin and Diversity of Cpt1 Genes in Vertebrate Species

The Carnitine palmitoyltransferase I (Cpt1) gene family plays a crucial role in energy homeostasis since it is required for the occurrence of fatty acid β-oxidation in the mitochondria. The exact gene repertoire in different vertebrate lineages is variable. Presently, four genes are documented: Cpt1a, also known as Cpt1a1, Cpt1a2; Cpt1b and Cpt1c. The later is considered a mammalian innovation resulting from a gene duplication event in the ancestor of mammals, after the divergence of sauropsids. In contrast, Cpt1a2 has been found exclusively in teleosts. Here, we reassess the overall evolutionary relationships of Cpt1 genes using a combination of approaches, including the survey of the gene repertoire in basal gnathostome lineages. Through molecular phylogenetics and synteny studies, we find that Cpt1c is most likely a rapidly evolving orthologue of Cpt1a2. Thus, Cpt1c is present in other lineages such as cartilaginous fish, reptiles, amphibians and the coelacanth. We show that genome duplications (2R) and variable rates of sequence evolution contribute to the history of Cpt1 genes in vertebrates. Finally, we propose that loss of Cpt1b is the likely cause for the unusual energy metabolism of elasmobranch.     

Expression of zebrafish aldh1a3 (raldh3) and absence of aldh1a1 in teleosts

The vitamin A-derived morphogen retinoic acid (RA) plays important roles during the development of chordate animals. The Aldh1a-family of RA-synthesizing enzymes consists of three members, Aldh1a1-3 (Raldh1-3), that are dynamically expressed throughout development. We have searched the known teleost genomes for the presence of Raldh family members and have found that teleost fish possess orthologs of Aldh1a2 and Aldh1a3 only. Here we describe the expression of aldh1a3 in the zebrafish, Danio rerio. Whole mount in situ hybridization shows that aldh1a3 is expressed during eye development in the retina flanking the optic stalks and later is expressed ventrally, opposite the expression domain of aldh1a2. During inner ear morphogenesis, aldh1a3 is expressed in developing sensory epithelia of the cristae and utricular macula and is specifically up-regulated in epithelial projections throughout the formation of the walls of the semicircular canals and endolymphatic duct. In contrast to the mouse inner ear, which expresses all three Raldhs, we find that only aldh1a3 is expressed in the zebrafish otocyst, while aldh1a2 is present in the periotic mesenchyme. During larval stages, additional expression domains of aldh1a3 appear in the anterior pituitary and the swim bladder. Our analyses provide a starting point for genetic studies to examine the role of RA in these organs and emphasize the suitability of the zebrafish inner ear in dissecting the contribution of RA signaling to the development of the vestibular system.     

Loss of protein interactions and regulatory divergence in yeast whole-genome duplicates

Whole-genome duplications are important for the growth of genome complexity. We investigated various factors involved in the evolution of yeast whole-genome duplicates (ohnologs) making emphasis on the analysis of protein interactions. We found that ohnologs have a lower number of protein interactions compared with small-scale duplicates and singletons (by about − 40%). The loss of interactions was proportional to their initial number and independent of ohnolog position in the protein interaction network. A faster evolving member of an ohnolog pair has a lower number of interactions compared to its counterpart. The Gene Ontology mapping of non-overlapping and overlapping interactants of paired ohnologs reveals a sharp asymmetry in GO terms related to regulation. The fraction of these terms is much higher in non-overlapping interactants (compared to overlapping interactants and total dataset). Network clustering coefficient is lower in ohnologs, yet they show an increased density of protein interactions restricted within the whole ohnologs set. These facts suggest that subfunctionalization (or subneofunctionalization) reflected in the loss of protein interactions was a prevailing process in the divergence of ohnologs, which distinguishes them from small-scale duplicates. The loss of protein interactions was associated with the regulatory divergence between the members of an ohnolog pair. A small-scale modularity (reflected in clustering coefficient) probably was not important for ohnologs retention, yet a larger-scale modularity could be involved in their evolution.