Linkage maps for the Pacific abalone (genus Haliotis) based on microsatellite DNA markers

Construction of a high-density linkage map and QTL detection for growth traits in South African abalone (Haliotis midae)

Haliotis midae is one of the most important molluscs in South African commercial aquaculture. In this study, a high-resolution integrated linkage map was constructed, and QTL identified using 2b-RADseq for genotyping SNPs in three families. The final integrated linkage map was composed by merging the individual family maps, resulting in 3290 informative SNPs mapping to 18 linkage groups, conforming to the known haploid chromosome number for H. midae. The total map spanned 1798.25 cM with an average marker interval of 0.55 cM, representing a genome coverage of 98.76%. QTL analysis, across all three families, resulted in a total of five QTL identified for growth-related traits, shell width, shell length, and total body weight. For shell width and total body weight, one QTL was identified for each trait respectively, whilst three QTL were identified for shell length. The identified QTL respectively explained between 7.20% and 11.40% of the observed phenotypic variance. All three traits were significantly correlated (r = 0.862-0.970; p < 0.01) and shared overlapping QTL. The QTL for growth traits were mapped back to the H. midae draft genome and BLAST searches revealed the identity of candidate genes, such as egf-1, megf10, megf6, tnx, sevp1, kcp, notch1, and scube2 with possible functional roles in H. midae growth. The constructed high-density linkage map and mapped QTL have given valuable insights regarding the genetic architecture of growth-related traits and will be important genetic resources for marker-assisted selection. It remains, however, important to validate causal variants through linkage disequilibrium fine mapping in future.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress

Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.

Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates

How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.

Genetic determination of sex and shell color in the Pacific abalone Haliotis discus hannai revealed by an integrated linkage map

Integrated linkage maps for each sex have been constructed for the Pacific abalone Haliotis discus hannai using three F₁ mapping families based on co-dominant markers. A total of 273 markers were placed on the female map, spanning 927.3 cM with an average interval of 3.64 cM, whereas 277 markers were mapped on the male map, covering 727.0 cM with an average spacing of 2.80 cM. Both female and male maps consisted of 18 linkage groups, corresponding well with the number of chromosomes. Furthermore, the sex-determining locus and the green/orange shell color controlling locus were mapped to the linkage group 3 (LG3) and LG9 respectively. A marker completely linked to phenotypic sex was identified, and the sex determination system was further concluded as paternal heterogametic (males XY and females XX). Based on the segregation ratio of the shell color in the progeny, a simple recessive model of epistasis was proposed to explain the distribution of different color morphs (green, orange and blue): the recessive allele determining orange type masks the effect of the locus controlling green and blue types, whereas the dominant allele at the green/orange locus permits the expression of green and blue types controlled by another locus. The current consensus map provides a useful framework for genetic studies in the Pacific abalone. Mapping of the sex-determining locus and the shell color-controlling locus leads to further understanding of the mechanisms underlying these important traits.

Alternative Splicing Profile and Sex-Preferential Gene Expression in the Female and Male Pacific Abalone Haliotis discus hannai

In order to characterize the female or male transcriptome of the Pacific abalone and further increase genomic resources, we sequenced the mRNA of full-length complementary DNA (cDNA) libraries derived from pooled tissues of female and male Haliotis discus hannai by employing the Iso-Seq protocol of the PacBio RSII platform. We successfully assembled whole full-length cDNA sequences and constructed a transcriptome database that included isoform information. After clustering, a total of 15,110 and 12,145 genes that coded for proteins were identified in female and male abalones, respectively. A total of 13,057 putative orthologs were retained from each transcriptome in abalones. Overall Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyzed in each database showed a similar composition between sexes. In addition, a total of 519 and 391 isoforms were genome-widely identified with at least two isoforms from female and male transcriptome databases. We found that the number of isoforms and their alternatively spliced patterns are variable and sex-dependent. This information represents the first significant contribution to sex-preferential genomic resources of the Pacific abalone. The availability of whole female and male transcriptome database and their isoform information will be useful to improve our understanding of molecular responses and also for the analysis of population dynamics in the Pacific abalone.

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.

Construction of the BAC Library of Small Abalone (Haliotis diversicolor) for Gene Screening and Genome Characterization

The small abalone (Haliotis diversicolor) is one of the most important aquaculture species in East Asia. To facilitate gene cloning and characterization, genome analysis, and genetic breeding of it, we constructed a large-insert bacterial artificial chromosome (BAC) library, which is an important genetic tool for advanced genetics and genomics research. The small abalone BAC library includes 92,610 clones with an average insert size of 120 Kb, equivalent to approximately 7.6× of the small abalone genome. We set up three-dimensional pools and super pools of 18,432 BAC clones for target gene screening using PCR method. To assess the approach, we screened 12 target genes in these 18,432 BAC clones and identified 16 positive BAC clones. Eight positive BAC clones were then sequenced and assembled with the next generation sequencing platform. The assembled contigs representing these 8 BAC clones spanned 928 Kb of the small abalone genome, providing the first batch of genome sequences for genome evaluation and characterization. The average GC content of small abalone genome was estimated as 40.33%. A total of 21 protein-coding genes, including 7 target genes, were annotated into the 8 BACs, which proved the feasibility of PCR screening approach with three-dimensional pools in small abalone BAC library. One hundred fifty microsatellite loci were also identified from the sequences for marker development in the future. The BAC library and clone pools provided valuable resources and tools for genetic breeding and conservation of H. diversicolor.

Cytogenetic analysis and chromosomal characteristics of the polymorphic 18S rDNA of Haliotis discus hannai from Fujian, China

We report on novel chromosomal characteristics of Haliotis discus hannai from a breeding population at Fujian, China. The karyotypes of H. discus hannai we obtained from an abalone farm include a common type 2n = 36 = 10M + 8SM (82%) and two rare types 2n = 36 = 11M + 7SM (14%) and 2n = 36 = 10M + 7SM + 1ST (4%). The results of silver staining showed that the NORs of H. discus hannai were usually located terminally on the long arms of chromosome pairs 14 and 17, NORs were also sometimes located terminally on the short arms of other chromosomes, either metacentric or submetacentric pairs. The number of Ag-nucleoli ranged from 2 to 8, and the mean number was 3.61 ± 0.93. Among the scored interphase cells, 41% had 3 detectable nucleoli and 37% had 4 nucleoli. The 18S rDNA FISH result is the first report of the location of 18S rDNA genes in H. discus hannai. The 18S rDNA locations were highly polymorphic in this species. Copies of the gene were observed in the terminal of long or/and short arms of submetacentric or/and metacentric chromosomes. Using FISH with probe for vertebrate-like telomeric sequences (CCCTAA)₃ displayed positive green FITC signals at telomere regions of all analyzed chromosome types. We found about 7% of chromosomes had breaks in prophase. A special form of nucleolus not previously described from H. discus hannai was observed in some interphase cells. It consists of many small silver-stained nucleoli gathered together to form a larger nucleolus and may correspond to prenucleolar bodies.

Genome wide SNP identification in chickpea for use in development of a high density genetic map and improvement of chickpea reference genome assembly

BACKGROUND

In the whole genome sequencing, genetic map provides an essential framework for accurate and efficient genome assembly and validation. The main objectives of this study were to develop a high-density genetic map using RAD-Seq (Restriction-site Associated DNA Sequencing) genotyping-by-sequencing (RAD-Seq GBS) and Illumina GoldenGate assays, and to examine the alignment of the current map with the kabuli chickpea genome assembly.

RESULTS

Genic single nucleotide polymorphisms (SNPs) totaling 51,632 SNPs were identified by 454 transcriptome sequencing of Cicer arietinum and Cicer reticulatum genotypes. Subsequently, an Illumina GoldenGate assay for 1,536 SNPs was developed. A total of 1,519 SNPs were successfully assayed across 92 recombinant inbred lines (RILs), of which 761 SNPs were polymorphic between the two parents. In addition, the next generation sequencing (NGS)-based GBS was applied to the same population generating 29,464 high quality SNPs. These SNPs were clustered into 626 recombination bins based on common segregation patterns. Data from the two approaches were used for the construction of a genetic map using a population derived from an intraspecific cross. The map consisted of 1,336 SNPs including 604 RAD recombination bins and 732 SNPs from Illumina GoldenGate assay. The map covered 653 cM of the chickpea genome with an average distance between adjacent markers of 0.5 cM. To date, this is the most extensive genetic map of chickpea using an intraspecific population. The alignment of the map with the CDC Frontier genome assembly revealed an overall conserved marker order; however, a few local inconsistencies within the Cicer arietinum pseudochromosome 1 (Ca1), Ca5 and Ca8 were detected. The map enabled the alignment of 215 unplaced scaffolds from the CDC Frontier draft genome assembly. The alignment also revealed varying degrees of recombination rates and hotspots across the chickpea genome.

CONCLUSIONS

A high-density genetic map using RAD-Seq GBS and Illumina GoldenGate assay was developed and aligned with the existing kabuli chickpea draft genome sequence. The analysis revealed an overall conserved marker order, although some localized inversions between draft genome assembly and the genetic map were detected. The current analysis provides an insight of the recombination rates and hotspots across the chickpea genome.

Amplified fragment length polymorphism analysis to assess crossover interference and homozygosity in gynogenetic diploid Pacific abalone (Haliotis discus hannai)

Recombination analysis in gynogenetic diploids is a powerful tool for assessing the degree of inbreeding, investigating crossover events and understanding chiasma interference during meiosis. To estimate the marker-centromere recombination rate, the inheritance pattern of 654 amplified fragment length polymorphism (AFLP) markers was examined in the 72-h veliger larvae of two meiogynogenetic diploid families in the Pacific abalone (Haliotis discus hannai). The second-division segregation frequency (y) of the AFLP loci ranged from 0.00 to 0.96, with 23.9% of loci showing y-values higher than 0.67, evidencing the existence of interference. The average recombination frequency across the 654 AFLP loci was 0.45, allowing estimation of the fixation index of 0.55, indicating that meiotic gynogenesis could provide an effective means of rapid inbreeding in the Pacific abalone. The AFLP loci have a small proportion (4.4%) of y-values greater than 0.90, suggesting that a relatively low or intermediate degree of chiasma interference occurred in the abalone chromosomes. The information obtained in this study will enhance our understanding of the abalone genome and will be useful for genetic studies in the species.

Genetic map construction and quantitative trait locus (QTL) detection of growth-related traits in Litopenaeus vannamei for selective breeding applications

Growth is a priority trait from the point of view of genetic improvement. Molecular markers linked to quantitative trait loci (QTL) have been regarded as useful for marker-assisted selection (MAS) in complex traits as growth. Using an intermediate F2 cross of slow and fast growth parents, a genetic linkage map of Pacific whiteleg shrimp, Litopenaeusvannamei, based on amplified fragment length polymorphisms (AFLP) and simple sequence repeats (SSR) markers was constructed. Meanwhile, QTL analysis was performed for growth-related traits. The linkage map consisted of 451 marker loci (429 AFLPs and 22 SSRs) which formed 49 linkage groups with an average marker space of 7.6 cM; they spanned a total length of 3627.6 cM, covering 79.50% of estimated genome size. 14 QTLs were identified for growth-related traits, including three QTLs for body weight (BW), total length (TL) and partial carapace length (PCL), two QTLs for body length (BL), one QTL for first abdominal segment depth (FASD), third abdominal segment depth (TASD) and first abdominal segment width (FASW), which explained 2.62 to 61.42% of phenotypic variation. Moreover, comparison of linkage maps between L. vannamei and Penaeusjaponicus was applied, providing a new insight into the genetic base of QTL affecting the growth-related traits. The new results will be useful for conducting MAS breeding schemes in L. vannamei .

An integrated genetic linkage map for white clover (Trifolium repens L.) with alignment to Medicago

BACKGROUND

White clover (Trifolium repens L.) is a temperate forage legume with an allotetraploid genome (2n=4×=32) estimated at 1093 Mb. Several linkage maps of various sizes, marker sources and completeness are available, however, no integrated map and marker set has explored consistency of linkage analysis among unrelated mapping populations. Such integrative analysis requires tools for homoeologue matching among populations. Development of these tools provides for a consistent framework map of the white clover genome, and facilitates in silico alignment with the model forage legume, Medicago truncatula.

RESULTS

This is the first report of integration of independent linkage maps in white clover, and adds to the literature on methyl filtered GeneThresher®-derived microsatellite (simple sequence repeat; SSR) markers for linkage mapping. Gene-targeted SSR markers were discovered in a GeneThresher® (TrGT) methyl-filtered database of 364,539 sequences, which yielded 15,647 SSR arrays. Primers were designed for 4,038 arrays and of these, 465 TrGT-SSR markers were used for parental consensus genetic linkage analysis in an F1 mapping population (MP2). This was merged with an EST-SSR consensus genetic map of an independent population (MP1), using markers to match homoeologues and develop a multi-population integrated map of the white clover genome. This integrated map (IM) includes 1109 loci based on 804 SSRs over 1274 cM, covering 97% of the genome at a moderate density of one locus per 1.2 cM. Eighteen candidate genes and one morphological marker were also placed on the IM. Despite being derived from disparate populations and marker sources, the component maps and the derived IM had consistent representations of the white clover genome for marker order and genetic length. In silico analysis at an E-value threshold of 1e-20 revealed substantial co-linearity with the Medicago truncatula genome, and indicates a translocation between T. repens groups 2 and 6 relative to M. truncatula.

CONCLUSIONS

This integrated genetic linkage analysis provides a consistent and comprehensive linkage analysis of the white clover genome, with alignment to a model forage legume. Associated marker locus information, particularly the homoeologue-specific markers, offers a new resource for forage legume research to enable genetic analysis and improvement of this forage and grassland species.

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.

An Expressed Sequence Tag (EST)-enriched genetic map of turbot (Scophthalmus maximus): a useful framework for comparative genomics across model and farmed teleosts

Background

The turbot (Scophthalmus maximus) is a relevant species in European aquaculture. The small turbot genome provides a source for genomics strategies to use in order to understand the genetic basis of productive traits, particularly those related to sex, growth and pathogen resistance. Genetic maps represent essential genomic screening tools allowing to localize quantitative trait loci (QTL) and to identify candidate genes through comparative mapping. This information is the backbone to develop marker-assisted selection (MAS) programs in aquaculture. Expressed sequenced tag (EST) resources have largely increased in turbot, thus supplying numerous type I markers suitable for extending the previous linkage map, which was mostly based on anonymous loci. The aim of this study was to construct a higher-resolution turbot genetic map using EST-linked markers, which will turn out to be useful for comparative mapping studies.

Results

A consensus gene-enriched genetic map of the turbot was constructed using 463 SNP and microsatellite markers in nine reference families. This map contains 438 markers, 180 EST-linked, clustered at 24 linkage groups. Linkage and comparative genomics evidences suggested additional linkage group fusions toward the consolidation of turbot map according to karyotype information. The linkage map showed a total length of 1402.7 cM with low average intermarker distance (3.7 cM; ~2 Mb). A global 1.6:1 female-to-male recombination frequency (RF) ratio was observed, although largely variable among linkage groups and chromosome regions. Comparative sequence analysis revealed large macrosyntenic patterns against model teleost genomes, significant hits decreasing from stickleback (54%) to zebrafish (20%). Comparative mapping supported particular chromosome rearrangements within Acanthopterygii and aided to assign unallocated markers to specific turbot linkage groups.

Conclusions

The new gene-enriched high-resolution turbot map represents a useful genomic tool for QTL identification, positional cloning strategies, and future genome assembling. This map showed large synteny conservation against model teleost genomes. Comparative genomics and data mining from landmarks will provide straightforward access to candidate genes, which will be the basis for genetic breeding programs and evolutionary studies in this species.

Centromere mapping in the Pacific abalone (Haliotis discus hannai) through half-tetrad analysis in gynogenetic diploid families

Centromere mapping is an essential prerequisite for our understanding of the composition and structure of genomes. For centromere mapping, in two meiogynogenetic families of the Pacific abalone (Haliotis discus hannai), we screened 97 microsatellite markers that cover all linkage groups from a currently available abalone linkage map. Microsatellite analysis showed that no unique paternal allele was found in all gynogenetic progeny, which confirmed 100% success of induction of gynogenesis. In the control crosses, all 97 microsatellite loci were compatible with Mendelian inheritance, while in meiogynogenetic progeny, 5.2% of the microsatellite loci showed segregation distortions from an expected 1:1 ratio of two homozygote classes. The second division segregation frequency of the microsatellites ranged from 0.037 to 0.950 with a mean of 0.399, indicating the existence of interference. Heterogeneity among linkage groups in the crossover distribution was observed. Centromere location was mostly in accordance with the abalone karyotype, but differences in marker order between linkage and centromere maps occurred. Information on the positions of centromeres in relation to the microsatellite loci will represent a contribution towards assembly of genetic maps in the commercially important abalone species.

Construction of an integrated map of Haliotis diversicolor using microsatellite markers

Small abalone, Haliotis diversicolor, is naturally distributed along the coastal waters of East Asia from Japan to the Philippines. It is an economically important maricultured species in southern China and Taiwan. Genetic linkage maps for small abalone were constructed using a total of 308 simple sequence repeat markers including 297 novel markers. Segregation data on 96 progeny were genotyped using a pseudo-testcross strategy. Sixteen linkage groups were identified in both female and male maps, consistent with the haploid chromosome number. The female linkage map covered 758.3 cM, with an average interval of 5.2 cM. The male linkage map spanned a total genetic distance of 676.2 cM, with an average interval of 4.5 cM. An integrated map was constructed by incorporating the homologous parental linkage groups, resulting in 16 linkage groups with a total of 762.1 cM. Genome coverage of the integrated linkage map was approximately 80.7%. The genetic linkage maps of small abalone may facilitate marker-assisted selection and quantitative trait loci mapping.

Construction of a genetic linkage map for silver carp (Hypophthalmichthys molitrix)

For silver carp (Hypophthalmichthys molitrix), a combined microsatellite (or simple sequence repeat) and amplified fragment length polymorphism (AFLP) sex average linkage map was constructed. A total of 483 markers (245 microsatellites and 238 AFLPs) were assigned to 33 linkage groups. The map spanned 1352.2 cM, covering 86.4% of the estimated genome size of silver carp. The maximum and average spaces between 420 loci were 21.5 cM and 3.2 cM, respectively. The length of linkage groups ranged from 3.6 cM to 98.5 cM with an average of 41.0 cM. The number of markers per group varied from 2 to 44 with an average of 14.6. The AFLP markers significantly improved the integrity of microsatellite-based linkage groups and increased the genome coverage and marker evenness. A genome-wide recombination suppression was observed in male. In an extreme case, six microsatellites co-segregated in male, but spanned a 45.1 cM region in female.

A second generation genetic linkage map of Japanese flounder (Paralichthys olivaceus)

BACKGROUND

Japanese flounder (Paralichthys olivaceus) is one of the most economically important marine species in Northeast Asia. Information on genetic markers associated with quantitative trait loci (QTL) can be used in breeding programs to identify and select individuals carrying desired traits. Commercial production of Japanese flounder could be increased by developing disease-resistant fish and improving commercially important traits. Previous maps have been constructed with AFLP markers and a limited number of microsatellite markers. In this study, improved genetic linkage maps are presented. In contrast with previous studies, these maps were built mainly with a large number of codominant markers so they can potentially be used to analyze different families and populations.

RESULTS

Sex-specific genetic linkage maps were constructed for the Japanese flounder including a total of 1,375 markers [1,268 microsatellites, 105 single nucleotide polymorphisms (SNPs) and two genes]; 1,167 markers are linked to the male map and 1,067 markers are linked to the female map. The lengths of the male and female maps are 1,147.7 cM and 833.8 cM, respectively. Based on estimations of map lengths, the female and male maps covered 79 and 82% of the genome, respectively. Recombination ratio in the new maps revealed F:M of 1:0.7. All linkage groups in the maps presented large differences in the location of sex-specific recombination hot-spots.

CONCLUSIONS

The improved genetic linkage maps are very useful for QTL analyses and marker-assisted selection (MAS) breeding programs for economically important traits in Japanese flounder. In addition, SNP flanking sequences were blasted against Tetraodon nigroviridis (puffer fish) and Danio rerio (zebrafish), and synteny analysis has been carried out. The ability to detect synteny among species or genera based on homology analysis of SNP flanking sequences may provide opportunities to complement initial QTL experiments with candidate gene approaches from homologous chromosomal locations identified in related model organisms.

Development of gene-targeted SNP markers for genomic mapping in Pacific abalone Haliotis discus hannai Ino

Useful and novel DNA markers are needed for aquaculture genetics and breeding. In this study, we report the discovery and development of gene-targeted single nucleotide polymorphisms (SNPs) for genomic mapping in the Pacific abalone Haliotis discus hannai Ino. Single EST or EST-contigs from 66 genes that had positive BLASTx matches (E-value ≤ 1e-8) were used for polymerase chain reaction (PCR) amplification. PCR products from the two parents of one mapping family were directly sequenced, and 83 SNP loci were found from 17 genes. Allele-specific PCR (AS-PCR) was developed and optimized for genotyping of 11 SNP loci in 120 progeny of the mapping family. Nine of the loci conformed to the expected Mendelian ratio of 1:1 based on the χ(2) test (P > 0.05) and could potentially be used for linkage map construction. Our data also indicate that the sequencing of two parents may be a practical strategy for the discovery of informative SNPs for linkage mapping in a particular mapping population.

Construction of a high-resolution genetic linkage map and comparative genome analysis for the reef-building coral Acropora millepora

A high-resolution genetic linkage map for the coral Acropora millepora is constructed and compared with other metazoan genomes, revealing syntenic blocks.

Background

Worldwide, coral reefs are in decline due to a range of anthropogenic disturbances, and are now also under threat from global climate change. Virtually nothing is currently known about the genetic factors that might determine whether corals adapt to the changing climate or continue to decline. Quantitative genetics studies aiming to identify the adaptively important genomic loci will require a high-resolution genetic linkage map. The phylogenetic position of corals also suggests important applications for a coral genetic map in studies of ancestral metazoan genome architecture.

Results

We constructed a high-resolution genetic linkage map for the reef-building coral Acropora millepora, the first genetic map reported for any coral, or any non-Bilaterian animal. More than 500 single nucleotide polymorphism (SNP) markers were developed, most of which are transferable in populations from Orpheus Island and Great Keppel Island. The map contains 429 markers (393 gene-based SNPs and 36 microsatellites) distributed in 14 linkage groups, and spans 1,493 cM with an average marker interval of 3.4 cM. Sex differences in recombination were observed in a few linkage groups, which may be caused by haploid selection. Comparison of the coral map with other metazoan genomes (human, nematode, fly, anemone and placozoan) revealed synteny regions.

Conclusions

Our study develops a framework that will be essential for future studies of adaptation in coral and it also provides an important resource for future genome sequence assembly and for comparative genomics studies on the evolution of metazoan genome structure.

An integrated genetic linkage map for silkworms with three parental combinations and its application to the mapping of single genes and QTL

Background

Bombyx mori, the domesticated silkworm, is a well-studied model insect with great economic and scientific significance. Although more than 400 mutations have been described in silkworms, most have not been identified, especially those affecting economically-important traits. Simple sequence repeats (SSRs) are effective and economical tools for mapping traits and genetic improvement. The current SSR linkage map is of low density and contains few polymorphisms. The purpose of this work was to develop a dense and informative linkage map that would assist in the preliminary mapping and dissection of quantitative trait loci (QTL) in a variety of silkworm strains.

Results

Through an analysis of > 50,000 genotypes across new mapping populations, we constructed two new linkage maps covering 27 assigned chromosomes and merged the data with previously reported data sets. The integrated consensus map contains 692 unique SSR sites, improving the density from 6.3 cM in the previous map to 4.8 cM. We also developed 497 confirmed neighboring markers for corresponding low-polymorphism sites, with 244 having polymorphisms. Large-scale statistics on the SSR type were suggestive of highly efficient markers, based upon which we searched 16,462 available genomic scaffolds for SSR loci. With the newly constructed map, we mapped single-gene traits, the QTL of filaments, and a number of ribosomal protein genes.

Conclusion

The integrated map produced in this study is a highly efficient genetic tool for the high-throughput mapping of single genes and QTL. Compared to previous maps, the current map offers a greater number of markers and polymorphisms; thus, it may be used as a resource for marker-assisted breeding.

Coevolution of interacting fertilization proteins

Reproductive proteins are among the fastest evolving in the proteome, often due to the consequences of positive selection, and their rapid evolution is frequently attributed to a coevolutionary process between interacting female and male proteins. Such a process could leave characteristic signatures at coevolving genes. One signature of coevolution, predicted by sexual selection theory, is an association of alleles between the two genes. Another predicted signature is a correlation of evolutionary rates during divergence due to compensatory evolution. We studied female-male coevolution in the abalone by resequencing sperm lysin and its interacting egg coat protein, VERL, in populations of two species. As predicted, we found intergenic linkage disequilibrium between lysin and VERL, despite our demonstration that they are not physically linked. This finding supports a central prediction of sexual selection using actual genotypes, that of an association between a male trait and its female preference locus. We also created a novel likelihood method to show that lysin and VERL have experienced correlated rates of evolution. These two signatures of coevolution can provide statistical rigor to hypotheses of coevolution and could be exploited for identifying coevolving proteins a priori. We also present polymorphism-based evidence for positive selection and implicate recent selective events at the specific structural regions of lysin and VERL responsible for their species-specific interaction. Finally, we observed deep subdivision between VERL alleles in one species, which matches a theoretical prediction of sexual conflict. Thus, abalone fertilization proteins illustrate how coevolution can lead to reproductive barriers and potentially drive speciation.

Construction of genetic linkage maps and comparative genome analysis of catfish using gene-associated markers

A genetic linkage map of the channel catfish genome (N=29) was constructed using EST-based microsatellite and single nucleotide polymorphism (SNP) markers in an interspecific reference family. A total of 413 microsatellites and 125 SNP markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 allowed mapping of 331 markers (259 microsatellites and 72 SNPs) to 29 linkage groups. Each linkage group contained 3-18 markers. The largest linkage group contained 18 markers and spanned 131.2 cM, while the smallest linkage group contained 14 markers and spanned only 7.9 cM. The linkage map covered a genetic distance of 1811 cM with an average marker interval of 6.0 cM. Sex-specific maps were also constructed; the recombination rate for females was 1.6 times higher than that for males. Putative conserved syntenies between catfish and zebrafish, medaka, and Tetraodon were established, but the overall levels of genome rearrangements were high among the teleost genomes. This study represents a first-generation linkage map constructed by using EST-derived microsatellites and SNPs, laying a framework for large-scale comparative genome analysis in catfish. The conserved syntenies identified here between the catfish and the three model fish species should facilitate structural genome analysis and evolutionary studies, but more importantly should facilitate functional inference of catfish genes. Given that determination of gene functions is difficult in nonmodel species such as catfish, functional genome analysis will have to rely heavily on the establishment of orthologies from model species.

A first-generation microsatellite-based genetic linkage map of the Siberian jay (Perisoreus infaustus): insights into avian genome evolution

BACKGROUND

Genomic resources for the majority of free-living vertebrates of ecological and evolutionary importance are scarce. Therefore, linkage maps with high-density genome coverage are needed for progress in genomics of wild species. The Siberian jay (Perisoreus infaustus; Corvidae) is a passerine bird which has been subject to lots of research in the areas of ecology and evolutionary biology. Knowledge of its genome structure and organization is required to advance our understanding of the genetic basis of ecologically important traits in this species, as well as to provide insights into avian genome evolution.

RESULTS

We describe the first genetic linkage map of Siberian jay constructed using 117 microsatellites and a mapping pedigree of 349 animals representing five families from a natural population breeding in western Finland from the years 1975 to 2006. Markers were resolved into nine autosomal and a Z-chromosome-specific linkage group, 10 markers remaining unlinked. The best-position map with the most likely positions of all significantly linked loci had a total sex-average size of 862.8 cM, with an average interval distance of 9.69 cM. The female map covered 988.4 cM, whereas the male map covered only 774 cM. The Z-chromosome linkage group comprised six markers, three pseudoautosomal and three sex-specific loci, and spanned 10.6 cM in females and 48.9 cM in males. Eighty-one of the mapped loci could be ordered on a framework map with odds of >1000:1 covering a total size of 809.6 cM in females and 694.2 cM in males. Significant sex specific distortions towards reduced male recombination rates were revealed in the entire best-position map as well as within two autosomal linkage groups. Comparative mapping between Siberian jay and chicken anchored 22 homologous loci on 6 different linkage groups corresponding to chicken chromosomes Gga1, 2, 3, 4, 5, and Z. Quite a few cases of intra-chromosomal rearrangements within the autosomes and three cases of inter-chromosomal rearrangement between the Siberian jay autosomal linkage groups (LG1, LG2 and LG3) and the chicken sex chromosome GgaZ were observed, suggesting a conserved synteny, but changes in marker order, within autosomes during about 100 million years of avian evolution.

CONCLUSION

The constructed linkage map represents a valuable resource for intraspecific genomics of Siberian jay, as well as for avian comparative genomic studies. Apart from providing novel insights into sex-specific recombination rates and patterns, the described maps - from a previously genomically uncharacterized superfamily (Corvidae) of passerine birds - provide new insights into avian genome evolution. In combination with high-resolution data on quantitative trait variability from the study population, they also provide a foundation for QTL-mapping studies.

Detection of QTL for growth rate in the blacklip abalone (Haliotis rubra Leach) using selective DNA pooling

The objective of this study was to identify QTL for growth rate in the blacklip abalone Haliotis rubra using selective DNA pooling. Three full-sibling families of H. rubra derived from crosses of wild broodstock were used. DNA was extracted from the largest and smallest 10% of progeny and combined into two pools for each phenotypic tail. The DNA pools were typed with 139 microsatellites, and markers showing significant differences between the peak height ratios of alleles inherited from the parents were individually genotyped and analysed by interval mapping. A strong correlation (r = 0.94, P < 0.001) was found between the t-values from the analysis of pools and the t-values from the analysis of individual genotypes. Based on the interval mapping analysis, QTL were detected on nine linkage groups at a chromosome-wide P < 0.01 and one linkage group at a chromosome-wide P < 0.05. The study demonstrated that selective DNA pooling is efficient and effective as a first-pass screen for the discovery of QTL in an aquaculture species.