Modeling Continuous Admixture Using Admixtureinduced Linkage Disequilibrium

• Journal List
• Front Genet
• PMC6539196

Front Genet. 2019; 10: 498.

Whole Genome Linkage Disequilibrium and Effective Population Size in a Coho Salmon (Oncorhynchus kisutch) Breeding Population Using a High-Density SNP Array

Agustín Barría, ¹ Kris A. Christensen, ² Grazyella Yoshida, ¹ Ana Jedlicki, ¹ Jong S. Leong, ² Eric B. Rondeau, ² Jean P. Lhorente, ³ Ben F. Koop, ² William S. Davidson, ⁴ and José M. Yáñez ^{1 , 5 , *}

Agustín Barría

¹ Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Santiago, Chile

Kris A. Christensen

² Department of Biology, Centre for Biomedical Research, University of Victoria, Victoria, BC, Canada

Grazyella Yoshida

¹ Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Santiago, Chile

Ana Jedlicki

¹ Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Santiago, Chile

Jong S. Leong

² Department of Biology, Centre for Biomedical Research, University of Victoria, Victoria, BC, Canada

Eric B. Rondeau

² Department of Biology, Centre for Biomedical Research, University of Victoria, Victoria, BC, Canada

Jean P. Lhorente

³ Benchmark Genetics Chile, Puerto Montt, Chile

Ben F. Koop

² Department of Biology, Centre for Biomedical Research, University of Victoria, Victoria, BC, Canada

William S. Davidson

⁴ Department of Molecular Biology and Biochemistry, Simon Fraser University, Burnaby, BC, Canada

José M. Yáñez

¹ Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Santiago, Chile

⁵ Nucleo Milenio INVASAL, Concepcion, Chile

Received 2018 May 13; Accepted 2019 May 7.

Abstract

The estimation of linkage disequilibrium between molecular markers within a population is critical when establishing the minimum number of markers required for association studies, genomic selection, and inferring historical events influencing different populations. This work aimed to evaluate the extent and decay of linkage disequilibrium in a coho salmon breeding population using a high-density SNP array. Linkage disequilibrium was estimated between a total of 93,502 SNPs found in 64 individuals (33 dams and 31 sires) from the breeding population. The markers encompass all 30 coho salmon chromosomes and comprise 1,684.62 Mb of the genome. The average density of markers per chromosome ranged from 48.31 to 66 per 1 Mb. The minor allele frequency averaged 0.26 (with a range from 0.22 to 0.27). The overall average linkage disequilibrium among SNPs pairs measured as r ² was 0.10. The Average r ² value decreased with increasing physical distance, with values ranging from 0.21 to 0.07 at a distance lower than 1 kb and up to 10 Mb, respectively. An r ² threshold of 0.2 was reached at distance of approximately 40 Kb. Chromosomes Okis05, Okis15 and Okis28 showed high levels of linkage disequilibrium (>0.20 at distances lower than 1 Mb). Average r ² values were lower than 0.15 for all chromosomes at distances greater than 4 Mb. An effective population size of 43 was estimated for the population 10 generations ago, and 325, for 139 generations ago. Based on the effective number of chromosome segments, we suggest that at least 74,000 SNPs would be necessary for an association mapping study and genomic predictions. Therefore, the SNP panel used allowed us to capture high-resolution information in the farmed coho salmon population. Furthermore, based on the contemporary N _e, a new mate allocation strategy is suggested to increase the effective population size.

Keywords: linkage disequilibrium, Oncorhynchus kisutch, selective breeding, GWAS, effective population size

Background

Coho salmon (Oncorhynchus kisutch) is one of the six Pacific salmon species found in North American and Asian watersheds (Groot and Margolis, 1991). This species was introduced into Chilean streams during the 1920s promoted by the Chilean Institute of Fisheries Department. Cultivation of coho salmon began in Chile at the end of the 1970s, when Chile imported almost 500,000 eggs from the Kitimat river (British Columbia) and Oregon, becoming the genetic basis of the broodstocks in Chile (Neira et al., 2014). Twenty years later, the production of the first eggs for commercial use was produced in Chile (SalmonChile, 2007). Currently, Chile is the main producer of farmed coho salmon, with the production of nearly 160,000 tons in 2014 (FAO, 2016). This represents more than 90% of the global farmed coho production (Canada and Japan are the other major coho salmon producers) (FAO, 2016). The temperature and the quality of the Chilean freshwater environments have reduced the coho reproductive cycle to only 2 years (Estay et al., 1997). To date, numerous genetic programs have been developed for coho salmon in Chile. These programs are mainly focused on growth, disease resistance, and flesh color (Neira et al., 2014).

With the advent of next generation sequencing (NGS) and high-throughput genotyping technologies, it has become possible to perform artificial selection through the use of genomic estimated breeding values (GEBVs). By using dense molecular markers from the whole genome, genomic selection (GS) can be used in broodstock enhancement (Bennewitz et al., 2009). This methodology makes it possible to estimate GEBVs with high accuracy, even with animals without recorded phenotypes (Meuwissen et al., 2001), which has improved the accuracy of selection in salmonid species (Ødegård et al., 2014; Tsai et al., 2016; Bangera et al., 2017; Correa et al., 2017; Yoshida et al., 2018; Barría et al., 2018a). Genome wide association studies (GWAs) and GS, exploit linkage disequilibrium (LD) between molecular markers. The amount of LD between loci is important in GWAs, as the extent of LD indicates the necessary number of SNPs to assure that causative mutations are in LD with genetic markers (Flint-Garcia et al., 2003). GWAs are key for mapping traits with commercial interest to specific variants in the genome. For GS, LD is related to the likelihood of successfully tagging the SNP effect in genomic breeding value prediction (Kemper and Goddard, 2012).

LD allows researchers to explore the genetic basis of traits influencing productivity. Through the comparison of the extent and pattern of LD, it is possible to elucidate the diversity among breeds with different phenotypic attributes, and even identify genomic regions subject to different selective pressures (McKay et al., 2007; López et al., 2015). The most common LD measurements are r ² and |D′|, both ranging from 0 to 1. When |D′| < 1, it indicates the occurrence of historical recombination between loci, while |D′| = 1 indicates no recombination. The r ² statistic represents the correlation between genotypes from molecular marker pairs. This latter parameter is preferred over |D′| because |D′| tends to be overestimated in small samples sizes and when low-frequency alleles are used (Teare et al., 2002). Moreover, in association studies, r ² is preferred due to the inverse relationship between its value and the sample size needed to detect a significant association between a causative variant and molecular markers (Wall and Pritchard, 2003).

Despite the many GWAs and GS analyses performed in Atlantic salmon (Gutierrez et al., 2015; Tsai et al., 2015, 2016; Bangera et al., 2017; Correa et al., 2017), rainbow trout (Vallejo et al., 2016, 2017) and coho salmon (Barría et al., 2018a), none of them have evaluated the LD in the studied populations. Further, most of the linkage disequilibrium studies have been focused on the extent and decay pattern of LD in livestock species, such as dairy (Sargolzaei et al., 2008; Bohmanova et al., 2010) and beef cattle (McKay et al., 2007; Makina et al., 2015), plants (Delourme et al., 2013; Porto-Neto et al., 2014), and pigs (Saura et al., 2015). Recently, LD has been evaluated in farmed rainbow trout (Oncorhynchus mykiss) (Rexroad and Vallejo, 2009) and in Atlantic salmon (Kijas et al., 2017; Barría et al., 2018b).

The first step to calculate the number of molecular markers necessary for genomic selection and association mapping is to estimate the extent and decline of LD within a population. To date, there have been no studies aimed to characterize the levels and extent of LD in coho salmon. The current work aimed to evaluate the effective population size and the extent of linkage disequilibrium, at the genomic and chromosome level, on a breeding coho salmon population using a 200K high-density SNP chip array.

Materials and Methods

Populations and Samples

The coho salmon samples were obtained from a breeding population belonging to a genetic improvement program established in 1998 run by Pesquera Antares (Puerto Montt, Chile). Using best linear unbiased prediction (BLUP), harvest weight had been selected over eight generations in this population. For LD estimations, a total of 64 animals (33 sires and 31 dams), corresponding to the parents of 33 families from a 2012-spawning year class, were selected. The mean relatedness among individuals (0.07) was estimated using Plink v1.09 (Purcell et al., 2007). For specific details about reproductive management, mating design, rearing conditions and inbreeding and breeding objectives of the genetic program for this population see Dufflocq et al. (2016) and Yáñez et al. (2014, 2016). Sampling protocols were approved by the Animal Bioethics Committee from Universidad de Chile (No. 08-2015).

Genotyping

Genomic DNA was extracted from fin clips from the 64 individuals. Genotyping was carried out using a 200K Affymetrix Axiom® myDesign Custom Array developed for coho salmon by the EPIC4 genome consortium¹ and constructed by ThermoFisher Scientific. Subsequent work will describe the design and performance of the array in greater detail, but in brief: This dense SNP array contains 203,077 polymorphic SNPs. Genotyping of the SNP array was performed by the McGill University and Genome Quebec Innovation Centre. Genotype calling was performed using Axiom Analysis Suite v3.1 (Thermo Scientific) following the Axiom Analysis user guide. A total of 22 SNPs were subsequently discarded prior to analysis due to unknown position on the coho salmon reference genome (GCF_002021735.1); a further 35,569 markers were discarded as they were identified as problematic (OTV, Call Rate Below Threshold, Other). This left 167,486 SNP markers for further consideration. The following parameters were used to exclude low-confidence SNPs using plink software: Hardy-Weinberg Equilibrium (HWE) p < 3.8e−7, Minor Allele Frequency (MAF) ≤ 0.05 and genotyping call rate < 0.95. Fish with genotyping call rates <0.95 were excluded from further analyses. Subsequent analyses were assessed using the SNPs markers which passed all quality control (QC) criteria.

LD Estimation

The LD between each pair of genetic markers was estimated using Pearson's squared correlation coefficient (r ²) statistic which is less sensitive to allelic frequencies (Ardlie et al., 2002), more suitable for biallelic markers (Zhao et al., 2005) and allows to compare estimations with previous studies. Pair-wise LD as r ² values, were estimated with Plink v1.09 (Purcell et al., 2007), based on the formula proposed by Hill and Robertson (1968). Genotypes were coded as 0, 1, and 2 relative to the number of non-reference alleles. The parameter -inter-chr, in conjunction with a ld-window-r ² set to zero, was used to obtain correlations between all the pairs of SNPs within each chromosome independently of their r ² value. Based on the physical distance, we created bins of 100 kb for each SNP pair. Decay and extend of the LD was estimated. LD decay curves for SNP pair were calculated as the average r ² within each bin, up to a distance of 10 Mb. Average r ² per chromosome was calculated sorting SNPs pairs into 10 bins according to an increasing average distance. The used distance was from 0 to 0.99 Mb and from 9 to 10 Mb (in the first and last bin, respectively), between SNPs pairs on each chromosome.

Effective Population Size and Number of Chromosome Segments

Contemporary effective population size (N _e) was estimated using NeEstimator v2.01 (Do et al., 2014). Estimation was assessed based on LD method (Waples and Do, 2008), setting a critical value of 0.05 and using a non-random mating model.

The historical effective population size was estimated using SNeP software v1.1 (Barbato et al., 2015). Based on estimated LD values, historical population size estimation was calculated with the following equation proposed by Corbin et al. (2012):

where N_t and c_t refers to the effective population size and recombination rate, t generations ago, respectively. Being the latter, proportional to the physical distance between the SNPs, r a d j 2 is the adjusted estimation of LD based on sample size, and α refers to the adjustment for mutation rate. Considering that mutation does occurs, we used an α = 2 (Tenesa et al., 2007; Vallejo et al., 2018). N _e estimation was calculated with a minimum and maximum distance between SNPs of 0 and 5 Mb, respectively. Data was arranged in 30 bins of 50 kb distance each. Thus, N _e was calculated from the r ² estimated for the average distance of each bin.

Effective number of chromosome segments (M _e) was estimated based on the following formula proposed by Goddard et al. (2011).

where mean(r ²) refers to the average linkage disequilibrium over all pairwise combination of SNPs within each chromosome. Thus, the total M_e within the population is the sum of each individual M_e .

Results

SNPs Quality Control

No individuals were removed after QC. From the 167,486 SNPs anchored to chromosomes, a total of 93,502 SNPs passed the quality control and were identified as segregating along the coho salmon genome. The MAF distribution of the identified SNPs was nearly uniform along the 30 chromosomes, with an average of 0.26 ± 0.01 (mean ± standard deviation), and a minimum and maximum value of 0.22 and 0.27, respectively (Table 1).

Table 1

Summary statistics for the evaluated SNPs and linkage disequilibrium values along coho salmon chromosomes.

Okis	Length (Mb)	Number of SNPs	SNP density (Mb)	Mean (r 2)	Median (r 2)	SD (r 2)	MAF
01	67.36	3,840	57.01	0.12	0.06	0.16	0.27
02	74.29	3,990	53.71	0.10	0.05	0.13	0.26
03	70.07	3,882	55.40	0.12	0.06	0.10	0.27
04	79.83	4,846	60.70	0.10	0.05	0.14	0.26
05	71.75	3,748	52.24	0.13	0.06	0.17	0.22
06	76.69	3,854	50.25	0.13	0.05	0.18	0.27
07	50.39	2,746	54.49	0.11	0.05	0.15	0.27
08	67.50	3,806	56.39	0.12	0.06	0.16	0.27
09	39.42	2,166	54.95	0.07	0.03	0.11	0.25
10	65.20	3,449	52.90	0.10	0.05	0.14	0.26
11	79.38	4,036	50.84	0.14	0.07	0.18	0.27
12	51.25	3,135	61.17	0.09	0.04	0.12	0.26
13	66.77	3,755	56.24	0.09	0.04	0.13	0.26
14	71.35	3,712	52.03	0.11	0.06	0.14	0.27
15	66.87	3,662	54.76	0.10	0.05	0.13	0.26
16	33.63	1,954	58.10	0.08	0.04	0.11	0.26
17	75.53	4,200	55.61	0.12	0.06	0.16	0.27
18	66.40	3,570	53.77	0.10	0.05	0.13	0.25
19	54.89	2,877	52.41	0.11	0.05	0.15	0.26
20	40.41	2,333	57.73	0.08	0.03	0.12	0.26
21	34.95	1,921	54.96	0.10	0.05	0.14	0.27
22	55.52	3,062	55.15	0.10	0.04	0.14	0.25
23	42.32	2,600	61.44	0.09	0.04	0.13	0.26
24	39.26	2,457	62.58	0.09	0.04	0.11	0.26
25	33.74	2,014	59.69	0.08	0.04	0.11	0.26
26	43.46	2,428	55.87	0.10	0.05	0.14	0.27
27	38.53	2,543	66.00	0.09	0.04	0.12	0.25
28	47.26	2,283	48.31	0.15	0.06	0.20	0.24
29	38.39	2,194	57.15	0.08	0.04	0.11	0.25
30	42.21	2,439	57.78	0.11	0.06	0.15	0.26
Mean	56.16	3,116	55.99	0.10	0.05	0.14	0.26

SNP, single-nucleotide polymorphism; MAF, minor allele frequency; SD, standard deviation.

Estimation of LD

Table 1 summarizes the mean, median and standard deviation of r ² values for each coho salmon chromosome. All of the 93,502 SNPs placed onto chromosomes and which passed quality control were included in this analysis. These markers encompassed 1,684.62 Mb of the genome, representing ~71.1% of the total genome size (assuming a genome size of 2,369 Mb based on the total sequence length of the final assembly GCF_002021735.1). The molecular marker density per chromosome per Mb, ranged from 48.31 to 66 with a mean of 55.99. In general, SNPs were uniformly distributed along the 30 chromosomes. The number of SNPs on each chromosome ranged from 1,954 on Okis16 to 4,846 on Okis04, which is in agreement with Okis16 and Okis04 being the shortest and longest chromosome, respectively. The overall mean linkage disequilibrium (measured as r ²) among SNP pairs was 0.10 ± 0.14. The global median was lower at 0.05. Low average LD among adjacent SNPs along the 30 chromosomes was observed in the current population, with values ranging from 0.07 to 0.15 (Table 1).

To estimate the decay of linkage disequilibrium as a function of physical distance, SNP pairs were sorted into bins of 100 kb, and mean values of r ² were calculated for each bin. As observed in other species (Lu et al., 2012; Kijas et al., 2017; Vos et al., 2017), LD declines smoothly as the physical distance increases between markers (Figure 1). A maximum average LD of 0.21 was estimated for SNPs less than 1 kb apart. This value declines quickly at marker distances up to 0.1 Mb, with a value of 0.16. From 1 Mb to 10 Mb LD range from 0.14 to 0.07. The latter value represents the lowest average LD estimated in the current data set. The r ² estimation drops below 0.2 at a distance of ~40 kb.

Decay of average LD (r ²) over distance among SNPs in coho salmon (Oncorhynchus kisutch) population. The blue line shows the mean LD in each 100 kb sliding window. Each black dot represents average r ² within each bin.

Comparison of average LD at different distance bins for each chromosome shows higher variation at closer distance bins (Figure 2). Suggesting that estimations of genome-wide linkage disequilibrium based on few chromosomes may be biased (Khatkar et al. 2008). Lower estimates of LD (<0.13) were found in Okis09, Okis12, Okis24, and Okis25, while higher levels of LD (>0.20) were estimated for Okis05, Okis15, and Okis28. When the distance among SNPs increased over 4 Mb, average LD values drops below 0.15 for all chromosomes. Furthermore, average r ² values <0.10 were estimated for all chromosomes, except for Okis11 and Okis28 at distances greater than 7 Mb.

Linkage disequilibrium estimations along the 30 chromosomes of coho salmon. Average values of LD measured as r ² per chromosome, according to distances between SNPs. Estimated values are shown from Okis01 to Okis015 (A) and from Okis16 to Okis30 (B).

Effective Population Size and Number of Chromosome Segments

Based on LD approach, estimated contemporary N _e reached up to 83.9. Figure 3 illustrates the estimated historical effective population size of the coho salmon, based on LD, from 10 to 241 generations ago. An increasing N _e as a function of the number of generation was observed, with a N _e of 43 estimated at 10 generations ago, and 543 for 241 generations ago.

Effective population size estimation in coho salmon population. Estimates of effective population size (N _e) over the past 241 generations based on the LD of an aquaculture strain of coho salmon. Black dots represent each N _e measurement.

Considering the estimated linkage disequilibrium (measured as r ²) among all SNP pairs within each chromosome, we estimated that at least 74K SNPs markers would be needed for an efficient GWAs or GS analysis in the current farmed coho salmon population.

Discussion

Understanding LD patterns enhances our knowledge of the demographic processes within the population. Biological factors such as recombination and mutation in conjunction with genetic drift, admixture and effective population size are important variables determining patterns of LD. For this reason, variation in LD among populations and genomic regions are of high interest and widely reported.

To our knowledge, this is the first study characterizing the whole-genome LD in a coho salmon population using a high-density SNP panel. The samples originated from the broodstock of a breeding program aimed at improving economically important traits for Chilean coho salmon aquaculture. Unrelated animals were chosen in order to avoid LD inflation that can occur when high kinship relationships are present in the sampled individuals (Gutierrez et al., 2015). Due to the increased bias of LD estimations, when estimating |D′| from small sample sizes (Bohmanova et al., 2010), we preferred to use the robust r ² statistic. Moreover, to predict the power of association mapping, r ² statistics is more useful. The minimum number of individuals necessary for an accurate r ² estimation has been suggested to range from 55 to 75 in cattle (Khatkar et al., 2008; Bohmanova et al., 2010). This range increases to 400 or more in case of |D′| (Khatkar et al., 2008). The number of individuals necessary to estimate LD depends on the demographic and genetic population history. Our sample size was within the range suggested above.

Sample sizes above 50 also provide accurate estimations of MAFs (>0.05) within a population, at a physical distance up to 10 Mb (Khatkar et al., 2008). Filtered markers showed an average MAF of 0.26 (ranging from 0.22 to 0.27) per chromosome (Table 1). A similar mean value was reported in Nellore cattle, ranging from 0.20 to 0.25 (Matukumalli et al., 2009; Espigolan et al., 2013) and from 0.28 to 0.30 in North American Holstein (Bohmanova et al., 2010). Despite that LD measurements of r ², tend to be less sensitive than |D′| to low MAF (Khatkar et al., 2008; Bohmanova et al., 2010; Kijas et al., 2017), low MAFs tends to underestimate LD measures (Espigolan et al., 2013). Thus, the high allele frequencies found in the current study suggest an unbiased estimation of our LD values.

Estimations of the extent and decay of linkage disequilibrium in the coho salmon breeding population provide insights into LD patterns in the coho salmon genome, which may have implications for GWAs, GS and for the design of SNP arrays. In terms of genomic predictions, the accuracy is influenced by multiple factors such as genetic architecture of the trait of interest, average relatedness between training and testing individuals, density of the SNP panel, effective population size, and effective number of chromosomes segments (M _e) (Daetwyler et al., 2010; Pszczola et al., 2012). Due to LD, loci do not segregate independently within a finite population. Thus, M _e can be considered as the number of chromosome segment which segregates independently within a population. When a lower number of segments are estimated, independent parameters are needed to be estimated from the same dataset, i.e., higher accuracy of genomic predictions (Yvonne et al., 2016; Lee et al., 2017). Based on the relatedness among pair of animals, Goddard (2009) suggested that the effective number of chromosomes is 2N _eL/Log(4N _eL). However, this approach may produce an underestimated M _e, leading to an overestimation on the EBVs accuracy (Goddard, 2009). Thus, we decided to used 1/mean(r ²) (Goddard et al., 2011) which is a much more direct approach. Even that the 74K SNPs predicted for this coho salmon population, is much higher than the estimated in a farmed rainbow trout population, in which authors suggested at least 20K markers for an efficient GWAs and GS study (Vallejo et al., 2018), the SNP array used in the current study allows coverage of the whole genome. However, further studies are needed to elucidate the minimum number of markers necessary to achieve a higher accuracy predictions compared to pedigree-based models, in the current coho salmon breeding population.

The variation in the average and standard deviation in the LD among chromosomes found in the current study (Table 1), is partly explained due to variation in recombination rates along different chromosomes (e.g., local hotspots for recombination), decreasing as a function of an increase in chromosome length (Arias et al., 2009; Espigolan et al., 2013), which in turn generates variation in LD along the different chromosomes, as observed in Figure 2. Therefore, inferences based on single or only on few chromosomes might be biased and inferences regarding LD would be best when using genome-wide data. LD information from the population may allow researchers to reduce the number of required SNPs for a genomic analysis by excluding redundant SNPs (Khatkar et al., 2008). This can be done by identifying tag SNPs, using information from haplotype block structure, as was previously done in Holstein-Friesian cattle (Khatkar et al., 2007).

Average r ² values estimated in our study were higher than those estimated in a wild Finnish Atlantic salmon population, with values ranging from 0.015 to 0.037 (Kijas et al., 2017). However, farmed Tasmanian Atlantic salmon showed mean LD (measured as r ²) values up to 0.67 for SNPs closer than 1 kb (Kijas et al., 2017), almost three times higher than in the current work (0.21). However, our estimation is similar as the one found in Chilean Atlantic salmon populations with European origin (Barría et al., 2018b). Some authors have found low linkage disequilibrium estimations in others Atlantic salmon populations, although these estimations were reported in units of recombination (Gutierrez et al., 2015) and using sliding windows of 20 SNPs (Johnston et al., 2014). The different estimation metrics make it difficult to compare directly with the current work. The origin of the current breeding coho population most likely involves two isolated wild populations (The Kitimat River and Oregon). The admixture of both founders populations could explain the induced long-range and reduced short-range LD (Pfaff et al., 2001) estimated in this Coho salmon breeding population. A similar pattern has also been suggested in a highly admixed Norwegian Atlantic salmon population (Ødegård et al., 2014) and in a recently admixed farmed rainbow trout population (Vallejo et al., 2018).

We found that historical effective population size shows a decline from 543 individuals, 241 generations ago, to 43 individuals 10 generations ago. Similar N _e pattern reduction has been observed in cattle populations (Villa-Angulo et al., 2009; Makina et al., 2015). Even though this is the first study aimed to estimate the effective population size of a coho salmon breeding population, caution must be taken when evaluating the estimations for the number of generations (Corbin et al., 2012). For recent generations, large c values are involved and do not necessarily fit the theoretical implications proposed by Hayes (Hayes et al., 2003) for N _e estimations. In the oldest generation, after 4N _e generations ago, none of the SNPs can be reliably sampled (Corbin et al., 2012). Therefore, N _e estimations after 4N _e generations ago may be questionable.

Contemporary effective population size calculated based on the LD is 83.9. This estimation is higher than the estimated for 10 generations ago (approximately 20 year ago, assuming an interval generation of 2 years), which agree exactly with the beginning of the breeding program of this farmed population (Neira et al., 2014). This value, is still below the minimum N _e value (500) suggested to retain the evolutionary potential and sustainability of a breeding program (Franklin and Frankham, 1998; Ponzoni et al., 2010). Thus, a modified mate allocation strategy that allows to increase the effective population size is suggested for the current breeding population.

Conclusions

In the current study, we performed an LD analysis with 64 coho salmon genotyped with 93,502 SNPs. We showed the feasibility to estimate LD and infer the effective population size based on the observed LD using data from a high-density SNP array. Furthermore, based on the current effective population size, a new mate allocation strategy that allows to increase it is suggested.

Ethics Statement

Coho salmon individuals and sampling procedures were approved by the Comité de Bioética Animal from the Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile (Certificate N08-2015).

Author Contributions

AB performed DNA extraction, LD and Ne analysis, and wrote the initial version of the manuscript. KC contributed on the data analysis and discussion. GY contributed with LD analysis and discussion. AJ performed DNA extraction. JPL contributed with study design. BK, JSL, and ER developed the chip array. WD contributed with writing and discussion. JMY conceived and designed the study, supervised work of AB, and contributed to the analysis, discussion, and writing. All authors have reviewed and approved the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

AB wants to acknowledge the National Commission of Scientific and Technologic Research (CONICYT) for the funding through the National PhD funding program and to the Government of Canada for the funding through the Canada-Chile Leadership Exchange Scholarship (ELAP).

Funding. This project was funded by the U-Inicia grant, from the Vicerrectoria de Investigación y Desarrollo, Universidad de Chile. This work has been conceived on the frame of the grant FONDEF NEWTON-PICARTE (IT14I10100), funded by CONICYT (Government of Chile). This work has been partially supported by Núcleo Milenio INVASAL from Iniciativa Científica Milenio (Ministerio de Economía, Fomento y Turismo, Gobierno de Chile). This research was carried out in conjunction with EPIC4 (Enhanced Production in Coho: Culture, Community, Catch), a project supported by the government of Canada through Genome Canada, Genome British Columbia, and Genome Quebec.

References

• Ardlie K. G., Kruglyak L., Seielstad M. (2002). Patterns of linkage disequilibrium in the human genome. Nat. Rev. Genet. 3, 299–309. 10.1038/nrg777, PMID: [PubMed] [CrossRef] [Google Scholar]
• Arias J. A., Keehan M., Fisher P., Coppieters W., Spelman R. (2009). A high density linkage map of the bovine genome. BMC Genet. 10, 1–12. 10.1186/1471-2156-10-18 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Bangera R., Correa K., Lhorente J. P., Figueroa R., Yáñez J. M. (2017). Genomic predictions can accelerate selection for resistance against Piscirickettsia salmonis in Atlantic salmon (Salmo salar). BMC Genomics 18:121. 10.1186/s12864-017-3487-y [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Barbato M., Orozco-terWengel P., Tapio M., Bruford M. W. (2015). SNeP: a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data. Front. Genet. 6, 1–6. 10.3389/fgene.2015.00109 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Barría A., Christensen K. A., Yoshida G. M., Correa K., Jedlicki A., Lhorente J. P., et al. (2018a). Genomic predictions and genome-wide association study of resistance against Piscirickettsia salmonis in coho salmon (Oncorhynchus kisutch) using ddRAD sequencing. G3 Genes Genomes Genet. 8, 1183–1194. 10.1534/g3.118.200053 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Barría A., López M. E., Yoshida G., Carvalheiro R., Yáñez J. M. (2018b). Population genomic structure and genome-wide linkage disequilibrium in farmed Atlantic salmon (Salmo salar L.) using dense SNP genotypes. Front. Genet. 9:649. 10.1534/g3.118.200053 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Bennewitz J., Solberg T., Meuwissen T. H. (2009). Genomic breeding value estimation using nonparametric additive regression models. Genet. Sel. Evol. 41:20. 10.1186/1297-9686-41-20, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Bohmanova J., Sargolzaei M., Schenkel F. S. (2010). Characteristics of linkage disequilibrium in north American Holsteins. BMC Genomics 11:421. 10.1186/1471-2164-11-421, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Corbin L. J., Liu A. Y. H., Bishop S. C., Woolliams J. A. (2012). Estimation of historical effective population size using linkage disequilibria with marker data. J. Anim. Breed. Genet. 129, 257–270. 10.1111/j.1439-0388.2012.01003.x, PMID: [PubMed] [CrossRef] [Google Scholar]
• Correa K., Bangera R., Figueroa R., Lhorente J. P., Yáñez J. M. (2017). The use of genomic information increases the accuracy of breeding value predictions for sea louse (Caligus rogercresseyi) resistance in Atlantic salmon (Salmo salar). Genet. Sel. Evol. 49:15. 10.1186/s12711-017-0291-8 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Daetwyler H. D., Pong-Wong R., Villanueva B., Woolliams J. A. (2010). The impact of genetic architecture on genome-wide evaluation methods. Genetics 185, 1021–1031. 10.1534/genetics.110.116855, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Delourme R., Falentin C., Fomeju B. F., Boillot M., Lassalle G., André I., et al. (2013). High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L. BMC Genomics 14, 1–18. 10.1186/1471-2164-14-120 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Do C., Waples R. S., Peel D., Macbeth G. M., Tillett B. J., Ovenden J. R. (2014). NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol. Ecol. Resour. 14, 209–214. 10.1111/1755-0998.12157, PMID: [PubMed] [CrossRef] [Google Scholar]
• Dufflocq P., Lhorente J. P., Bangera R., Neira R., Newman S., Yáñez J. M. (2016). Correlated response of flesh color to selection for harvest weight in coho salmon (Oncorhynchus kisutch). Aquaculture 472, 6–11. 10.1016/j.aquaculture.2016.08.037 [CrossRef] [Google Scholar]
• Espigolan R., Baldi F., Boligon A. A., Souza F. R. P., Gordo D. G. M., Tonussi R. L., et al.. (2013). Study of whole genome linkage disequilibrium in Nellore cattle. BMC Genomics 14:305. 10.1186/1471-2164-14-305, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Estay F., Diaz N., Neira R., Garcia X. (1997). Reproductive performance of cultured female coho salmon in Chile. Prog. Fish-Cult 59, 36–40. 10.1577/1548-8640(1997)059<0036:RPOCFC>2.3.CO;2 [CrossRef] [Google Scholar]
• FAO (2016). Fisheries and aquaculture information and statistical branch. (Rome: FAO; ). [Google Scholar]
• Flint-Garcia S. A., Thornsberry J. M., Buckler E. S. (2003). Structure of linkage disequilibrium in plants. Annu. Rev. Plant Biol. 54, 357–374. 10.1146/annurev.arplant.54.031902.134907, PMID: [PubMed] [CrossRef] [Google Scholar]
• Franklin I. R., Frankham R. (1998). How large must populations be to retain evolutionary potential? Anim. Conserv. 1, 69–73. 10.1111/j.1469-1795.1998.tb00228.x [CrossRef] [Google Scholar]
• Goddard M. (2009). Genomic selection: prediction of accuracy and maximisation of long term response. Genetica 136, 245–257. 10.1007/s10709-008-9308-0, PMID: [PubMed] [CrossRef] [Google Scholar]
• Goddard M. E., Hayes B. J., Meuwissen T. H. E. (2011). Using the genomic relationship matrix to predict the accuracy of genomic selection. J. Anim. Breed. Genet. 128, 409–421. 10.1111/j.1439-0388.2011.00964.x, PMID: [PubMed] [CrossRef] [Google Scholar]
• Groot C., Margolis L. (eds) (1991). Pacific salmon life histories. (British Columbia, Vancouver, Canada: University of British Columbia Press; ). [Google Scholar]
• Gutierrez A. P., Yáñez J. M., Fukui S., Swift B., Davidson W. S. (2015). Genome-wide association study (GWAS) for growth rate and age at sexual maturation in Atlantic salmon (Salmo salar). PLoS One 10:e0119730. 10.1371/journal.pone.0119730, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Hayes B. J., Visscher P. M., Mcpartlan H. C., Goddard M. E. (2003). Novel multilocus measure of linkage disequilibrium to estimate past effective population size. Genome Res. 13, 635–643. 10.1101/gr.387103, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Hill W., Robertson A. (1968). Linkage disequilibrium in finite populations. Theor. Appl. Genet. 38, 226–231. 10.1007/BF01245622, PMID: [PubMed] [CrossRef] [Google Scholar]
• Johnston S. E., Orell P., Pritchard V. L., Kent M. P., Lien S., Niemelä E., et al.. (2014). Genome-wide SNP analysis reveals a genetic basis for sea-age variation in a wild population of Atlantic salmon (Salmo salar). Mol. Ecol. 23, 3452–3468. 10.1111/mec.12832, PMID: [PubMed] [CrossRef] [Google Scholar]
• Kemper K. E., Goddard M. E. (2012). Understanding and predicting complex traits: knowledge from cattle. Hum. Mol. Genet. 21, 45–51. 10.1093/hmg/dds332 [PubMed] [CrossRef] [Google Scholar]
• Khatkar M. S., Nicholas F. W., Collins A. R., Zenger K. R., Cavanagh J. A. L., Barris W., et al.. (2008). Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel. BMC Genomics 9:187. 10.1186/1471-2164-9-187, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Khatkar M. S., Zenger K. R., Hobbs M., Hawken R. J., Cavanagh J. A. L., Barris W., et al.. (2007). A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in Holstein-Friesian cattle. Genetics 176, 763–772. 10.1534/genetics.106.069369, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Kijas J., Elliot N., Kube P., Evans B., Botwright N., King H., et al.. (2017). Diversity and linkage disequilibrium in farmed Tasmanian Atlantic salmon. Anim. Genet. 48, 237–241. 10.1111/age.12513, PMID: [PubMed] [CrossRef] [Google Scholar]
• Lee S. H., Weerasinghe W. M. S. P., Wray N. R., Goddard M. E., Van Der Werf J. H. J. (2017). Using information of relatives in genomic prediction to apply effective stratified medicine. Sci. Rep. 7, 1–13. 10.1038/srep42091 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• López M. E., Neira R., Yáñez J. M. (2015). Applications in the search for genomic selection signatures in fish. Front. Genet. 5, 1–12. 10.3389/fgene.2014.00458 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Lu D., Sargolzaei M., Kelly M., Li C., Voort G. V., Wang Z., et al. (2012). Linkage disequilibrium in Angus, Charolais, and crossbred beef cattle. Front. Genet. 3, 1–10. 10.3389/fgene.2012.00152 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Makina S. O., Taylor J. F., Van Marle-Köster E., Muchadeyi F. C., Makgahlela M. L., MacNeil M. D., et al. (2015). Extent of linkage disequilibrium and effective population size in four south African sanga cattle breeds. Front. Genet. 6, 1–12. 10.3389/fgene.2015.00337 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Matukumalli L. K., Lawley C. T., Schnabel R. D., Taylor J. F., Allan M. F., Heaton M. P., et al.. (2009). Development and characterization of a high density SNP genotyping assay for cattle. PLoS One 4:e5350. 10.1371/journal.pone.0005350, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• McKay S. D., Schnabel R. D., Murdoch B. M., Matukumalli L. K., Aerts J., Coppieters W., et al.. (2007). Whole genome linkage disequilibrium maps in cattle. BMC Genet. 8:74. 10.1186/1471-2156-8-74, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Meuwissen T. H., Hayes B. J., Goddard M. E. (2001). Prediction of total genetic value using genome-wide dense marker maps. Genetics 157, 1819–1829. PMID: [PMC free article] [PubMed] [Google Scholar]
• Neira R., Lhorente J., Newman S., Yáñez J., Araneda M., Filp M. (2014). "Evolution of coho salmon (Oncorhynchus kisutch) breeding programs" in 10th world congress of genetics applied to livestock production, Vol. 472. [Google Scholar]
• Ødegård J., Moen T., Santi N., Korsvoll S. A., Kjøglum S., Meuwissen T. H. E. (2014). Genomic prediction in an admixed population of Atlantic salmon (Salmo salar). Front. Genet. 5, 1–8. 10.3389/fgene.2014.00402 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Pfaff C. L., Parra E. J., Bonilla C., Hiester K., McKeigue P. M., Kamboh M. I., et al.. (2001). Population structure in admixed populations: effect of admixture dynamics on the pattern of linkage disequilibrium. Am. J. Hum. Genet. 68, 198–207. 10.1086/316935, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Ponzoni R., Ling H., Nguyen N., Hamzah A. (2010). Inbreeding and effective population size in the Malaysian nucleus of the GIFT strain of Nile tilapia (Oreochromis niloticus). Aquaculture 302, 42–48. 10.1016/j.aquaculture.2010.02.009 [CrossRef] [Google Scholar]
• Porto-Neto L. R., Kijas J. W., Reverter A. (2014). The extent of linkage disequilibrium in beef cattle breeds using high-density SNP genotypes. Genet. Sel. Evol. 46, 1–5. 10.1186/1297-9686-46-22 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Pszczola M., Strabel T., Mulder H. A., Calus M. P. (2012). Reliability of direct genomic values for animals with different relationships within and to the reference population. J. Dairy Sci. 95, 389–400. 10.3168/jds.2011-4338, PMID: [PubMed] [CrossRef] [Google Scholar]
• Purcell S., Neale B., Todd-Brown K., Thomas L., Ferreira M. A. R., Bender D., et al.. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. 10.1086/519795, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Rexroad C. E., Vallejo R. L. (2009). Estimates of linkage disequilibrium and effective population size in rainbow trout. BMC Genet. 10:83. 10.1186/1471-2156-10-83 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• SalmonChile (2007). Historia en Chile. (Chile: ). Available at: http://www.salmonchile.cl/es/historia-en-chile.php#1990-2007 (Accessed February 20, 2017). [Google Scholar]
• Sargolzaei M., Schenkel F. S., Jansen G. B., Schaeffer L. R. (2008). Extent of linkage disequilibrium in Holstein cattle in North America. J. Dairy Sci. 91, 2106–2117. 10.3168/jds.2007-0553, PMID: [PubMed] [CrossRef] [Google Scholar]
• Saura M., Tenesa A., Woolliams J. A., Fernández A., Villanueva B. (2015). Evaluation of the linkage-disequilibrium method for the estimation of effective population size when generations overlap: an empirical case. BMC Genomics 16, 1–10. 10.1186/s12864-015-2167-z [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Teare M. D., Dunning A. M., Durocher F., Rennart G., Easton D. F. (2002). Sampling distribution of summary linkage disequilibrium measures. Ann. Hum. Genet. 66, 223–233. 10.1017/S0003480002001082, PMID: [PubMed] [CrossRef] [Google Scholar]
• Tenesa A., Navarro P., Hayes B. J., Duffy D. L., Clarke G. M., Goddard M. E., et al.. (2007). Recent human effective population size estimated from linkage disequilibrium. Genome Res. 17, 520–526. 10.1101/gr.6023607, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Tsai H. Y., Hamilton A., Tinch A. E., Guy D. R., Bron J. E., Taggart J. B., et al. (2016). Genomic prediction of host resistance to sea lice in farmed Atlantic salmon populations. Genet. Sel. Evol. 48, 1–11. 10.1186/s12711-016-0226-9 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Tsai H. Y., Hamilton A., Tinch A. E., Guy D. R., Gharbi K., Stear M. J., et al. (2015). Genome wide association and genomic prediction for growth traits in juvenile farmed Atlantic salmon using a high density SNP array. BMC Genomics 16d, 1–9. 10.1186/s12864-015-2117-9 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Vallejo R. L., Leeds T. D., Fragomeni B. O., Gao G., Hernandez A. G., Misztal I., et al. (2016). Evaluation of genome-enabled selection for bacterial cold water disease resistance using progeny performance data in rainbow trout: insights on genotyping methods and genomic prediction models. Front. Genet. 7, 1–13. 10.3389/fgene.2016.00096 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Vallejo R. L., Leeds T. D., Gao G., Parsons J. E., Martin K. E., Evenhuis J. P., et al. (2017). Genomic selection models double the accuracy of predicted breeding values for bacterial cold water disease resistance compared to a traditional pedigree-based model in rainbow trout aquaculture. Genet. Sel. Evol. 49:17. 10.1186/s12711-017-0293-6 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Vallejo R. L., Silva R. M. O., Evenhuis J. P., Gao G., Sixin L., Parsons J. E., et al.. (2018). Accurate genomic predictions for BCWD resistance in rainbow trout are achieved using low-density SNP panels: evidence that long-range LD is a major contributing factor. J. Anim. Breed. Genet. 135, 263–274. 10.1111/jbg.12335, PMID: [PubMed] [CrossRef] [Google Scholar]
• Villa-Angulo R., Matukumalli L. K., Gill C. A., Choi J., Van Tassell C. P., Grefenstette J. J. (2009). High-resolution haplotype block structure in the cattle genome. BMC Genet. 10, 1–13. 10.1186/1471-2156-10-19 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Vos P. G., Paulo M. J., Voorrips R. E., Visser R. G. F., van Eck H. J., van Eeuwijk F. A. (2017). Evaluation of LD decay and various LD-decay estimators in simulated and SNP-array data of tetraploid potato. Theor. Appl. Genet. 130, 123–135. 10.1007/s00122-016-2798-8, PMID: [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Wall J. D., Pritchard J. K. (2003). Haplotype blocks and linkage disequilibrium in the human genome. Nat. Rev. Genet. 4, 587–597. 10.1038/nrg1123 [PubMed] [CrossRef] [Google Scholar]
• Waples R. S., Do C. (2008). LDNE: a program for estimating effective population size from data on linkage disequilibrium. Mol. Ecol. Resour. 8, 753–756. 10.1111/j.1755-0998.2007.02061.x, PMID: [PubMed] [CrossRef] [Google Scholar]
• Yáñez J. M., Bangera R., Lhorente J. P., Barría A., Oyarzún M., Neira R., et al. (2016). Negative genetic correlation between resistance against Piscirickettsia salmonis and harvest weight in coho salmon (Oncorhynchus kisutch). Aquaculture 459, 8–13. 10.1016/j.aquaculture.2016.03.020 [CrossRef] [Google Scholar]
• Yáñez J. M., Bassini L. N., Filp M., Lhorente J. P., Ponzoni R. W., Neira R. (2014). Inbreeding and effective population size in a coho salmon (Oncorhynchus kisutch) breeding nucleus in Chile. Aquaculture 420–421, S15–S19. 10.1016/j.aquaculture.2013.05.028 [CrossRef] [Google Scholar]
• Yoshida G. M., Bangera R., Carvalheiro R., Correa K., Figueroa R., Lhorente J. P., et al. (2018). Genomic prediction accuracy for resistance against Piscirickettsia salmonis in farmed rainbow trout. G3 Genes Genomes Genet. 8, 719–726. 10.1534/g3.117.300499 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Yvonne W., Piter B., Roel V., Mario C. (2016). An equation to predict the accuracy of genomic values by combining data from multiple traits, populations, or environments. Genetics 202, 799–823. 10.1534/genetics.115.183269 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
• Zhao H., Nettleton D., Zoller M., Dekkers J. C. M. (2005). Evaluation of linkage disequilibrium measures between multi-allelic markers as predictors of linkage disequilibrium between markers and QTL. Genet. Res. 86, 77–87. 10.1017/S001667230500769X, PMID: [PubMed] [CrossRef] [Google Scholar]

---

Articles from Frontiers in Genetics are provided here courtesy of Frontiers Media SA

---

wadestalwasoming.blogspot.com

Source: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6539196/

Share this post