Genomic Selection Strategies Research Articles

PSIV-37 Development of a genomic selection strategy to include meat quality traits in turkeys (Meleagris gallopavo).

PSIV-37 Development of a genomic selection strategy to include meat quality traits in turkeys (Meleagris gallopavo).

Stable Quantitative Resistance Loci to Blackleg Disease in Canola (Brassica napus L.) Over Continents.

The hemibiotrophic fungus, Leptosphaeria maculans is the most devastating pathogen, causing blackleg disease in canola (Brassica napus L). To study the genomic regions involved in quantitative resistance (QR), 259–276 DH lines from Darmor-bzh/Yudal (DYDH) population were assessed for resistance to blackleg under shade house and field conditions across 3 years. In different experiments, the broad sense heritability varied from 43 to 95%. A total of 27 significant quantitative trait loci (QTL) for QR were detected on 12 chromosomes and explained between 2.14 and 10.13% of the genotypic variance. Of the significant QTL, at least seven were repeatedly detected across different experiments on chromosomes A02, A07, A09, A10, C01, and C09. Resistance alleles were mainly contributed by ‘Darmor-bzh’ but ‘Yudal’ also contributed few of them. Our results suggest that plant maturity and plant height may have a pleiotropic effect on QR in our conditions. We confirmed that Rlm9 which is present in ‘Darmor-bzh’ is not effective to confer resistance in our Australian field conditions. Comparative mapping showed that several R genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors map in close proximity (within 200 Kb) of the significant trait-marker associations on the reference ‘Darmor-bzh’ genome assembly. More importantly, eight significant QTL regions were detected across diverse growing environments: Australia, France, and United Kingdom. These stable QTL identified herein can be utilized for enhancing QR in elite canola germplasm via marker- assisted or genomic selection strategies.

Karyomapping for simultaneous genomic evaluation and aneuploidy screening of preimplantation bovine embryos: The first live-born calves

In cattle breeding, the development of genomic selection strategies based on single nucleotide polymorphism (SNP) interrogation has led to improved rates of genetic gain. Additionally, the application of genomic selection to in-vitro produced (IVP) embryos is expected to bring further benefits thanks to the ability to test a greater number of individuals before establishing a pregnancy and to ensure only carriers of desirable traits are born. However, aneuploidy, a leading cause of developmental arrest, is known to be common in IVP embryos. Karyomapping is a comprehensive screening test based on SNP typing that can be used for simultaneous genomic selection and aneuploidy detection, offering the potential to maximize pregnancy rates. Moreover, Karyomapping can be used to characterize the frequency and parental origin of aneuploidy in bovine IVP embryos, which have remained underexplored to date. Here, we report the use of Karyomapping to characterize the frequency and parental origin of aneuploidy in IVP bovine embryos in order to establish an estimate of total aneuploidy rates in each parental germline. We report an estimate of genome wide recombination rate in cattle and demonstrate, for the first time, a proof of principle for the application of Karyomapping to cattle breeding, with the birth of five calves after screening. This combined genomic selection and aneuploidy screening approach was highly reliable, with calves showing 98% concordance with their respective embryo biopsies for SNP typing and 100% concordance with their respective biopsies for aneuploidy screening. This approach has the potential to simultaneously improve pregnancy rates following embryo transfer and the rate of genetic gain in cattle breeding, and is applicable to basic research to investigate meiosis and aneuploidy.

Effects of genomic selection for intramuscular fat content in breast muscle in Chinese local chickens.

Improvements in living standards have resulted in consumers having higher expectations for chicken meat quality. This is particularly true in Asia, where there is high consumer preference for local breeds. Nothing is presently known about the effectiveness of using genomic selection (GS) strategies in chickens to genetically improve meat quality traits that cannot be measured in living potential parents. In this study, 724 Beijing-You chickens were used as a training population; all were genotyped using Illumina 60K SNP chips, and intramuscular fat content in breast muscle (IMFbr ) was measured. Birds in the GS line were selected based on genomic estimated breeding values, IMFbr being the sole trait. Genetic progress in one generation was compared to that from conventional family-based selection, and both were evaluated against random-bred controls. Results showed that relative to the random-bred controls, IMF percentage was improved 9.62% using GS, comparable to the 10.38% improvement using family-based selection. We quantified the effectiveness of GS when applied to a meat quality trait with low heritability in chickens. We plan to introduce custom SNP chips, appropriate for native chicken breeds in China, to assist in applying GS in local breeding and accelerate genetic gain.

Benefits of testing in both bio-secure and production environments in genomic selection breeding programs for commercial broiler chicken

BackgroundA breeding program for commercial broiler chicken that is carried out under strict biosecure conditions can show reduced genetic gain due to genotype by environment interactions (G × E) between bio-secure (B) and commercial production (C) environments. Accuracy of phenotype-based best linear unbiased prediction of breeding values of selection candidates using sib-testing in C is low. Genomic prediction based on dense genetic markers may improve accuracy of selection. Stochastic simulation was used to explore the benefits of genomic selection in breeding schemes for broiler chicken that include birds in both B and C for assessment of phenotype.ResultsWhen genetic correlations (r_{g}) between traits measured in B and C were equal to 0.5 and 0.7, breeding schemes with 15, 30 and 45% of birds assessed in C resulted in higher genetic gain for performance in C compared to those without birds in C. The optimal proportion of birds phenotyped in C for genetic gain was 30%. When the proportion of birds in C was optimal and genotyping effort was limited, allocating 30% of the genotyping effort to birds in C was also the optimal genotyping strategy for genetic gain. When r_{g} was equal to 0.9, genetic gain for performance in C was not improved with birds in C compared to schemes without birds in C. Increasing the heritability of traits assessed in C increased genetic gain significantly. Rates of inbreeding decreased when the proportion of birds in C increased because of a lower selection intensity among birds retained in B and a reduction in the probability of co-selecting close relatives.ConclusionsIf G × E interactions (r_{g} of 0.5 and 0.7) are strong, a genomic selection scheme in which 30% of the birds hatched are phenotyped in C has larger genetic gain for performance in C compared to phenotyping all birds in B. Rates of inbreeding decreased as the proportion of birds moved to C increased from 15 to 45%.

Genomic Selection for Growth Traits in Pacific Oyster (Crassostrea gigas): Potential of Low-Density Marker Panels for Breeding Value Prediction.

Pacific oysters are a key aquaculture species globally, and genetic improvement via selective breeding is a major target. Genomic selection has the potential to expedite genetic gain for key target traits of a breeding program, but has not yet been evaluated in oyster. The recent development of SNP arrays for Pacific oyster (Crassostrea gigas) raises the opportunity to test genomic selection strategies for polygenic traits. In this study, a population of 820 oysters (comprising 23 full-sibling families) were genotyped using a medium density SNP array (23 K informative SNPs), and the genetic architecture of growth-related traits [shell height (SH), shell length (SL), and wet weight (WW)] was evaluated. Heritability was estimated to be moderate for the three traits (0.26 ± 0.06 for SH, 0.23 ± 0.06 for SL and 0.35 ± 0.05 for WW), and results of a GWAS indicated that the underlying genetic architecture was polygenic. Genomic prediction approaches were used to estimate breeding values for growth, and compared to pedigree based approaches. The accuracy of the genomic prediction models (GBLUP) outperformed the traditional pedigree approach (PBLUP) by ∼25% for SL and WW, and ∼30% for SH. Further, reduction in SNP marker density had little impact on prediction accuracy, even when density was reduced to a few hundred SNPs. These results suggest that the use of genomic selection in oyster breeding could offer benefits for the selection of breeding candidates to improve complex economic traits at relatively modest cost.

Phenotypic Data from Inbred Parents Can Improve Genomic Prediction in Pearl Millet Hybrids.

Pearl millet is a non-model grain and fodder crop adapted to extremely hot and dry environments globally. In India, a great deal of public and private sectors’ investment has focused on developing pearl millet single cross hybrids based on the cytoplasmic-genetic male sterility (CMS) system, while in Africa most pearl millet production relies on open pollinated varieties. Pearl millet lines were phenotyped for both the inbred parents and hybrids stage. Many breeding efforts focus on phenotypic selection of inbred parents to generate improved parental lines and hybrids. This study evaluated two genotyping techniques and four genomic selection schemes in pearl millet. Despite the fact that 6× more sequencing data were generated per sample for RAD-seq than for tGBS, tGBS yielded more than 2× as many informative SNPs (defined as those having MAF > 0.05) than RAD-seq. A genomic prediction scheme utilizing only data from hybrids generated prediction accuracies (median) ranging from 0.73-0.74 (1000-grain weight), 0.87-0.89 (days to flowering time), 0.48-0.51 (grain yield) and 0.72-0.73 (plant height). For traits with little to no heterosis, hybrid only and hybrid/inbred prediction schemes performed almost equivalently. For traits with significant mid-parent heterosis, the direct inclusion of phenotypic data from inbred lines significantly (P < 0.05) reduced prediction accuracy when all lines were analyzed together. However, when inbreds and hybrid trait values were both scored relative to the mean trait values for the respective populations, the inclusion of inbred phenotypic datasets moderately improved genomic predictions of the hybrid genomic estimated breeding values. Here we show that modern approaches to genotyping by sequencing can enable genomic selection in pearl millet. While historical pearl millet breeding records include a wealth of phenotypic data from inbred lines, we demonstrate that the naive incorporation of this data into a hybrid breeding program can reduce prediction accuracy, while controlling for the effects of heterosis per se allowed inbred genotype and trait data to improve the accuracy of genomic estimated breeding values for pearl millet hybrids.

Training genomic selection models across several breeding cycles increases genetic gain in oil palm in silico study

Genomic selection (GS) is expected to increase the rate of genetic gain in oil palm. In a GS scheme, breeding cycles with progeny tests (phenotypic selection, PS) used to calibrate the GS predictive model and for selection alternate with GS cycles, making it possible to train the GS model with aggregated data from several cycles. To evaluate this possibility, we simulated four cycles of hybrid breeding for bunch production and compared two methods of calibrating the GS model, one using aggregated data from the two most recent cycles (Tr2Gen), the other using data from the last cycle (Tr1Gen). We also compared a GS scheme with two PS cycles and two GS cycles (2PT-2noPT), and a scheme with PS every other cycle and GS otherwise (PT-noPT). We showed that Tr2Gen had a 10.7% higher genetic gain per cycle than Tr1Gen, mostly due to increased selection accuracy, particularly in across-cycle selection, despite the decreased relationship between training individuals and selection candidates. After four cycles, Tr2Gen had a 5% higher cumulative genetic gain than Tr1Gen, with a lower coefficient of variation. PT-noPT benefited more from the advantages offered by Tr2Gen than 2PT-2noPT. Over four breeding cycles, combining PT-noPT and Tr2Gen largely outperformed conventional reciprocal recurrent selection (RRS), with an increase in annual genetic gain ranging from 37.6 to 57.5%, depending on the number of GS candidates. This study confirms the advantages of GS over RRS and indicated that oil palm breeders should train GS models using all data from past breeding cycles.

Genomic relationships reveal significant dominance effects for growth in hybrid Eucalyptus

Non-additive genetic effects can be effectively exploited in control-pollinated families with the availability of genome-wide markers. We used 41,304 SNP markers and compared pedigree vs. marker-based genetic models by analysing height, diameter, basic density and pulp yield for Eucalyptus urophylla × E.grandis control-pollinated families represented by 949 informative individuals. We evaluated models accounting for additive, dominance, and first-order epistatic interactions (additive by additive, dominance by dominance, and additive by dominance). We showed that the models can capture a large proportion of the genetic variance from dominance and epistasis for growth traits as those components are typically not independent. We also showed that we could partition genetic variances more precisely when using relationship matrices derived from markers compared to using only pedigree information. In addition, phenotypic prediction accuracies were only slightly increased by including dominance effects for growth traits since estimates of non-additive variances yielded rather high standard errors. This novel result improves our current understanding of the architecture of quantitative traits and recommends accounting for dominance variance when developing genomic selection strategies in hybrid Eucalyptus.

Two‐Part Genomic Selection Strategy

CSA NewsVolume 62, Issue 10 p. 12-12 Society Science Two-Part Genomic Selection Strategy First published: 05 October 2017 https://doi.org/10.2134/csa2017.62.1009Read the full textAboutPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onFacebookTwitterLinkedInRedditWechat No abstract is available for this article. Volume62, Issue10October 2017Pages 12-12 RelatedInformation

Exploring and Harnessing Haplotype Diversity to Improve Yield Stability in Crops.

In order to meet future food, feed, fiber, and bioenergy demands, global yields of all major crops need to be increased significantly. At the same time, the increasing frequency of extreme weather events such as heat and drought necessitates improvements in the environmental resilience of modern crop cultivars. Achieving sustainably increase yields implies rapid improvement of quantitative traits with a very complex genetic architecture and strong environmental interaction. Latest advances in genome analysis technologies today provide molecular information at an ultrahigh resolution, revolutionizing crop genomic research, and paving the way for advanced quantitative genetic approaches. These include highly detailed assessment of population structure and genotypic diversity, facilitating the identification of selective sweeps and signatures of directional selection, dissection of genetic variants that underlie important agronomic traits, and genomic selection (GS) strategies that not only consider major-effect genes. Single-nucleotide polymorphism (SNP) markers today represent the genotyping system of choice for crop genetic studies because they occur abundantly in plant genomes and are easy to detect. SNPs are typically biallelic, however, hence their information content compared to multiallelic markers is low, limiting the resolution at which SNP–trait relationships can be delineated. An efficient way to overcome this limitation is to construct haplotypes based on linkage disequilibrium, one of the most important features influencing genetic analyses of crop genomes. Here, we give an overview of the latest advances in genomics-based haplotype analyses in crops, highlighting their importance in the context of polyploidy and genome evolution, linkage drag, and co-selection. We provide examples of how haplotype analyses can complement well-established quantitative genetics frameworks, such as quantitative trait analysis and GS, ultimately providing an effective tool to equip modern crops with environment-tailored characteristics.

An innovative procedure of genome-wide association analysis fits studies on germplasm population and plant breeding.

The innovative RTM-GWAS procedure provides a relatively thorough detection of QTL and their multiple alleles for germplasm population characterization, gene network identification, and genomic selection strategy innovation in plant breeding. The previous genome-wide association studies (GWAS) have been concentrated on finding a handful of major quantitative trait loci (QTL), but plant breeders are interested in revealing the whole-genome QTL-allele constitution in breeding materials/germplasm (in which tremendous historical allelic variation has been accumulated) for genome-wide improvement. To match this requirement, two innovations were suggested for GWAS: first grouping tightly linked sequential SNPs into linkage disequilibrium blocks (SNPLDBs) to form markers with multi-allelic haplotypes, and second utilizing two-stage association analysis for QTL identification, where the markers were preselected by single-locus model followed by multi-locus multi-allele model stepwise regression. Our proposed GWAS procedure is characterized as a novel restricted two-stage multi-locus multi-allele GWAS (RTM-GWAS, https://github.com/njau-sri/rtm-gwas ). The Chinese soybean germplasm population (CSGP) composed of 1024 accessions with 36,952 SNPLDBs (generated from 145,558 SNPs, with reduced linkage disequilibrium decay distance) was used to demonstrate the power and efficiency of RTM-GWAS. Using the CSGP marker information, simulation studies demonstrated that RTM-GWAS achieved the highest QTL detection power and efficiency compared with the previous procedures, especially under large sample size and high trait heritability conditions. A relatively thorough detection of QTL with their multiple alleles was achieved by RTM-GWAS compared with the linear mixed model method on 100-seed weight in CSGP. A QTL-allele matrix (402 alleles of 139 QTL×1024 accessions) was established as a compact form of the population genetic constitution. The 100-seed weight QTL-allele matrix was used for genetic characterization, candidate gene prediction, and genomic selection for optimal crosses in the germplasm population.

A Two‐Part Strategy for Using Genomic Selection to Develop Inbred Lines

We propose a strategy for implementing genomic selection in plant breeding programs for developing inbred lines that reorganizes traditional breeding programs into two distinct components. These components are: (i) a population improvement component to develop improved germplasm through rapid recurrent selection and (ii) a product development component to identify new inbred varieties or parents for hybrids using traditional breeding program designs. Stochastic simulations of entire breeding programs over 40 yr were used to evaluate the effectiveness of this strategy relative to a conventional program without genomic selection and programs using three standard strategies of implementing genomic selection. Cost effectiveness was measured by constraining all programs to approximately equal annual operating costs and directly comparing each program's overall performance. Programs using the two‐part strategy generated between 2.36 and 2.47 times more genetic gain than the conventional program and between 1.31 and 1.46 times more genetic gain than the best performing standard genomic selection strategy. These results indicate that the two‐part strategy is a cost‐effective strategy for implementing genomic selection in plant breeding programs.

Accuracy of genomic breeding values revisited: Assessment of two established approaches and a novel one to determine the accuracy in two-step genomic prediction.

Selection decisions in genomic selection schemes are made based on genomic breeding values (GBV) of candidates. Thus, the accuracy of GBV is a relevant parameter, as it reflects the stability of prediction and the possibility that the GBV might change when more information becomes available. Accuracy of genomic prediction defined as the correlation between GBV and true breeding values (TBV), however, is difficult to assess, considering TBV of the candidates are not available in reality. In previous studies, several methods were proposed to assess the accuracy of GBV including methods using population parameters or parameters inferred from mixed-model equations. In practice, most approaches tended to overestimate the accuracy of genomic prediction. We thus tested approaches used in previous studies in order to assess the magnitude of bias. Analyses were performed based on simulated data under a variety of scenarios mimicking different livestock breeding programmes. Furthermore, we proposed a novel method and tested it both with simulated data and in a real Holstein data set. The new method provided a better prediction for the accuracy of GBV in the simulated scenarios.

Detection and assessment of copy number variation using PacBio long-read and Illumina sequencing in New Zealand dairy cattle

Single nucleotide polymorphisms have been the DNA variant of choice for genomic prediction, largely because of the ease of single nucleotide polymorphism genotype collection. In contrast, structural variants (SV), which include copy number variants (CNV), translocations, insertions, and inversions, have eluded easy detection and characterization, particularly in nonhuman species. However, evidence increasingly shows that SV not only contribute a substantial proportion of genetic variation but also have significant influence on phenotypes. Here we present the discovery of CNV in a prominent New Zealand dairy bull using long-read PacBio (Pacific Biosciences, Menlo Park, CA) sequencing technology and the Sniffles SV discovery tool (version 0.0.1; https://github.com/fritzsedlazeck/Sniffles). The CNV identified from long reads were compared with CNV discovered in the same bull from Illumina sequencing using CNVnator (read depth-based tool; Illumina Inc., San Diego, CA) as a means of validation. Subsequently, further validation was undertaken using whole-genome Illumina sequencing of 556 cattle representing the wider New Zealand dairy cattle population. Very limited overlap was observed in CNV discovered from the 2 sequencing platforms, in part because of the differences in size of CNV detected. Only a few CNV were therefore able to be validated using this approach. However, the ability to use CNVnator to genotype the 557 cattle for copy number across all regions identified as putative CNV allowed a genome-wide assessment of transmission level of copy number based on pedigree. The more highly transmissible a putative CNV region was observed to be, the more likely the distribution of copy number was multimodal across the 557 sequenced animals. Furthermore, visual assessment of highly transmissible CNV regions provided evidence supporting the presence of CNV across the sequenced animals. This transmission-based approach was able to confirm a subset of CNV that segregates in the New Zealand dairy cattle population. Genome-wide identification and validation of CNV is an important step toward their inclusion in genomic selection strategies.

Mitigation of inbreeding while preserving genetic gain in genomic breeding programs for outbred plants.

Heuristic genomic inbreeding controls reduce inbreeding in genomic breeding schemes without reducing genetic gain. Genomic selection is increasingly being implemented in plant breeding programs to accelerate genetic gain of economically important traits. However, it may cause significant loss of genetic diversity when compared with traditional schemes using phenotypic selection. We propose heuristic strategies to control the rate of inbreeding in outbred plants, which can be categorised into three types: controls during mate allocation, during selection, and simultaneous selection and mate allocation. The proposed mate allocation measure GminF allocates two or more parents for mating in mating groups that minimise coancestry using a genomic relationship matrix. Two types of relationship-adjusted genomic breeding values for parent selection candidates ([Formula: see text]) and potential offspring ([Formula: see text]) are devised to control inbreeding during selection and even enabling simultaneous selection and mate allocation. These strategies were tested in a case study using a simulated perennial ryegrass breeding scheme. As compared to the genomic selection scheme without controls, all proposed strategies could significantly decrease inbreeding while achieving comparable genetic gain. In particular, the scenario using [Formula: see text] in simultaneous selection and mate allocation reduced inbreeding to one-third of the original genomic selection scheme. The proposed strategies are readily applicable in any outbred plant breeding program.

Economic evaluation of progeny-testing and genomic selection schemes for small-sized nucleus dairy cattle breeding programs in developing countries

In developing countries minimal and erratic performance and pedigree recording impede implementation of large-sized breeding programs. Small-sized nucleus programs offer an alternative but rely on their economic performance for their viability. We investigated the economic performance of 2 alternative small-sized dairy nucleus programs [i.e., progeny testing (PT) and genomic selection (GS)] over a 20-yr investment period. The nucleus was made up of 453 male and 360 female animals distributed in 8 non-overlapping age classes. Each year 10 active sires and 100 elite dams were selected. Populations of commercial recorded cows (CRC) of sizes 12,592 and 25,184 were used to produce test daughters in PT or to create a reference population in GS, respectively. Economic performance was defined as gross margins, calculated as discounted revenues minus discounted costs following a single generation of selection. Revenues were calculated as cumulative discounted expressions (CDE, kg) × 0.32 (€/kg of milk) × 100,000 (size commercial population). Genetic superiorities, deterministically simulated using pseudo-BLUP index and CDE, were determined using gene flow. Costs were for one generation of selection. Results show that GS schemes had higher cumulated genetic gain in the commercial cow population and higher gross margins compared with PT schemes. Gross margins were between 3.2- and 5.2-fold higher for GS, depending on size of the CRC population. The increase in gross margin was mostly due to a decreased generation interval and lower running costs in GS schemes. In PT schemes many bulls are culled before selection. We therefore also compared 2 schemes in which semen was stored instead of keeping live bulls. As expected, semen storage resulted in an increase in gross margins in PT schemes, but gross margins remained lower than those of GS schemes. We conclude that implementation of small-sized GS breeding schemes can be economically viable for developing countries.

Application of a Bayesian dominance model improves power in quantitative trait genome-wide association analysis

BackgroundMulti-marker methods, which fit all markers simultaneously, were originally tailored for genomic selection purposes, but have proven to be useful also in association analyses, especially the so-called BayesC Bayesian methods. In a recent study, BayesD extended BayesC towards accounting for dominance effects and improved prediction accuracy and persistence in genomic selection. The current study investigated the power and precision of BayesC and BayesD in genome-wide association studies by means of stochastic simulations and applied these methods to a dairy cattle dataset.MethodsThe simulation protocol was designed to mimic the genetic architecture of quantitative traits as realistically as possible. Special emphasis was put on the joint distribution of the additive and dominance effects of causative mutations. Additive marker effects were estimated by BayesC and additive and dominance effects by BayesD. The dependencies between additive and dominance effects were modelled in BayesD by choosing appropriate priors. A sliding-window approach was used. For each window, the R. Fernando window posterior probability of association was calculated and this was used for inference purpose. The power to map segregating causal effects and the mapping precision were assessed for various marker densities up to full sequence information and various window sizes.ResultsPower to map a QTL increased with higher marker densities and larger window sizes. This held true for both methods. Method BayesD had improved power compared to BayesC. The increase in power was between −2 and 8% for causative genes that explained more than 2.5% of the genetic variance. In addition, inspection of the estimates of genomic window dominance variance allowed for inference about the magnitude of dominance at significant associations, which remains hidden in BayesC analysis. Mapping precision was not substantially improved by BayesD.ConclusionsBayesD improved power, but precision only slightly. Application of BayesD needs large datasets with genotypes and own performance records as phenotypes. Given the current efforts to establish cow reference populations in dairy cattle genomic selection schemes, such datasets are expected to be soon available, which will enable the application of BayesD for association mapping and genomic prediction purposes.

Single Nucleotide Polymorphism Discovery in Bovine Pituitary Gland Using RNA-Seq Technology.

Examination of bovine pituitary gland transcriptome by strand-specific RNA-seq allows detection of putative single nucleotide polymorphisms (SNPs) within potential candidate genes (CGs) or QTLs regions as well as to understand the genomics variations that contribute to economic trait. Here we report a breed-specific model to successfully perform the detection of SNPs in the pituitary gland of young growing bulls representing Polish Holstein-Friesian (HF), Polish Red, and Hereford breeds at three developmental ages viz., six months, nine months, and twelve months. A total of 18 bovine pituitary gland polyA transcriptome libraries were prepared and sequenced using the Illumina NextSeq 500 platform. Sequenced FastQ databases of all 18 young bulls were submitted to NCBI-SRA database with NCBI-SRA accession numbers SRS1296732. For the investigated young bulls, a total of 113,882,3098 raw paired-end reads with a length of 156 bases were obtained, resulting in an approximately 63 million paired-end reads per library. Breed-wise, a total of 515.38, 215.39, and 408.04 million paired-end reads were obtained for Polish HF, Polish Red, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed 93.04%, 94.39%, and 83.46% of the mapped sequencing reads were properly paired to the Polish HF, Polish Red, and Hereford breeds, respectively. Constructed breed-specific SNP-db of three cattle breeds yielded at 13,775,885 SNPs. On an average 765,326 breed-specific SNPs per young bull were identified. Using two stringent filtering parameters, i.e., a minimum 10 SNP reads per base with an accuracy ≥ 90% and a minimum 10 SNP reads per base with an accuracy = 100%, SNP-db records were trimmed to construct a highly reliable SNP-db. This resulted in a reduction of 95,7% and 96,4% cut-off mark of constructed raw SNP-db. Finally, SNP discoveries using RNA-Seq data were validated by KASP™ SNP genotyping assay. The comprehensive QTLs/CGs analysis of 76 QTLs/CGs with RNA-seq data identified KCNIP4, CCSER1, DPP6, MAP3K5 and GHR CGs with highest SNPs hit loci in all three breeds and developmental ages. However, CAST CG with more than 100 SNPs hits were observed only in Polish HF and Hereford breeds.These findings are important for identification and construction of novel tissue specific SNP-db and breed specific SNP-db dataset by screening of putative SNPs according to QTL db and candidate genes for bovine growth and reproduction traits, one can develop genomic selection strategies for growth and reproductive traits.

Genetic effects of FASN, PPARGC1A, ABCG2 and IGF1 revealing the association with milk fatty acids in a Chinese Holstein cattle population based on a post genome-wide association study.

BackgroundA previous genome-wide association study deduced that one (ARS-BFGL-NGS-39328), two (Hapmap26001-BTC-038813 and Hapmap31284-BTC-039204), two (Hapmap26001-BTC-038813 and BTB-00246150), and one (Hapmap50366-BTA-46960) genome-wide significant single nucleotide polymorphisms (SNPs) associated with milk fatty acids were close to or within the fatty acid synthase (FASN), peroxisome proliferator-activated receptor gamma, coactivator 1 alpha (PPARGC1A), ATP-binding cassette, sub-family G, member 2 (ABCG2) and insulin-like growth factor 1 (IGF1) genes. To further confirm the linkage and reveal the genetic effects of these four candidate genes on milk fatty acid composition, genetic polymorphisms were identified and genotype-phenotype associations were performed in a Chinese Holstein cattle population.ResultsNine SNPs were identified in FASN, among which SNP rs41919985 was predicted to result in an amino acid substitution from threonine (ACC) to alanine (GCC), five SNPs (rs136947640, rs134340637, rs41919992, rs41919984 and rs41919986) were synonymous mutations, and the remaining three (rs41919999, rs132865003 and rs133498277) were found in FASN introns. Only one SNP each was identified for PPARGC1A, ABCG2 and IGF1.Association studies revealed that FASN, PPARGC1A, ABCG2 and IGF1 were mainly associated with medium-chain saturated fatty acids and long-chain unsaturated fatty acids, especially FASN for C10:0, C12:0 and C14:0. Strong linkage disequilibrium was observed among ARS-BFGL-NGS-39328 and rs132865003 and rs134340637 in FASN (D´ > 0.9), and among Hapmap26001-BTC-038813 and Hapmap31284-BTC-039204 and rs109579682 in PPARGC1A (D´ > 0.9). Subsequently, haplotype-based analysis revealed significant associations of the haplotypes encompassing eight FASN SNPs (rs41919999, rs132865003, rs134340637, rs41919992, rs133498277, rs41919984, rs41919985 and rs41919986) with C10:0, C12:0, C14:0, C18:1n9c, saturated fatty acids (SFA) and unsaturated fatty acids (UFA) (P = 0.0204 to P < 0.0001).ConclusionOur study confirmed the linkage between the significant SNPs in our previous genome-wide association study and variants in FASN and PPARGC1A. SNPs within FASN, PPARGC1A, ABCG2 and IGF1 showed significant genetic effects on milk fatty acid composition in dairy cattle, indicating their potential functions in milk fatty acids synthesis and metabolism. The findings presented here provide evidence for the selection of dairy cows with healthier milk fatty acid composition by marker-assisted breeding or genomic selection schemes, as well as furthering our understanding of technological processing aspects of cows’ milk.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-016-0418-x) contains supplementary material, which is available to authorized users.