A type I and type II microsatellite linkage map of rainbow trout (Oncorhynchus mykiss) with presumptive coverage of all chromosome arms

APIS: an updated parentage assignment software managing triploids induced from diploid parents

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields, and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximizes the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96-192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1,232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low false discovery rate (<0.01), and excellent specificity (0.96-0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents.

Developmental competence of eggs produced by rainbow trout Doubled Haploids (DHs) and generation of the clonal lines

Poor quality eggs produced by the fully homozygous doubled haploids (DHs) may impair generation of clonal lines in fish species. In the present research, gynogenetic development of rainbow trout (Oncorhynchus mykiss) was induced in eggs originated from the DH females. Eggs were activated with the UV-irradiated grayling (Thymallus thymallus) spermatozoa and subjected to the high hydrostatic pressure (HHP) shock to provide diploid clonal individuals. Only two of four DH females produced eggs that were successfully activated by the irradiated spermatozoa and subsequently developed into the gynogenetic embryos. Survival rates of rainbow trout from the clonal lines equalled 21.5% and 19.8% during embryogenesis and decreased after hatching to 18.6% and 14.9%, respectively. Some of the dead rainbow trout clones collected between hatching and swim-up stage were emaciated and exhibited spinal deformities including scoliosis. Provided results confirmed limited developmental competences of eggs produced by rainbow trout DH females. Clonal rainbow trout developing in such eggs exhibited reduced survival and increased frequency of the body abnormalities.

Influence of epistatic segregation distortion loci on genetic marker linkages in Japanese flounder

For genetic linkage analysis of Japanese flounder, 160 doubled haploids (DH) were artificially produced using mitotic gynogenesis and were genotyped for 458 simple sequence repeat (SSR) markers, 101 of which show distortional segregation. The genetic linkage map was constructed by modifying recombination fractions between the distorted markers. Between the corrected and uncorrected genetic maps, there were considerable differences in genetic distance, but not in relative locations among markers. Using a liability model, a segregation distortion locus (SDL), with an additive genetic effect of 1.772, was mapped between markers BDHYP387 and Poli56TUF of chromosome 24 in the corrected genetic map. Additionally, six pairs of epistatic SDLs were identified on chromosomes 1, 5, 8, 9, 23, and 24. Changes in genetic distances between markers did not occur on chromosome regions with main effect SDLs. However, most chromosome regions where genetic distances changed covered the detected epistatic SDLs. This study concluded that epistatic SDLs decrease linkages between markers and lengthen genetic distances in Japanese flounder. This finding has been partially validated in other DH populations derived from three female Japanese flounders.

Gene-centromere mapping in meiotic gynogenetic European seabass

Background

Fully isogenic lines in fish can be developed using “mitotic” gynogenesis (suppression of first zygotic mitosis following inactivation of the sperm genome). However, genome-wide verification of the steps in this process has seldom been applied. We used ddRADseq to generate SNP markers in a meiotic gynogenetic family of European seabass (Dicentrarchus labrax): (i) to verify the lack of paternal contribution in a meiotic gynogenetic family; (ii) to generate a gene-centromere map from this family; (iii) to identify telomeric markers that could distinguish mitotic gynogenetics from meiotic gynogenetics, which sometimes arise spontaneously in mitotic gynogenetic families.

Results

From a single meiotic gynogenetic family consisting of 79 progeny, 42 million sequencing reads (Illumina, trimmed to 148 bases) resolved 6866 unique RAD-tags. The 340 male-informative SNP markers that were identified confirmed the lack of paternal contribution. A gene-centromere map was constructed based on 804 female-informative SNPs in 24 linkage groups (2n = 48) with a total length of 1251.02 cM (initial LG assignment was based on the seabass genome assembly, dicLab v1). Chromosome arm structure could be clearly discerned from the pattern of heterozygosity in each linkage group in 18 out of 24 LGs: the other six showed anomalies that appeared to be related to issues in the genome assembly.

Conclusion

Genome-wide screening enabled substantive verification of the production of the gynogenetic family used in this study. The large number of telomeric and subtelomeric markers with high heterozygosity values in the meiotic gynogenetic family indicate that such markers could be used to clearly distinguish between meiotic and mitotic gynogenetics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3826-z) contains supplementary material, which is available to authorized users.

Construction of a High-Density Genetic Map and Quantitative Trait Locus Mapping in the Manila clam Ruditapes philippinarum

Genetic linkage maps are indispensable tools in a wide range of genetic and genomic research. With the advancement of genotyping-by-sequencing (GBS) methods, the construction of a high-density linkage maps has become achievable in marine organisms lacking sufficient genomic resources, such as mollusks. In this study, high-density linkage map was constructed for an ecologically and commercially important clam species, Ruditapes philippinarum. For the consensus linkage map, a total of 9658 markers spanning 1926.98 cM were mapped to 18 sex-averaged linkage groups, with an average marker distance of 0.42 cM. Based on the high-density linkage map, ten QTLs for growth-related traits and shell color were detected. The coverage and density of the current map are sufficient for us to effectively detect QTL for segregating traits, and two QTL positions were all coincident with the closest markers. This high-density genetic linkage map reveals basic genomic architecture and will be useful for comparative genomics research, genome assembly and genetic improvement of R. philippinarum and other bivalve molluscan species.

Preliminary genetic linkage map of Indian major carp, Labeo rohita (Hamilton 1822) based on microsatellite markers

Linkage map with wide marker coverage is an essential resource for genetic improvement study for any species. Sex-averaged genetic linkage map of Labeo rohita, popularly known as 'rohu', widely cultured in the Indian subcontinent, was developed by placing 68 microsatellite markers generated by a simplified method. The parents and their F1 progeny (92 individuals) were used as segregating populations. The genetic linkage map spans a sex-averaged total length of 1462.2 cM, in 25 linkage groups. The genome length of rohu was estimated to be 3087.9 cM. This genetic linkage map may facilitate systematic searches of the genome to identify genes associated with commercially important characters and marker-assisted selection programmes of this species.

Genome survey and high-density genetic map construction provide genomic and genetic resources for the Pacific White Shrimp Litopenaeus vannamei

The Pacific white shrimp Litopenaeus vannamei is the dominant crustacean species in global seafood mariculture. Understanding the genome and genetic architecture is useful for deciphering complex traits and accelerating the breeding program in shrimp. In this study, a genome survey was conducted and a high-density linkage map was constructed using a next-generation sequencing approach. The genome survey was used to identify preliminary genome characteristics and to generate a rough reference for linkage map construction. De novo SNP discovery resulted in 25,140 polymorphic markers. A total of 6,359 high-quality markers were selected for linkage map construction based on marker coverage among individuals and read depths. For the linkage map, a total of 6,146 markers spanning 4,271.43 cM were mapped to 44 sex-averaged linkage groups, with an average marker distance of 0.7 cM. An integration analysis linked 5,885 genome scaffolds and 1,504 BAC clones to the linkage map. Based on the high-density linkage map, several QTLs for body weight and body length were detected. This high-density genetic linkage map reveals basic genomic architecture and will be useful for comparative genomics research, genome assembly and genetic improvement of L. vannamei and other penaeid shrimp species.

Effects of crossovers between homeologs on inheritance and population genomics in polyploid-derived salmonid fishes

A whole genome duplication occurred in the ancestor of all salmonid fishes some 50-100 million years ago. Early inheritance studies with allozymes indicated that loci in the salmonid genome are inherited disomically in females. However, some pairs of duplicated loci showed patterns of inheritance in males indicating pairing and recombination between homeologous chromosomes. Nearly 20% of loci in the salmonid genome are duplicated and share the same alleles (isoloci), apparently due to homeologous recombination. Half-tetrad analysis revealed that isoloci tend to be telomeric. These results suggested that residual tetrasomic inheritance of isoloci results from homeologous recombination near chromosome ends and that continued disomic inheritance resulted from homologous pairing of centromeric regions. Many current genetic maps of salmonids are based on single nucleotide polymorphisms and microsatellites that are no longer duplicated. Therefore, long sections of chromosomes on these maps are poorly represented, especially telomeric regions. In addition, preferential multivalent pairing of homeologs from the same species in F1 hybrids results in an excess of nonparental gametes (so-called pseudolinkage). We consider how not including duplicated loci has affected our understanding of population and evolutionary genetics of salmonids, and we discuss how incorporating these loci will benefit our understanding of population genomics.

Genetic mapping of centromeres in the nine Citrus clementina chromosomes using half-tetrad analysis and recombination patterns in unreduced and haploid gametes

BACKGROUND

Mapping centromere locations in plant species provides essential information for the analysis of genetic structures and population dynamics. The centromere's position affects the distribution of crossovers along a chromosome and the parental heterozygosity restitution by 2n gametes is a direct function of the genetic distance to the centromere. Sexual polyploidisation is relatively frequent in Citrus species and is widely used to develop new seedless triploid cultivars. The study's objectives were to (i) map the positions of the centromeres of the nine Citrus clementina chromosomes; (ii) analyse the crossover interference in unreduced gametes; and (iii) establish the pattern of genetic recombination in haploid clementine gametes along each chromosome and its relationship with the centromere location and distribution of genic sequences.

RESULTS

Triploid progenies were derived from unreduced megagametophytes produced by second-division restitution. Centromere positions were mapped genetically for all linkage groups using half-tetrad analysis. Inference of the physical locations of centromeres revealed one acrocentric, four metacentric and four submetacentric chromosomes. Crossover interference was observed in unreduced gametes, with variation seen between chromosome arms. For haploid gametes, a strong decrease in the recombination rate occurred in centromeric and pericentromeric regions, which contained a low density of genic sequences. In chromosomes VIII and IX, these low recombination rates extended beyond the pericentromeric regions. The genomic region corresponding to a genetic distance < 5cM from a centromere represented 47% of the genome and 23% of the genic sequences.

CONCLUSIONS

The centromere positions of the nine citrus chromosomes were genetically mapped. Their physical locations, inferred from the genetic ones, were consistent with the sequence constitution and recombination pattern along each chromosome. However, regions with low recombination rates extended beyond the pericentromeric regions of some chromosomes into areas richer in genic sequences. The persistence of strong linkage disequilibrium between large numbers of genes promotes the stability of epistatic interactions and multilocus-controlled traits over successive generations but also maintains multi-trait associations. Identification of the centromere positions will allow the development of simple methods to analyse unreduced gamete formation mechanisms in a large range of genotypes and further modelling of genetic inheritance in sexual polyploidisation breeding schemes.

First haploid genetic map based on microsatellite markers in Senegalese sole (Solea senegalensis, Kaup 1858)

The Senegalese sole (Solea senegalensis, Kaup 1858) is a flatfish species of great value for aquaculture. In this study, we develop the first linkage map in this species based on microsatellite markers characterized from genomic DNA libraries and EST databases of Senegalese sole and from other flatfish species. Three reference gynogenetic families were obtained by chromosome-manipulation techniques: two haploid gynogenetics, used to assign and order microsatellites to linkage groups and another diploid gynogenetic family, used for estimating marker-centromere distances. The consensus map consists of 129 microsatellites distributed in 27 linkage groups (LG), with an average density of 4.7 markers per LG and comprising 1,004 centimorgans (cM). Additionally, 15 markers remained unlinked. Through half-tetrad analysis, we were able to estimate the centromere distance for 81 markers belonging to 24 LG, representing an average of 3 markers per LG. Comparative mapping was performed between flatfish species LG and model fish species chromosomes (stickleback, Tetraodon, medaka, fugu and zebrafish). The usefulness of microsatellite markers and the genetic map as tools for comparative mapping and evolution studies is discussed.

Divergence in expression of candidate genes for the smoltification process between juvenile resident rainbow and anadromous steelhead trout

Rainbow and steelhead trout (Oncorhynchus mykiss), among other salmonid fishes, exhibit tremendous life history diversity, foremost of which is variation in migratory propensity. While some individuals possess the ability to undertake an anadromous marine migration, others remain resident in freshwater throughout their life cycle. Those that will migrate undergo tremendous physiological, morphological, and behavioral transformations in a process called smoltification which transitions freshwater-adapted parr to marine-adapted smolts. While the behavior, ecology, and physiology of smoltification are well described, our understanding of the proximate genetic mechanisms that trigger the process are not well known. Quantitative genetic analyses have identified several genomic regions associated with smoltification and migration-related traits within this species. Here we investigate the divergence in gene expression of 18 functional and positional candidate genes for the smoltification process in the brain, gill, and liver tissues of migratory smolts, resident parr, and precocious mature male trout at the developmental stage of out-migration. Our analysis reveals several genes differentially expressed between life history classes and validates the candidate nature of several genes in the parr-smolt transformation including Clock1α, FSHβ, GR, GH2, GHR1, GHR2, NDK7, p53, SC6a7, Taldo1, THRα, THRβ, and Vdac2.

Genetic mapping and QTL analysis of growth traits in the large yellow croaker Larimichthys crocea

Large yellow croaker (Larimichthys crocea) is an important maricultured species in China. A genetic linkage map of the large yellow croaker was constructed using type II microsatellites and expressed sequence tag (EST)-derived microsatellites in two half-sib families (two females and one male). A total of 289 microsatellite markers (contained 93 EST-SSRs) were integrated into 24 linkage groups, which agreed with the haploid chromosome number. The map spanned a length of 1,430.8 cm with an average interval of 5.4 cm, covering 83.9 % of the estimated genome size (1,704.8 cm). A total of seven quantitative trait locis (QTLs) were detected for growth traits on five linkage groups, including two 1 % and five 5 % chromosome-wide significant QTLs, and explained from 2.33 to 5.31 % of the trait variation. The identified QTLs can be applied in marker-assisted selection programs to improve the growth traits.

Genetics and genomics of disease resistance in salmonid species

Infectious and parasitic diseases generate large economic losses in salmon farming. A feasible and sustainable alternative to prevent disease outbreaks may be represented by genetic improvement for disease resistance. To include disease resistance into the breeding goal, prior knowledge of the levels of genetic variation for these traits is required. Furthermore, the information from the genetic architecture and molecular factors involved in resistance against diseases may be used to accelerate the genetic progress for these traits. In this regard, marker assisted selection and genomic selection are approaches which incorporate molecular information to increase the accuracy when predicting the genetic merit of selection candidates. In this article we review and discuss key aspects related to disease resistance in salmonid species, from both a genetic and genomic perspective, with emphasis in the applicability of disease resistance traits into breeding programs in salmonids.

Microsatellite-centromere mapping in common carp through half-tetrad analysis in diploid meiogynogenetic families

Gene-centromere (G-C) mapping provides insights into the understanding of the composition, structure, and evolution of vertebrate genomes. Common carp (Cyprinus carpio) is an important aquaculture fish and has been proposed to undertake tetraploidization. In this study, we selected 214 informative microsatellite markers across 50 linkage groups of a common carp genetic map to perform gene-centromere mapping using half-tetrad analysis. A total of 199 microsatellites were segregated under the Mendelian expectations in at least one of the three gynogenetic families and were used for G-C distance estimation. The G-C recombination frequency (y) ranged from 0 to 0.99 (0.43 on average), corresponding to a fixation index (F) of 0.57 after one generation of gynogenesis. Large y values for some loci together with significant correlation between G-C distances and genetic linkage map distances suggested the presence of high interference in common carp. Under the assumption of complete interference, 50 centromeres were localized onto corresponding linkage groups (LGs) of common carp, with G-C distances of centromere-linked markers per LG ranging from 0 to 10.3 cM (2.9 cM on average). Based on the information for centromere positions, we proposed a chromosome formula of 2n = 100 = 58 m/sm + 42 t/st with 158 chromosome arms for common carp, which was similar to a study observed by cytogenetic method. The examination of crossover distributions along 10 LGs revealed that the proportion of crossover chromatids was overall higher than that of non-crossover chromatids in gynogenetic progenies, indicating high recombination levels across most LGs. Comparative genomics analyses suggested that the chromosomes of common carp have undergone extensive rearrangement after genome duplication. This study would be valuable to elucidate the mechanism of genome evolution and integrate physical and genetic maps in common carp.

Molecular markers for variation in spawning date in a hatchery population of rainbow trout (Oncorhynchus mykiss)

We examined the distribution of alleles at 63 microsatellite loci distributed across 29 linkage groups in broodstock females from a commercial population of rainbow trout spawning on different dates throughout the season (August to January). A total of 368 females, 184 and 117 females from each of the tail-ends of the spawning distribution and a subsample of 67 females spawning in the middle, were used to detect marker-trait associations. Twenty-one loci in a subset of genomic regions (RT-5, 7, 8, 10, 12, 14, 15, 22, 23, 24, 25, 29, 30, and 31) were significantly associated with variation in spawning date. Many of these markers localize to regions with known spawning date quantitative trait loci based on previous studies. An individual assignment analysis was used to test how well the molecular data could be used to assign individuals to their correct spawning group, and markers were given a ranking reflecting their contribution to the accuracy of assignment. The top 15 ranked markers were successful at assigning the majority of females to the correct spawning group based on genotype with an average accuracy of 76 %. The most likely genes that could contribute to these differences in spawning date are discussed. Together, these data indicate that the loci could be incorporated into a selection index with phenotype data to increase the accuracy of selection for spawning date.

Mapping and Expression of Candidate Genes for Development Rate in Rainbow Trout (Oncorhynchus mykiss)

Development rate has important implications for individual fitness and physiology. In salmonid fishes, development rate correlates with many traits later in life, including life-history diversity, growth, and age and size at sexual maturation. In rainbow trout (Oncorhynchus mykiss), a quantitative trait locus for embryonic development rate has been detected on chromosome 5 across populations. However, few candidate genes have been identified within this region. In this study, we use gene mapping, gene expression, and quantitative genetic methods to further identify the genetic basis of embryonic developmental rate in O. mykiss Among the genes located in the region of the major development rate quantitative trait locus (GHR1, Clock1a, Myd118-1, and their paralogs), all were expressed early in embryonic development (fertilization through hatch), but none were differentially expressed between individuals with the fast- or slow-developing alleles for a major embryonic development rate quantitative trait locus. In a follow-up study of migratory and resident rainbow trout from natural populations in Alaska, we found significant additive variation in development rate and, moreover, found associations between development rate and allelic variation in all 3 candidate genes within the quantitative trait locus for embryonic development. The mapping of these genes to this region and associations in multiple populations provide positional candidates for further study of their roles in growth, development, and life-history diversity in this model salmonid.

A dense linkage map for Chinook salmon (Oncorhynchus tshawytscha) reveals variable chromosomal divergence after an ancestral whole genome duplication event

Comparisons between the genomes of salmon species reveal that they underwent extensive chromosomal rearrangements following whole genome duplication that occurred in their lineage 58−63 million years ago. Extant salmonids are diploid, but occasional pairing between homeologous chromosomes exists in males. The consequences of re-diploidization can be characterized by mapping the position of duplicated loci in such species. Linkage maps are also a valuable tool for genome-wide applications such as genome-wide association studies, quantitative trait loci mapping or genome scans. Here, we investigated chromosomal evolution in Chinook salmon (Oncorhynchus tshawytscha) after genome duplication by mapping 7146 restriction-site associated DNA loci in gynogenetic haploid, gynogenetic diploid, and diploid crosses. In the process, we developed a reference database of restriction-site associated DNA loci for Chinook salmon comprising 48528 non-duplicated loci and 6409 known duplicated loci, which will facilitate locus identification and data sharing. We created a very dense linkage map anchored to all 34 chromosomes for the species, and all arms were identified through centromere mapping. The map positions of 799 duplicated loci revealed that homeologous pairs have diverged at different rates following whole genome duplication, and that degree of differentiation along arms was variable. Many of the homeologous pairs with high numbers of duplicated markers appear conserved with other salmon species, suggesting that retention of conserved homeologous pairing in some arms preceded species divergence. As chromosome arms are highly conserved across species, the major resources developed for Chinook salmon in this study are also relevant for other related species.

Centromere localization for Bighead Carp (Aristichthys nobilis) through half-tetrad analysis in diploid gynogenetic families

Gene-centromere (G-C) mapping provides insights into structural and behavioural properties of chromosomes. In this study, G-C mapping using microsatellite markers and meiogynogenetic (meiotic gynogenetic) families were performed in bighead carp (Aristichthys nobilis, 2N = 48), which belongs to Cyprinidae. A total of 218 microsatellites were selected across 24 linkage groups (LGs) of a recently well-defined genetic linkage map for bighead carp, with 151 being heterozygous in at least one of six dams in diploid meiogynogenetic families. After tests for Mendelian segregation in two diploid control families, 103 microsatellites were used for G-C distance calculation in 383 gynogens. The second division segregation frequency (y) was computed through half-tetrad analyses, and the values ranged from 0 to 0.97 (mean 0.40). High G-C recombination frequencies (over 0.667) were observed in 18 (17.5%) of the loci examined, which revealed a low level of chiasma interferences compared with other fishes studied previously. Distribution of G-C distances across LGs ranged from 0 cM to 48.5 cM (mean 20 cM) under the assumption of complete interference. All 24 centromeres were localized according to their closest-related microsatellites at 95% confident intervals. The average distance between centromeres and their closest-linked markers was 6.1 cM with 15 out of 24 LGs having a distance below 5 cM. Based on the centromere positions in this study, we proposed a formula of 24 m/sm+24 t/st chromosomes with 92 arms for bighead carp, which was mostly in accordance with a previously reported karyotype for bighead carp (24 m/sm+24 st). These results of centromere localization provide a basic framework and important resources for genetics and comparative genomics studies in bighead carp and its closely-related cyprinid species.

Comparative genome mapping between Chinook salmon (Oncorhynchus tshawytscha) and rainbow trout (O. mykiss) based on homologous microsatellite loci

Comparative genome mapping can rapidly facilitate the transfer of DNA sequence information from a well-characterized species to one that is less described. Chromosome arm numbers are conserved between members of the teleost family Salmonidae, order Salmoniformes, permitting rapid alignment of large syntenic blocks of DNA between members of the group. However, extensive Robertsonian rearrangements after an ancestral whole-genome duplication event has resulted in different chromosome numbers across Salmonid taxa. In anticipation of the rapid application of genomic data across members of the Pacific salmon genus Oncorhynchus, we mapped the genome of Chinook salmon (O. tshawytscha) by using 361 microsatellite loci and compared linkage groups to those already derived for a well-characterized species rainbow trout (O. mykiss). The Chinook salmon female map length was 1526 cM, the male map 733 cM, and the consensus map between the two sexes was 2206 cM. The average female to male recombination ratio was 5.43 (range 1-42.8 across all pairwise marker comparisons). We detected 34 linkage groups that corresponded with all chromosome arms mapped with homologous loci in rainbow trout and inferred that 16 represented metacentric chromosomes and 18 represented acrocentric chromosomes. Up to 13 chromosomes were conserved between the two species, suggesting that their structure precedes the divergence between Chinook salmon and rainbow trout. However, marker order differed in one of these linkage groups. The remaining linkage group structures reflected independent Robertsonian chromosomal arrangements, possibly after divergence. The putative linkage group homologies presented here are expected to facilitate future DNA sequencing efforts in Chinook salmon.

Assignment of Chinook salmon (Oncorhynchus tshawytscha) linkage groups to specific chromosomes reveals a karyotype with multiple rearrangements of the chromosome arms of rainbow trout (Oncorhynchus mykiss)

The Chinook salmon genetic linkage groups have been assigned to specific chromosomes using fluorescence in situ hybridization with bacterial artificial chromosome probes containing genetic markers mapped to each linkage group in Chinook salmon and rainbow trout. Comparison of the Chinook salmon chromosome map with that of rainbow trout provides strong evidence for conservation of large syntenic blocks in these species, corresponding to entire chromosome arms in the rainbow trout as expected. In almost every case, the markers were found at approximately the same location on the chromosome arm in each species, suggesting conservation of marker order on the chromosome arms of the two species in most cases. Although theoretically a few centric fissions could convert the karyotype of rainbow trout (2N = 58–64) into that of Chinook salmon (2N = 68) or vice versa, our data suggest that chromosome arms underwent multiple centric fissions and subsequent new centric fusions to form the current karyotypes. The morphology of only approximately one-third of the chromosome pairs have been conserved between the two species.

Whole genome comparative analysis of channel catfish (Ictalurus punctatus) with four model fish species

BACKGROUND

Comparative mapping is a powerful tool to study evolution of genomes. It allows transfer of genome information from the well-studied model species to non-model species. Catfish is an economically important aquaculture species in United States. A large amount of genome resources have been developed from catfish including genetic linkage maps, physical maps, BAC end sequences (BES), integrated linkage and physical maps using BES-derived markers, physical map contig-specific sequences, and draft genome sequences. Application of such genome resources should allow comparative analysis at the genome scale with several other model fish species.

RESULTS

In this study, we conducted whole genome comparative analysis between channel catfish and four model fish species with fully sequenced genomes, zebrafish, medaka, stickleback and Tetraodon. A total of 517 Mb draft genome sequences of catfish were anchored to its genetic linkage map, which accounted for 62% of the total draft genome sequences. Based on the location of homologous genes, homologous chromosomes were determined among catfish and the four model fish species. A large number of conserved syntenic blocks were identified. Analysis of the syntenic relationships between catfish and the four model fishes supported that the catfish genome is most similar to the genome of zebrafish.

CONCLUSION

The organization of the catfish genome is similar to that of the four teleost species, zebrafish, medaka, stickleback, and Tetraodon such that homologous chromosomes can be identified. Within each chromosome, extended syntenic blocks were evident, but the conserved syntenies at the chromosome level involve extensive inter-chromosomal and intra-chromosomal rearrangements. This whole genome comparative map should facilitate the whole genome assembly and annotation in catfish, and will be useful for genomic studies of various other fish species.

Chromosome rearrangements, recombination suppression, and limited segregation distortion in hybrids between Yellowstone cutthroat trout (Oncorhynchus clarkii bouvieri) and rainbow trout (O. mykiss)

Background

Introgressive hybridization is an important evolutionary process that can lead to the creation of novel genome structures and thus potentially new genetic variation for selection to act upon. On the other hand, hybridization with introduced species can threaten native species, such as cutthroat trout (Oncorhynchus clarkii) following the introduction of rainbow trout (O. mykiss). Neither the evolutionary consequences nor conservation implications of rainbow trout introgression in cutthroat trout is well understood. Therefore, we generated a genetic linkage map for rainbow-Yellowstone cutthroat trout (O. clarkii bouvieri) hybrids to evaluate genome processes that may help explain how introgression affects hybrid genome evolution.

Results

The hybrid map closely aligned with the rainbow trout map (a cutthroat trout map does not exist), sharing all but one linkage group. This linkage group (RYHyb20) represented a fusion between an acrocentric (Omy28) and a metacentric chromosome (Omy20) in rainbow trout. Additional mapping in Yellowstone cutthroat trout indicated the two rainbow trout homologues were fused in the Yellowstone genome. Variation in the number of hybrid linkage groups (28 or 29) likely depended on a Robertsonian rearrangement polymorphism within the rainbow trout stock. Comparison between the female-merged F₁ map and a female consensus rainbow trout map revealed that introgression suppressed recombination across large genomic regions in 5 hybrid linkage groups. Two of these linkage groups (RYHyb20 and RYHyb25_29) contained confirmed chromosome rearrangements between rainbow and Yellowstone cutthroat trout indicating that rearrangements may suppress recombination. The frequency of allelic and genotypic segregation distortion varied among parents and families, suggesting few incompatibilities exist between rainbow and Yellowstone cutthroat trout genomes.

Conclusions

Chromosome rearrangements suppressed recombination in the hybrids. This result supports several previous findings demonstrating that recombination suppression restricts gene flow between chromosomes that differ by arrangement. Conservation of synteny and map order between the hybrid and rainbow trout maps and minimal segregation distortion in the hybrids suggest rainbow and Yellowstone cutthroat trout genomes freely introgress across chromosomes with similar arrangement. Taken together, these results suggest that rearrangements impede introgression. Recombination suppression across rearrangements could enable large portions of non-recombined chromosomes to persist within admixed populations.

Comparative genomic analysis of catfish linkage group 8 reveals two homologous chromosomes in zebrafish and other teleosts with extensive inter-chromosomal rearrangements

Background

Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study.

Results

With existing genome resources, 287 unique genes were identified in LG8. Comparative genome analysis indicated that most of these 287 genes on catfish LG8 are located on two homologous chromosomes of zebrafish, medaka, stickleback, and three chromosomes of green-spotted pufferfish. Large numbers of conserved syntenies were identified. Detailed analysis of the conserved syntenies in relation to chromosome level similarities revealed extensive inter-chromosomal and intra-chromosomal rearrangements during evolution. Of the 287 genes, 35 genes were found to be duplicated in the catfish genome, with the vast majority of the duplications being interchromosomal.

Conclusions

Comparative genome analysis is a powerful tool even in the absence of a well-assembled whole genome sequence. In spite of sequence stacking due to low resolution of the linkage and physical maps, conserved syntenies can be identified although the exact gene order and orientation are unknown at present. Through chromosome-level comparative analysis, homologous chromosomes among teleosts can be identified. Syntenic analysis should facilitate annotation of the catfish genome, which in turn, should facilitate functional inference of genes based on their orthology.

The sexually dimorphic on the Y-chromosome gene (sdY) is a conserved male-specific Y-chromosome sequence in many salmonids

All salmonid species investigated to date have been characterized with a male heterogametic sex-determination system. However, as these species do not share any Y-chromosome conserved synteny, there remains a debate on whether they share a common master sex-determining gene. In this study, we investigated the extent of conservation and evolution of the rainbow trout (Oncorhynchus mykiss) master sex-determining gene, sdY (sexually dimorphic on the Y-chromosome), in 15 different species of salmonids. We found that the sdY sequence is highly conserved in all salmonids and that sdY is a male-specific Y-chromosome gene in the majority of these species. These findings demonstrate that most salmonids share a conserved sex-determining locus and also strongly suggest that sdY may be this conserved master sex-determining gene. However, in two whitefish species (subfamily Coregoninae), sdY was found both in males and females, suggesting that alternative sex-determination systems may have also evolved in this family. Based on the wide conservation of sdY as a male-specific Y-chromosome gene, efficient and easy molecular sexing techniques can now be developed that will be of great interest for studying these economically and environmentally important species.

Resistance to a rhabdovirus (VHSV) in rainbow trout: identification of a major QTL related to innate mechanisms

Health control is a major issue in animal breeding and a better knowledge of the genetic bases of resistance to diseases is needed in farm animals including fish. The detection of quantitative trait loci (QTL) will help uncovering the genetic architecture of important traits and understanding the mechanisms involved in resistance to pathogens. We report here the detection of QTL for resistance to Viral Haemorrhagic Septicaemia Virus (VHSV), a major threat for European aquaculture industry. Two induced mitogynogenetic doubled haploid F2 rainbow trout (Oncorhynchus mykiss) families were used. These families combined the genome of susceptible and resistant F0 breeders and contained only fully homozygous individuals. For phenotyping, fish survival after an immersion challenge with the virus was recorded, as well as in vitro virus replication on fin explants. A bidirectional selective genotyping strategy identified seven QTL associated to survival. One of those QTL was significant at the genome-wide level and largely explained both survival and viral replication in fin explants in the different families of the design (up to 65% and 49% of phenotypic variance explained respectively). These results evidence the key role of innate defence in resistance to the virus and pave the way for the identification of the gene(s) responsible for resistance. The identification of a major QTL also opens appealing perspectives for selective breeding of fish with improved resistance.

Genetic positioning of centromeres through half-tetrad analysis in gynogenetic diploid families of the Zhikong scallop (Chlamys farreri)

Centromere mapping is a powerful tool for improving linkage maps, investigating crossover events, and understanding chiasma interference during meiosis. Ninety microsatellite markers selected across all linkage groups (LGs) from a previous Chlamys farreri genetic map were studied in three artificially induced meiogynogenetic families for centromere mapping by half-tetrad analysis. Inheritance analyses showed that all 90 microsatellite loci conformed to Mendelian inheritance in the control crosses, while 4.4 % of the microsatellite loci showed segregation departures from an expected 1:1 ratio of two homozygote classes in meiogynogenetic progeny. The second division segregation frequency (y) of the microsatellites ranged from 0.033 to 0.778 with a mean of 0.332, confirming the occurrence of partial chiasma interference in this species. Heterogeneity of y is observed in one of 42 cases in which markers were typed in more than one family, suggesting variation in gene-centromere recombination among families. Centromere location was mostly in accordance with the C. farreri karyotype, but differences in marker order between linkage and centromere maps occurred. Overall, this study makes the genetic linkage map a more complete and informative tool for genomic studies and it will also facilitate future research of the structure and function of the scallop centromeres.

A microsatellite-based linkage map of salt tolerant tilapia (Oreochromis mossambicus x Oreochromis spp.) and mapping of sex-determining loci

Background

Tilapia is the common name for a group of cichlid fishes and is one of the most important aquacultured freshwater food fish. Mozambique tilapia and its hybrids, including red tilapia are main representatives of salt tolerant tilapias. A linkage map is an essential framework for mapping QTL for important traits, positional cloning of genes and understanding of genome evolution.

Results

We constructed a consensus linkage map of Mozambique tilapia and red tilapia using 95 individuals from two F₁ families and 401 microsatellites including 282 EST-derived markers. In addition, we conducted comparative mapping and searched for sex-determining loci on the whole genome. These 401 microsatellites were assigned to 22 linkage groups. The map spanned 1067.6 cM with an average inter-marker distance of 3.3 cM. Comparative mapping between tilapia and stickleback, medaka, pufferfish and zebrafish revealed clear homologous relationships between chromosomes from different species. We found evidence for the fusion of two sets of two independent chromosomes forming two new chromosome pairs, leading to a reduction of 24 chromosome pairs in their ancestor to 22 pairs in tilapias. The XY sex determination locus in Mozambique tilapia was mapped on LG1, and verified in five families containing 549 individuals. The major XY sex determination locus in red tilapia was located on LG22, and verified in two families containing 275 individuals.

Conclusions

A first-generation linkage map of salt tolerant tilapia was constructed using 401 microsatellites. Two separate fusions of two sets of two independent chromosomes may lead to a reduction of 24 chromosome pairs in their ancestor to 22 pairs in tilapias. The XY sex-determining loci from Mozambique tilapia and red tilapia were mapped on LG1 and LG22, respectively. This map provides a useful resource for QTL mapping for important traits and comparative genome studies. The DNA markers linked to the sex-determining loci could be used in the selection of YY males for breeding all-male populations of salt tolerant tilapia, as well as in studies on mechanisms of sex determination in fish.

Genetic basis of differences in myxospore count between whirling disease-resistant and susceptible strains of rainbow trout

We used a quantitative genetics approach and estimated broad sense heritability (h2b) of myxospore count and the number of genes involved in myxospore formation to gain a better understanding of how resistance to Myxobolus cerebralis, the parasite responsible for whirling disease, is inherited in rainbow trout Oncorhynchus mykiss. An M. cerebralis-resistant strain of rainbow trout, the German Rainbow (GR), and a wild, susceptible strain of rainbow trout, the Colorado River Rainbow (CRR), were spawned to create 3 intermediate crossed populations (an F1 cross, F2 intercross, and a B2 backcross between the F1 and the CRR). Within each strain or cross, h2b was estimated from the between-family variance of myxospore counts using full-sibling families. Estimates of h2b and average myxospore counts were lowest in the GR strain, F1 cross, and F2 intercross (h2b = 0.34, 0.42, and 0.34; myxospores fish-1 = 275, 9566, and 45780, respectively), and highest in the B2 backcross and CRR strain (h2b = 0.93 and 0.89; myxospores fish-1 = 97865 and 187595, respectively). Comparison of means and a joint-scaling test suggest that resistance alleles arising from the GR strain are dominant to susceptible alleles from the CRR strain. Resistance was retained in the intermediate crosses but decreased as filial generation number increased (F2) or backcrossing occurred (B2). The estimated number of segregating loci responsible for differences in myxospore count in the parental strains was 9 ± 5. Our results indicate that resistance to M. cerebralis is a heritable trait within these populations and would respond to either artificial selection in hatcheries or natural selection in the wild.

Construction of a high-density microsatellite genetic linkage map and mapping of sexual and growth-related traits in half-smooth tongue sole (Cynoglossus semilaevis)

High-density genetic linkage maps of half-smooth tongue sole were developed with 1007 microsatellite markers, two SCAR markers and an F1 family containing 94. The female map was composed of 828 markers in 21 linkage groups, covering a total of 1447.3 cM, with an average interval 1.83 cM between markers. The male map consisted of 794 markers in 21 linkage groups, spanning 1497.5 cM, with an average interval of 1.96 cM. The female and male maps had 812 and 785 unique positions, respectively. The genome length of half-smooth tongue sole was estimated to be 1527.7 cM for the females and 1582.1 cM for the males. Based on estimations of the map lengths, the female and male maps covered 94.74 and 94.65% of the genome, respectively. The consensus map was composed of 1007 microsatellite markers and two SCAR markers in 21 linkage groups, covering a total of 1624 cM with an average interval of 1.67 cM. Furthermore, 159 sex-linked SSR markers were identified. Five sex-linked microsatellite markers were confirmed in their association with sex in a large number of individuals selected from different families. These sex-linked markers were mapped on the female map LG1f with zero recombination. Two QTLs that were identified for body weight, designated as We-1 and We-2, accounted for 26.39% and 10.60% of the phenotypic variation. Two QTLs for body width, designated Wi-1 and Wi-2, were mapped in LG4f and accounted for 14.33% and 12.83% of the phenotypic variation, respectively. Seven sex-related loci were mapped in LG1f, LG14f and LG1m by CIM, accounting for 12.5–25.2% of the trait variation. The results should prove to be very useful for improving growth traits using molecular MAS.

Construction of high-density genetic linkage maps and mapping of growth-related quantitative trail loci in the Japanese flounder (Paralichthys olivaceus)

High-density genetic linkage maps were constructed for the Japanese flounder (Paralichthys olivaceus). A total of 1624 microsatellite markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 resulted in the mapping of 1487 markers to 24 linkage groups, a result which was consistent with the 24 chromosomes seen in chromosome spreads. The female map was composed of 1257 markers, covering a total of 1663.8 cM with an average interval 1.35 cM between markers. The male map consisted of 1224 markers, spanning 1726.5 cM, with an average interval of 1.44 cM. The genome length in the Japanese flounder was estimated to be 1730.3 cM for the females and 1798.0 cM for the males, a coverage of 96.2% for the female and 96.0% for the male map. The mean recombination at common intervals throughout the genome revealed a slight difference between sexes, i.e. 1.07 times higher in the male than female. High-density genetic linkage maps are very useful for marker-assisted selection (MAS) programs for economically valuable traits in this species and for further evolutionary studies in flatfish and vertebrate species. Furthermore, four quantiative trait loci (QTL) associated with growth traits were mapped on the genetic map. One QTL was identified for body weight on LG 14 f, which explained 14.85% of the total variation of the body weight. Three QTL were identified for body width on LG14f and LG14m, accounting for 16.75%, 13.62% and 13.65% of the total variation in body width, respectively. The additive effects were evident as negative values. There were four QTL for growth traits clustered on LG14, which should prove to be very useful for improving growth traits using molecular MAS.

Genetic architecture of migration-related traits in rainbow and steelhead trout, Oncorhynchus mykiss

Although migration plays a critical role in the evolution and diversification of species, relatively little is known of the genetic architecture underlying this life history in any species. Rainbow and steelhead trout (Oncorhynchus mykiss) naturally segregate for both resident and migratory life-history types, respectively, as do other members of the salmonid family of fishes. Using an experimental cross derived from wild resident rainbow and wild migratory steelhead trout from Southeast Alaska and high throughput restriction-site associated DNA (RAD) tag sequencing, we perform a quantitative trait locus (QTL) analysis to identify the number, position, and relative contribution of genetic effects on a suite of 27 physiological and morphological traits associated with the migratory life history in this species. In total, 37 QTL are localized to 19 unique QTL positions, explaining 4-13.63% of the variation for 19 of the 27 migration-related traits measured. Two chromosomal positions, one on chromosome Omy12 and the other on Omy14 each harbor 7 QTL for migration-related traits, suggesting that these regions could harbor master genetic controls for the migratory life-history tactic in this species. Another QTL region on Omy5 has been implicated in several studies of adaptive life histories within this species and could represent another important locus underlying the migratory life history. We also evaluate whether loci identified in this out-crossed QTL study colocalize to genomic positions previously identified for associations with migration-related traits in a doubled haploid mapping family.

Comparative mapping reveals quantitative trait loci that affect spawning time in coho salmon (Oncorhynchus kisutch)

Spawning time in salmonids is a sex-limited quantitative trait that can be modified by selection. In rainbow trout (Oncorhynchus mykiss), various quantitative trait loci (QTL) that affect the expression of this trait have been discovered. In this study, we describe four microsatellite loci associated with two possible spawning time QTL regions in coho salmon (Oncorhynchus kisutch). The four loci were identified in females from two populations (early and late spawners) produced by divergent selection from the same base population. Three of the loci (OmyFGT34TUF, One2ASC and One19ASC) that were strongly associated with spawning time in coho salmon (p < 0.0002) were previously associated with QTL for the same trait in rainbow trout; a fourth loci (Oki10) with a suggestive association (p = 0.00035) mapped 10 cM from locus OmyFGT34TUF in rainbow trout. The changes in allelic frequency observed after three generations of selection were greater than expected because of genetic drift. This work shows that comparing information from closely-related species is a valid strategy for identifying QTLs for marker-assisted selection in species whose genomes are poorly characterized or lack a saturated genetic map.

An immune-related gene evolved into the master sex-determining gene in rainbow trout, Oncorhynchus mykiss

Since the discovery of Sry in mammals [1, 2], few other master sex-determining genes have been identified in vertebrates [3-7]. To date, all of these genes have been characterized as well-known factors in the sex differentiation pathway, suggesting that the same subset of genes have been repeatedly and independently selected throughout evolution as master sex determinants [8, 9]. Here, we characterized in rainbow trout an unknown gene expressed only in the testis, with a predominant expression during testicular differentiation. This gene is a male-specific genomic sequence that is colocalized along with the sex-determining locus. This gene, named sdY for sexually dimorphic on the Y chromosome, encodes a protein that displays similarity to the C-terminal domain of interferon regulatory factor 9. The targeted inactivation of sdY in males using zinc-finger nuclease induces ovarian differentiation, and the overexpression of sdY in females using additive transgenesis induces testicular differentiation. Together, these results demonstrate that sdY is a novel vertebrate master sex-determining gene not related to any known sex-differentiating gene. These findings highlight an unexpected evolutionary plasticity in vertebrate sex determination through the demonstration that master sex determinants can arise from the de novo evolution of genes that have not been previously implicated in sex differentiation.

Development and characterisation of an expressed sequence tags (EST)-derived single nucleotide polymorphisms (SNPs) resource in rainbow trout

Background

There is considerable interest in developing high-throughput genotyping with single nucleotide polymorphisms (SNPs) for the identification of genes affecting important ecological or economical traits. SNPs are evenly distributed throughout the genome and are likely to be functionally relevant. In rainbow trout, in silico screening of EST databases represents an attractive approach for de novo SNP identification. Nevertheless, EST sequencing errors and assembly of EST paralogous sequences can lead to the identification of false positive SNPs which renders the reliability of EST-derived SNPs relatively low. Further validation of EST-derived SNPs is therefore required. The objective of this work was to assess the quality of and to validate a large number of rainbow trout EST-derived SNPs.

Results

A panel of 1,152 EST-derived SNPs was selected from the INRA Sigenae SNP database and was genotyped in standard and double haploid individuals from several populations using the Illumina GoldenGate BeadXpress assay. High-quality genotyping data were obtained for 958 SNPs representing a genotyping success rate of 83.2 %, out of which, 350 SNPs (36.5 %) were polymorphic in at least one population and were designated as true SNPs. They also proved to be a potential tool to investigate genetic diversity of the species, as the set of SNP successfully sorted individuals into three main groups using STRUCTURE software. Functional annotations revealed 28 non-synonymous SNPs, out of which four substitutions were predicted to affect protein functions. A subset of 223 true SNPs were polymorphic in the two INRA mapping reference families and were integrated into the INRA microsatellite-based linkage map.

Conclusions

Our results represent the first study of EST-derived SNPs validation in rainbow trout, a species whose genome sequences is not yet available. We designed several specific filters in order to improve the genotyping yield. Nevertheless, our selection criteria should be further improved in order to reduce the observed high rate of false positive SNPs which results from the occurrence of whole genome duplications.

A second generation integrated map of the rainbow trout (Oncorhynchus mykiss) genome: analysis of conserved synteny with model fish genomes

DNA fingerprints and end sequences from bacterial artificial chromosomes (BACs) from two new libraries were generated to improve the first generation integrated physical and genetic map of the rainbow trout (Oncorhynchus mykiss) genome. The current version of the physical map is composed of 167,989 clones of which 158,670 are assembled into contigs and 9,319 are singletons. The number of contigs was reduced from 4,173 to 3,220. End sequencing of clones from the new libraries generated a total of 11,958 high quality sequence reads. The end sequences were used to develop 238 new microsatellites of which 42 were added to the genetic map. Conserved synteny between the rainbow trout genome and model fish genomes was analyzed using 188,443 BAC end sequence (BES) reads. The fractions of BES reads with significant BLASTN hits against the zebrafish, medaka, and stickleback genomes were 8.8%, 9.7%, and 10.5%, respectively, while the fractions of significant BLASTX hits against the zebrafish, medaka, and stickleback protein databases were 6.2%, 5.8%, and 5.5%, respectively. The overall number of unique regions of conserved synteny identified through grouping of the rainbow trout BES into fingerprinting contigs was 2,259, 2,229, and 2,203 for stickleback, medaka, and zebrafish, respectively. These numbers are approximately three to five times greater than those we have previously identified using BAC paired ends. Clustering of the conserved synteny analysis results by linkage groups as derived from the integrated physical and genetic map revealed that despite the low sequence homology, large blocks of macrosynteny are conserved between chromosome arms of rainbow trout and the model fish species.

A microsatellite linkage map of striped bass (Morone saxatilis) reveals conserved synteny with the three-spined stickleback (Gasterosteus aculeatus)

The striped bass (Morone saxatilis) and its relatives (genus Morone) are of great importance to fisheries and aquaculture in North America. As part of a collaborative effort to employ molecular genetics technologies in striped bass breeding programs, we previously developed nearly 500 microsatellite markers. The objectives of this study were to construct a microsatellite linkage map of striped bass and to examine conserved synteny between striped bass and three-spined stickleback (Gasterosteus aculeatus). Of 480 microsatellite markers screened for polymorphism, 289 informative markers were identified and used to genotype two half-sib mapping families. Twenty-six linkage groups were assembled, and only two markers remain unlinked. The sex-averaged map spans 1,623.8 cM with an average marker density of 5.78 cM per marker. Among 287 striped bass microsatellite markers assigned to linkage groups, 169 (58.9%) showed homology to sequences on stickleback chromosomes or scaffolds. Comparison between the stickleback genome and the striped bass linkage map revealed conserved synteny between these two species. This is the first linkage map for any of the Morone species. This map will be useful for molecular mapping and marker-assisted selection of genes of interest in striped bass breeding programs. The conserved synteny between striped bass and stickleback will facilitate fine mapping of genome regions of interest and will serve as a new resource for comparative mapping with other Perciform fishes such as European sea bass (Dicentrarchus labrax), gilthead sea bream (Sparus aurata), and tilapia (Oreochromis ssp.).

A synthetic rainbow trout linkage map provides new insights into the salmonid whole genome duplication and the conservation of synteny among teleosts

Background

Rainbow trout is an economically important fish and a suitable experimental organism in many fields of biology including genome evolution, owing to the occurrence of a salmonid specific whole-genome duplication (4^th WGD). Rainbow trout is among some of the most studied teleosts and has benefited from substantial efforts to develop genomic resources (e.g., linkage maps. Here, we first generated a synthetic map by merging segregation data files derived from three independent linkage maps. Then, we used it to evaluate genome conservation between rainbow trout and three teleost models, medaka, stickleback and zebrafish and to further investigate the extent of the 4^th WGD in trout genome.

Results

The INRA linkage map was updated by adding 211 new markers. After standardization of marker names, consistency of marker assignment to linkage groups and marker orders was checked across the three different data sets and only loci showing consistent location over all or almost all of the data sets were kept. This resulted in a synthetic map consisting of 2226 markers and 29 linkage groups spanning over 3600 cM. Blastn searches against medaka, stickleback, and zebrafish genomic databases resulted in 778, 824 and 730 significant hits respectively while blastx searches yielded 505, 513 and 510 significant hits. Homology search results revealed that, for most rainbow trout chromosomes, large syntenic regions encompassing nearly whole chromosome arms have been conserved between rainbow trout and its closest models, medaka and stickleback. Large conserved syntenies were also found between the genomes of rainbow trout and the reconstructed teleost ancestor. These syntenies consolidated the known homeologous affinities between rainbow trout chromosomes due to the 4^th WGD and suggested new ones.

Conclusions

The synthetic map constructed herein further highlights the stability of the teleost genome over long evolutionary time scales. This map can be easily extended by incorporating new data sets and should help future rainbow trout whole genome sequence assembly. Finally, the persistence of large conserved syntenies across teleosts should facilitate the identification of candidate genes through comparative mapping, even if the occurrence of intra-chromosomal micro-rearrangement may hinder the accurate prediction their genomic location.

Isolation and characterization of polymorphic microsatellite loci from Metapenaeopsis barbata using PCR-Based Isolation of Microsatellite Arrays (PIMA)

The red-spot prawn, Metapenaeopsis barbata, is a commercially important, widely distributed demersal species in the Indo-West Pacific Ocean. Overfishing has made its populations decline in the past decade. To study conservation genetics, eight polymorphic microsatellite loci were isolated. Genetic characteristics of the SSR (simple sequence repeat) fingerprints were estimated in 61 individuals from adjacent seas of Taiwan and China. The number of alleles, ranging from 2 to 4, as well as observed and expected heterozygosities in populations, ranging from 0.048 to 0.538, and 0.048 and 0.654, respectively, were detected. No deviation from Hardy–Weinberg expectations was detected at either locus. No significant linkage disequilibrium was detected in locus pairs. The polymorphic microsatellite loci will be useful for investigations of the genetic variation, population structure, and conservation genetics of this species.

A dense SNP-based linkage map for Atlantic salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns

Background

The Atlantic salmon genome is in the process of returning to a diploid state after undergoing a whole genome duplication (WGD) event between 25 and100 million years ago. Existing data on the proportion of paralogous sequence variants (PSVs), multisite variants (MSVs) and other types of complex sequence variation suggest that the rediplodization phase is far from over. The aims of this study were to construct a high density linkage map for Atlantic salmon, to characterize the extent of rediploidization and to improve our understanding of genetic differences between sexes in this species.

Results

A linkage map for Atlantic salmon comprising 29 chromosomes and 5650 single nucleotide polymorphisms (SNPs) was constructed using genotyping data from 3297 fish belonging to 143 families. Of these, 2696 SNPs were generated from ESTs or other gene associated sequences. Homeologous chromosomal regions were identified through the mapping of duplicated SNPs and through the investigation of syntenic relationships between Atlantic salmon and the reference genome sequence of the threespine stickleback (Gasterosteus aculeatus). The sex-specific linkage maps spanned a total of 2402.3 cM in females and 1746.2 cM in males, highlighting a difference in sex specific recombination rate (1.38:1) which is much lower than previously reported in Atlantic salmon. The sexes, however, displayed striking differences in the distribution of recombination sites within linkage groups, with males showing recombination strongly localized to telomeres.

Conclusion

The map presented here represents a valuable resource for addressing important questions of interest to evolution (the process of re-diploidization), aquaculture and salmonid life history biology and not least as a resource to aid the assembly of the forthcoming Atlantic salmon reference genome sequence.

Analysis of BAC-end sequences in rainbow trout: content characterization and assessment of synteny between trout and other fish genomes

BACKGROUND

Rainbow trout (Oncorhynchus mykiss) are cultivated worldwide for aquaculture production and are widely used as a model species to gain knowledge of many aspects of fish biology. The common ancestor of the salmonids experienced a whole genome duplication event, making extant salmonids such as the rainbow trout an excellent model for studying the evolution of tetraploidization and re-diploidization in vertebrates. However, the lack of a reference genome sequence hampers research progress for both academic and applied purposes. In order to enrich the genomic tools already available in this species and provide further insight on the complexity of its genome, we sequenced a large number of rainbow trout BAC-end sequences (BES) and characterized their contents.

RESULTS

A total of 176,485 high quality BES, were generated, representing approximately 4% of the trout genome. BES analyses identified 6,848 simple sequence repeats (SSRs), of which 3,854 had high quality flanking sequences for PCR primers design. The first rainbow trout repeat elements database (INRA RT rep1.0) containing 735 putative repeat elements was developed, and identified almost 59.5% of the BES database in base-pairs as repetitive sequence. Approximately 55% of the BES reads (97,846) had more than 100 base pairs of contiguous non-repetitive sequences. The fractions of the 97,846 non-repetitive trout BES reads that had significant BLASTN hits against the zebrafish, medaka and stickleback genome databases were 15%, 16.2% and 17.9%, respectively, while the fractions of the non-repetitive BES reads that had significant BLASTX hits against the zebrafish, medaka, and stickleback protein databases were 10.7%, 9.5% and 9.5%, respectively. Comparative genomics using paired BAC-ends revealed several regions of conserved synteny across all the fish species analyzed in this study.

CONCLUSIONS

The characterization of BES provided insights on the rainbow trout genome. The discovery of specific repeat elements will facilitate analyses of sequence content (e.g. for SNPs discovery and for transcriptome characterization) and future genome sequence assemblies. The numerous microsatellites will facilitate integration of the linkage and physical maps and serve as valuable resource for fine mapping QTL and positional cloning of genes affecting aquaculture production traits. Furthermore, comparative genomics through BES can be used for identifying positional candidate genes from QTL mapping studies, aid in future assembly of a reference genome sequence and elucidating sequence content and complexity in the rainbow trout genome.

A major effect quantitative trait locus for whirling disease resistance identified in rainbow trout (Oncorhynchus mykiss)

Whirling disease, caused by the pathogen Myxobolus cerebralis, leads to skeletal deformation, neurological impairment and under certain conditions, mortality of juvenile salmonid fishes. The disease has impacted the propagation and survival of many salmonid species over six continents, with particularly negative consequences for rainbow trout. To assess the genetic basis of whirling disease resistance in rainbow trout, genome-wide mapping was initiated using a large outbred F(2) rainbow trout family (n=480) and results were confirmed in three additional outbred F(2) families (n=96 per family). A single quantitative trait locus (QTL) region on chromosome Omy9 was identified in the large mapping family and confirmed in all additional families. This region explains 50-86% of the phenotypic variance across families. Therefore, these data establish that a single QTL region is capable of explaining a large percentage of the phenotypic variance contributing to whirling disease resistance. This is the first genetic region discovered that contributes directly to the whirling disease phenotype and the finding moves the field closer to a mechanistic understanding of resistance to this important disease of salmonid fish.

Detection of QTL with effects on osmoregulation capacities in the rainbow trout (Oncorhynchus mykiss)

Background

There is increasing evidence that the ability to adapt to seawater in teleost fish is modulated by genetic factors. Most studies have involved the comparison of species or strains and little is known about the genetic architecture of the trait. To address this question, we searched for QTL affecting osmoregulation capacities after transfer to saline water in a nonmigratory captive-bred population of rainbow trout.

Results

A QTL design (5 full-sib families, about 200 F2 progeny each) was produced from a cross between F0 grand-parents previously selected during two generations for a high or a low cortisol response after a standardized confinement stress. When fish were about 18 months old (near 204 g body weight), individual progeny were submitted to two successive hyper-osmotic challenges (30 ppt salinity) 14 days apart. Plasma chloride and sodium concentrations were recorded 24 h after each transfer. After the second challenge, fish were sacrificed and a gill index (weight of total gill arches corrected for body weight) was recorded. The genome scan was performed with 196 microsatellites and 85 SNP markers. Unitrait and multiple-trait QTL analyses were carried out on the whole dataset (5 families) through interval mapping methods with the QTLMap software. For post-challenge plasma ion concentrations, significant QTL (P < 0.05) were found on six different linkage groups and highly suggestive ones (P < 0.10) on two additional linkage groups. Most QTL affected concentrations of both chloride and sodium during both challenges, but some were specific to either chloride (2 QTL) or sodium (1 QTL) concentrations. Six QTL (4 significant, 2 suggestive) affecting gill index were discovered. Two were specific to the trait, while the others were also identified as QTL for post-challenge ion concentrations. Altogether, allelic effects were consistent for QTL affecting chloride and sodium concentrations but inconsistent for QTL affecting ion concentrations and gill morphology. There was no systematic lineage effect (grand-parental origin of QTL alleles) on the recorded traits.

Conclusions

For the first time, genomic loci associated with effects on major physiological components of osmotic adaptation to seawater in a nonmigratory fish were revealed. The results pave the way for further deciphering of the complex regulatory mechanisms underlying seawater adaptation and genes involved in osmoregulatory physiology in rainbow trout and other euryhaline fishes.

A first generation integrated map of the rainbow trout genome

Background

Rainbow trout (Oncorhynchus mykiss) are the most-widely cultivated cold freshwater fish in the world and an important model species for many research areas. Coupling great interest in this species as a research model with the need for genetic improvement of aquaculture production efficiency traits justifies the continued development of genomics research resources. Many quantitative trait loci (QTL) have been identified for production and life-history traits in rainbow trout. An integrated physical and genetic map is needed to facilitate fine mapping of QTL and the selection of positional candidate genes for incorporation in marker-assisted selection (MAS) programs for improving rainbow trout aquaculture production.

Results

The first generation integrated map of the rainbow trout genome is composed of 238 BAC contigs anchored to chromosomes of the genetic map. It covers more than 10% of the genome across segments from all 29 chromosomes. Anchoring of 203 contigs to chromosomes of the National Center for Cool and Cold Water Aquaculture (NCCCWA) genetic map was achieved through mapping of 288 genetic markers derived from BAC end sequences (BES), screening of the BAC library with previously mapped markers and matching of SNPs with BES reads. In addition, 35 contigs were anchored to linkage groups of the INRA (French National Institute of Agricultural Research) genetic map through markers that were not informative for linkage analysis in the NCCCWA mapping panel. The ratio of physical to genetic linkage distances varied substantially among chromosomes and BAC contigs with an average of 3,033 Kb/cM.

Conclusions

The integrated map described here provides a framework for a robust composite genome map for rainbow trout. This resource is needed for genomic analyses in this research model and economically important species and will facilitate comparative genome mapping with other salmonids and with model fish species. This resource will also facilitate efforts to assemble a whole-genome reference sequence for rainbow trout.