Building in vitro tools for livestock genomics: chromosomal variation within the PK15 cell line.

BACKGROUND: Cultured porcine cell lines are powerful tools for functional genomics and in vitro phenotypic testing of candidate causal variants. However, to be utilised for genomic or variant interrogation assays, the genome sequence and structure of cultured cell lines must be realised. In this work, we called variants and used read coverage in combination with within-sample allele frequency to detect potential aneuploidy in two immortalised porcine kidney epithelial (PK15) cell lines and in a pig embryonic fibroblast line. RESULTS: We compared two PK15 cultured cells samples: a new American Type Culture Collection (ATCC) sample and one that has been utilised and passaged within the laboratory for an extended period (> 10 years). Read coverage and within-sample allele frequencies showed that several chromosomes are fully or partially aneuploid in both PK15 lines, including potential trisomy of chromosome 4 and tetrasomy of chromosome 17. The older PK15 line showed evidence of additional structural variation and potentially clonal variation. By comparison, the pig embryonic fibroblast line was free from the gross aneuploidies seen in the PK15s. CONCLUSIONS: Our results show that the PK15 cell lines examined have aneuploidies and complex structural variants in their genomes. We propose that screening for aneuploidy should be considered for cell lines, and discuss implications for livestock genomics.

Efficient use of genomic information for sustainable genetic improvement in small cattle populations.

In this study, we compared genetic gain, genetic variation, and the efficiency of converting variation into gain under different genomic selection scenarios with truncation or optimum contribution selection in a small dairy population by simulation. Breeding programs have to maximize genetic gain but also ensure sustainability by maintaining genetic variation. Numerous studies have shown that genomic selection increases genetic gain. Although genomic selection is a well-established method, small populations still struggle with choosing the most sustainable strategy to adopt this type of selection. We developed a simulator of a dairy population and simulated a model after the Slovenian Brown Swiss population with ∼10,500 cows. We compared different truncation selection scenarios by varying (1) the method of sire selection and their use on cows or bull-dams, and (2) selection intensity and the number of years a sire is in use. Furthermore, we compared different optimum contribution selection scenarios with optimization of sire selection and their usage. We compared scenarios in terms of genetic gain, selection accuracy, generation interval, genetic and genic variance, rate of coancestry, effective population size, and conversion efficiency. The results showed that early use of genomically tested sires increased genetic gain compared with progeny testing, as expected from changes in selection accuracy and generation interval. A faster turnover of sires from year to year and higher intensity increased the genetic gain even further but increased the loss of genetic variation per year. Although maximizing intensity gave the lowest conversion efficiency, faster turnover of sires gave an intermediate conversion efficiency. The largest conversion efficiency was achieved with the simultaneous use of genomically and progeny-tested sires that were used over several years. Compared with truncation selection, optimizing sire selection and their usage increased the conversion efficiency by achieving either comparable genetic gain for a smaller loss of genetic variation or higher genetic gain for a comparable loss of genetic variation. Our results will help breeding organizations implement sustainable genomic selection.

Accuracy of genotype imputation in Labrador Retrievers.

The dog is a valuable model species for the genetic analysis of complex traits, and the use of genotype imputation in dogs will be an important tool for future studies. It is of particular interest to analyse the effect of factors like single nucleotide polymorphism (SNP) density of genotyping arrays and relatedness between dogs on imputation accuracy due to the acknowledged genetic and pedigree structure of dog breeds. In this study, we simulated different genotyping strategies based on data from 1179 Labrador Retriever dogs. The study involved 5826 SNPs on chromosome 1 representing the high density (HighD) array; the low-density (LowD) array was simulated by masking different proportions of SNPs on the HighD array. The correlations between true and imputed genotypes for a realistic masking level of 87.5% ranged from 0.92 to 0.97, depending on the scenario used. A correlation of 0.92 was found for a likely scenario (10% of dogs genotyped using HighD, 87.5% of HighD SNPs masked in the LowD array), which indicates that genotype imputation in Labrador Retrievers can be a valuable tool to reduce experimental costs while increasing sample size. Furthermore, we show that genotype imputation can be performed successfully even without pedigree information and with low relatedness between dogs in the reference and validation sets. Based on these results, the impact of genotype imputation was evaluated in a genome-wide association analysis and genomic prediction in Labrador Retrievers.

Genomic evaluations using similarity between haplotypes.

Long-range phasing and haplotype library imputation methodologies are accurate and efficient methods to provide haplotype information that could be used in prediction of breeding value or phenotype. Modelling long haplotypes as independent effects in genomic prediction would be inefficient due to the many effects that need to be estimated and phasing errors, even if relatively low in frequency, exacerbate this problem. One approach to overcome this is to use similarity between haplotypes to model covariance of genomic effects by region or of animal breeding values. We developed a simple method to do this and tested impact on genomic prediction by simulation. Results show that the diagonal and off-diagonal elements of a genomic relationship matrix constructed using the haplotype similarity method had higher correlations with the true relationship between pairs of individuals than genomic relationship matrices built using unphased genotypes or assumed unrelated haplotypes. However, the prediction accuracy of such haplotype-based prediction methods was not higher than those based on unphased genotype information.

Heterogeneity of genetic parameters for calving difficulty in Holstein heifers in Ireland.

Calving difficulty is a trait that greatly affects animal welfare, herd profitability, and the amount of labor required by cattle farmers. It is influenced by direct and maternal genetic components. Selection and breeding strategies can optimize the accuracy of genetic evaluations and correctly emphasize calving difficulty in multiple-trait indices provided there are accurate estimates of genetic parameters. In Ireland, large differences exist in the age at which heifers first give birth to calves. The objective of this study was to estimate genetic parameters for calving difficulty in first-parity Holsteins and to determine whether these differed with age of the heifer at calving. Transformed calving difficulty records for 18,798 Holstein heifers, which calved between January 2002 and May 2006, were analyzed using univariate, multitrait, and random regression linear sire-maternal grandsire models. The model that 1) fitted a second-order random regression of dam age at first parity for the direct component, 2) treated the maternal component as a single trait regardless of dam age, and 3) fitted a single residual variance component was optimal. Heritabilities for direct (0.13) and maternal (0.04) calving difficulty were significantly different from zero. These 2 components were moderately negatively correlated (-0.47). Estimates of direct genetic variance and heritability were heterogeneous along the dam age trajectory, decreasing initially with dam age before subsequently increasing. Heritability estimates ranged between 0.11 and 0.37 and were higher for records with younger and older dams at parturition. Genetic correlations between the direct components of calving difficulty decreased from unity to 0.5 with increasing distance between dam ages at parturition. The results of this study indicated that heterogeneity of direct genetic variance existed for calving difficulty, depending on dam age at first parturition.

The role of correctional officers in multidisciplinary mental health care in prisons.

Prisons have become the homes of thousands of inmates who have mental disorders. The stress of incarceration can cause morbidity among these individuals, resulting in more severe symptoms and more disruptive behavior. Effective treatment for such inmates often involves services provided by a multidisciplinary treatment team that includes correctional officers. Correctional officers can assist in observations and interventions, and they play a unique role on specialized housing units. Successful collaboration between correctional officers and treatment teams requires a foundation of mutual respect, shared training, and ongoing communication and cooperation. With these elements in place, correctional officers can assist the treatment team and make important and constructive contributions to the assessment and management of offenders who have mental disorders.

Evaluation of measures of correctness of genotype imputation in the context of genomic prediction: a review of livestock applications.

In livestock, many studies have reported the results of imputation to 50k single nucleotide polymorphism (SNP) genotypes for animals that are genotyped with low-density SNP panels. The objective of this paper is to review different measures of correctness of imputation, and to evaluate their utility depending on the purpose of the imputed genotypes. Across studies, imputation accuracy, computed as the correlation between true and imputed genotypes, and imputation error rates, that counts the number of incorrectly imputed alleles, are commonly used measures of imputation correctness. Based on the nature of both measures and results reported in the literature, imputation accuracy appears to be a more useful measure of the correctness of imputation than imputation error rates, because imputation accuracy does not depend on minor allele frequency (MAF), whereas imputation error rate depends on MAF. Therefore imputation accuracy can be better compared across loci with different MAF. Imputation accuracy depends on the ability of identifying the correct haplotype of a SNP, but many other factors have been identified as well, including the number of genotyped immediate ancestors, the number of animals with genotypes at the high-density panel, the SNP density on the low- and high-density panel, the MAF of the imputed SNP and whether imputed SNP are located at the end of a chromosome or not. Some of these factors directly contribute to the linkage disequilibrium between imputed SNP and SNP on the low-density panel. When imputation accuracy is assessed as a predictor for the accuracy of subsequent genomic prediction, we recommend that: (1) individual-specific imputation accuracies should be used that are computed after centring and scaling both true and imputed genotypes; and (2) imputation of gene dosage is preferred over imputation of the most likely genotype, as this increases accuracy and reduces bias of the imputed genotypes and the subsequent genomic predictions.

Practical implementation of cost-effective genomic selection in commercial pig breeding using imputation.

Genomic selection can be implemented in pig breeding at a reduced cost using genotype imputation. Accuracy of imputation and the impact on resulting genomic breeding values (gEBV) was investigated. High-density genotype data was available for 4,763 animals from a single pig line. Three low-density genotype panels were constructed with SNP densities of 450 (L450), 3,071 (L3k) and 5,963 (L6k). Accuracy of imputation was determined using 184 test individuals with no genotyped descendants in the data but with parents and grandparents genotyped using the Illumina PorcineSNP60 Beadchip. Alternative genotyping scenarios were created in which parents, grandparents, and individuals that were not direct ancestors of test animals (Other) were genotyped at high density (S1), grandparents were not genotyped (S2), dams and granddams were not genotyped (S3), and dams and granddams were genotyped at low density (S4). Four additional scenarios were created by excluding Other animal genotypes. Test individuals were always genotyped at low density. Imputation was performed with AlphaImpute. Genomic breeding values were calculated using the single-step genomic evaluation. Test animals were evaluated for the information retained in the gEBV, calculated as the correlation between gEBV using imputed genotypes and gEBV using true genotypes. Accuracy of imputation was high for all scenarios but decreased with fewer SNP on the low-density panel (0.995 to 0.965 for S1) and with reduced genotyping of ancestors, where the largest changes were for L450 (0.965 in S1 to 0.914 in S3). Exclusion of genotypes for Other animals resulted in only small accuracy decreases. Imputation accuracy was not consistent across the genome. Information retained in the gEBV was related to genotyping scenario and thus to imputation accuracy. Reducing the number of SNP on the low-density panel reduced the information retained in the gEBV, with the largest decrease observed from L3k to L450. Excluding Other animal genotypes had little impact on imputation accuracy but caused large decreases in the information retained in the gEBV. These results indicate that accuracy of gEBV from imputed genotypes depends on the level of genotyping in close relatives and the size of the genotyped dataset. Fewer high-density genotyped individuals are needed to obtain accurate imputation than are needed to obtain accurate gEBV. Strategies to optimize development of low-density panels can improve both imputation and gEBV accuracy.

Estimation of accuracy and bias in genetic evaluations with genetic groups using sampling.

Accuracy and bias of EBV are important measures of the quality of genetic evaluations. A sampling method that accounts for the uncertainty in the estimation of genetic group effects was used to calculate accuracy and bias of estimated effects. The method works by repeatedly simulating phenotypes for multiple traits for a defined data and pedigree structure. These simulated values are analyzed using BLUP with genetic groups in the relationship matrix. Accuracies and biases are then calculated as correlations among and differences between true and estimated values across all replicates, respectively. The method was applied to the Irish beef production data set for 15 traits and with 15 genetic groups to account for differences in breed means. Accuracy and bias of estimated genetic group effects, estimated comparisons between genetic group effects, EBV within genetic group, and EBV across genetic group were calculated. Small biases were detected for most estimated genetic group effects and most estimated comparisons between genetic group effects. Most of these were not important relative to the phenotypic SD of the traits involved. For example, a bias of 0.78% of the phenotypic SD was detected for carcass conformation in Aberdeen Angus. However, one trait, calf quality, which had few performance records in the data set, displayed larger bias, ranging from -10.31 to 5.85% of the phenotypic SD across the different estimated genetic group effects. Large differences were observed in the accuracies of genetic group effects, ranging from 0.02 for feed intake in Holstein, which had no data recorded, to >0.97 for carcass conformation, a trait with large amounts of data recorded in the different genetic groups. Large differences were also observed in the accuracies of the comparisons among genetic group effects. The accuracies of the EBV within genetic group and EBV across genetic group were sometimes different; for example, carcass conformation in Belgian Blue had an average accuracy within genetic group of 0.69 compared with an average accuracy across genetic group of 0.89. This suggests that the accuracy of genetic groups should be taken into account when publishing EBV across genetic groups.

Genetic parameters for EUROP carcass traits within different groups of cattle in Ireland.

The first objective of this study was to test the ability of systems of weighing and classifying bovine carcasses used in commercial abattoirs in Ireland to provide information that can be used for the purposes of genetic evaluation of carcass weight, carcass fatness class, and carcass conformation class. Secondly, the study aimed to test whether genetic and phenotypic variances differed by breed of sire. Variance components for carcass traits were estimated for crosses between dairy cows and 8 breeds of sire commonly found in the Irish cattle population. These 8 breeds were Aberdeen Angus, Belgian Blue, Charolais, Friesian, Hereford, Holstein, Limousin, and Simmental. A multivariate animal model was used to estimate genetic parameters within the Holstein sire breed group. Univariate analyses were used to estimate variance components for the remaining 7 sire breed groups. Multivariate sire models were used to formally test differences in genetic variances in sire breed groups. Field data on 64,443 animals, which were slaughtered in commercial abattoirs between the ages of 300 and 875 d, were analyzed in 8 analyses. Carcass fat class and carcass conformation class were measured using the European Union beef carcass classification system (EUROP) scale. For all 3 traits, the sire breed group with the greatest genetic variance had a value of more than 8 times the sire breed group with least genetic variance. Heritabilities ranged from zero to moderate for carcass fatness class (0.00 to 0.40), from low to moderate for carcass conformation class (0.04 to 0.36), and from low to high for carcass weight (0.06 to 0.65). Carcass weight was the most heritable (0.26) of the 3 traits. Carcass conformation class and carcass fatness class were equally heritable (0.17). Genetic and phenotypic correlations were all positive in the Holstein sire breed group. The genetic correlations varied from 0.11 for the relationship between carcass weight and carcass fatness class to 0.44 for the relationship between carcass conformation class and carcass fatness class. Carcass weight and classification data collected in Irish abattoirs are useful for the purposes of genetic evaluation for beef traits of Irish cattle. There were significantly different variance components across the sire breed groups.