Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider.

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.

Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication.

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid Evolution.

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes.

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

A Functional Aqp1 Gene Product Localizes on The Contractile Vacuole Complex in Paramecium multimicronucleatum.

In a ciliate Paramecium, the presence of water channels on the membrane of contractile vacuole has long been predicted by both morphological and physiological data, however, to date either the biochemical or the molecular biological data have not been provided. In the present study, to examine the presence of aquaporin in Paramecium, we carried out RT-PCR with degenerated primers designed based on the ParameciumDB, and an aquaporin cDNA (aquaporin 1, aqp1) with a full-length ORF encoding 251 amino acids was obtained from Paramecium multimicronucleatum by using RACE. The deduced amino acid sequence of AQP1 had NPA-NPG motifs, and the prediction of protein secondary structure by CNR5000 and hydropathy plot showed the presence of six putative transmembrane domains and five connecting loops. Phylogenetic analysis results showed that the amino acid sequence of AQP1 was close to that of the Super-aquaporin group. The AQP1-GFP fusion protein clearly demonstrated the subcellular localization of AQP1 on the contractile vacuole complex, except for the decorated spongiome membrane. The functional analyses of aqp1 were done by RNA interference-based gene silencing, using an established feeding method. The aqp1 was found to be crucial for the total fluid output of the cell, the function of contractile vacuole membranes.

Evolution of vertebrate nicotinic acetylcholine receptors.

BACKGROUND: Many physiological processes are influenced by nicotinic acetylcholine receptors (nAChR), ranging from neuromuscular and parasympathetic signaling to modulation of the reward system and long-term memory. Due to the complexity of the nAChR family and variable evolutionary rates among its members, their evolution in vertebrates has been difficult to resolve. In order to understand how and when the nAChR genes arose, we have used a broad approach of analyses combining sequence-based phylogeny, chromosomal synteny and intron positions. RESULTS: Our analyses suggest that there were ten subunit genes present in the vertebrate predecessor. The two basal vertebrate tetraploidizations (1R and 2R) then expanded this set to 19 genes. Three of these have been lost in mammals, resulting in 16 members today. None of the ten ancestral genes have kept all four copies after 2R. Following 2R, two of the ancestral genes became triplicates, five of them became pairs, and three seem to have remained single genes. One triplet consists of CHRNA7, CHRNA8 and the previously undescribed CHRNA11, of which the two latter have been lost in mammals but are still present in lizards and ray-finned fishes. The other triplet consists of CHRNB2, CHRNB4 and CHRNB5, the latter of which has also been lost in mammals. In ray-finned fish the neuromuscular subunit gene CHRNB1 underwent a local gene duplication generating CHRNB1.2. The third tetraploidization in the predecessor of teleosts (3R) expanded the repertoire to a total of 31 genes, of which 27 remain in zebrafish. These evolutionary relationships are supported by the exon-intron organization of the genes. CONCLUSION: The tetraploidizations explain all gene duplication events in vertebrates except two. This indicates that the genome doublings have had a substantial impact on the complexity of this gene family leading to a very large number of members that have existed for hundreds of millions of years.

Inference of Causative Genes for Alzheimer's Disease Due to Dosage Imbalance.

Copy number variations (CNVs) have recently drawn attention as an important genetic factor for diseases, especially common neuropsychiatric disorders including Alzheimer's disease (AD). Because most of the pathogenic CNV regions overlap with multiple genes, it has been challenging to identify the true disease-causing genes amongst them. Notably, a recent study reported that CNV regions containing ohnologs, which are dosage-sensitive genes, are likely to be deleterious. Utilizing the unique feature of ohnologs could be useful for identifying causative genes with pathogenic CNVs, however its effectiveness is still unclear. Although it has been reported that AD is strongly affected by CNVs, most of AD-causing genes with pathogenic CNVs have not been identified yet. Here, we show that dosage-sensitive ohnologs within CNV regions reported in patients with AD are related to the nervous system and are highly expressed in the brain, similar to other known susceptible genes for AD. We found that CNV regions in patients with AD contained dosage-sensitive genes, which are ohnologs not overlapping with control CNV regions, frequently. Furthermore, these dosage-sensitive genes in pathogenic CNV regions had a strong enrichment in the nervous system for mouse knockout phenotype and high expression in the brain similar to the known susceptible genes for AD. Our results demonstrated that selecting dosage-sensitive ohnologs out of multiple genes with pathogenic CNVs is effective in identifying the causative genes for AD. This methodology can be applied to other diseases caused by dosage imbalance and might help to establish the medical diagnosis by analysis of CNVs.

Gene Function is a Driver of Activin Signaling Pathway Evolution Following Whole-Genome Duplication in Rainbow Trout (Oncorhynchus mykiss).

The genomes of plant and animal species are influenced by ancestral whole-genome duplication (WGD) events, which have profound impacts on the regulation and function of gene networks. To gain insight into the consequences of WGD events, we characterized the sequence conservation and expression patterns of ohnologs in the highly duplicated activin receptor signaling pathway in rainbow trout (RBT). The RBT activin receptor signaling pathway is defined by tissue-specific expression of inhibitors and ligands and broad expression of receptors and Co-Smad signaling molecules. Signaling pathway ligands exhibited shared expression, while inhibitors and Smad signaling molecules primarily express a single dominant ohnolog. Our findings suggest that gene function influences ohnolog evolution following duplication of the activin signaling pathway in RBT.

A homology guide for Pacific salmon genus Oncorhynchus resolves patterns of ohnolog retention, resolution and local adaptation following the salmonid-specific whole-genome duplication event.

Salmonid fishes have emerged as a tractable model to study whole-genome duplications (WGDs) as this group has undergone four rounds of WGDs. While most of the salmonid genome has returned to a diploid state, a significant proportion of genes are maintained as duplicates and are referred to as ohnologs. The fact that much of the modern salmonid gene repertoire is comprised of ohnologs, while other genes have returned to their singleton state creates complications for genetic studies by obscuring homology relationships. The difficulty this creates is particularly prominent in Pacific salmonids belonging to genus Oncorhynchus who are the focus of intense genetics-based conservation and management efforts owing to the important ecological and cultural roles these fish play. To address this gap, we generated a homology guide for six species of Oncorhynchus with available genomes and used this guide to describe patterns of ohnolog retention and resolution. Overall, we find that ohnologs comprise approximately half of each species modern gene repertoires, which are functionally enriched for genes involved in DNA binding, while the less numerous singleton genes are heavily enriched in dosage-sensitive processes such as mitochondrial metabolism. Additionally, by reanalyzing published expression data from locally adapted strains of O. mykiss, we show that numerous ohnologs exhibit adaptive expression profiles; however, ohnologs are not more likely to display adaptive signatures than either paralogs or singletons. Finally, we demonstrate the utility of our homology guide by investigating the evolutionary relationship among genes highlighted as playing a role in salmonid life-history traits or gene editing targets.

Gene Duplication, Shifting Selection, and Dosage Balance of Silicon Transporter Proteins in Marine and Freshwater Diatoms.

Numerous factors shape the evolution of protein-coding genes, including shifts in the strength or type of selection following gene duplications or changes in the environment. Diatoms and other silicifying organisms use a family of silicon transporters (SITs) to import dissolved silicon from the environment. Freshwaters contain higher silicon levels than oceans, and marine diatoms have more efficient uptake kinetics and less silicon in their cell walls, making them better competitors for a scarce resource. We compiled SITs from 37 diatom genomes to characterize shifts in selection following gene duplications and marine-freshwater transitions. A deep gene duplication, which coincided with a whole-genome duplication, gave rise to two gene lineages. One of them (SIT1-2) is present in multiple copies in most species and is known to actively import silicon. These SITs have evolved under strong purifying selection that was relaxed in freshwater taxa. Episodic diversifying selection was detected but not associated with gene duplications or habitat shifts. In contrast, genes in the second SIT lineage (SIT3) were present in just half the species, the result of multiple losses. Despite conservation of SIT3 in some lineages for the past 90-100 million years, repeated losses, relaxed selection, and low expression highlighted the dispensability of SIT3, consistent with a model of deterioration and eventual loss due to relaxed selection on SIT3 expression. The extensive but relatively balanced history of duplications and losses, together with paralog-specific expression patterns, suggest diatoms continuously balance gene dosage and expression dynamics to optimize silicon transport across major environmental gradients.

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.

The divergence of alternative splicing between ohnologs in teleost fishes.

BACKGROUND: Gene duplication and alternative splicing (AS) are two distinct mechanisms generating new materials for genetic innovations. The evolutionary link between gene duplication and AS is still controversial, due to utilizing duplicates from inconsistent ages of duplication events in earlier studies. With the aid of RNA-seq data, we explored evolutionary scenario of AS divergence between duplicates with ohnologs that resulted from the teleost genome duplication event in zebrafish, medaka, and stickleback. RESULTS: Ohnologs in zebrafish have fewer AS forms compared to their singleton orthologs, supporting the function-sharing model of AS divergence between duplicates. Ohnologs in stickleback have more AS forms compared to their singleton orthologs, which supports the accelerated model of AS divergence between duplicates. The evolution of AS in ohnologs in medaka supports a combined scenario of the function-sharing and the accelerated model of AS divergence between duplicates. We also found a small number of ohnolog pairs in each of the three teleosts showed significantly asymmetric AS divergence. For example, the well-known ovary-factor gene cyp19a1a has no AS form but its ohnolog cyp19a1b has multiple AS forms in medaka, suggesting that functional divergence between duplicates might have result from AS divergence. CONCLUSIONS: We found that a combined scenario of function-sharing and accelerated models for AS evolution in ohnologs in teleosts and rule out the independent model that assumes a lack of correlation between gene duplication and AS. Our study thus provided insights into the link between gene duplication and AS in general and ohnolog divergence in teleosts from AS perspective in particular.

Independent pseudogenizations and losses of sox15 during amniote diversification following asymmetric ohnolog evolution.

BACKGROUND: Four ohnologous genes (sox1, sox2, sox3, and sox15) were generated by two rounds of whole-genome duplication in a vertebrate ancestor. In eutherian mammals, Sox1, Sox2, and Sox3 participate in central nervous system (CNS) development. Sox15 has a function in skeletal muscle regeneration and has little functional overlap with the other three ohnologs. In contrast, the frog Xenopus laevis and zebrafish orthologs of sox15 as well as sox1-3 function in CNS development. We previously reported that Sox15 is involved in mouse placental development as neofunctionalization, but is pseudogenized in the marsupial opossum. These findings suggest that sox15 might have evolved with divergent gene fates during vertebrate evolution. However, knowledge concerning sox15 in other vertebrate lineages than therian mammals, anuran amphibians, and teleost fish is scarce. Our purpose in this study was to clarify the fate and molecular evolution of sox15 during vertebrate evolution. RESULTS: We searched for sox15 orthologs in all vertebrate classes from agnathans to mammals by significant sequence similarity and synteny analyses using vertebrate genome databases. Interestingly, sox15 was independently pseudogenized at least twice during diversification of the marsupial mammals. Moreover, we observed independent gene loss of sox15 at least twice during reptile evolution in squamates and crocodile-bird diversification. Codon-based phylogenetic tree and selective analyses revealed an increased dN/dS ratio for sox15 compared to the other three ohnologs during jawed vertebrate evolution. CONCLUSIONS: The findings revealed an asymmetric evolution of sox15 among the four ohnologs during vertebrate evolution, which was supported by the increased dN/dS values in cartilaginous fishes, anuran amphibians, and amniotes. The increased dN/dS value of sox15 may have been caused mainly by relaxed selection. Notably, independent pseudogenizations and losses of sox15 were observed during marsupial and reptile evolution, respectively. Both might have been caused by strong relaxed selection. The drastic gene fates of sox15, including neofunctionalization and pseudogenizations/losses during amniote diversification, might be caused by a release from evolutionary constraints.

Addressing incomplete lineage sorting and paralogy in the inference of uncertain salmonid phylogenetic relationships.

Recent and continued progress in the scale and sophistication of phylogenetic research has yielded substantial advances in knowledge of the tree of life; however, segments of that tree remain unresolved and continue to produce contradicting or unstable results. These poorly resolved relationships may be the product of methodological shortcomings or of an evolutionary history that did not generate the signal traits needed for its eventual reconstruction. Relationships within the euteleost fish family Salmonidae have proven challenging to resolve in molecular phylogenetics studies in part due to ancestral autopolyploidy contributing to conflicting gene trees. We examine a sequence capture dataset from salmonids and use alternative strategies to accommodate the effects of gene tree conflict based on aspects of salmonid genome history and the multispecies coalescent. We investigate in detail three uncertain relationships: (1) subfamily branching, (2) monophyly of Coregonus and (3) placement of Parahucho. Coregoninae and Thymallinae are resolved as sister taxa, although conflicting topologies are found across analytical strategies. We find inconsistent and generally low support for the monophyly of Coregonus, including in results of analyses with the most extensive dataset and complex model. The most consistent placement of Parahucho is as sister lineage of Salmo.

Sensitivity to gene dosage and gene expression affects genes with copy number variants observed among neuropsychiatric diseases.

BACKGROUND: Copy number variants (CNVs) have been reported to be associated with diseases, traits, and evolution. However, it is hard to determine which gene should have priority as a target for further functional experiments if a CNV is rare or a singleton. In this study, we attempted to overcome this issue by using two approaches: by assessing the influences of gene dosage sensitivity and gene expression sensitivity. Dosage sensitive genes derived from two-round whole-genome duplication in previous studies. In addition, we proposed a cross-sectional omics approach that utilizes open data from GTEx to assess the effect of whole-genome CNVs on gene expression. METHODS: Affymetrix Genome-Wide SNP Array 6.0 was used to detect CNVs by PennCNV and CNV Workshop. After quality controls for population stratification, family relationship and CNV detection, 287 patients with narcolepsy, 133 patients with essential hypersomnia, 380 patients with panic disorders, 164 patients with autism, 784 patients with Alzheimer disease and 1280 healthy individuals remained for the enrichment analysis. RESULTS: Overall, significant enrichment of dosage sensitive genes was found across patients with narcolepsy, panic disorders and autism. Particularly, significant enrichment of dosage-sensitive genes in duplications was observed across all diseases except for Alzheimer disease. For deletions, less or no enrichment of dosage-sensitive genes with deletions was seen in the patients when compared to the healthy individuals. Interestingly, significant enrichments of genes with expression sensitivity in brain were observed in patients with panic disorder and autism. While duplications presented a higher burden, deletions did not cause significant differences when compared to the healthy individuals. When we assess the effect of sensitivity to genome dosage and gene expression at the same time, the highest ratio of enrichment was observed in the group including dosage-sensitive genes and genes with expression sensitivity only in brain. In addition, shared CNV regions among the five neuropsychiatric diseases were also investigated. CONCLUSIONS: This study contributed the evidence that dosage-sensitive genes are associated with CNVs among neuropsychiatric diseases. In addition, we utilized open data from GTEx to assess the effect of whole-genome CNVs on gene expression. We also investigated shared CNV region among neuropsychiatric diseases.

Evolutionary Origin and Nomenclature of Vertebrate Wnt11-Family Genes.

To adequately connect zebrafish medical models to human biology, it is essential that gene nomenclature reflects gene orthology. Analysis of gene phylogenies and conserved syntenies shows that the zebrafish gene currently called wnt11 (ENSDARG00000004256, ZFIN ID: ZDB-GENE-990603-12) is not the ortholog of the human gene called WNT11 (ENSG00000085741); instead, the gene currently called wnt11r (ENSDARG00000014796, ZFIN ID: ZDB-GENE-980526-249) is the zebrafish ortholog of human WNT11. Genomic analysis of Wnt11-family genes suggests a model for the birth of Wnt11-family gene ohnologs in genome duplication events, provides a mechanism for the death of a Wnt11-family ohnolog in mammals after they diverged from birds, and suggests revised nomenclature to better connect teleost disease models to human biology.

Characterization of Human Dosage-Sensitive Transcription Factor Genes.

Copy number changes in protein-coding genes are detrimental if the consequent changes in protein concentrations disrupt essential cellular functions. The dosage sensitivity of transcription factor (TF) genes is particularly interesting because their products are essential in regulating the expression of genetic information. From four recently curated data sets of dosage-sensitive genes (genes with conserved copy numbers across mammals, ohnologs, and two data sets of haploinsufficient genes), we compiled a data set of the most reliable dosage-sensitive (MRDS) genes and a data set of the most reliable dosage-insensitive (MRDIS) genes. The MRDS genes were those present in all four data sets, while the MRDIS genes were those absent from any one of the four data sets and with the probability of being loss of function-intolerant (pLI) values < 0.5 in both of the haploinsufficient gene data sets. Enrichment analysis of TF genes among the MRDS and MRDIS gene data sets showed that TF genes are more likely to be dosage-sensitive than other genes in the human genome. The nuclear receptor family was the most enriched TF family among the dosage-sensitive genes. TF families with very few members were also deemed more likely to be dosage-sensitive than TF families with more members. In addition, we found a certain number of dosage-insensitive TFs. The most typical were the Krüppel-associated box domain-containing zinc-finger proteins (KZFPs). Gene ontology (GO) enrichment analysis showed that the MRDS TFs were enriched for many more terms than the MRDIS TFs; however, the proteins interacting with these two groups of TFs did not show such sharp differences. Furthermore, we found that the MRDIS KZFPs were not significantly enriched for any GO terms, whereas their interacting proteins were significantly enriched for thousands of GO terms. Further characterizations revealed significant differences between MRDS TFs and MRDIS TFs in the lengths and nucleotide compositions of DNA-binding sites as well as in expression level, protein size, and selective force.

Structure, evolution and expression of collagen XXVIII: Lessons from the zebrafish.

Collagen XXVIII is the last discovered member of the collagen superfamily and thus has been only sparsely investigated. We studied collagen XXVIII in zebrafish to gain insight into its structure, evolution and expression. In contrast to human and mouse, the zebrafish genome contains four collagen XXVIII genes, col28a1a and -b, and col28a2a and -b. Genomic context and phylogenetic analysis revealed that the a2 branch was lost during evolution of mammals, whereas the duplication of the a1 and a2 branches results from the whole genome duplication in the teleost lineage. Sequence analysis revealed conservation of domain structure and the unique imperfections in the triple helical domain. Two major forms of collagen XXVIII were identified, Col28a1b in adult and Col28a2a in 3-5dpf zebrafish. Composite agarose/polyacrylamide gel electrophoresis revealed that both these chains mainly form dimers of trimers, although Col28a1b appears to be more polydisperse. Homodimers are abundant, although it is possible that complexes consisting of Col28a2a and Col28a1a or -a2b occur. Peptide mass fingerprint analysis revealed that the C-terminal Kunitz domain is often proteolytically processed. In contrast to murine collagen XXVIII, the zebrafish orthologs are widely expressed and not only present in the nervous system. They are differentially expressed in the liver, thymus, muscle, intestine and skin. Altogether our results point to a unique nature of collagen XXVIII within the collagen family.

Connectivity of vertebrate genomes: Paired-related homeobox (Prrx) genes in spotted gar, basal teleosts, and tetrapods.

Teleost fish are important models for human biology, health, and disease. Because genome duplication in a teleost ancestor (TGD) impacts the evolution of teleost genome structure and gene repertoires, we must discriminate gene functions that are shared and ancestral from those that are lineage-specific in teleosts or tetrapods to accurately apply inferences from teleost disease models to human health. Generalizations must account both for the TGD and for divergent evolution between teleosts and tetrapods after the likely two rounds of genome duplication shared by all vertebrates. Progress in sequencing techniques provides new opportunities to generate genomic and transcriptomic information from a broad range of phylogenetically informative taxa that facilitate detailed understanding of gene family and gene function evolution. We illustrate here the use of new sequence resources from spotted gar (Lepisosteus oculatus), a rayfin fish that diverged from teleosts before the TGD, as well as RNA-Seq data from gar and multiple teleost lineages to reconstruct the evolution of the Paired-related homeobox (Prrx) transcription factor gene family, which is involved in the development of mesoderm and neural crest-derived mesenchyme. We show that for Prrx genes, the spotted gar genome and gene expression patterns mimic mammals better than teleosts do. Analyses force the seemingly paradoxical conclusion that regulatory mechanisms for the limb expression domains of Prrx genes existed before the evolution of paired appendages. Detailed evolutionary analyses like those reported here are required to identify fish species most similar to the human genome to optimally connect fish models to human gene functions in health and disease.