An update on sheep and goat linkage maps and other genomic resources

Identification of genomic regions that exhibit sexual dimorphism for size and muscularity in cattle

Sexual dimorphism, the phenomenon whereby males and females of the same species are distinctive in some aspect of appearance or size, has previously been documented in cattle for traits such as growth rate and carcass merit using a quantitative genetics approach. No previous study in cattle has attempted to document sexual dimorphism at a genome level; therefore, the objective of the present study was to determine whether genomic regions associated with size and muscularity in cattle exhibited signs of sexual dimorphism. Analyses were undertaken on 10 linear-type traits that describe the muscular and skeletal characteristics of both males and females of five beef cattle breeds: 1,444 Angus (AA), 6,433 Charolais (CH), 1,129 Hereford, 8,745 Limousin (LM), and 1,698 Simmental. Genome-wide association analyses were undertaken using imputed whole-genome sequence data for each sex separately by breed. For each single-nucleotide polymorphism (SNP) that was segregating in both sexes, the difference between the allele substitution effect sizes for each sex, in each breed separately, was calculated. Suggestively (P ≤ 1 × 10-5) sexually dimorphic SNPs that were segregating in both males and females were detected for all traits in all breeds, although the location of these SNPs differed by both trait and breed. Significantly (P ≤ 1 × 10-8) dimorphic SNPs were detected in just three traits in the AA, seven traits in the CH, and three traits in the LM. The vast majority of all segregating autosomal SNPs (86% in AA to 94% in LM) had the same minor allele in both males and females. Differences (P ≤ 0.05) in allele frequencies between the sexes were observed for between 36% (LM) and 66% (AA) of the total autosomal SNPs that were segregating in both sexes. Dimorphic SNPs were located within a number of genes related to muscularity and/or size including the NAB1, COL5A2, and IWS1 genes on BTA2 that are located close to, and thought to be co-inherited with, the MSTN gene. Overall, sexual dimorphism exists in cattle at the genome level, but it is not consistent by either trait or breed.

Cattle Sex-Specific Recombination and Genetic Control from a Large Pedigree Analysis

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.

Detection of chromosomal segments underlying scrotal circumference in ram lambs and age at onset of puberty in ewe lambs

Scrotal circumference (SC) is highly correlated with total sperm production and quality of the produced semen. In ewe lambs, puberty is an important reproductive trait and extreme delay in reaching puberty will have a negative effect on breeding efficiency. To identify genomic regions (QTL) underlying variation in SC and age at onset of puberty in ewe lambs (AP), a whole genome scan was performed with 169 microsatellites covering the ovine autosomes. Progeny (360 animals) from six half-sib families in a population of Kermani sheep, an indigenous fat tailed sheep breed in south-east of Iran, were recorded for SC and AP. Linkage analysis using the interval mapping approach indicated a total of 21 QTL for the studied traits on nine chromosomes (OAR1, 2, 3, 5, 6, 9, 11, 18 and 24), including a 5% genome wide significant QTL for SC on OAR5. There was evidence for the overlap of QTL affecting SC and AP on OAR5 and OAR9. The size of the QTL effect ranged from 0.8 to 1.6 in units of the phenotypic standard deviation. Although QTL for testis size in ram lambs have been mapped, this is the first report of QTL for age at onset of puberty in ewe lambs. Pursuing the underlying causes of these QTL may shed light on the molecular aspects of male and female reproductive development in sheep.

Strong artificial selection in domestic mammals did not result in an increased recombination rate

Recombination rates vary in intensity and location at the species, individual, sex and chromosome levels. Despite the fundamental biological importance of this process, the selective forces that operate to shape recombination rate and patterns are unclear. Domestication offers a unique opportunity to study the interplay between recombination and selection. In domesticates, intense selection for particular traits is imposed on small populations over many generations, resulting in organisms that differ, sometimes dramatically, in morphology and physiology from their wild ancestor. Although earlier studies suggested increased recombination rate in domesticates, a formal comparison of recombination rates between domestic mammals and their wild congeners was missing. In order to determine broad-scale recombination rate, we used immunolabeling detection of MLH1 foci as crossover markers in spermatocytes in three pairs of closely related wild and domestic species (dog and wolf, goat and ibex, and sheep and mouflon). In the three pairs, and contrary to previous suggestions, our data show that contemporary recombination rate is higher in the wild species. Subsequently, we inferred recombination breakpoints in sequence data for 16 genomic regions in dogs and wolves, each containing a locus associated with a dog phenotype potentially under selection during domestication. No difference in the number and distribution of recombination breakpoints was found between dogs and wolves. We conclude that our data indicate that strong directional selection did not result in changes in recombination in domestic mammals, and that both upper and lower bounds for crossover rates may be tightly regulated.

The sex-specific impact of meiotic recombination on nucleotide composition

Meiotic recombination is an important evolutionary force shaping the nucleotide landscape of genomes. For most vertebrates, the frequency of recombination varies slightly or considerably between the sexes (heterochiasmy). In humans, male, rather than female, recombination rate has been found to be more highly correlated with the guanine and cytosine (GC) content across the genome. In the present study, we review the results in human and extend the examination of the evolutionary impact of heterochiasmy beyond primates to include four additional eutherian mammals (mouse, dog, pig, and sheep), a metatherian mammal (opossum), and a bird (chicken). Specifically, we compared sex-specific recombination rates (RRs) with nucleotide substitution patterns evaluated in transposable elements. Our results, based on a comparative approach, reveal a great diversity in the relationship between heterochiasmy and nucleotide composition. We find that the stronger male impact on this relationship is a conserved feature of human, mouse, dog, and sheep. In contrast, variation in genomic GC content in pig and opossum is more strongly correlated with female, rather than male, RR. Moreover, we show that the sex-differential impact of recombination is mainly driven by the chromosomal localization of recombination events. Independent of sex, the higher the RR in a genomic region and the longer this recombination activity is conserved in time, the stronger the bias in nucleotide substitution pattern, through such mechanisms as biased gene conversion. Over time, this bias will increase the local GC content of the region.

Extended cytogenetic maps of sheep chromosome 1 and their cattle and river buffalo homoeologues: comparison with the OAR1 RH map and human chromosomes 2, 3, 21 and 1q

Cytogenetic maps are useful tools for several applications, such as the physical anchoring of linkage and RH maps or genome sequence contigs to specific chromosome regions or the analysis of chromosome rearrangements. Recently, a detailed RH map was reported in OAR1. In the present study, we selected 38 markers equally distributed in this RH map for identification of ovine genomic DNA clones within the ovine BAC library CHORI-243 using the virtual sheep genome browser and performed FISH mapping for both comparison of OAR1 and homoeologous chromosomes BBU1q-BBU6 and BTA1-BTA3 and considerably extending the cytogenetic maps of the involved species-specific chromosomes. Comparison of the resulting maps with human-identified homology with HSA2q, HSA3, HSA21 and HSA1q reveals complex chromosome rearrangements differentiating human and bovid chromosomes. In addition, we identified 2 new small human segments from HSA2q and HSA3q conserved in the telomeric regions of OAR1p and homoeologous chromosome regions of BTA3 and BBU6, and OAR1q, respectively. Evaluation of the present OAR1 cytogenetic map and the OAR1 RH map supports previous RH assignments with 2 main exceptions. The 2 loci BMS4011 and CL638002 occupy inverted positions in these 2 maps.

Genetic mapping of quantitative trait loci for aseasonal reproduction in sheep

The productivity and economic prosperity of sheep farming could benefit greatly from more effective methods of selection for year-round lambing. Identification of QTL for aseasonal reproduction in sheep could lead to more accurate selection and faster genetic improvement. One hundred and twenty microsatellite markers were genotyped on 159 backcross ewes from a Dorset × East Friesian crossbred pedigree. Interval mapping was undertaken to map the QTL underlying several traits describing aseasonal reproduction including the number of oestrous cycles, maximum level of progesterone prior to breeding, pregnancy status determined by progesterone level, pregnancy status determined by ultrasound, lambing status and number of lambs born. Seven chromosomes (1, 3, 12, 17, 19, 20 and 24) were identified to harbour putative QTL for one or more component traits used to describe aseasonal reproduction. Ovine chromosomes 12, 17, 19 and 24 harbour QTL significant at the 5% chromosome-wide level, chromosomes 3 and 20 harbour QTL that exceeded the threshold at the 1% chromosome-wide level, while the QTL identified on chromosome 1 exceeded the 1% experiment-wide significance level. These results are a first step towards understanding the genetic mechanism of this complex trait and show that variation in aseasonal reproduction is associated with multiple chromosomal regions.

Chromosomal localisation and genetic variation of the SLC11A1 gene in goats (Capra hircus)

The solute carrier family 11 member A1 (SLC11A1) gene is associated with resistance to infectious diseases. Chromosomal localisation, genomic regions corresponding to functional domains and the genetic variability of microsatellites in the 3' untranslated region (3'-UTR) of this gene were investigated in 427 goats (Capra hircus) of six breeds. Using dual colour fluorescence in situ hybridisation, SLC11A1 was localised to goat chromosome 2. Single strand conformation polymorphism was used to screen for polymorphisms in SLC11A1 exons 2, 10 and 15. There was no variation among goat breeds in the sarcoma homology 3 (SH3) binding motif, the protein kinase C phosphorylation site or the two N-linked glycosylation sites. Exon 15 exhibited variability due to the presence of two polymorphic microsatellites. Genotyping of the upstream guanine-thymine repeat (GTn) at 3'-UTR revealed eight alleles (GT11, GT12, GT14-GT19) in goats, whereas GT13 (present in cattle) was absent. Most goats carried the GT16 allele and no allele was found to be exclusive to only one breed. The coefficient of genetic differentiation value (G(ST)) was 0.084. This microsatellite appears to be an informative DNA marker for genetic linkage analysis in goats.

Mapping quantitative trait loci (QTL) in sheep. II. Meta-assembly and identification of novel QTL for milk production traits in sheep

An (Awassi × Merino) × Merino backcross family of 172 ewes was used to map quantitative trait loci (QTL) for different milk production traits on a framework map of 200 loci across all autosomes. From five previously proposed mathematical models describing lactation curves, the Wood model was considered the most appropriate due to its simplicity and its ability to determine ovine lactation curve characteristics. Derived milk traits for milk, fat, protein and lactose yield, as well as percentage composition and somatic cell score were used for single and two-QTL approaches using maximum likelihood estimation and regression analysis. A total of 15 significant (P < 0.01) and additional 25 suggestive (P < 0.05) QTL were detected across both single QTL methods and all traits. In preparation of a meta-analysis, all QTL results were compared with a meta-assembly of QTL for milk production traits in dairy ewes from various public domain sources and can be found on the ReproGen ovine gbrowser http://crcidp.vetsci.usyd.edu.au/cgi-bin/gbrowse/oaries_genome/. Many of the QTL for milk production traits have been reported on chromosomes 1, 3, 6, 16 and 20. Those on chromosomes 3 and 20 are in strong agreement with the results reported here. In addition, novel QTL were found on chromosomes 7, 8, 9, 14, 22 and 24. In a cross-species comparison, we extended the meta-assembly by comparing QTL regions of sheep and cattle, which provided strong evidence for synteny conservation of QTL regions for milk, fat, protein and somatic cell score data between cattle and sheep.

Mapping quantitative trait loci (QTL) in sheep. I. A new male framework linkage map and QTL for growth rate and body weight

A male sheep linkage map comprising 191 microsatellites was generated from a single family of 510 Awassi-Merino backcross progeny. Except for ovine chromosomes 1, 2, 10 and 17, all other chromosomes yielded a LOD score difference greater than 3.0 between the best and second-best map order. The map is on average 11% longer than the Sheep Linkage Map v4.7 male-specific map. This map was employed in quantitative trait loci (QTL) analyses on body-weight and growth-rate traits between birth and 98 weeks of age. A custom maximum likelihood program was developed to map QTL in half-sib families for non-inbred strains (QTL-MLE) and is freely available on request. The new analysis package offers the advantage of enabling QTL × fixed effect interactions to be included in the model. Fifty-four putative QTL were identified on nine chromosomes. Significant QTL with sex-specific effects (i.e. QTL × sex interaction) in the range of 0.4 to 0.7 SD were found on ovine chromosomes 1, 3, 6, 11, 21, 23, 24 and 26.

Genome-wide cross-amplification of domestic sheep microsatellites in bighorn sheep and mountain goats

We tested for cross-species amplification of microsatellite loci located throughout the domestic sheep (Ovis aries) genome in two north American mountain ungulates (bighorn sheep, Ovis canadensis, and mountain goats, Oreamnos americanus). We identified 247 new polymorphic markers in bighorn sheep (≥ 3 alleles in one of two study populations) and 149 in mountain goats (≥ 2 alleles in a single study population) using 648 and 576 primer pairs, respectively. Our efforts increased the number of available polymorphic microsatellite markers to 327 for bighorn sheep and 180 for mountain goats. The average distance between successive polymorphic bighorn sheep and mountain goat markers inferred from the Australian domestic sheep genome linkage map (mean ± 1 SD) was 11.9 ± 9.2 and 15.8 ± 13.8 centimorgans, respectively. The development of genomic resources in these wildlife species enables future studies of the genetic architecture of trait variation.

A high-resolution comparative radiation hybrid map of ovine chromosomal regions that are homologous to human chromosome 6 (HSA6)

In this study, we constructed high-resolution radiation hybrid (RH) and comparative maps of ovine chromosomes or chromosomal segments that are homologous to human chromosome 6 (HSA6). A total of 251 markers were successfully genotyped across the recently developed USUoRH5000 whole-genome panel; 208 of these markers were assigned to five RH linkage groups distributed on three ovine chromosomes (OAR8, 9 and 20). The RH maps have good correspondence with previous chromosome painting data, although a small centromeric region on OAR9 that is homologous to HSA6 had not been previously detected using human chromosome paints on ovine chromosomal spreads. High percentages of the ovine markers were identified as orthologues in the bovine (86.3%), dog (85.8%), horse (69.3%) and human (88.7%) genomes. These maps contribute to investigations in mammalian chromosome evolution and the search for economic trait loci in sheep.

Biomedical applications of sheep models: from asthma to vaccines

Although rodent models are very popular for scientific studies, it is becoming more evident that large animal models can provide unique opportunities for biomedical research. Sheep are docile in nature and large in size, which facilitates surgical manipulation, and their physiology is similar to humans. As a result, for decades they have been chosen for several models and continue to be used to study an ever-increasing array of applications. Despite this, their full potential has not been exploited. Here, we review the use of sheep as an animal model for human vaccine development, asthma pathogenesis and treatment, the study of neonatal development, and the optimization of drug delivery and surgical techniques.