Analyses of Signatures of Selection in the Bovine Genome

This thesis presents a novel method called meta-selection-scores (MSS) to construct a meta-assembly of genome-wide signatures of selection within breeds (n = 24) and across geographical archetypes (European, African, Zebu and composite). The meta-assembly and review of published studies highlight the historical selection events, the role of underlying genes for various traits, limitations of the available bovine genomic resources and implications of using different methodologies. This study also developed another new method called composite selection signals (CSS) – to improve the power of selection signature scans – in which multiple pieces of evidence for selection derived from the rank distribution of individual tests are combined in a single score. CSS has shown improved power by localizing the genomic regions under selection for major traits in multi-breed panels of cattle and sheep. CSS was also used to detect regions associated with the complex trait of bovine stature and implicated 12 (nine were novel) regions harbouring multiple candidate genes in contrasting cohorts within the European and African cattle. Without phenotypic records on individual animals, the CSS method provides a firsthand scan of the genome to detect putative regions associated with complex traits. The genome-wide CSS scans for individual breeds (n = 60) of cattle using previously available and new genotypes from 50K SNP chip assay identified 177 genomic regions under selection. Finally, an ultra-high density dataset (1.6 million SNPs) was investigated with CSS for Angus and Holstein. A high (8 of 9) reproducibility of the 50K SNPs based CSS regions and up to 14 fold additional genomic regions under selection were detected using 1.6 million SNPs. Overall, this study provides a detailed investigation about the core traits influenced by the historical selection events in worldwide cattle breeds and presents novel insights about the hotspots of positive selection in the bovine genome

Comparison of Gene Editing Versus Conventional Breeding to Introgress the POLLED Allele Into the Tropically Adapted Australian Beef Cattle Population

Dehorning is the process of physically removing horns to protect animals and humans from injury, but the process is costly, unpleasant, and faces increasing public scrutiny. Genetic selection for polled (hornless), which is genetically dominant to horned, is a long-term solution to eliminate the need for dehorning. However, due to the limited number of polled Australian Brahman bulls, the northern Australian beef cattle population remains predominantly horned. The potential to use gene editing to produce high-genetic-merit polled cattle was recently demonstrated. To further explore the concept, this study simulated introgression of the POLLED allele into a tropically adapted Australian beef cattle population via conventional breeding or gene editing (top 1% or 10% of seedstock bulls/year) for 3 polled mating schemes and compared results to baseline selection on genetic merit (Japan Ox selection index, $JapOx) alone, over the course of 20 years. The baseline scenario did not significantly decrease the 20-year HORNED allele frequency (80$8.00/year). Compared to the baseline, the conventional breeding scenarios where polled bulls were preferentially used for breeding, regardless of their genetic merit, significantly decreased the 20-year HORNED allele frequency (30%), but resulted in a significantly slower rate of genetic gain ($6.70/year, P ≤ 0.05). The mating scheme that required the exclusive use of homozygous polled bulls, resulted in the lowest 20-year HORNED allele frequency (8$5.50/year). The addition of gene editing the top 1% or 10% of seedstock bull calves/year to each conventional breeding scenario resulted in significantly faster rates of genetic gain (up to $8.10/year, P ≤ 0.05). Overall, our study demonstrates that, due to the limited number of polled Australian Brahman bulls, strong selection pressure on polled will be necessary to meaningfully increase the number of polled animals in this population. Moreover, these scenarios illustrate how gene editing could be a tool for accelerating the development of high-genetic-merit homozygous polled sires to mitigate the current trade-off of slower genetic gain associated with decreasing HORNED allele frequency in the Australian Brahman population

Composite selection signals can localize the trait specific genomic regions in multi-breed populations of cattle and sheep

Background: Discerning 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.Results: The 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.Conclusions: CSS 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

A meta-assembly of selection signatures in cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed "selection signatures". This article reviews 64 recent (2009-2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species

Composite selection signals for complex traits exemplified through bovine stature using multibreed cohorts of European and African Bos taurus

Understanding the evolution and molecular architecture of complex traits is important in domestic animals. Due to phenotypic selection, genomic regions develop unique patterns of genetic diversity called signatures of selection, which are challenging to detect, especially for complex polygenic traits. In this study, we applied the composite selection signals (CSS) method to investigate evidence of positive selection in a complex polygenic trait by examining stature in phenotypically diverse cattle comprising 47 European and 8 African Bos taurus breeds, utilizing a panel of 38,033 SNPs genotyped on 1106 animals. CSS were computed for phenotypic contrasts between multibreed cohorts of cattle by classifying the breeds according to their documented wither height to detect the candidate regions under selection. Using the CSS method, clusters of signatures of selection were detected at 26 regions (9 in European and 17 in African cohorts) on 13 bovine autosomes. Using comparative mapping information on human height, 30 candidate genes mapped at 12 selection regions (on 8 autosomes) could be linked to bovine stature diversity. Of these 12 candidate gene regions, three contained known genes (i.e., NCAPGLCORL, FBP2-PTCH1, and PLAG1-CHCHD7) related to bovine stature, and nine were not previously described in cattle (five in European and four in African cohorts). Overall, this study demonstrates the utility of CSS coupled with strategies of combining multibreed datasets in the identification and discovery of genomic regions underlying complex traits. Characterization of multiple signatures of selection and their underlying candidate genes will elucidate the polygenic nature of stature across cattle breeds

A list of selection tests used in published studies on Bovine selection signatures.

<p>A list of selection tests used in published studies on Bovine selection signatures.</p

Meta-assembly of selection signatures in four groups within European, Zebu, African and Composite breeds of cattle.

<p>Meta-assembly of selection signatures in four groups within European, Zebu, African and Composite breeds of cattle.</p

Comparison of Gene Editing Versus Conventional Breeding to Introgress the POLLED Allele Into the Tropically Adapted Australian Beef Cattle Population

Dehorning is the process of physically removing horns to protect animals and humans from injury, but the process is costly, unpleasant, and faces increasing public scrutiny. Genetic selection for polled (hornless), which is genetically dominant to horned, is a long-term solution to eliminate the need for dehorning. However, due to the limited number of polled Australian Brahman bulls, the northern Australian beef cattle population remains predominantly horned. The potential to use gene editing to produce high-genetic-merit polled cattle was recently demonstrated. To further explore the concept, this study simulated introgression of the POLLED allele into a tropically adapted Australian beef cattle population via conventional breeding or gene editing (top 1% or 10% of seedstock bulls/year) for 3 polled mating schemes and compared results to baseline selection on genetic merit (Japan Ox selection index, $JapOx) alone, over the course of 20 years. The baseline scenario did not significantly decrease the 20-year HORNED allele frequency (80$8.00/year). Compared to the baseline, the conventional breeding scenarios where polled bulls were preferentially used for breeding, regardless of their genetic merit, significantly decreased the 20-year HORNED allele frequency (30%), but resulted in a significantly slower rate of genetic gain ($6.70/year, P ≤ 0.05). The mating scheme that required the exclusive use of homozygous polled bulls, resulted in the lowest 20-year HORNED allele frequency (8$5.50/year). The addition of gene editing the top 1% or 10% of seedstock bull calves/year to each conventional breeding scenario resulted in significantly faster rates of genetic gain (up to $8.10/year, P ≤ 0.05). Overall, our study demonstrates that, due to the limited number of polled Australian Brahman bulls, strong selection pressure on polled will be necessary to meaningfully increase the number of polled animals in this population. Moreover, these scenarios illustrate how gene editing could be a tool for accelerating the development of high-genetic-merit homozygous polled sires to mitigate the current trade-off of slower genetic gain associated with decreasing HORNED allele frequency in the Australian Brahman population