Metodologias e estratégias de imputação de marcadores genéticos em bovinos da raça Cachim

High-density panels (HD) have strong level linkage disequilibrium among genetic markers (i.e. single nucleotide polymorphism - SNP), which allows better predictions of genomic breeding values. However, HD genotyping still expensive and became a limitation for the quantity of candidate animals used in genomic studies. As an alternative to decrease costs, imputation methods are powerful tools to infer missing marker genotypes from low-density (LD) panels to HD. Imputation uses information from a reference population of animals genotyped with a HD panel to impute variants that are not directly genotyped in LD panels. The objective of this study was to compare different scenarios and methodologies of imputation for the Canchim cattle. Data set was provided by Embrapa Pecuária Sudeste and comprised 285 Canchim animals, 114 MA genetic group animals, and 1 ancestor Charolais bull. Animals born between 1999 and 2005 were genotyped with the Illumina BovineHD panel (786,799 SNP). Data editing was performed in the R software and in C ++ language. Multiple scenarios combining different minor allele frequencies (MAF) thresholds for SNPs were tested: no MAF filter (QC1), and exclusion of SNPs with MAF lower than 0.0025 (QC2) and MAF lower than 0.10 (QC3). LD panels were created by masking SNPs originally present in the HD panel, and then assigning markers into the Illumina Bovine3K (3K), Illumina BovineLD (6K), Beef LD GeneSeek Genomic Profiler (9K), Indicus LD GeneSeek Genomic Profiler (20K), Illumina BovineSNP50 (50K), GeneSeek Genomic Profiler Beef HD (80K) and GeneSeek Genomic Profiler Indicus HD (90K) panels. Reference and target populations were defined as scenario 1 (C1), reference animals were born up until 2004 and target animals were born in 2005; scenario 2A (C2A), reference animals from Canchim breed and target animals from MA genetic group; scenario 2B (C2B), reference animals from MA genetic ...Painéis de marcadores genéticos de alta densidade (HD) possuem forte desequilíbrio de ligação, que permite melhores predições de valores genômicos. Entretanto, genotipar animais com estes painéis apresenta custo elevado, tornando-se uma limitação para a genotipagem de todos os candidatos à seleção. Uma alternativa para a redução desses custos é utilizar imputação de genótipos. A imputação é um método em que marcadores de uma população genotipada com painéis de baixa densidade (LD) são inferidos utilizando informações provenientes de uma população referência genotipada com painéis HD. O objetivo deste trabalho foi comparar em diferentes cenários metodologias de imputação de marcadores moleculares de polimorfismos de nucleotídeos únicos (SNP) em bovinos de corte da raça Canchim. Foram utilizadas informações de 285 animais da raça Canchim, 114 do grupo genético “MA” e 1 touro da raça Charolês genotipados com painel Illumina BovineHD BeadChip (786.799 SNP), nascidos entre 1999 e 2005 e provenientes da base de dados genômicos da Embrapa Pecuária Sudeste, São Carlos, SP. A edição dos dados foi realizada no software R e em linguagem C++. Para a frequência mínima de alelos (MAF) foram aplicados 3 diferentes critérios: sem remover MAF (QC1); SNP com MAF menor que 0,0025 (QC2) e menor que 0,10 (QC3) foram excluídos. O painel HD original foi reduzido para painéis de baixa densidade (LD) 3K, 6K, 9K, 50K, 20K, 80K e 90K, selecionando os marcadores em comum entre o painel HD original e os painéis comerciais Illumina Bovine3K (3K), BovineLD (6K), GeneSeek Genomic Profiler (GGP) Beef LD (9K), BovineSNP50 (50K), GGP Indicus LD (20K) ,GGP Beef HD (80K) e GGP Indicus HD (90K). Os animais foram divididos em diferentes cenários, denominados de população referência e imputação, sendo o cenário 1 (C1): População referência formada por animais nascidos ..

Genetic parameter estimates and principal component analysis on performance and carcass traits of a terminal pig sire line

This study aimed to estimate genetic parameters for five traits related to performance and carcass in a terminal pig sire line using multitrait (MUV) and reduced rank animal models, and to analyze the relationship between the breeding values (EBVs) of the traits using genetic principal component analysis (PCA) applied directly on the mixed model equations (one-step PCA) and on the EBVs obtained using the MUV (two-step PCA). The heritabilities were very similar using both procedures ranged from 0.032 +/- 0.014 for feed conversion to 0.102 +/- 0.029 for lean meat. Except for some differences in the heritability for feed conversion ratio, the heritability estimates were remarkably similar using One-step PCA or MUV. The genetic correlations obtained using One-Step PCA were more reliable due to their lower standard errors compared to those obtained using MUV. The one-step PCA had better fit than MUV and could reduce the time computing due to the reduction of the additive genetic matrix.Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq

Genetic analysis on accumulated productivity and calving intervals in Nelore cattle

The aim of this study was to estimate genetic parameters for accumulated productivity (ACP), first calving interval (CI1), second calving interval (CI2), and mean calving interval (MCI) in Nelore beef cattle. The ACP trait is a reproduction index and comprises the total number of calves born per dam, weight of weaned calves, and age of the dam at calving. Genetic parameters were estimated by the average information restricted maximum likelihood method in two-trait analyses. The average heritability estimate for ACP was 0.17 (0.03). For CI1, CI2, and MCI, the heritability estimates were 0.02, 0.02, and 0.06, respectively. Genetic correlations between ACP with CI1, CI2, and MCI were -0.16 ± 0.47, -0.29 ± 0.53, and -0.40 ± 0.27, respectively. Despite of the low heritability estimates obtained in our study, reproduction traits should be further studied and their inclusion in the selection criteria must be evaluated in order to improve the performance of females because these traits are of great economic importance in beef cattle. Accumulated productivity could contribute to decrease the mean calving interval in Nelore cattle.Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq

Genomic structure of a crossbred Landrace pig population.

Single nucleotide polymorphism (SNP) markers are used to study population structure and conservation genetics, which permits assessing similarities regarding the linkage disequilibrium and information about the relationship among individuals. To investigate the population genomic structure of 300 females and 25 males from a commercial maternal pig line we analyzed linkage disequilibrium extent, inbreeding coefficients using genomic and conventional pedigree data, and population stratification. The average linkage disequilibrium (r2) was 0.291 ± 0.312 for all adjacent SNPs, distancing less than 100 Kb (kilobase) between markers. The average inbreeding coefficients obtained from runs of homozygosity (ROH) and pedigree analyses were 0.119 and 0.0001, respectively. Low correlation was observed between the inbreeding coefficients possibly as a result of genetic recombination effect accounted for the ROH estimates or caused by pedigree identification errors. A large number of long ROHs might indicate recent inbreeding events in the studied population. A total of 36 homozygous segments were found in more than 30% of the population and these ROH harbor genes associated with reproductive traits. The population stratification analysis indicated that this population was possibly originated from two distinct populations, which is a result from crossings between the eastern and western breeds used in the formation of the line. Our findings provide support to understand the genetic structure of swine populations and may assist breeding companies to avoid a high level of inbreeding coefficients to maintain genetic diversity, showing the effectiveness of using genome-wide SNP information for quantifying inbreeding when the pedigree was incomplete or incorrect

Strategies for genotype imputation in composite beef cattle

Genotype imputation has been used to increase genomic information, allow more animals in genome-wide analyses, and reduce genotyping costs. In Brazilian beef cattle production, many animals are resulting from crossbreeding and such an event may alter linkage disequilibrium patterns. Thus, the challenge is to obtain accurately imputed genotypes in crossbred animals. The objective of this study was to evaluate the best fitting and most accurate imputation strategy on the MA genetic group (the progeny of a Charolais sire mated with crossbred Canchim X Zebu cows) and Canchim cattle. The data set contained 400 animals (born between 1999 and 2005) genotyped with the Illumina BovineHD panel. Imputation accuracy of genotypes from the Illumina-Bovine3K (3K), Illumina-BovineLD (6K), GeneSeek-Genomic-Profiler (GGP) BeefLD (GGP9K), GGP-IndicusLD (GGP20Ki), Illumina-BovineSNP50 (50K), GGP-IndicusHD (GGP75Ki), and GGP-BeefHD (GGP80K) to Illumina-BovineHD (HD) SNP panels were investigated. Seven scenarios for reference and target populations were tested; the animals were grouped according with birth year (S1), genetic groups (S2 and S3), genetic groups and birth year (S4 and S5), gender (S6), and gender and birth year (S7). Analyses were performed using FImpute and BEAGLE software and computation run-time was recorded. Genotype imputation accuracy was measured by concordance rate (CR) and allelic R square (R(2)). The highest imputation accuracy scenario consisted of a reference population with males and females and a target population with young females. Among the SNP panels in the tested scenarios, from the 50K, GGP75Ki and GGP80K were the most adequate to impute to HD in Canchim cattle. FImpute reduced computation run-time to impute genotypes from 20 to 100 times when compared to BEAGLE. The genotyping panels possessing at least 50 thousands markers are suitable for genotype imputation to HD with acceptable accuracy. The FImpute algorithm demonstrated a higher efficiency of imputed markers, especially in lower density panels. These considerations may assist to increase genotypic information, reduce genotyping costs, and aid in genomic selection evaluations in crossbred animals

Comparison of BovineHD BeadChip array genotypes to sequencing SNVs.

<p>Comparison of BovineHD BeadChip array genotypes to sequencing SNVs.</p

Single nucleotide variants and InDels identified from whole-genome re-sequencing of Guzerat, Gyr, Girolando and Holstein cattle breeds

<div><p>Whole-genome re-sequencing, alignment and annotation analyses were undertaken for 12 sires representing four important cattle breeds in Brazil: Guzerat (multi-purpose), Gyr, Girolando and Holstein (dairy production). A total of approximately 4.3 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.7 to 16.4-fold genome coverage. A total of 27,441,279 single nucleotide variations (SNVs) and 3,828,041 insertions/deletions (InDels) were detected in the samples, of which 2,557,670 SNVs and 883,219 InDels were novel. The submission of these genetic variants to the dbSNP database significantly increased the number of known variants, particularly for the indicine genome. The concordance rate between genotypes obtained using the Bovine HD BeadChip array and the same variants identified by sequencing was about 99.05%. The annotation of variants identified numerous non-synonymous SNVs and frameshift InDels which could affect phenotypic variation. Functional enrichment analysis was performed and revealed that variants in the olfactory transduction pathway was over represented in all four cattle breeds, while the ECM-receptor interaction pathway was over represented in Girolando and Guzerat breeds, the ABC transporters pathway was over represented only in Holstein breed, and the metabolic pathways was over represented only in Gyr breed. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Gyr, Girolando, Guzerat and Holstein breeding programs.</p></div

Description of insertion-to-deletion, heterozygous-to-homozygous (Het/Hom) and transition-to-transversion (Ti/Tv) ratios for InDels and SNVs.

<p>Description of insertion-to-deletion, heterozygous-to-homozygous (Het/Hom) and transition-to-transversion (Ti/Tv) ratios for InDels and SNVs.</p

Gene Ontology (GO) terms and KEGG pathways enriched (FDR<0.10) in Gyr, Girolando, Guzerat and Holstein breeds.

<p>Gene Ontology (GO) terms and KEGG pathways enriched (FDR<0.10) in Gyr, Girolando, Guzerat and Holstein breeds.</p