Construction of a high-density linkage map and fine mapping of QTL for growth in Asian seabass

A Dense Brown Trout (Salmo trutta) Linkage Map Reveals Recent Chromosomal Rearrangements in the Salmo Genus and the Impact of Selection on Linked Neutral Diversity

High-density linkage maps are valuable tools for conservation and eco-evolutionary issues. In salmonids, a complex rediploidization process consecutive to an ancient whole genome duplication event makes linkage maps of prime importance for investigating the evolutionary history of chromosome rearrangements. Here, we developed a high-density consensus linkage map for the brown trout (Salmo trutta), a socioeconomically important species heavily impacted by human activities. A total of 3977 ddRAD markers were mapped and ordered in 40 linkage groups using sex- and lineage-averaged recombination distances obtained from two family crosses. Performing map comparison between S. trutta and its sister species, S. salar, revealed extensive chromosomal rearrangements. Strikingly, all of the fusion and fission events that occurred after the S. salar/S. trutta speciation happened in the Atlantic salmon branch, whereas the brown trout remained closer to the ancestral chromosome structure. Using the strongly conserved synteny within chromosome arms, we aligned the brown trout linkage map to the Atlantic salmon genome sequence to estimate the local recombination rate in S. trutta at 3721 loci. A significant positive correlation between recombination rate and within-population nucleotide diversity (π) was found, indicating that selection constrains variation at linked neutral sites in brown trout. This new high-density linkage map provides a useful genomic resource for future aquaculture, conservation, and eco-evolutionary studies in brown trout.

Development of a 690 K SNP array in catfish and its application for genetic mapping and validation of the reference genome sequence

Single nucleotide polymorphisms (SNPs) are capable of providing the highest level of genome coverage for genomic and genetic analysis because of their abundance and relatively even distribution in the genome. Such a capacity, however, cannot be achieved without an efficient genotyping platform such as SNP arrays. In this work, we developed a high-density SNP array with 690,662 unique SNPs (herein 690 K array) that were relatively evenly distributed across the entire genome, and covered 98.6% of the reference genome sequence. Here we also report linkage mapping using the 690 K array, which allowed mapping of over 250,000 SNPs on the linkage map, the highest marker density among all the constructed linkage maps. These markers were mapped to 29 linkage groups (LGs) with 30,591 unique marker positions. This linkage map anchored 1,602 scaffolds of the reference genome sequence to LGs, accounting for over 97% of the total genome assembly. A total of 1,007 previously unmapped scaffolds were placed to LGs, allowing validation and in few instances correction of the reference genome sequence assembly. This linkage map should serve as a valuable resource for various genetic and genomic analyses, especially for GWAS and QTL mapping for genes associated with economically important traits.

Construction of high-resolution recombination maps in Asian seabass

Background

A high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing.

Results

A high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males.

Conclusions

These high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass.

Electronic supplementary material

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

Genetic load in marine animals: a review

Marine invertebrates and fish are well known for their remarkable genetic diversity, which is commonly explained by large population size and the characteristic dispersive nature of their early, planktonic life history. Other potential sources of diversity in marine animals, such as a higher mutation rate, have been much less considered, though evidence for a high genetic load in marine bivalves has been accumulating for nearly half a century. In this review, I examine evidence for a higher genetic load in marine animals from studies of molecular marker segregation and linkage over the last 40 years, and survey recent work examining mutational load with molecular evolution approaches. Overall, marine animals appear to have higher genetic load than terrestrial animals (higher dn/ds ratios, inbreeding load, and segregation dis`tortion), though results are mixed for marine fish and data are lacking for many marine animal groups. Bivalves (oysters) have the highest loads observed among marine animals, comparable only to long-lived plants; however, more data is needed from other bivalves and more marine invertebrate taxa generally. For oysters, a higher load may be related to a chronically lower effective population size that, in concert with a higher mutational rate, elevate the number of deleterious mutations observed. I suggest that future studies use high-throughput sequencing approaches to examine (1) polymorphism in genome-scale datasets across a wider range of marine animals at the population level and (2) intergenerational mutational changes between parents and offspring in crosses of aquaculture species to quantify mutation rates.

Mapping QTL for Sex and Growth Traits in Salt-Tolerant Tilapia (Oreochromis spp. X O. mossambicus)

In aquaculture, growth and sex are economically important traits. To accelerate genetic improvement in increasing growth in salt-tolerant tilapia, we conducted QTL mapping for growth traits and sex with an F2 family, including 522 offspring and two parents. We used 144 polymorphic microsatellites evenly covering the genome of tilapia to genotype the family. QTL analyses were carried out using interval mapping for all individuals, males and females in the family, respectively. Using all individuals, three suggestive QTL for body weight, body length and body thickness respectively were detected in LG20, LG22 and LG12 and explained 2.4% to 3.1% of phenotypic variance (PV). When considering only males, five QTL for body weight were detected on five LGs, and explained 4.1 to 6.3% of PV. Using only females from the F2 family, three QTL for body weight were detected on LG1, LG6 and LG8, and explained 7.9-14.3% of PV. The QTL for body weight in males and females were located in different LGs, suggesting that in salt-tolerant tilapia, different set of genes 'switches' control the growth in males and females. QTL for sex were mapped on LG1 and LG22, indicating multigene sex determination in the salt-tolerant tilapia. This study provides new insights on the locations and effects of QTL for growth traits and sex, and sets the foundation for fine mapping for future marker-assisted selection for growth and sex in salt-tolerant tilapia aquaculture.

Construction of a high-density genetic map and QTL mapping for pearl quality-related traits in Hyriopsis cumingii

A high-density genetic map is essential for quantitative trait locus (QTL) fine mapping. In this study, 4,508 effective single nucleotide polymorphism markers (detected using specific-locus amplified fragment sequencing) and 475 microsatellites were mapped to 19 linkage groups (LGs) using a family with 157 individuals. The map spanned 2,713 cM, with an average of 259 markers and 79 loci per LG and an average inter-marker distance of 1.81 cM. To identify QTLs for pearl quality traits, 26 putatively significant QTLs were detected for 10 traits, including, three for shell width, seven for body weight, two for shell weight, two for margin mantle weight, five for inner mantle weight, and seven for shell nacre colour. Among them, five QTLs associated with shell nacre colour were mapped to LG17 and explained 19.7% to 22.8% of the trait variation; this suggests that some important genes or loci determine shell nacre colour in LG17. The linkage map and mapped QTLs for shell nacre colour would be useful for improving the quality of Hyriopsis cumingii via marker-assisted selection.

Fine mapping QTL for resistance to VNN disease using a high-density linkage map in Asian seabass

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

Genetic variability, local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass

Genomewide analysis of genetic divergence is critically important in understanding the genetic processes of allopatric speciation. We sequenced RAD tags of 131 Asian seabass individuals of six populations from South-East Asia and Australia/Papua New Guinea. Using 32 433 SNPs, we examined the genetic diversity and patterns of population differentiation across all the populations. We found significant evidence of genetic heterogeneity between South-East Asian and Australian/Papua New Guinean populations. The Australian/Papua New Guinean populations showed a rather lower level of genetic diversity. FST and principal components analysis revealed striking divergence between South-East Asian and Australian/Papua New Guinean populations. Interestingly, no evidence of contemporary gene flow was observed. The demographic history was further tested based on the folded joint site frequency spectrum. The scenario of ancient migration with historical population size changes was suggested to be the best fit model to explain the genetic divergence of Asian seabass between South-East Asia and Australia/Papua New Guinea. This scenario also revealed that Australian/Papua New Guinean populations were founded by ancestors from South-East Asia during mid-Pleistocene and were completely isolated from the ancestral population after the last glacial retreat. We also detected footprints of local selection, which might be related to differential ecological adaptation. The ancient gene flow was examined and deemed likely insufficient to counteract the genetic differentiation caused by genetic drift. The observed genomic pattern of divergence conflicted with the 'genomic islands' scenario. Altogether, Asian seabass have likely been evolving towards allopatric speciation since the split from the ancestral population during mid-Pleistocene.

Construction of high-density genetic linkage map and identification of flowering-time QTLs in orchardgrass using SSRs and SLAF-seq

Orchardgrass (Dactylis glomerata L.) is one of the most economically important perennial, cool-season forage species grown and pastured worldwide. High-density genetic linkage mapping is a valuable and effective method for exploring complex quantitative traits. In this study, we developed 447,177 markers based on SLAF-seq and used them to perform a comparative genomics analysis. Perennial ryegrass sequences were the most similar (5.02%) to orchardgrass sequences. A high-density linkage map of orchardgrass was constructed using 2,467 SLAF markers and 43 SSRs, which were distributed on seven linkage groups spanning 715.77 cM. The average distance between adjacent markers was 0.37 cM. Based on phenotyping in four environments, 11 potentially significant quantitative trait loci (QTLs) for two target traits–heading date (HD) and flowering time (FT)–were identified and positioned on linkage groups LG1, LG3, and LG5. Significant QTLs explained 8.20–27.00% of the total phenotypic variation, with the LOD ranging from 3.85–12.21. Marker167780 and Marker139469 were associated with FT and HD at the same location (Ya’an) over two different years. The utility of SLAF markers for rapid generation of genetic maps and QTL analysis has been demonstrated for heading date and flowering time in a global forage grass.

Construction and Annotation of a High Density SNP Linkage Map of the Atlantic Salmon (Salmo salar) Genome

High density linkage maps are useful tools for fine-scale mapping of quantitative trait loci, and characterization of the recombination landscape of a species’ genome. Genomic resources for Atlantic salmon (Salmo salar) include a well-assembled reference genome, and high density single nucleotide polymorphism (SNP) arrays. Our aim was to create a high density linkage map, and to align it with the reference genome assembly. Over 96,000 SNPs were mapped and ordered on the 29 salmon linkage groups using a pedigreed population comprising 622 fish from 60 nuclear families, all genotyped with the ‘ssalar01’ high density SNP array. The number of SNPs per group showed a high positive correlation with physical chromosome length (r = 0.95). While the order of markers on the genetic and physical maps was generally consistent, areas of discrepancy were identified. Approximately 6.5% of the previously unmapped reference genome sequence was assigned to chromosomes using the linkage map. Male recombination rate was lower than females across the vast majority of the genome, but with a notable peak in subtelomeric regions. Finally, using RNA-Seq data to annotate the reference genome, the mapped SNPs were categorized according to their predicted function, including annotation of ∼2500 putative nonsynonymous variants. The highest density SNP linkage map for any salmonid species has been created, annotated, and integrated with the Atlantic salmon reference genome assembly. This map highlights the marked heterochiasmy of salmon, and provides a useful resource for salmonid genetics and genomics research.

Transcriptome analysis of genes responding to NNV infection in Asian seabass epithelial cells

Asian seabass is an important food fish in Southeast Asia. Viral nervous necrosis (VNN) disease, triggered by nervous necrosis virus (NNV) infection, has caused mass mortality of Asian seabass larvae, resulting in enormous economic losses in the Asian seabass industry. In order to better understand the complex molecular interaction between Asian seabass and NNV, we investigated the transcriptome profiles of Asian seabass epithelial cells, which play an essential role in immune regulation, after NNV infection. Using the next generation sequencing (NGS) technology, we sequenced mRNA from eight samples (6, 12, 24, 48 h post-inoculation) of mock and NNV-infected Asian seabass epithelial cell line, respectively. Clean reads were de novo assembled into a transcriptome consisting of 89026 transcripts with a N50 of 2617 bp. Furthermore, 251 differentially expressed genes (DEGs) in response to NNV infection were identified. Top DEGs include protein asteroid homolog 1-like (ASTE1), receptor-transporting protein 3 (RTP3), heat shock proteins 30 (HSP30) and 70 (HSP70), Viperin, interferon regulatory factor 3 (IRF3) and other genes related to innate immunity. Our data suggest that abundant and diverse genes corresponding to NNV infection. The results of this study could also offer vital information not only for identification of novel genes involved in Asian seabass-NNV interaction, but also for our understanding of the molecular mechanism of Asian seabass' response to viral infection. In addition, 24807 simple sequence repeats (SSRs) were detected in the assembled transcriptome, providing valuable resources for studying genetic variations and accelerating quantitative trait loci (QTL) mapping for disease resistance in Asian seabass in the future.

A high-density genetic map and growth related QTL mapping in bighead carp (Hypophthalmichthys nobilis)

Growth related traits in fish are controlled by quantitative trait loci (QTL), but no QTL for growth have been detected in bighead carp (Hypophthalmichthys nobilis) due to the lack of high-density genetic map. In this study, an ultra-high density genetic map was constructed with 3,121 SNP markers by sequencing 117 individuals in a F₁ family using 2b-RAD technology. The total length of the map was 2341.27 cM, with an average marker interval of 0.75 cM. A high level of genomic synteny between our map and zebrafish was detected. Based on this genetic map, one genome-wide significant and 37 suggestive QTL for five growth-related traits were identified in 6 linkage groups (i.e. LG3, LG11, LG15, LG18, LG19, LG22). The phenotypic variance explained (PVE) by these QTL varied from 15.4% to 38.2%. Marker within the significant QTL region was surrounded by CRP1 and CRP2, which played an important role in muscle cell division. These high-density map and QTL information provided a solid base for QTL fine mapping and comparative genomics in bighead carp.

An ultra-high density linkage map and QTL mapping for sex and growth-related traits of common carp (Cyprinus carpio)

High density genetic linkage maps are essential for QTL fine mapping, comparative genomics and high quality genome sequence assembly. In this study, we constructed a high-density and high-resolution genetic linkage map with 28,194 SNP markers on 14,146 distinct loci for common carp based on high-throughput genotyping with the carp 250 K single nucleotide polymorphism (SNP) array in a mapping family. The genetic length of the consensus map was 10,595.94 cM with an average locus interval of 0.75 cM and an average marker interval of 0.38 cM. Comparative genomic analysis revealed high level of conserved syntenies between common carp and the closely related model species zebrafish and medaka. The genome scaffolds were anchored to the high-density linkage map, spanning 1,357 Mb of common carp reference genome. QTL mapping and association analysis identified 22 QTLs for growth-related traits and 7 QTLs for sex dimorphism. Candidate genes underlying growth-related traits were identified, including important regulators such as KISS2, IGF1, SMTLB, NPFFR1 and CPE. Candidate genes associated with sex dimorphism were also identified including 3KSR and DMRT2b. The high-density and high-resolution genetic linkage map provides an important tool for QTL fine mapping and positional cloning of economically important traits, and improving common carp genome assembly.