Russell McCulloch

Genome-Wide Analysis of the World's Sheep Breeds Reveals High Levels of Historic Mixture and Strong Recent Selection

Genomic structure in a global collection of domesticated sheep reveals a history of artificial selection for horn loss and traits relating to pigmentation, reproduction, and body size.

A Genome Wide Survey of SNP Variation Reveals the Genetic Structure of Sheep Breeds

The genetic structure of sheep reflects their domestication and subsequent formation into discrete breeds. Understanding genetic structure is essential for achieving genetic improvement through genome-wide association studies, genomic selection and the dissection of quantitative traits. After identifying the first genome-wide set of SNP for sheep, we report on levels of genetic variability both within and between a diverse sample of ovine populations. Then, using cluster analysis and the partitioning of genetic variation, we demonstrate sheep are characterised by weak phylogeographic structure, overlapping genetic similarity and generally low differentiation which is consistent with their short evolutionary history. The degree of population substructure was, however, sufficient to cluster individuals based on geographic origin and known breed history. Specifically, African and Asian populations clustered separately from breeds of European origin sampled from Australia, New Zealand, Europe and North America. Furthermore, we demonstrate the presence of stratification within some, but not all, ovine breeds. The results emphasize that careful documentation of genetic structure will be an essential prerequisite when mapping the genetic basis of complex traits. Furthermore, the identification of a subset of SNP able to assign individuals into broad groupings demonstrates even a small panel of markers may be suitable for applications such as traceability.

Comprehensive mapping of the bull sperm surface proteome

While the mechanisms that underpin maturation, capacitation, and sperm–egg interactions remain elusive it is known that these essential fertilisation events are driven by the protein complement of the sperm surface. Understanding these processes is critical to the regulation of animal reproduction, but few studies have attempted to define the full repertoire of sperm surface proteins in animals of agricultural importance. Recent developments in proteomics technologies, subcellular fractionation, and optimised solubilisation strategies have enhanced the potential for the comprehensive characterisation of the sperm surface proteome. Here we report the identification of 419 proteins from a mature bull sperm plasma membrane fraction. Protein domain enrichment analyses indicate that 67% of all the proteins identified may be membrane associated. A large number of the proteins identified are conserved between mammalian species and are reported to play key roles in sperm–egg communication, capacitation and fertility. The major functional pathways identified were related to protein catabolism (26S proteasome complex), chaperonin‐containing TCP‐1 (CCT) complex and fundamental metabolic processes such as glycolysis and energy production. We have also identified 118 predicted transmembrane proteins, some of which are implicated in cell adhesion, acrosomal exocytosis, vesicle transport and immunity and fertilisation events, while others have not been reported in mammalian LC‐MS‐derived sperm proteomes to date. Comparative proteomics and functional network analyses of these proteins expand our systems level of understanding of the bull sperm proteome and provide important clues toward finding the essential conserved function of these proteins.

Genome-wide association study for inhibin, luteinizing hormone, insulin-like growth factor 1, testicular size and semen traits in bovine species.

The fertility of young bulls impacts on reproduction rates, farm profit and the rate of genetic progress in beef herds. Cattle researchers and industry therefore routinely collect data on the reproductive performance of bulls. Genome‐wide association studies were carried out to identify genomic regions and genes associated with reproductive traits measured during the pubertal development of Tropical Composite bulls, from 4 to 24 months of age. Data from 1 085 bulls were collected for seven traits: blood hormone levels of inhibin at 4 months (IN), luteinizing hormone following a gonadotropin releasing hormone challenge at 4 months (LH), insulin‐like growth factor 1 at 6 months (IGF1), scrotal circumference at 12 months (SC), sperm motility at 18 months (MOT), percentage of normal spermatozoa at 24 months (PNS) and age at a scrotal circumference of 26 cm (AGE26, or pubertal age). Data from 729 068 single‐nucleotide polymorphisms were used in the association analysis. Significant polymorphism associations were discovered for IN, IGF1, SC, AGE26 and PNS. Based on these associations, INHBE, INHBC and HELB are proposed as candidate genes for IN regulation. Polymorphisms associated with IGF1 mapped to the PLAG1 gene region, validating a reported quantitative trait locus on chromosome 14 for IGF1. The X chromosome contained most of the significant associations found for SC, AGE26 and PNS. These findings will contribute to the identification of diagnostic genetic markers and informed genomic selection strategies to assist breeding of cattle with improved fertility. Furthermore, this work provides evidence contributing to gene function annotation in the context of male fertility.

Genetic diversity and investigation of polledness in divergent goat populations using 52 088 SNPs

The recent availability of a genome-wide SNP array for the goat genome dramatically increases the power to investigate aspects of genetic diversity and to conduct genome-wide association studies in this important domestic species. We collected and analysed genotypes from 52 088 SNPs in Boer, Cashmere and Rangeland goats that had both polled and horned individuals. Principal components analysis revealed a clear genetic division between animals for each population, and model-based clustering successfully detected evidence of admixture that matched aspects of their recorded history. For example, shared co-ancestry was detected, suggesting Boer goats have been introgressed into the Rangeland population. Further, allele frequency data successfully tracked the altered genetic profile that has taken place after 40 years of breeding Australian Cashmere goats using the Rangeland animals as the founding population. Genome-wide association mapping of the POLL locus revealed a strong signal on goat chromosome 1. The 769-kb critical interval contained the polled intersex syndrome locus, confirming the genetic basis in non-European animals is the same as identified previously in Saanen goats. Interestingly, analysis of the haplotypes carried by a small set of sex-reversed animals, known to be associated with polledness, revealed some animals carried the wild-type chromosome associated with the presence of horns. This suggests a more complex basis for the relationship between polledness and the intersex condition than initially thought while validating the application of the goat SNP50 BeadChip for fine-mapping traits in goat.

Linkage disequilibrium over short physical distances measured in sheep using a high‐density SNP chip

The extent of linkage disequilibrium (LD) between genetic loci has implications for both association studies and the accuracy of genomic prediction. To characterise the persistence of LD in diverse sheep breeds, two SNP genotyping platforms were used. First, existing SNP genotypes from 63 breeds obtained using the ovine SNP50 BeadChip (49,034 loci) were used to estimate LD decay in populations with contrasting levels of genetic diversity. Given the paucity of marker pairs separated by short physical distances on the SNP50 BeadChip, genotyping was subsequently performed for four breeds using the recently developed ovine HD BeadChip that assays approximately 600,000 SNPs with an average genomic spacing of 5 kb. This facilitated a highly accurate estimate of LD over short genomic distances (<30 kb) and revealed LD varies considerably between sheep breeds. Further, sheep appear to contain generally lower levels of LD than do other domestic species, likely a reflection of aspects of their past population history.

Reference genome of wild goat (capra aegagrus) and sequencing of goat breeds provide insight into genic basis of goat domestication

BackgroundDomestic goats (Capra hircus) have been selected to play an essential role in agricultural production systems, since being domesticated from their wild progenitor, bezoar (Capra aegagrus). A detailed understanding of the genetic consequences imparted by the domestication process remains a key goal of evolutionary genomics.ResultsWe constructed the reference genome of bezoar and sequenced representative breeds of domestic goats to search for genomic changes that likely have accompanied goat domestication and breed formation. Thirteen copy number variation genes associated with coat color were identified in domestic goats, among which ASIP gene duplication contributes to the generation of light coat-color phenotype in domestic goats. Analysis of rapidly evolving genes identified genic changes underlying behavior-related traits, immune response and production-related traits.ConclusionBased on the comparison studies of copy number variation genes and rapidly evolving genes between wild and domestic goat, our findings and methodology shed light on the genetic mechanism of animal domestication and will facilitate future goat breeding.

Photochemically crosslinked matrices of gelatin and fibrinogen promote rapid cell proliferation

Here we report the use of a facile photochemical crosslinking method to fabricate stable polymer matrices from unmodified gelatin and fibrinogen. Gels were produced by covalent crosslinking of the proteins in a rapid photo‐oxidative process, catalysed by a ruthenium metal complex and irradiation with visible light. For generation of macroporous, spongy matrices, the proteins and crosslinking reagents were mixed with catalase and hydrogen peroxide to achieve a foaming reaction, producing a stable, foamed matrix that was subsequently photo‐crosslinked. C2C12 cells were either seeded onto the matrices after photo‐curing or embedded in the protein matrix prior to foaming and crosslinking. Cells seeded onto scaffolds post‐curing showed high cell viability and rapid proliferation in vitro. For cells embedded in the matrix prior to crosslinking there was some loss of initial viability, but surviving cells were able to proliferate after a period of in vitro cultivation. The matrices were shown to be biocompatible when implanted into nude mice, with evidence of proliferation and differentiation of cells seeded into the scaffolds. The results are promising for further development of tissue‐engineering scaffolds based on this ruthenium‐catalysed photo‐crosslinking method. Copyright

East Friesian sheep carry a Myostatin allele known to cause muscle hypertrophy in other breeds

The East Friesian breed of sheep was developed in northern Germany and the Netherlands, and has become one of the world’s most productive dairy sheep. It is likely to have contributed to the foundation of other breeds, such as the Texel, which originated in the Netherland’s chain of West Friesian islands. The Texel is a meat breed that displays a muscle hypertrophy phenotype caused by a G to A substitution (g.+6723G>A) in the Myostatin gene.1 Given the likelihood of a common population history linking Texel and East Frisian, we sought to determine if the latter also carries the mutant g.+6723A Myostatin allele despite the divergent production profiles of the two breeds.

Numerical analysis of intensity signals resulting from genotyping pooled DNA samples in beef cattle and broiler chicken1

Pooled genomic DNA has been proposed as a cost-effective approach in genomewide association studies (GWAS). However, algorithms for genotype calling of biallelic SNP are not adequate with pooled DNA samples because they assume the presence of 2 fluorescent signals, 1 for each allele, and operate under the expectation that at most 2 copies of the variant allele can be found for any given SNP and DNA sample. We adapt analytical methodology from 2-channel gene expression microarray technology to SNP genotyping of pooled DNA samples. Using 5 datasets from beef cattle and broiler chicken of varying degrees of complexity in terms of design and phenotype, continuous and dichotomous, we show that both differential hybridization (M = green minus red intensity signal) and abundance (A = average of red and green intensities) provide useful information in the prediction of SNP allele frequencies. This is predominantly true when making inference about extreme SNP that are either nearly fixed or highly polymorphic. We propose the use of model-based clustering via mixtures of bivariate normal distributions as an optimal framework to capture the relationship between hybridization intensity and allele frequency from pooled DNA samples. The range of M and A values observed here are in agreement with those reported within the context of gene expression microarray and also with those from SNP array data within the context of analytical methodology for the identification of copy number variants. In particular, we confirm that highly polymorphic SNP yield a strong signal from both channels (red and green) while lowly or nonpolymorphic SNP yield a strong signal from 1 channel only. We further confirm that when the SNP allele frequencies are known, either because the individuals in the pools or from a closely related population are themselves genotyped, a multiple regression model with linear and quadratic components can be developed with high prediction accuracy. We conclude that when these approaches are applied to the estimation of allele frequencies, the resulting estimates allow for the development of cost-effective and reliable GWAS.