Dataset Information

Dynamics of genomic architecture during composite breed development in cattle.

ABSTRACT: Some livestock breeds face the challenge of reduced genetic variation, increased inbreeding depression owing to genetic drift and selection. Hybridization can reverse these processes and increase levels of productivity and adaptation to various environmental stressors. Samples from American Brangus were used to evaluate the indicine/taurine composition through nine generations (~45 years) after the hybridization process was completed. The purpose was to determine how hybridization alters allelic combinations of a breed over time when genetic factors such as selection and drift are operating. Furthermore, we explored genomic regions with deviations from the expected composition from the progenitor breeds and related these regions to traits under selection. The Brangus composition deviated from the theoretical expectation, defined by the breed association, of 62.5% taurine, showing taurine composition to be 70.4 ± 0.6%. Taurine and indicine proportion were not consistent across chromosomes. Furthermore, these non-uniform areas were found to be associated with traits that were probably under selection such as intermuscular fat and average daily gain. Interestingly, the sex chromosomes were predominantly taurine, which could be due to the composite being formed particularly in the final cross that resulted in progeny designated as purebred Brangus. This work demonstrated the process of new breed formation on a genomic level. It suggests that factors like genetic drift, selection and complementarity shift the genetic architecture into a uniquely different population. These findings are important to better understand how hybridization and crossbreeding systems shape the genetic architecture of composite populations.

SUBMITTER: Paim TDP

PROVIDER: S-EPMC7065137 | biostudies-literature | 2020 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Dynamics of genomic architecture during composite breed development in cattle.

Paim T do P TDP Hay E H A EHA Wilson C C Thomas M G MG Kuehn L A LA Paiva S R SR McManus C C Blackburn H D HD

Animal genetics 20200121 2

Some livestock breeds face the challenge of reduced genetic variation, increased inbreeding depression owing to genetic drift and selection. Hybridization can reverse these processes and increase levels of productivity and adaptation to various environmental stressors. Samples from American Brangus were used to evaluate the indicine/taurine composition through nine generations (~45 years) after the hybridization process was completed. The purpose was to determine how hybridization alters allelic ...[more]

PMID: 31961956

Similar Datasets

Project description:BackgroundDiscerning the traits evolving under neutral conditions from those traits evolving rapidly because of various selection pressures is a great challenge. We propose a new method, composite selection signals (CSS), which unifies the multiple pieces of selection evidence from the rank distribution of its diverse constituent tests. The extreme CSS scores capture highly differentiated loci and underlying common variants hauling excess haplotype homozygosity in the samples of a target population.ResultsThe data on high-density genotypes were analyzed for evidence of an association with either polledness or double muscling in various cohorts of cattle and sheep. In cattle, extreme CSS scores were found in the candidate regions on autosome BTA-1 and BTA-2, flanking the POLL locus and MSTN gene, for polledness and double muscling, respectively. In sheep, the regions with extreme scores were localized on autosome OAR-2 harbouring the MSTN gene for double muscling and on OAR-10 harbouring the RXFP2 gene for polledness. In comparison to the constituent tests, there was a partial agreement between the signals at the four candidate loci; however, they consistently identified additional genomic regions harbouring no known genes. Persuasively, our list of all the additional significant CSS regions contains genes that have been successfully implicated to secondary phenotypic diversity among several subpopulations in our data. For example, the method identified a strong selection signature for stature in cattle capturing selective sweeps harbouring UQCC-GDF5 and PLAG1-CHCHD7 gene regions on BTA-13 and BTA-14, respectively. Both gene pairs have been previously associated with height in humans, while PLAG1-CHCHD7 has also been reported for stature in cattle. In the additional analysis, CSS identified significant regions harbouring multiple genes for various traits under selection in European cattle including polledness, adaptation, metabolism, growth rate, stature, immunity, reproduction traits and some other candidate genes for dairy and beef production.ConclusionsCSS successfully localized the candidate regions in validation datasets as well as identified previously known and novel regions for various traits experiencing selection pressure. Together, the results demonstrate the utility of CSS by its improved power, reduced false positives and high-resolution of selection signals as compared to individual constituent tests.

Project description:Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from -0.009 to 0.029. A low (-0.12-0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (<4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5-1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.

Project description:Genomic breed composition (GBC) of an individual animal refers to the partition of its genome according to the inheritance from its ancestors or ancestral breeds. For crossbred or composite animals, knowing their GBC is useful to estimate heterosis, to characterize their actual inheritance from foundation breeds, and to make management decisions for crossbreeding programs. Various statistical approaches have been proposed to estimate GBC in animals, but the interpretations of estimates have varied with these methods. In the present study, we proposed a causality interpretation of GBC based on path analysis. We applied this method to estimating GBC in two composite breeds of beef cattle, namely Brangus and Beefmaster. Three SNP panels were used to estimate GBC: a 10K SNP panel consisting of 10,226 common SNPs across three GeneSeek Genomic Profiler (GGP) bovine SNP arrays (GGP 30K, GGP 40K, and GGP 50K), and two subsets (1K and 5K) of uniformly distributed SNPs. The path analysis decomposed the relationships between the ancestors and the composite animals into direct and indirect path effects, and GBC was measured by the relative ratio of the coefficients of direct (D-GBC) and combined (C-GBC) effects from each ancestral breed to the progeny, respectively. Estimated GBC varied only slightly between different genotyping platforms and between the three SNP panels. In the Brangus cattle, because the two ancestral breeds had a very distant relationship, the estimated D-GBC and C-GBC were comparable to each other in the path analysis, and they corresponded roughly to the estimated GBC from the linear regression and the admixture model. In the Beefmaster, however, the strong relationship in allelic frequencies between Hereford and Shorthorn imposed a challenge for the linear regression and the admixture model to estimated GBC reliably. Instead, D-GBC by the path analysis included only direct ancestral effects, and it was robust to bias due to high genomic correlations between reference (ancestral) breeds.

Project description:Reducing the incidence of both the degree of assistance required at calving, as well as the extent of perinatal mortality (PM) has both economic and societal benefits. The existence of heritable genetic variability in both traits signifies the presence of underlying genomic variability. The objective of the present study was to locate regions of the genome, and by extension putative genes and mutations, that are likely to be underpinning the genetic variability in direct calving difficulty (DCD), maternal calving difficulty (MCD), and PM. Imputed whole-genome single-nucleotide polymorphism (SNP) data on up to 8,304 Angus (AA), 17,175 Charolais (CH), 16,794 Limousin (LM), and 18,474 Holstein-Friesian (HF) sires representing 5,866,712 calving events from descendants were used. Several putative quantitative trait loci (QTL) regions associated with calving performance both within and across dairy and beef breeds were identified, although the majority were both breed- and trait-specific. QTL surrounding and encompassing the myostatin (MSTN) gene were associated (P < 5 × 10-8) with DCD and PM in both the CH and LM populations. The well-known Q204X mutation was the fifth strongest association with DCD in the CH population and accounted for 5.09% of the genetic variance in DCD. In contrast, none of the 259 segregating variants in MSTN were associated (P > × 10-6) with DCD in the LM population but a genomic region 617 kb downstream of MSTN was associated (P < 5 × 10-8). The genetic architecture for DCD differed in the HF population relative to the CH and LM, where two QTL encompassing ZNF613 on Bos taurus autosome (BTA)18 and PLAG1 on BTA14 were identified in the former. Pleiotropic SNP associated with all three calving performance traits were also identified in the three beef breeds; 5 SNP were pleiotropic in AA, 116 in LM, and 882 in CH but no SNP was associated with more than one trait within the HF population. The majority of these pleiotropic SNP were on BTA2 surrounding MSTN and were associated with both DCD and PM. Multiple previously reported, but also novel QTL, associated with calving performance were detected in this large study. These also included QTL regions harboring SNP with the same direction of allele substitution effect for both DCD and MCD thus contributing to a more effective simultaneous selection for both traits.

Project description:Crossbreeding is a common strategy used in tropical and subtropical regions to enhance beef production, and having accurate knowledge of breed composition is essential for the success of a crossbreeding program. Although pedigree records have been traditionally used to obtain the breed composition of crossbred cattle, the accuracy of pedigree-based breed composition can be reduced by inaccurate and/or incomplete records and Mendelian sampling. Breed composition estimation from genomic data has multiple advantages including higher accuracy without being affected by missing, incomplete, or inaccurate records and the ability to be used as independent authentication of breed in breed-labeled beef products. The present study was conducted with 676 Angus-Brahman crossbred cattle with genotype and pedigree information to evaluate the feasibility and accuracy of using genomic data to determine breed composition. We used genomic data in parametric and non-parametric methods to detect population structure due to differences in breed composition while accounting for the confounding effect of close familial relationships. By applying principal component analysis (PCA) and the maximum likelihood method of ADMIXTURE to genomic data, it was possible to successfully characterize population structure resulting from heterogeneous breed ancestry, while accounting for close familial relationships. PCA results offered additional insight into the different hierarchies of genetic variation structuring. The first principal component was strongly correlated with Angus-Brahman proportions, and the second represented variation within animals that have a relatively more extended Brangus lineage-indicating the presence of a distinct pattern of genetic variation in these cattle. Although there was strong agreement between breed proportions estimated from pedigree and genetic information, there were significant discrepancies between these two methods for certain animals. This was most likely due to inaccuracies in the pedigree-based estimation of breed composition, which supported the case for using genomic information to complement and/or replace pedigree information when estimating breed composition. Comparison with a supervised analysis where purebreds are used as the training set suggest that accurate predictions can be achieved even in the absence of purebred population information.

Project description:BackgroundReference-guided read alignment and variant genotyping are prone to reference allele bias, particularly for samples that are greatly divergent from the reference genome. A Hereford-based assembly is the widely accepted bovine reference genome. Haplotype-resolved genomes that exceed the current bovine reference genome in quality and continuity have been assembled for different breeds of cattle. Using whole genome sequencing data of 161 Brown Swiss cattle, we compared the accuracy of read mapping and sequence variant genotyping as well as downstream genomic analyses between the bovine reference genome (ARS-UCD1.2) and a highly continuous Angus-based assembly (UOA_Angus_1).ResultsRead mapping accuracy did not differ notably between the ARS-UCD1.2 and UOA_Angus_1 assemblies. We discovered 22,744,517 and 22,559,675 high-quality variants from ARS-UCD1.2 and UOA_Angus_1, respectively. The concordance between sequence- and array-called genotypes was high and the number of variants deviating from Hardy-Weinberg proportions was low at segregating sites for both assemblies. More artefactual INDELs were genotyped from UOA_Angus_1 than ARS-UCD1.2 alignments. Using the composite likelihood ratio test, we detected 40 and 33 signatures of selection from ARS-UCD1.2 and UOA_Angus_1, respectively, but the overlap between both assemblies was low. Using the 161 sequenced Brown Swiss cattle as a reference panel, we imputed sequence variant genotypes into a mapping cohort of 30,499 cattle that had microarray-derived genotypes using a two-step imputation approach. The accuracy of imputation (Beagle R2) was very high (0.87) for both assemblies. Genome-wide association studies between imputed sequence variant genotypes and six dairy traits as well as stature produced almost identical results from both assemblies.ConclusionsThe ARS-UCD1.2 and UOA_Angus_1 assemblies are suitable for reference-guided genome analyses in Brown Swiss cattle. Although differences in read mapping and genotyping accuracy between both assemblies are negligible, the choice of the reference genome has a large impact on detecting signatures of selection that already reached fixation using the composite likelihood ratio test. We developed a workflow that can be adapted and reused to compare the impact of reference genomes on genome analyses in various breeds, populations and species.

Dataset Information

Dynamics of genomic architecture during composite breed development in cattle.

Publications

Dynamics of genomic architecture during composite breed development in cattle.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets