An extended anchored linkage map and virtual mapping for the American mink genome based on homology to human and dog

Population genomics of American mink using genotype data

Understanding the genetic structure of the target population is critically important to develop an efficient genomic selection program in domestic animals. In this study, 2,973 American mink of six color types from two farms (Canadian Centre for Fur Animal Research (CCFAR), Truro, NS and Millbank Fur Farm (MFF), Rockwood, ON) were genotyped with the Affymetrix Mink 70K panel to compute their linkage disequilibrium (LD) patterns, effective population size (Ne), genetic diversity, genetic distances, and population differentiation and structure. The LD pattern represented by average r ², decreased to <0.2 when the inter-marker interval reached larger than 350 kb and 650 kb for CCFAR and MFF, respectively, and suggested at least 7,700 and 4,200 single nucleotide polymorphisms (SNPs) be used to obtain adequate accuracy for genomic selection programs in CCFAR and MFF respectively. The Ne for five generations ago was estimated to be 76 and 91 respectively. Our results from genetic distance and diversity analyses showed that American mink of the various color types had a close genetic relationship and low genetic diversity, with most of the genetic variation occurring within rather than between color types. Three ancestral genetic groups was considered the most appropriate number to delineate the genetic structure of these populations. Black (in both CCFAR and MFF) and pastel color types had their own ancestral clusters, while demi, mahogany, and stardust color types were admixed with the three ancestral genetic groups. This study provided essential information to utilize the first medium-density SNP panel for American mink in their genomic studies.

A chromosome-level genome assembly reveals genomic characteristics of the American mink (Neogale vison)

Availability of a contiguous chromosome-level genome assembly is the foundational step to develop genome-based studies in American mink (Neogale vison). The main objective of this study was to provide a high quality chromosome-level genome assembly for American mink. An initial draft of the genome assembly was generated using 2,884,047 PacBio long reads. Integration of Hi-C data into the initial draft led to an assembly with 183 scaffolds and scaffold N50 of 220 Mb. This gap-free genome assembly of American mink (ASM_NN_V1) had a length of 2.68 Gb in which about 98.6% of the whole genome was covered by 15 chromosomes. In total, 25,377 genes were predicted across the American mink genome using the NCBI Eukaryotic Genome Annotation Pipeline. In addition, gene orthology, demographic history, synteny blocks, and phylogenetic relationships were studied in connection with the genomes of other related Carnivora. Furthermore, population-based statistics of 100 sequenced mink were presented using the newly assembled genome. Remarkable improvements were observed in genome contiguity, the number of scaffolds, and annotation compared to the first draft of mink genome assembly (NNQGG.v01). This high-quality genome assembly will support the development of efficient breeding strategies as well as conservation programs for American mink.

Genetic and Phenotypic Parameters for Pelt Quality and Body Length and Weight Traits in American Mink

Understanding the genetics of fur characteristics and skin size is important for developing effective breeding programs in the mink industry. Therefore, the objectives of this study were to estimate the genetic and phenotypic parameters for pelt quality traits including live grading overall quality (LQU), live grading nap size (LNAP), dried pelt size (DPS), dried pelt nap size (DNAP) and overall quality of dried pelt (DQU), and body length and weight traits, including November body weight (Nov_BW), November body length (Nov_BL), harvest weight (HW) and harvest length (HL) in American mink. Dried pelt quality traits on 1195 mink and pelt quality traits on live animals on 1680 were collected from mink raised at two farms, in Nova Scotia and Ontario. A series of univariate analyses were implemented in ASReml 4.1 software to identify the significance (p < 0.05) of random effects (maternal genetic effects, and common litter effects) and fixed effects (farm, sex, color type, year, and age) for each trait. Subsequently, bivariate models were used to estimate the genetic and phenotypic parameters using ASReml 4.1. Heritability (±SE) estimates were 0.41 ± 0.06 for DPS, 0.23 ± 0.10 for DNAP, 0.12 ± 0.04 for DQU, 0.28 ± 0.06 for LQU, 0.44 ± 0.07 for LNAP, 0.29 ± 0.10 for Nov_BW, 0.28 ± 0.09 for Nov_BL, 0.41 ± 0.07 for HW and 0.31 ± 0.06 for HL. DPS had high positive genetic correlations (±SE) with Nov_BW (0.89 ± 0.10), Nov_BL (0.81 ± 0.07), HW (0.85 ± 0.05) and HL (0.85 ± 0.06). These results suggested that body weight and length measured on live animals in November of the first year were reliable indicators of dried pelt size. DQU had favorable genetic correlations with Nov_BL (0.55 ± 0.24) and HL (0.46 ± 0.20), and nonsignificant genetic correlations with DNAP (0.13 ± 0.25), Nov_BW (0.25 ± 0.25) and HW (0.06 ± 0.20), which made body length traits an appealing trait for selection for increased pelt size. High positive genetic correlation (±SE) was observed between LNAP and DNAP (0.82 ± 0.22), which revealed that nap size measurement on live animals is a reliable indicator trait for dried pelt nap size. However, nonsignificant (p > 0.05) low genetic correlation (±SE) was obtained between LQU and DQU (0.08 ± 0.45), showing that indirect selection based on live grading might not lead to the satisfactory improvement of dried pelt overall quality. The estimated genetic parameters for live grading, dried pelt quality, and body weight and body length traits may be incorporated into breeding programs to improve fur characteristics in Canadian mink populations.

Genetic and phenotypic parameters for feed efficiency and component traits in American mink

Feed cost is the largest expense of mink production systems, and therefore, improvement of feed efficiency (FE) through selection for high feed efficient mink is a practical way to increase the mink industry's sustainability. In this study, we estimated the heritability, phenotypic and genetic correlations for different FE measures and component traits, including harvest weight (HW), harvest length (HL), final body length (FBL), final body weight (FBW), average daily gain (ADG), daily feed intake (DFI), feed conversion ratio (FCR), residual feed intake (RFI), residual gain (RG), residual intake and gain (RIG), and Kleiber ratio (KR), using data from 2,288 American mink (for HW and HL), and 1,038-1,906 American mink (for other traits). Significance (P < 0.05) of fixed effects (farm, sex, and color-type), a covariate (age of animal), and random effects (additive genetic, maternal, and common litter) were evaluated through univariate models implemented in ASReml-R version 4. Genetic parameters were estimated via fitting a set of bivariate models using ASReml-R version 4. Estimates of heritabilities (±SE) were 0.28±0.06, 0.23±0.06, 0.28±0.10, 0.27±0.11, 0.25±0.09, 0.26±0.09, 0.20±0.09, 0.23±0.09, 0.21±0.10, 0.25±0.10, and 0.26±0.10 for HW, HL, FBL, FBW, ADG, DFI, FCR, RFI, RG, RIG, and KR, respectively. RIG had favorable genetic correlations with DFI (-0.62±0.24) and ADG (0.58±0.21), and non-significant (P > 0.05) genetic correlations with FBW (0.14±0.31) and FBL (-0.15±0.31). These results revealed that RIG might be superior trait as it guarantees reduced feed intake with faster-growing mink yet with no negative impacts on body weight and length. In addition, the strong positive genetic correlations (±SE) between KR with component traits (0.88±0.11 with FBW; 0.68±0.17 with FBL; and 0.97±0.02 with ADG) suggested KR as an applicable indirect measure of FE for improvement of component traits as it did not require the individual feed intake to be measured. Overall, our results confirmed the possibility of including FE traits in mink breeding programs to effectively select feed-efficient animals.

SNP markers associated with body size and pelt length in American mink (Neovison vison)

Background

Identification of genes underlying production traits is a key aim of the mink research community. Recent availability of genomic tools have opened the possibility for faster genetic progress in mink breeding. Availability of mink genome assembly allows genome-wide association studies in mink.

Results

In this study, we used genotyping-by-sequencing to obtain single nucleotide polymorphism (SNP) genotypes of 2496 mink. After multiple rounds of filtering, we retained 28,336 high quality SNPs and 2352 individuals for a genome-wide association study (GWAS). We performed the first GWAS for body weight, behavior, along with 10 traits related to fur quality in mink.

Conclusions

Combining association results with existing functional information of genes and mammalian phenotype databases, we proposed WWC3, MAP2K4, SLC7A1 and USP22 as candidate genes for body weight and pelt length in mink.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0688-6) contains supplementary material, which is available to authorized users.

The first draft reference genome of the American mink (Neovison vison)

The American mink (Neovison vison) is a semiaquatic species of mustelid native to North America. It's an important animal for the fur industry. Many efforts have been made to locate genes influencing fur quality and color, but this search has been impeded by the lack of a reference genome. Here we present the first draft genome of mink. In our study, two mink individuals were sequenced by Illumina sequencing with 797 Gb sequence generated. Assembly yielded 7,175 scaffolds with an N50 of 6.3 Mb and length of 2.4 Gb including gaps. Repeat sequences constitute around 31% of the genome, which is lower than for dog and cat genomes. The alignments of mink, ferret and dog genomes help to illustrate the chromosomes rearrangement. Gene annotation identified 21,053 protein-coding sequences present in mink genome. The reference genome's structure is consistent with the microsatellite-based genetic map. Mapping of well-studied genes known to be involved in coat quality and coat color, and previously located fur quality QTL provide new knowledge about putative candidate genes for fur traits. The draft genome shows great potential to facilitate genomic research towards improved breeding for high fur quality animals and strengthen our understanding on evolution of Carnivora.

Identifying QTL and genetic correlations between fur quality traits in mink (Neovison vison)

Mapping of QTL affecting fur quality traits (guard hair length, guard hair thickness, density of wool, surface of the fur and quality) and skin length was performed in a three-generation mink population (F2 design). In the parental generation, Nordic Brown mink were crossed reciprocally with American Black short nap mink. In all, 1082 mink encompassing three generations were used for the analyses. The mink were genotyped for 104 microsatellites covering all 14 autosomes. The QTL analyses were performed by least-square regression implemented in gridqtl software. Genetic and phenotypic correlations and heritabilities were estimated using the average information-restricted maximum-likelihood method. Evidence was found for QTL affecting fur quality traits on nine autosomes. QTL were detected for guard hair thickness on chromosomes 1, 2, 3, 6 and 13; for guard hair length on chromosomes 2, 3 and 6; for wool density on chromosomes 6 and 13; for surface on chromosomes 7, 12 and 13; for quality on chromosomes 6, 7, 11 and 13; and for skin length on chromosomes 7 and 9. Proximity of locations of QTL for guard hair length, guard hair thickness and for wool density and quality suggests that some of the traits are in part under the influence of the same genes. Traits under the influence of QTL at close or identical positions also were traits that were strongly genotypically correlated. Based on the results of correlation analyses, the most important single traits influencing the quality were found to be density of wool, guard hair thickness and appearance of the surface.

A frameshift mutation in the LYST gene is responsible for the Aleutian color and the associated Chédiak-Higashi syndrome in American mink

One of the colors of mink is Aleutian (aa)-a specific gun-metal gray pigmentation of the fur-commonly used in combination with other color loci to generate popular colors such as Violet (aammpp) and Sapphire (aapp). The Aleutian color allele is a manifestation of mink Chédiak-Higashi syndrome (CHS), which has been described in humans and several other species. As with forms of CHS in other species, we report that the mink CHS is linked to the lysosomal trafficking regulator ( LYST ) gene. Furthermore, we have identified a base deletion (c.9468delC) in exon 40 of LYST, which causes a frameshift and virtually terminates the LYST product prematurely (p.Leu3156Phefs*37). We investigated the blood parameters of three wild-type mink and three CHS mink. No difference in the platelet number between the two groups was observed, but an accumulation of platelets between the groups appears different when collagen is used as a coagulant. Microscopic analysis of peripheral blood indicates giant inclusions in the neutrophils of the Aleutian mink types. Molecular findings at the LYST locus enable the development of genetic tests for analyzing the color selection in American mink.

Construction of an American mink bacterial artificial chromosome (BAC) library and sequencing candidate genes important for the fur industry

Background

Bacterial artificial chromosome (BAC) libraries continue to be invaluable tools for the genomic analysis of complex organisms. Complemented by the newly and fast growing deep sequencing technologies, they provide an excellent source of information in genomics projects.

Results

Here, we report the construction and characterization of the CHORI-231 BAC library constructed from a Danish-farmed, male American mink (Neovison vison). The library contains approximately 165,888 clones with an average insert size of 170 kb, representing approximately 10-fold coverage. High-density filters, each consisting of 18,432 clones spotted in duplicate, have been produced for hybridization screening and are publicly available. Overgo probes derived from expressed sequence tags (ESTs), representing 21 candidate genes for traits important for the mink industry, were used to screen the BAC library. These included candidate genes for coat coloring, hair growth and length, coarseness, and some receptors potentially involved in viral diseases in mink. The extensive screening yielded positive results for 19 of these genes. Thirty-five clones corresponding to 19 genes were sequenced using 454 Roche, and large contigs (184 kb in average) were assembled. Knowing the complete sequences of these candidate genes will enable confirmation of the association with a phenotype and the finding of causative mutations for the targeted phenotypes.

Additionally, 1577 BAC clones were end sequenced; 2505 BAC end sequences (80% of BACs) were obtained. An excess of 2 Mb has been analyzed, thus giving a snapshot of the mink genome.

Conclusions

The availability of the CHORI-321 American mink BAC library will aid in identification of genes and genomic regions of interest. We have demonstrated how the library can be used to identify specific genes of interest, develop genetic markers, and for BAC end sequencing and deep sequencing of selected clones. To our knowledge, this is the first report of 454 sequencing of selected BAC clones in mammals and re-assures the suitability of this technique for obtaining the sequence information of genes of interest in small genomics projects. The BAC end sequences described in this paper have been deposited in the GenBank data library [HN339419-HN341884, HN604664-HN604702]. The 454 produced contigs derived from selected clones are deposited with reference numbers [GenBank: JF288166-JF288183 &JF310744].

A genetic linkage map and comparative mapping of the prairie vole (Microtus ochrogaster) genome

BACKGROUND

The prairie vole (Microtus ochrogaster) is an emerging rodent model for investigating the genetics, evolution and molecular mechanisms of social behavior. Though a karyotype for the prairie vole has been reported and low-resolution comparative cytogenetic analyses have been done in this species, other basic genetic resources for this species, such as a genetic linkage map, are lacking.

RESULTS

Here we report the construction of a genome-wide linkage map of the prairie vole. The linkage map consists of 406 markers that are spaced on average every 7 Mb and span an estimated ~90% of the genome. The sex average length of the linkage map is 1707 cM, which, like other Muroid rodent linkage maps, is on the lower end of the length distribution of linkage maps reported to date for placental mammals. Linkage groups were assigned to 19 out of the 26 prairie vole autosomes as well as the X chromosome. Comparative analyses of the prairie vole linkage map based on the location of 387 Type I markers identified 61 large blocks of synteny with the mouse genome. In addition, the results of the comparative analyses revealed a potential elevated rate of inversions in the prairie vole lineage compared to the laboratory mouse and rat.

CONCLUSIONS

A genetic linkage map of the prairie vole has been constructed and represents the fourth genome-wide high-resolution linkage map reported for Muroid rodents and the first for a member of the Arvicolinae sub-family. This resource will advance studies designed to dissect the genetic basis of a variety of social behaviors and other traits in the prairie vole as well as our understanding of genome evolution in the genus Microtus.