Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4× coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element–like and long interspersed element–like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.     

A compact cartilaginous fish model genome

The genomes of several vertebrates, including six mammals, the chicken, Xenopus and four ray-finned fishes have been sequenced or are currently being sequenced to provide a better understanding of the human genome through comparative analysis. However, this list does not include cartilaginous fishes, which are the most basal living jawed vertebrates [1]. The genomes of the current ‘popular’ cartilaginous fishes such as the nurse shark, dogfish, and horn shark are larger than the human genome (∼3800 Mb to 7000 Mb) [2], and are not attractive for whole-genome sequencing. Here, we report the characterization of the relatively small genome (1200 Mb) of a cartilaginous fish, the elephant fish (Callorhinchus milii), and propose it as a model for whole-genome sequencing.     

Ancient vertebrate conserved noncoding elements have been evolving rapidly in teleost fishes

Vertebrate genomes contain thousands of conserved noncoding elements (CNEs) that often function as tissue-specific enhancers. In this study, we have identified CNEs in human, dog, chicken, Xenopus, and four teleost fishes (zebrafish, stickleback, medaka, and fugu) using elephant shark, a cartilaginous vertebrate, as the base genome and investigated the evolution of these ancient vertebrate CNEs (aCNEs) in bony vertebrate lineages. Our analysis shows that aCNEs have been evolving at different rates in different bony vertebrate lineages. Although 78-83% of CNEs have diverged beyond recognition ("lost") in different teleost fishes, only 24% and 40% have been lost in the chicken and mammalian lineages, respectively. Relative rate tests of substitution rates in CNEs revealed that the teleost fish CNEs have been evolving at a significantly higher rate than those in other bony vertebrates. In the ray-finned fish lineage, 68% of aCNEs were lost before the divergence of the four teleosts. This implicates the "fish-specific" whole-genome duplication in the accelerated evolution and the loss of a large number of both copies of duplicated CNEs in teleost fishes. The aCNEs are rich in tissue-specific enhancers and thus many of them are likely to be evolutionarily constrained cis-regulatory elements. The rapid evolution of aCNEs might have affected the expression patterns driven by them. Transgenic zebrafish assay of some human CNE enhancers that have been lost in teleosts has indicated instances of conservation or changes in trans-acting factors between mammals and fishes.     

Divergent Toll-like receptors in catfish (Ictalurus punctatus): TLR5S, TLR20, TLR21

Toll-like receptors (TLR) mediate pathogen recognition in vertebrate species through detection of conserved microbial ligands. Families of TLR molecules have been described from the genomes of the teleost fish model species zebrafish and Takifugu, but much research remains to characterize the full length sequences and pathogen specificities of individual TLR members in fish. While the majority of these pathogen receptors are conserved among vertebrate species with clear orthologues present in fish for most mammalian TLRs, several interesting differences are present in the TLR repertoire of teleost fish when compared to that of mammals. A soluble form of TLR5 has been reported from salmonid fish and Takifugu rubripes which is not present in mammals, and a large group of TLRs (arbitrarily numbered 19-23) was identified from teleost genomes with no easily discernible orthologues in mammals. To better understand these teleost adaptations to the TLR family, we have isolated, sequenced, and characterized the full-length cDNA and gene sequences of TLR5S, TLR20, and TLR21 from catfish as well as studied their expression pattern in tissues. We also mapped these genes to bacterial artificial chromosome (BAC) clones for genome analysis. While TLR5S appeared to be common in teleost fish, and TLR21 is common to birds, amphibians and fish, TLR20 has only been identified in zebrafish and catfish. Phylogenetic analysis of catfish TLR20 indicated that it is closely related to murine TLR11 and TLR12, two divergent TLRs about which little is known. All three genes appear to exist in catfish as single copy genes.     

Towards the delineation of the ancestral eutherian genome organization: comparative genome maps of human and the African elephant (Loxodonta africana) generated by chromosome painting

This study presents a whole-genome comparison of human and a representative of the Afrotherian clade, the African elephant, generated by reciprocal Zoo-FISH. An analysis of Afrotheria genomes is of special interest, because recent DNA sequence comparisons identify them as the oldest placental mammalian clade. Complete sets of whole-chromosome specific painting probes for the African elephant and human were constructed by degenerate oligonucleotide-primed PCR amplification of flow-sorted chromosomes. Comparative genome maps are presented based on their hybridization patterns. These maps show that the elephant has a moderately rearranged chromosome complement when compared to humans. The human paint probes identified 53 evolutionary conserved segments on the 27 autosomal elephant chromosomes and the X chromosome. Reciprocal experiments with elephant probes delineated 68 conserved segments in the human genome. The comparison with a recent aardvark and elephant Zoo-FISH study delineates new chromosomal traits which link the two Afrotherian species phylogenetically. In the absence of any morphological evidence the chromosome painting data offer the first non-DNA sequence support for an Afrotherian clade. The comparative human and elephant genome maps provide new insights into the karyotype organization of the proto-afrotherian, the ancestor of extant placental mammals, which most probably consisted of 2n=46 chromosomes.     

Subfunction partitioning, the teleost radiation and the annotation of the human genome

Half of all vertebrate species are teleost fish. What accounts for this explosion of biodiversity? Recent evidence and advances in evolutionary theory suggest that genomic features could have played a significant role in the teleost radiation. This review examines evidence for an ancient whole-genome duplication (tetraploidization) event that probably occurred just before the teleost radiation. The partitioning of ancestral subfunctions between gene copies arising from this duplication could have contributed to the genetic isolation of populations, to lineage-specific diversification of developmental programs, and ultimately to phenotypic variation among teleost fish. Beyond its importance for understanding mechanisms that generate biodiversity, the partitioning of subfunctions between teleost co-orthologs of human genes can facilitate the identification of tissue-specific conserved noncoding regions and can simplify the analysis of ancestral gene functions obscured by pleiotropy or haploinsufficiency. Applying these principles on a genomic scale can accelerate the functional annotation of the human genome and understanding of the roles of human genes in health and disease.     

The Evolution of the Secreted Regulatory Protein Progranulin

Progranulin is a secreted growth factor that is active in tumorigenesis, wound repair, and inflammation. Haploinsufficiency of the human progranulin gene, GRN, causes frontotemporal dementia. Progranulins are composed of chains of cysteine-rich granulin modules. Modules may be released from progranulin by proteolysis as 6kDa granulin polypeptides. Both intact progranulin and some of the granulin polypeptides are biologically active. The granulin module occurs in certain plant proteases and progranulins are present in early diverging metazoan clades such as the sponges, indicating their ancient evolutionary origin. There is only one Grn gene in mammalian genomes. More gene-rich Grn families occur in teleost fish with between 3 and 6 members per species including short-form Grns that have no tetrapod counterparts. Our goals are to elucidate progranulin and granulin module evolution by investigating (i): the origins of metazoan progranulins (ii): the evolutionary relationships between the single Grn of tetrapods and the multiple Grn genes of fish (iii): the evolution of granulin module architectures of vertebrate progranulins (iv): the conservation of mammalian granulin polypeptide sequences and how the conserved granulin amino acid sequences map to the known three dimensional structures of granulin modules. We report that progranulin-like proteins are present in unicellular eukaryotes that are closely related to metazoa suggesting that progranulin is among the earliest extracellular regulatory proteins still employed by multicellular animals. From the genomes of the elephant shark and coelacanth we identified contemporary representatives of a precursor for short-from Grn genes of ray-finned fish that is lost in tetrapods. In vertebrate Grns pathways of exon duplication resulted in a conserved module architecture at the amino-terminus that is frequently accompanied by an unusual pattern of tandem nearly identical module repeats near the carboxyl-terminus. Polypeptide sequence conservation of mammalian granulin modules identified potential structure-activity relationships that may be informative in designing progranulin based therapeutics.     

Rearrangement rate following the whole-genome duplication in teleosts

It is now clear that a whole-genome duplication (WGD) occurred at the base of the teleost fish lineage. Like the other anciently polyploid genomes investigated so far, teleost genomes now behave like diploids with chromosomes forming pairs at meiosis. The diploidization process is currently poorly understood. It is associated with many gene deletions, such that one of the duplicates is lost at most loci and has also been proposed to coincide with an increase in genomic instability. Here we ask whether WGD is a determinant of the genomic rearrangement rate in teleosts. We study variability of the rates of rearrangement along a vertebrate phylogenetic tree, composed of 3 tetrapods (human, chicken, and mouse) and 3 teleost fishes (zebrafish, Tetraodon, and Takifugu), whose complete genome sequences are available. We devise a simple parsimony method for counting rearrangements, which takes into account various methodological complications caused by the WGD and the subsequent gene losses. We show that there does appear to be an increase in rearrangement rate after WGD, but that there is also a great deal of additional variability in rearrangement rates across species.     

The DNA sequence of medaka chromosome LG22

We report the genomic DNA sequence of a single chromosome (linkage group 22; LG22) of the small teleost fish medaka (Oryzias latipes) as a first whole chromosome sequence from a non-mammalian vertebrate. The order and orientation of 633 protein-coding genes were deduced from 18,803,338 bp of DNA sequence, providing the opportunity to analyze chromosome evolution of vertebrate genomes by direct comparison with the human genome. The average number of genes in the "conserved gene cluster" (CGC), a strict definition of "synteny" at the sequence basis, between medaka and human was 1.6. These and other data suggest that approximately 38.8% of pair-wise gene relationships would have been broken from their common ancestor in the human and medaka lineages and further imply that approx 20,000 (15,520-23,280) breaks would have occurred from the entire genome of the common ancestor. These breaks were generated mainly by intra-chromosomal shufflings at a specific era in the vertebrate lineage. These precise comparative genomics allowed us to identify the pieces of ancient chromosomes of the common vertebrate ancestor and estimate chromosomal evolution in the vertebrate lineage.     

A model of the statistical power of comparative genome sequence analysis

Comparative genome sequence analysis is powerful, but sequencing genomes is expensive. It is desirable to be able to predict how many genomes are needed for comparative genomics, and at what evolutionary distances. Here I describe a simple mathematical model for the common problem of identifying conserved sequences. The model leads to some useful rules of thumb. For a given evolutionary distance, the number of comparative genomes needed for a constant level of statistical stringency in identifying conserved regions scales inversely with the size of the conserved feature to be detected. At short evolutionary distances, the number of comparative genomes required also scales inversely with distance. These scaling behaviors provide some intuition for future comparative genome sequencing needs, such as the proposed use of “phylogenetic shadowing” methods using closely related comparative genomes, and the feasibility of high-resolution detection of small conserved features.     

High-resolution comparative mapping among man,cattle and mouse suggests a role for repeat sequences in mammalian genome evolution

Background

Sharks are members of the taxonomic class Chondrichthyes, the oldest living jawed vertebrates. Genomic studies of this group, in comparison to representative species in other vertebrate taxa, will allow us to theorize about the fundamental genetic, developmental, and functional characteristics in the common ancestor of all jawed vertebrates.

Aims

In order to obtain mapping and sequencing data for comparative genomics, we constructed a bacterial artificial chromosome (BAC) library for the nurse shark, Ginglymostoma cirratum.

Results

The BAC library consists of 313,344 clones with an average insert size of 144 kb, covering ~4.5 × 10¹⁰ bp and thus providing an 11-fold coverage of the haploid genome. BAC end sequence analyses revealed, in addition to LINEs and SINEs commonly found in other animal and plant genomes, two new groups of nurse shark-specific repetitive elements, NSRE1 and NSRE2 that seem to be major components of the nurse shark genome. Screening the library with single-copy or multi-copy gene probes showed 6–28 primary positive clones per probe of which 50–90% were true positives, demonstrating that the BAC library is representative of the different regions of the nurse shark genome. Furthermore, some BAC clones contained multiple genes, making physical mapping feasible.

Conclusion

We have constructed a deep-coverage, high-quality, large insert, and publicly available BAC library for a cartilaginous fish. It will be very useful to the scientific community interested in shark genomic structure, comparative genomics, and functional studies. We found two new groups of repetitive elements specific to the nurse shark genome, which may contribute to the architecture and evolution of the nurse shark genome.     

Hox Gene Clusters of Early Vertebrates: Do They Serve as Reliable Markers for Genome Evolution?

Hox genes,responsible for regional specification along the anteroposterior axis in embryogenesis,are found as clusters in most eumetazoan genomes sequenced to date.Invertebrates possess a single Hox gene cluster with some exceptions of secondary cluster breakages, while osteichthyans (bony vertebrates) have multiple Hox clusters. In tetrapods, four Hox clusters,derived from the so-called two-round whole genome duplications (2R-WGDs),are observed.Overall,the number of Hox gene clusters has been regarded as a reliable marker of ploidy levels in animal genomes. In fact, this scheme also fits the situations in teleost fishes that experienced an additional WGD. In this review, I focus on cyclostomes and cartilaginous fishes as lineages that would fill the gap between invertebrates and osteichthyans.A recent study highlighted a possible loss of the HoxC cluster in the galeomorph shark lineage, while other aspects of cartilaginous fish Hox clusters usually mark their conserved nature. In contrast,existing resources suggest that the cyclostomes exhibit a different mode of Hox cluster organization.For this group of species,whose genomes could have differently responded to the 2R-WGDs from jawed vertebrates,therefore the number of Hox clusters may not serve as a good indicator of their ploidy level.     

Multi-species comparative mapping in silico using the COMPASS strategy

MOTIVATION: The completion of human and mouse genome sequences provides a valuable resource for decoding other mammalian genomes. The comparative mapping by annotation and sequence similarity (COMPASS) strategy takes advantage of the resource and has been used in several genome-mapping projects. It uses existing comparative genome maps based on conserved regions to predict map locations of a sequence. An automated multiple-species COMPASS tool can facilitate in the genome sequencing effort and comparative genomics study of other mammalian species. RESULTS: The prerequisite of COMPASS is a comparative map table between the reference genome and the predicting genome. We have built and collected comparative maps among five species including human, cattle, pig, mouse and rat. Cattle-human and pig-human comparative maps were built based on the positions of orthologous markers and the conserved synteny groups between human and cattle and human and pig genomes, respectively. Mouse-human and rat-human comparative maps were based on the conserved sequence segments between the two genomes. With a match to human genome sequences, the approximate location of a query sequence can be predicted in cattle, pig, mouse and rat genomes based on the position of the match relatively to the orthologous markers or the conserved segments. AVAILABILITY: The COMPASS-tool and databases are available at http://titan.biotec.uiuc.edu/COMPASS/     

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.