Viral nervous necrosis resistance in gilthead sea bream (Sparus aurata) at the larval stage: heritability and accuracy of genomic prediction with different training and testing settings.

BACKGROUND: The gilthead sea bream (Sparus aurata) has long been considered resistant to viral nervous necrosis (VNN), until recently, when significant mortalities caused by a reassortant nervous necrosis virus (NNV) strain were reported. Selective breeding to enhance resistance against NNV might be a preventive action. In this study, 972 sea bream larvae were subjected to a NNV challenge test and the symptomatology was recorded. All the experimental fish and their parents were genotyped using a genome-wide single nucleotide polymorphism (SNP) array consisting of over 26,000 markers. RESULTS: Estimates of pedigree-based and genomic heritabilities of VNN symptomatology were consistent with each other (0.21, highest posterior density interval at 95% (HPD95%): 0.1-0.4; 0.19, HPD95%: 0.1-0.3, respectively). The genome-wide association study suggested one genomic region, i.e., in linkage group (LG) 23 that might be involved in sea bream VNN resistance, although it was far from the genome-wide significance threshold. The accuracies (r) of the predicted estimated breeding values (EBV) provided by three Bayesian genomic regression models (Bayes B, Bayes C, and Ridge Regression) were consistent and on average were equal to 0.90 when assessed in a set of cross-validation (CV) procedures. When genomic relationships between training and testing sets were minimized, accuracy decreased greatly (r = 0.53 for a validation based on genomic clustering, r = 0.12 for a validation based on a leave-one-family-out approach focused on the parents of the challenged fish). Classification of the phenotype using the genomic predictions of the phenotype or using the genomic predictions of the pedigree-based, all data included, EBV as classifiers was moderately accurate (area under the ROC curve 0.60 and 0.66, respectively). CONCLUSIONS: The estimate of the heritability for VNN symptomatology indicates that it is feasible to implement selective breeding programs for increased resistance to VNN of sea bream larvae/juveniles. Exploiting genomic information offers the opportunity of developing prediction tools for VNN resistance, and genomic models can be trained on EBV using all data or phenotypes, with minimal differences in classification performance of the trait phenotype. In a long-term view, the weakening of the genomic ties between animals in the training and test sets leads to decreased genomic prediction accuracies, thus periodical update of the reference population with new data is mandatory.

Resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax L.): heritability and relationships with body weight, cortisol concentration, and antibody titer.

BACKGROUND: Susceptibility of European sea bass (Dicentrarchus labrax L.) to viral nervous necrosis (VNN) is well-known. Interest towards selective breeding as a tool to enhance genetic resistance in this species has increased sharply due to the major threat represented by VNN for farmed sea bass and limitations concerning specific therapeutical measures. A sea bass experimental population (N = 650) was challenged with nervous necrosis virus (NNV) to investigate genetic variation in VNN mortality. In addition, relationships of this trait with serum cortisol concentration after stress exposure, antibody titer against NNV antigens, and body weight at a fixed age were studied to identify potential indicator traits of VNN resistance. RESULTS: The estimate of heritability for VNN mortality was moderate and ranged from 0.15 (HPD95%, 95% highest posterior density interval: 0.02, 0.31) to 0.23 (HPD95%: 0.06, 0.47). Heritability estimates for cortisol concentration, antibody titer, and body weight were 0.19 (HPD95%: 0.07, 0.34), 0.36 (HPD95%: 0.16, 0.59) and 0.57 (HPD95%: 0.33, 0.84), respectively. Phenotypic relationships between traits were trivial and not statistically significant, except for the estimated correlation between antibody titer and body weight (0.24). Genetic correlations of mortality with body weight or antibody titer (- 0.39) exhibited a 0.89 probability of being negative. A negligible genetic correlation between mortality and cortisol concentration was detected. Antibody titer was estimated to be positively correlated with body weight (0.49). CONCLUSIONS: Antibody titer against NNV offers the opportunity to use indirect selection to enhance resistance, while the use of cortisol concentration as an indicator trait in breeding programs for VNN resistance is questionable. The estimate of heritability for VNN mortality indicates the feasibility of selective breeding to enhance resistance to NNV and raises attention to the development of genomic prediction tools to simplify testing procedures for selection candidates.

Helping decision making for reliable and cost-effective 2b-RAD sequencing and genotyping analyses in non-model species.

High-throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b-RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b-RAD protocols on non-model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b-RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade-off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective-base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b-RAD protocols on non-model organisms with different genome sizes, helping decision-making for a reliable and cost-effective genotyping.

Genomic analysis of Sparus aurata reveals the evolutionary dynamics of sex-biased genes in a sequential hermaphrodite fish.

Sexual dimorphism is a fascinating subject in evolutionary biology and mostly results from sex-biased expression of genes, which have been shown to evolve faster in gonochoristic species. We report here genome and sex-specific transcriptome sequencing of Sparus aurata, a sequential hermaphrodite fish. Evolutionary comparative analysis reveals that sex-biased genes in S. aurata are similar in number and function, but evolved following strikingly divergent patterns compared with gonochoristic species, showing overall slower rates because of stronger functional constraints. Fast evolution is observed only for highly ovary-biased genes due to female-specific patterns of selection that are related to the peculiar reproduction mode of S. aurata, first maturing as male, then as female. To our knowledge, these findings represent the first genome-wide analysis on sex-biased loci in a hermaphrodite vertebrate species, demonstrating how having two sexes in the same individual profoundly affects the fate of a large set of evolutionarily relevant genes.

The population genomics of yellowfin tuna (Thunnus albacares) at global geographic scale challenges current stock delineation.

Yellowfin tuna, Thunnus albacares, is one of the most important seafood commodities in the world. Despite its great biological and economic importance, conflicting evidence arises from classical genetic and tagging studies concerning the yellowfin tuna population structure at local and global oceanic scales. Access to more powerful and cost effective genetic tools would represent the first step towards resolving the population structure of yellowfin tuna across its distribution range. Using a panel of 939 neutral Single Nucleotide Polymorphisms (SNPs), and the most comprehensive data set of yellowfin samples available so far, we found genetic differentiation among the Atlantic, Indian and Pacific oceans. The genetic stock structure analysis carried out with 33 outlier SNPs, putatively under selection, identified discrete populations within the Pacific Ocean and, for the first time, also within the Atlantic Ocean. Stock assessment approaches that consider genetic differences at neutral and adaptive genomic loci should be routinely implemented to check the status of the yellowfin tuna, prevent illegal trade, and develop more sustainable management measures.

Genetic differentiation and phylogeography of Mediterranean-North Eastern Atlantic blue shark (Prionace glauca, L. 1758) using mitochondrial DNA: panmixia or complex stock structure?

BACKGROUND: The blue shark (Prionace glauca, Linnaeus 1758) is one of the most abundant epipelagic shark inhabiting all the oceans except the poles, including the Mediterranean Sea, but its genetic structure has not been confirmed at basin and interoceanic distances. Past tagging programs in the Atlantic Ocean failed to find evidence of migration of blue sharks between the Mediterranean and the adjacent Atlantic, despite the extreme vagility of the species. Although the high rate of by-catch in the Mediterranean basin, to date no genetic study on Mediterranean blue shark was carried out, which constitutes a significant knowledge gap, considering that this population is classified as "Critically Endangered", unlike its open-ocean counterpart. METHODS: Blue shark phylogeography and demography in the Mediterranean Sea and North-Eastern Atlantic Ocean were inferred using two mitochondrial genes (Cytb and control region) amplified from 207 and 170 individuals respectively, collected from six localities across the Mediterranean and two from the North-Eastern Atlantic. RESULTS: Although no obvious pattern of geographical differentiation was apparent from the haplotype network, Φst analyses indicated significant genetic structure among four geographical groups. Demographic analyses suggest that these populations have experienced a constant population expansion in the last 0.4-0.1 million of years. DISCUSSION: The weak, but significant, differences in Mediterranean and adjacent North-eastern Atlantic blue sharks revealed a complex phylogeographic structure, which appears to reject the assumption of panmixia across the study area, but also supports a certain degree of population connectivity across the Strait of Gibraltar, despite the lack of evidence of migratory movements observed by tagging data. Analyses of spatial genetic structure in relation to sex-ratio and size could indicate some level of sex/stage biased migratory behaviour.

An integrated genomic approach for the study of mandibular prognathism in the European seabass (Dicentrarchus labrax).

Skeletal anomalies in farmed fish are a relevant issue affecting animal welfare and health and causing significant economic losses. Here, a high-density genetic map of European seabass for QTL mapping of jaw deformity was constructed and a genome-wide association study (GWAS) was carried out on a total of 298 juveniles, 148 of which belonged to four full-sib families. Out of 298 fish, 107 were affected by mandibular prognathism (MP). Three significant QTLs and two candidate SNPs associated with MP were identified. The two GWAS candidate markers were located on ChrX and Chr17, both in close proximity with the peaks of the two most significant QTLs. Notably, the SNP marker on Chr17 was positioned within the Sobp gene coding region, which plays a pivotal role in craniofacial development. The analysis of differentially expressed genes in jaw-deformed animals highlighted the "nervous system development" as a crucial pathway in MP. In particular, Zic2, a key gene for craniofacial morphogenesis in model species, was significantly down-regulated in MP-affected animals. Gene expression data revealed also a significant down-regulation of Sobp in deformed larvae. Our analyses, integrating transcriptomic and GWA methods, provide evidence for putative mechanisms underlying seabass jaw deformity.

Genomic Prediction of Resistance to Pasteurellosis in Gilthead Sea Bream (Sparus aurata) Using 2b-RAD Sequencing.

Gilthead sea bream (Sparus aurata) is a species of paramount importance to the Mediterranean aquaculture industry, with an annual production exceeding 140,000 metric tons. Pasteurellosis due to the Gram-negative bacterium Photobacterium damselae subsp. piscicida (Phdp) causes significant mortality, especially during larval and juvenile stages, and poses a serious threat to bream production. Selective breeding for improved resistance to pasteurellosis is a promising avenue for disease control, and the use of genetic markers to predict breeding values can improve the accuracy of selection, and allow accurate calculation of estimated breeding values of nonchallenged animals. In the current study, a population of 825 sea bream juveniles, originating from a factorial cross between 67 broodfish (32 sires, 35 dams), were challenged by 30 min immersion with 1 × 105 CFU virulent Phdp. Mortalities and survivors were recorded and sampled for genotyping by sequencing. The restriction-site associated DNA sequencing approach, 2b-RAD, was used to generate genome-wide single nucleotide polymorphism (SNP) genotypes for all samples. A high-density linkage map containing 12,085 SNPs grouped into 24 linkage groups (consistent with the karyotype) was constructed. The heritability of surviving days (censored data) was 0.22 (95% highest density interval: 0.11-0.36) and 0.28 (95% highest density interval: 0.17-0.4) using the pedigree and the genomic relationship matrix respectively. A genome-wide association study did not reveal individual SNPs significantly associated with resistance at a genome-wide significance level. Genomic prediction approaches were tested to investigate the potential of the SNPs obtained by 2b-RAD for estimating breeding values for resistance. The accuracy of the genomic prediction models (r = 0.38-0.46) outperformed the traditional BLUP approach based on pedigree records (r = 0.30). Overall results suggest that major quantitative trait loci affecting resistance to pasteurellosis were not present in this population, but highlight the effectiveness of 2b-RAD genotyping by sequencing for genomic selection in a mass spawning fish species.

RAD SNP markers as a tool for conservation of dolphinfish Coryphaena hippurus in the Mediterranean Sea: Identification of subtle genetic structure and assessment of populations sex-ratios.

Dolphinfish is an important fish species for both commercial and sport fishing, but so far limited information is available on genetic variability and pattern of differentiation of dolphinfish populations in the Mediterranean basin. Recently developed techniques allow genome-wide identification of genetic markers for better understanding of population structure in species with limited genome information. Using restriction-site associated DNA analysis we successfully genotyped 140 individuals of dolphinfish from eight locations in the Mediterranean Sea at 3324 SNP loci. We identified 311 sex-related loci that were used to assess sex-ratio in dolphinfish populations. In addition, we identified a weak signature of genetic differentiation of the population closer to Gibraltar Strait in comparison to other Mediterranean populations, which might be related to introgression of individuals from Atlantic. No further genetic differentiation could be detected in the other populations sampled, as expected considering the known highly mobility of the species. The results obtained improve our knowledge of the species and can help managing dolphinfish stock in the future.

Methodological assessment of 2b-RAD genotyping technique for population structure inferences in yellowfin tuna (Thunnus albacares).

Global population genetic structure of yellowfin tuna (Thunnus albacares) is still poorly understood despite its relevance for the tuna fishery industry. Low levels of genetic differentiation among oceans speak in favour of the existence of a single panmictic population worldwide of this highly migratory fish. However, recent studies indicated genetic structuring at a much smaller geographic scales than previously considered, pointing out that YFT population genetic structure has not been properly assessed so far. In this study, we demonstrated for the first time, the utility of 2b-RAD genotyping technique for investigating population genetic diversity and differentiation in high gene-flow species. Running de novo pipeline in Stacks, a total of 6772 high-quality genome-wide SNPs were identified across Atlantic, Indian and Pacific population samples representing all major distribution areas. Preliminary analyses showed shallow but significant population structure among oceans (FST=0.0273; P-value<0.01). Discriminant Analysis of Principal Components endorsed the presence of genetically discrete yellowfin tuna populations among three oceanic pools. Although such evidence needs to be corroborated by increasing sample size, these results showed the efficiency of this genotyping technique in assessing genetic divergence in a marine fish with high dispersal potential.

Second generation genetic linkage map for the gilthead sea bream Sparus aurata L.

An updated second linkage map was constructed for the gilthead sea bream, Sparus aurata L., a fish species of great economic importance for the Mediterranean aquaculture industry. In contrast to the first linkage map which mainly consisted of genomic microsatellites (SSRs), the new linkage map is highly enriched with SSRs found in Expressed Sequence Tags (EST-SSRs), which greatly facilitates comparative mapping with other teleosts. The new map consists of 321 genetic markers in 27 linkage groups (LGs): 232 genomic microsatellites, 85 EST-SSRs and 4 SNPs; of those, 13 markers were linked to LGs but were not ordered. Eleven markers (5 SSRs, 5 EST-SSRs and 1 SNP) are not assigned to any LG. The total length of the sex-averaged map is 1769.7cM, 42% longer than the previously published one, and the number of markers in each LG ranges from 2 to 30. The inter-marker distance varies from 0 to 75.6cM, with an average of 5.75cM. The male and female maps have a length of 1349.2 and 2172.1cM, respectively, and the average distance between markers is 4.38 and 7.05cM, respectively. Comparative mapping with the three-spined stickleback (Gasterosteus acuulatus) chromosomes and scaffolds showed conserved synteny with 132 S. aurata markers (42.9% of those mapped) having a hit on the stickleback genome.

High mortality of juvenile gilthead sea bream (Sparus aurata) from photobacteriosis is associated with alternative macrophage activation and anti-inflammatory response: results of gene expression profiling of early responses in the head kidney.

The halophilic bacterium Photobacterium damselae subsp. piscicida (Phdp) represents a substantial health problem for several fish species in aquaculture. Bacteria that reside free and inside phagocytes cause acute and chronic forms of photobacteriosis. Infections of juveniles rapidly kill up to 90-100% fish. Factors underlying failure of the immune protection against bacteria remain largely unknown. The reported study used a transcriptomic approach to address this issue. Juvenile sea breams (0.5 g) were challenged by immersion in salt water containing 2.89 × 10(8) CFU of a virulent Phdp and the head kidney was sampled after 24- and 48-h. Analyses were performed using the second version of a 44 k oligonucleotide DNA microarray that represents 19,734 sea bream unique transcripts and covers diverse immune pathways. Expression changes of selected immune genes were validated with qPCR. Results suggested rapid recognition of the pathogen, as testified by up-regulation of lectins and antibacterial proteins (bactericidal permeability-increasing protein lectins, lysozyme, intracellular and extracellular proteases), chemokines and chemokine receptors. Increased expression of proteins involved in iron and heme metabolism also could be a response against bacteria that are dependent on iron. However, negative regulators of immune/inflammatory response were preponderant among the up-regulated genes. A remarkable finding was the increased expression of IL-10 in concert with up-regulation of arginase I and II and proteins of the polyamine biosynthesis pathway that diverts the arginine flux from the production of reactive nitrogen species. Such expression changes are characteristic for alternatively activated macrophages that do not develop acute inflammatory responses. Immune suppression can be induced by the host to reduce tissue damages or by the pathogen to evade host response.

A gene-based radiation hybrid map of the gilthead sea bream Sparus aurata refines and exploits conserved synteny with Tetraodon nigroviridis.

BACKGROUND: Comparative teleost studies are of great interest since they are important in aquaculture and in evolutionary issues. Comparing genomes of fully sequenced model fish species with those of farmed fish species through comparative mapping offers shortcuts for quantitative trait loci (QTL) detections and for studying genome evolution through the identification of regions of conserved synteny in teleosts. Here a comparative mapping study is presented by radiation hybrid (RH) mapping genes of the gilthead sea bream Sparus aurata, a non-model teleost fish of commercial and evolutionary interest, as it represents the worldwide distributed species-rich family of Sparidae. RESULTS: An additional 74 microsatellite markers and 428 gene-based markers appropriate for comparative mapping studies were mapped on the existing RH map of Sparus aurata. The anchoring of the RH map to the genetic linkage map resulted in 24 groups matching the karyotype of Sparus aurata. Homologous sequences to Tetraodon were identified for 301 of the gene-based markers positioned on the RH map of Sparus aurata. Comparison between Sparus aurata RH groups and Tetraodon chromosomes (karyotype of Tetraodon consists of 21 chromosomes) in this study reveals an unambiguous one-to-one relationship suggesting that three Tetraodon chromosomes correspond to six Sparus aurata radiation hybrid groups. The exploitation of this conserved synteny relationship is furthermore demonstrated by in silico mapping of gilthead sea bream expressed sequence tags (EST) that give a significant similarity hit to Tetraodon. CONCLUSION: The addition of primarily gene-based markers increased substantially the density of the existing RH map and facilitated comparative analysis. The anchoring of this gene-based radiation hybrid map to the genome maps of model species broadened the pool of candidate genes that mainly control growth, disease resistance, sex determination and reversal, reproduction as well as environmental tolerance in this species, all traits of great importance for QTL mapping and marker assisted selection. Furthermore this comparative mapping approach will facilitate to give insights into chromosome evolution and into the genetic make up of the gilthead sea bream.

A genetic linkage map of the hermaphrodite teleost fish Sparus aurata L.

The gilthead sea bream (Sparus aurata L.) is a marine fish of great importance for fisheries and aquaculture. It has also a peculiar sex-determination system, being a protandrous hermaphrodite. Here we report the construction of a first-generation genetic linkage map for S. aurata, based on 204 microsatellite markers. Twenty-six linkage groups (LG) were found. The total map length was 1241.9 cM. The ratio between sex-specific map lengths was 1:1.2 (male:female). Comparison with a preliminary radiation hybrid (RH) map reveals a good concordance, as all markers located in a single LG are located in a single RH group, except for Ad-25 and CId-31. Comparison with the Tetraodon nigroviridis genome revealed a considerable number of evolutionary conserved regions (ECRs) between the two species. The mean size of ECRs was 182 bp (sequence identity 60-90%). Forty-one ECRs have a known chromosomal location in the pufferfish genome. Despite the limited number of anchoring points, significant syntenic relationships were found. The linkage map presented here provides a robust comparative framework for QTL analysis in S. aurata and is a step toward the identification of genetic loci involved both in the determination of economically important traits and in the individual timing of sex reversal.

Differential availability/processing of decorin precursor in arterial and venous smooth muscle cells.

The existence of specific differentiation markers for arterial smooth muscle (SM) cells is still a matter of debate. A clone named MM1 was isolated from a library of monoclonal antibodies to adult porcine aorta, which in vivo binds to arterial but not venous SM cells, except for the pulmonary vein. MM1 immunoreactivity in Western blotting involved bands in the range of M(r) 33-226 kDa, in both arterial and venous SM tissues. However, immunoprecipitation experiments revealed that MM1 bound to a 100-kDa polypeptide that was present only in the arterial SM extract. By mass spectrometry analysis of tryptic digests from MM1-positive 130- and 120-kDa polypeptides of aorta SM extract, the antigen recognized by the antibody was identified as a decorin precursor. Using a crude decorin preparation from this tissue MM1 reacted strongly with the 33-kDa polypeptide and this pattern did not change after chondroitinase ABC treatment. In vitro, decorin immunoreactivity was found in secreted grainy material produced by confluent arterial SM cells, although lesser amounts were also seen in venous SM cells. Western blotting of extracts from these cultures showed the presence of the 33-kDa band but not of the high-molecular-weight components, except for the 100-kDa monomer. The 100/33-kDa combination was more abundant in arterial SM cells than in the venous counterpart. In the early phase of neointima formation, induced by endothelial injury of the carotid artery or vein-to-artery transposition, the decorin precursor was not expressed, but it was up-regulated in the SM cells of the media underlying the neointima in both models. Collectively, these data suggest a different processing/utilization of the 100-kDa monomer of proteoglycan decorin in arterial and venous SM cells, which is abolished after vein injury.

Full-length sequence and expression analysis of Toll-like receptor 9 in the gilthead seabream (Sparus aurata L.).

Toll-Like Receptors (TLRs) have recently emerged as key sensors of invading microbes, acting through recognition of pathogen-associated molecular patterns. It has been demonstrated that TLR9 is involved in the recognition of unmethylated CpG motifs in mice, humans, and pigs. We report here the full-length sequence of TLR9 cDNA in the gilthead sea bream (Sparus aurata L.). The predicted protein (1063 amino acids) was similar to mammalian TLR9s, showing 21 leucine-rich repeats in the extracellular region and a typical Toll/IL-1R (TIR) domain in the intracellular region. Comparative analysis of TLR9 sequences indicated that critical residues for ligand-binding are conserved across vertebrate lineages, although evidence of functional divergence was observed. Analysis of the genomic structure of sea bream TLR9 gene revealed the presence of two intervening sequences. Retention of the second intron produced an alternatively spliced mRNA (TLR9B) showing differential expression among tissues or developmental stages compared to the wild-type isoform (TLR9A). RT-PCR analysis indicated a broad expression of TLR9A, especially in immune-related organs (spleen, head-kidney) and mucosal-epithelial barriers (gills, gut, skin). Using quantitative Real-Time RT-PCR, no statistically significant variation was observed for TLR9 mRNAs expression in the spleen of experimentally infected animals compared to healthy controls. Comparing sequence and expression profile of sea bream TLR9 with mammalian TLR9s suggested that the main function of TLR9 might be conserved across vertebrates, although species-specific features are present (modulation of ligand-binding specificity, alternative splicing).

The first radiation hybrid map of a perch-like fish: the gilthead seabream (Sparus aurata L).

Among Teleosts, Perciformes are the largest order of fishes and include numerous species of commercial importance. Perciformes also comprise species of primary interest for evolutionary studies and analysis of the sex determination systems and sex chromosome plasticity. Unfortunately, genomics tools and resources for Perciformes remain to be developed. Here, we report the production of a seabream whole-genome radiation hybrid (RH) panel in which quality was ascertained by the construction of a 2-Mb-resolution RH map. The map encompasses 440 markers (288 microsatellites, 82 gene-based markers, and 70 STS) suitable for linkage analysis and comparative mapping studies. Achievement of a RH panel and a whole-genome RH map should contribute to establishing seabream as a fish model among the Perciformes and should be of importance in aquaculture for marker-assisted selection, improvement of growth performance, and disease management. Development of RH maps in a cost-effective manner for other fishes with the described methodology will offer a powerful approach in aquaculture and will provide extended capabilities for comparing vertebrate genome evolution.