The application of genome-wide SNP genotyping methods in studies on livestock genomes

Comparative assessment of the effective population size and linkage disequilibrium of Karan Fries cattle revealed viable population dynamics

OBJECTIVE

Karan Fries (KF), a high-producing composite cattle was developed through crossing indicine Tharparkar cows with taurine bulls (Holstein Friesian, Brown Swiss, and Jersey), to increase the milk yield across India. This composite cattle population must maintain sufficient genetic diversity for long-term development and breed improvement in the coming years. The level of linkage disequilibrium (LD) measures the influence of population genetic forces on the genomic structure and provides insights into the evolutionary history of populations, while the decay of LD is important in understanding the limits of genome-wide association studies for a population. Effective population size (Ne) which is genomically based on LD accumulated over the course of previous generations, is a valuable tool for e valuation of the genetic diversity and level of inbreeding. The present study was undertaken to understand KF population dynamics through the estimation of Ne and LD for the longterm sustainability of these breeds.

METHODS

The present study included 96 KF samples genotyped using Illumina HDBovine array to estimate the effective population and examine the LD pattern. The genotype data were also obtained for other crossbreds (Santa Gertrudis, Brangus, and Beefmaster) and Holstein Friesian cattle for comparison purposes.

RESULTS

The average LD between single nucleotide polymorphisms (SNPs) was r2 = 0.13 in the present study. LD decay (r2 = 0.2) was observed at 40 kb inter-marker distance, indicating a panel with 62,765 SNPs was sufficient for genomic breeding value estimation in KF cattle. The pedigree-based Ne of KF was determined to be 78, while the Ne estimates obtained using LD-based methods were 52 (SNeP) and 219 (genetic optimization for Ne estimation), respectively.

CONCLUSION

KF cattle have an Ne exceeding the FAO's minimum recommended level of 50, which was desirable. The study also revealed significant population dynamics of KF cattle and increased our understanding of devising suitable breeding strategies for longterm sustainable development.

Assessment of runs of homozygosity, heterozygosity-rich regions and genomic inbreeding estimates in a subpopulation of Guzerá (Bos indicus) dual-purpose cattle

For decades, inbreeding in cattle has been evaluated using pedigree information. Nowadays, inbreeding coefficients can be obtained using genomic information such as runs of homozygosity (ROH). The aims of this study were to quantify ROH and heterozygosity-rich regions (HRR) in a subpopulation of Guzerá dual-purpose cattle, to examine ROH and HRR islands, and to compare inbreeding coefficients obtained by ROH with alternative genomic inbreeding coefficients. A subpopulation of 1733 Guzerá animals genotyped for 50k SNPs was used to obtain the ROH and HRR segments. Inbreeding coefficients by ROH (FROH ), by genomic relationship matrix based on VanRaden's method 1 using reference allele frequency in the population (FGRM ), by genomic relationship matrix based on VanRaden's method 1 using allele frequency fixed in 0.5 (FGRM_0.5 ), and by the proportion of homozygous loci (FHOM ) were calculated. A total of 15,660 ROH were identified, and the chromosome with the highest number of ROH was BTA6. A total of 4843 HRRs were identified, and the chromosome with the highest number of HRRs was BTA23. No ROH and HRR islands were identified according to established criteria, but the regions closest to the definition of an island were examined from 64 to 67 Mb of BTA6, from 36 to 37 Mb of BTA2 and from 0.50 to 1.25 Mb of BTA23. The genes identified in ROH islands have previously been associated with dairy and beef traits, while genes identified on HRR islands have previously been associated with reproductive traits and disease resistance. FROH was equal to 0.095 ± 0.084, and its Spearman correlation with FGRM was low (0.44) and moderate-high with FHOM (0.79) and with FGRM_0.5 (0.80). The inbreeding coefficients determined by ROH were higher than other cattle breeds' and higher than pedigree-based inbreeding in the Guzerá breed obtained in previous studies. It is recommended that future studies investigate the effects of inbreeding determined by ROH on the traits under selection in the subpopulation studied.

Estimation of Linkage Disequilibrium, Effective Population Size, and Genetic Parameters of Phenotypic Traits in Dabieshan Cattle

Dabieshan cattle (DBSC) are a valuable genetic resource for indigenous cattle breeds in China. It is a small to medium-sized breed with slower growth, but with good meat quality and fat deposition. Genetic markers could be used for the estimation of population genetic structure and genetic parameters. In this work, we genotyped the DBSC breeding population (n = 235) with the GeneSeek Genomic Profiler (GGP) 100 k density genomic chip. Genotype data of 222 individuals and 81,579 SNPs were retained after quality control. The average minor allele frequency (MAF) was 0.20 and the average linkage disequilibrium (LD) level (r²) was 0.67 at a distance of 0-50 Kb. The estimated relationship coefficient and effective population size (Ne) were 0.023 and 86 for the current generation. In addition, we used genotype data to estimate the genetic parameters of the population's phenotypic traits. Among them, height at hip cross (HHC) and shin circumference (SC) were rather high heritability traits, with heritability of 0.41 and 0.54, respectively. The results reflected the current cattle population's extent of inbreeding and history. Through the principal breeding parameters, genomic breeding would significantly improve the genetic progress of breeding.

Screening Discriminating SNPs for Chinese Indigenous Pig Breeds Identification Using a Random Forests Algorithm

Chinese indigenous pig breeds have unique genetic characteristics and a rich diversity; however, effective breed identification methods have not yet been well established. In this study, a genotype file of 62,822 single-nucleotide polymorphisms (SNPs), which were obtained from 1059 individuals of 18 Chinese indigenous pig breeds and 5 cosmopolitan breeds, were used to screen the discriminating SNPs for pig breed identification. After linkage disequilibrium (LD) pruning filtering, this study excluded 396 SNPs on non-constant chromosomes and retained 20.92~-27.84% of SNPs for each of the 18 autosomes, leaving a total of 14,823 SNPs. The principal component analysis (PCA) showed the largest differences between cosmopolitan and Chinese pig breeds (PC1 = 10.452%), while relatively small differences were found among the 18 indigenous pig breeds from the Yangtze River Delta region of China. Next, a random forest (RF) algorithm was used to filter these SNPs and obtain the optimal number of decision trees (ntree = 1000) using corresponding out-of-bag (OOB) error rates. By comparing two different SNP ranking methods in the RF analysis, the mean decreasing accuracy (MDA) and mean decreasing Gini index (MDG), the effects of panels with different numbers of SNPs on the assignment accuracy, and the statistics of SNP distribution on each chromosome in the panels, a panel of 1000 of the most breed-discriminative tagged SNPs were finally selected based on the MDA screening method. A high accuracy (>99.3%) was obtained by the breed prediction of 318 samples in the RF test set; thus, a machine learning classification method was established for the multi-breed identification of Chinese indigenous pigs based on a low-density panel of SNPs.

Feedomics provides bidirectional omics strategies between genetics and nutrition for improved production in cattle

Increasing the efficiency and sustainability of cattle production is an effective way to produce valuable animal proteins for a growing human population. Genetics and nutrition are the 2 major research topics in selecting cattle with beneficial phenotypes and developing genetic potentials for improved performance. There is an inextricable link between genetics and nutrition, which urgently requires researchers to uncover the underlying molecular mechanisms to optimize cattle production. Feedomics integrates a range of omic techniques to reveal the mechanisms at different molecular levels related to animal production and health, which can provide novel insights into the relationships of genes and nutrition/nutrients. In this review, we summarized the applications of feedomics techniques to reveal the effect of genetic elements on the response to nutrition and investigate how nutrients affect the functional genome of cattle from the perspective of both nutrigenetics and nutrigenomics. We highlighted the roles of rumen microbiome in the interactions between host genes and nutrition. Herein, we discuss the importance of feedomics in cattle nutrition research, with a view to ensure that cattle exhibit the best production traits for human consumption from both genetic and nutritional aspects.

Targeted genome-wide SNP genotyping in feral horses using non-invasive fecal swabs

The development of high-throughput sequencing has prompted a transition in wildlife genetics from using microsatellites toward sets of single nucleotide polymorphisms (SNPs). However, genotyping large numbers of targeted SNPs using non-invasive samples remains challenging due to relatively large DNA input requirements. Recently, target enrichment has emerged as a promising approach requiring little template DNA. We assessed the efficacy of Tecan Genomics’ Allegro Targeted Genotyping (ATG) for generating genome-wide SNP data in feral horses using DNA isolated from fecal swabs. Total and host-specific DNA were quantified for 989 samples collected as part of a long-term individual-based study of feral horses on Sable Island, Nova Scotia, Canada, using dsDNA fluorescence and a host-specific qPCR assay, respectively. Forty-eight samples representing 44 individuals containing at least 10 ng of host DNA (ATG’s recommended minimum input) were genotyped using a custom multiplex panel targeting 279 SNPs. Genotyping accuracy and consistency were assessed by contrasting ATG genotypes with those obtained from the same individuals with SNP microarrays, and from multiple samples from the same horse, respectively. 62% of swabs yielded the minimum recommended amount of host DNA for ATG. Ignoring samples that failed to amplify, ATG recovered an average of 88.8% targeted sites per sample, while genotype concordance between ATG and SNP microarrays was 98.5%. The repeatability of genotypes from the same individual approached unity with an average of 99.9%. This study demonstrates the suitability of ATG for genome-wide, non-invasive targeted SNP genotyping, and will facilitate further ecological and conservation genetics research in equids and related species.

Genetic structure analysis in several populations of cattle using SNP genotypes

Parameters such as effective population size (Ne), runs of homozygosity (ROH), and inbreeding based on ROH (FROH) can give new insight into the level of genetic diversity for the population under selection. This research aimed to measure the extent of linkage disequilibrium (LD), effective population size (Ne), Haplotype Block Structure, and runs of homozygosity (ROHs) in several populations of cattle using SNP genotypes. In this study, that the average r² decreased with the increasing distance of SNP pairs. A general decrease in Ne can be seen for all four populations, indicating a loss of genetic diversity. The Iranian Holstein had the lowest level of genomic inbreeding at an ROH of 1, 5, 10 Mb, while the French Holstein had the highest. The maximum number of ROH is seen at a distance of less than 1 Mb, and the lowest number of ROH is seen at a distance of 10 Mb. The number of ROH decreases with increasing distance due to the increased recombination rate. This is a concern as an increase in inbreeding leads to a reduction in the effective population size, which was also evident in the study populations.

Genome-wide analysis of homozygosity regions in european simmental bulls

The study of Runs of Homozygosity (ROH) is a useful approach for the characterization of the genome of livestock populations. Due to their high relationship with autozygosity, ROH allow to make inference about population genetic history, to estimate the level of inbreeding, to assess within breed heterogeneity and to detect the footprints of selection on livestock genomes. Aim of this study was to investigate the distribution of runs of homozygosity in bulls belonging to five European Simmental populations and to assess the relationship between three production traits (milk yield, fat and protein contents) and autozygosity. ROH count, distribution and ROH-based coefficient of inbreeding (F_ROH ) were calculated for 3,845 Simmental bulls of five different European countries: Austria (AT), Switzerland (CH), Czech Republic (CZ), Germany (DE) and Italy (IT). Average values of ROH number per animal, and total genome length covered by ROH were 77.8 ± 20.7 and 205 ± 74.4 Mb, respectively. Bulls from AT, DE and IT exhibited similar ROH characteristics. Swiss animals showed the highest (12.6%), while CZ the lowest (4.6%) F_ROH coefficient. The relationship between ROH occurrence and milk production traits was investigated through a genome-wide ROH-traits association analysis (GWRA). A total of 34 regions previously associated with milk traits (yield and/or composition) were identified by GWRA. Results of the present research highlight a mixed genetic background in the 5 European Simmental populations, with the possible presence of three subgroups. Moreover, a strong relationship between autozygosity and production traits has been detected.

Reduced representation approach for identification of genome-wide SNPs and their annotation for economically important traits in Indian Tharparkar cattle

The present study was carried out in Tharparkar cattle for identification of genome-wide SNPs and microsatellites, and then annotate the identified high-quality SNPs to milk production, fertility, carcass, adaptability and immune response of economically important traits. A total of 146,011 SNPs were identified with respect to Bos taurus reference genome which are indicus specific, out of which 10,519 SNPs were found to be novel. Similarly, a total of 87,047 SNPs were identified with respect to Bos indicus reference genome. After final annotation of SNPs identified with respect to Bos indicus reference genome, 2871 SNPs were found to be associated in 383 candidate genes having to do with milk production, fertility, carcass, immune response and adaptability traits. Following that, 2571 microsatellites were identified. The information mined from the data might be of importance for the future breed improvement programs, conservation efforts and for enhancing the SNPs density of the existing bovine SNP chips.

Haplotype-Based Genome-Wide Association Study and Identification of Candidate Genes Associated with Carcass Traits in Hanwoo Cattle

Hanwoo, is the most popular native beef cattle in South Korea. Due to its extensive popularity, research is ongoing to enhance its carcass quality and marbling traits. In this study we conducted a haplotype-based genome-wide association study (GWAS) by constructing haplotype blocks by three methods: number of single nucleotide polymorphisms (SNPs) in a haplotype block (nsnp), length of genomic region in kb (Len) and linkage disequilibrium (LD). Significant haplotype blocks and genes associated with them were identified for carcass traits such as BFT (back fat thickness), EMA (eye Muscle area), CWT (carcass weight) and MS (marbling score). Gene-set enrichment analysis and functional annotation of genes in the significantly-associated loci revealed candidate genes, including PLCB1 and PLCB4 present on BTA13, coding for phospholipases, which might be important candidates for increasing fat deposition due to their role in lipid metabolism and adipogenesis. CEL (carboxyl ester lipase), a bile-salt activated lipase, responsible for lipid catabolic process was also identified within the significantly-associated haplotype block on BTA11. The results were validated in a different Hanwoo population. The genes and pathways identified in this study may serve as good candidates for improving carcass traits in Hanwoo cattle.

A detailed characteristics of bias associated with long runs of homozygosity identification based on medium density SNP microarrays

In the present study, runs of homozygosity (ROH) detected with the use of a standard bovine 54k single nucleotide polymorphism (SNP) genotyping assay and two different ROH detection approaches, based on 50 (M1) or 15 (M2) consecutive SNPs, were compared with results of whole genome sequencing. Both microarray-based methods accurately recognised medium-sized ROH, however, it was found that M2 method seemed to better than M1 identify short ROH, but highly overestimated their number, leading to numerous false positive calls. Moreover, long ROH identified with microarray data tended to break into shorter segments in sequencing data because of the presence of regions with high heterozygosity within the ROH sequences. This may indicate, that these long ROH are formed by closely positioned shorter homozygous segments that may be of older origin or may be created by two similar but not identical haplotypes, showing minor internal recombination signs. Such finding also suggests that at least some of the results of previous studies in regard to long ROH may be biased leading to inaccurate estimations of genomes autozygosity via ROH classification into length categories.

Evaluation of Linkage Disequilibrium, Effective Population Size and Haplotype Block Structure in Chinese Cattle

Simple Summary

Evaluation of the population structure and linkage disequilibrium can offer important insights to fully understand the genetic diversity and population history of cattle, which can enable us to appropriately design and implement GWAS and GS in cattle. In this study, we characterized the extent of genome-wide LD and the haplotype block structure, and estimated the persistence of phase of Chinese indigenous cattle with Illumina BovineHD BeadChip. According to our study, 58K, 87K, 95K, 52K, and 52K markers would be necessary for SCHC, NCC, SWC, SIM, and WAG, respectively, in the implementation of GWAS and GS and combining a multipopulation with high persistence of phase is feasible for the implication of genomic selection for Chinese beef cattle.

Abstract

Understanding the linkage disequilibrium (LD) across the genome, haplotype structure, and persistence of phase between breeds can enable us to appropriately design and implement the genome-wide association (GWAS) and genomic selection (GS) in beef cattle. We estimated the extent of genome-wide LD, haplotype block structure, and the persistence of phase in 10 Chinese cattle population using high density BovinHD BeadChip. The overall LD measured by r² between adjacent SNPs were 0.60, 0.67, 0.58, 0.73, and 0.71 for South Chinese cattle (SCHC), North Chinese cattle (NCC), Southwest Chinese cattle (SWC), Simmental (SIM), and Wagyu (WAG). The highest correlation (0.53) for persistence of phase across groups was observed for SCHC vs. SWC at distances of 0–50 kb, while the lowest correlation was 0.13 for SIM vs. SCHC at the same distances. In addition, the estimated current effective population sizes were 27, 14, 31, 34, and 43 for SCHC, NCC, SWC, SIM, and WAG, respectively. Our result showed that 58K, 87K, 95K, 52K, and 52K markers were required for implementation of GWAS and GS in SCHC, NCC, SWC, SIM, and WAG, respectively. Also, our findings suggested that the implication of genomic selection for multipopulation with high persistence of phase is feasible for Chinese cattle.

Breeding for resistance to gastrointestinal nematodes - the potential in low-input/output small ruminant production systems

The control of gastrointestinal nematodes (GIN) is mainly based on the use of drugs, grazing management, use of copper oxide wire particles and bioactive forages. Resistance to anthelmintic drugs in small ruminants is documented worldwide. Host genetic resistance to parasites, has been increasingly used as a complementary control strategy, along with the conventional intervention methods mentioned above. Genetic diversity in resistance to GIN has been well studied in experimental and commercial flocks in temperate climates and more developed economies. However, there are very few report outputs from the more extensive low-input/output smallholder systems in developing and emerging countries. Furthermore, results on quantitative trait loci (QTL) associated with nematode resistance from various studies have not always been consistent, mainly due to the different nematodes studied, different host breeds, ages, climates, natural infections versus artificial challenges, infection level at sampling periods, among others. The increasing use of genetic markers (Single Nucleotide Polymorphisms, SNPs) in GWAS or the use of whole genome sequence data and a plethora of analytic methods offer the potential to identify loci or regions associated nematode resistance. Genomic selection as a genome-wide level method overcomes the need to identify candidate genes. Benefits in genomic selection are now being realised in dairy cattle and sheep under commercial settings in the more advanced countries. However, despite the commercial benefits of using these tools, there are practical problems associated with incorporating the use of marker-assisted selection or genomic selection in low-input/output smallholder farming systems breeding schemes. Unlike anthelmintic resistance, there is no empirical evidence suggesting that nematodes will evolve rapidly in response to resistant hosts. The strategy of nematode control has evolved to a more practical manipulation of host-parasite equilibrium in grazing systems by implementation of various strategies, in which improvement of genetic resistance of small ruminant should be included. Therefore, selection for resistant hosts can be considered as one of the sustainable control strategy, although it will be most effective when used to complement other control strategies such as grazing management and improving efficiency of anthelmintics currently.

Concepts and relevance of genome-wide association studies

Genome-wide association study (GWAS) technology has been a primary method for identifying the genes responsible for diseases and other traits for the past 10 years. Over 2,000 human GWAS reports now appear in the scientific journals. The technology is continuing to improve, and has recently become accessible to researchers studying a wide variety of animals, plants and model organisms. Here, we present an overview of GWAS concepts: the underlying biology, the origins of the method, and the primary components of a GWAS experiment.

Identification of genome-wide selection signatures in the Limousin beef cattle breed

The study is aimed at identifying selection footprints within the genome of Limousin cattle. With the use of Extended Haplotype Homozygosity test, supplemented with correction for variation in recombination rates across the genome, we created map of selection footprints and detected 173 significant (p < 0.01) core haplotypes being potentially under positive selection. Within these regions, a number of candidate genes associated inter alia with skeletal muscle growth (GDF15, BMP7, BMP4 and TGFB3) or postmortem proteolysis and meat maturation (CAPN1 and CAPN5) were annotated. Noticeable clusters of selection footprints were detected on chromosomes 1, 4, 8 and 14, which are known to carry several quantitative trait loci for growth traits and meat quality. The study provides information about the genes and metabolic pathways potentially modified under the influence of directional selection, aimed at improving beef production characteristics in Limousin cattle.

Copy number variations in high and low fertility breeding boars

Background

In this study we applied the extreme groups/selective genotyping approach for identifying copy number variations in high and low fertility breeding boars. The fertility indicator was the calculated Direct Boar Effect on litter size (DBE) that was obtained as a by-product of the national genetic evaluation for litter size (BLUP). The two groups of animals had DBE values at the upper (high fertility) and lower (low fertility) end of the distribution from a population of more than 38,000 boars. Animals from these two diverse phenotypes were genotyped with the Porcine SNP60K chip and compared by several approaches in order to prove the feasibility of our CNV analysis and to identify putative markers of fertility.

Results

We have identified 35 CNVRs covering 36.5 Mb or ~1.3% of the porcine genome. Among these 35 CNVRs, 14 were specific to the high fertility group, while 19 CNVRs were specific to the low fertility group which overlap with 137 QTLs of various reproductive traits. The identified 35 CNVRs encompassed 50 genes, among them 40 were specific to the low fertility group, seven to the high fertility group, while three were found in regions that were present in both groups but with opposite gain/loss status. A functional analysis of several databases revealed that the genes found in CNVRs from the low fertility group have been significantly enriched in members of the innate immune system, Toll-like receptor and RIG-I-like receptor signaling and fatty acid oxidation pathways.

Conclusions

We have demonstrated that our analysis pipeline could identify putative CNV markers of fertility, especially in case of low fertility boars.

Electronic supplementary material

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