SNP芯片数据估计动物个体基因组品种构成的方法及应用

自然和人工选择、地理隔离和遗传漂移等原因使动物基因组中许多位点的等位基因频率在群体间会产生差异。源于不同品种(祖先)杂交(交配)的动物个体,其基因组与这些品种(祖先)的基因频率(基因型)会存在一定的相关性。因此采用合适的统计模型和分析方法,可以估计出每个品种(祖先)对于个体基因组的遗传贡献比例,又称为个体的基因组品种构成(genomic breed composition, GBC)。本文介绍了利用SNP芯片数据估计动物个体GBC的原理、方法及步骤,并且通过对198头待鉴定的日本红毛和牛GBC的评估,演示了用回归模型和混合分布模型估计动物个体GBC的具体步骤,其中包括SNP子集的筛选、参考群体中动物个体选择以及待测定动物GBC的计算。参考动物群体选自日本红毛和牛(Akaushi)、安格斯牛(Angus)、海福特牛(Hereford)、荷斯坦牛(Holstein)和娟珊牛(Jersey) 5个品种共36 574头,每个个体有40K或50K芯片数据。本文在现有商用 SNP芯片基础上筛选用于品种鉴定和估计动物个体GBC的SNP子集,是对现有SNP芯片功能的拓展和深入开发利用。此外,在基因组选择中如何利用SNP基因型估计动物个体GBC的结果,提高纯种和杂种动物的预测准确度,也是值得深入研究的领域。     

Accuracies of genomically estimated breeding values from pure-breed and across-breed predictions in Australian beef cattle

Background

The major obstacles for the implementation of genomic selection in Australian beef cattle are the variety of breeds and in general, small numbers of genotyped and phenotyped individuals per breed. The Australian Beef Cooperative Research Center (Beef CRC) investigated these issues by deriving genomic prediction equations (PE) from a training set of animals that covers a range of breeds and crosses including Angus, Murray Grey, Shorthorn, Hereford, Brahman, Belmont Red, Santa Gertrudis and Tropical Composite. This paper presents accuracies of genomically estimated breeding values (GEBV) that were calculated from these PE in the commercial pure-breed beef cattle seed stock sector.

Methods

PE derived by the Beef CRC from multi-breed and pure-breed training populations were applied to genotyped Angus, Limousin and Brahman sires and young animals, but with no pure-breed Limousin in the training population. The accuracy of the resulting GEBV was assessed by their genetic correlation to their phenotypic target trait in a bi-variate REML approach that models GEBV as trait observations.

Results

Accuracies of most GEBV for Angus and Brahman were between 0.1 and 0.4, with accuracies for abattoir carcass traits generally greater than for live animal body composition traits and reproduction traits. Estimated accuracies greater than 0.5 were only observed for Brahman abattoir carcass traits and for Angus carcass rib fat. Averaged across traits within breeds, accuracies of GEBV were highest when PE from the pooled across-breed training population were used. However, for the Angus and Brahman breeds the difference in accuracy from using pure-breed PE was small. For the Limousin breed no reasonable results could be achieved for any trait.

Conclusion

Although accuracies were generally low compared to published accuracies estimated within breeds, they are in line with those derived in other multi-breed populations. Thus PE developed by the Beef CRC can contribute to the implementation of genomic selection in Australian beef cattle breeding.     

Evidence of positive selection towards Zebuine haplotypes in the BoLA region of Brangus cattle

The Brangus breed was developed to combine the superior characteristics of both of its founder breeds, Angus and Brahman. It combines the high adaptability to tropical and subtropical environments, disease resistance, and overall hardiness of Zebu cattle with the reproductive potential and carcass quality of Angus. It is known that the major histocompatibility complex (MHC, also known as bovine leucocyte antigen: BoLA), located on chromosome 23, encodes several genes involved in the adaptive immune response and may be responsible for adaptation to harsh environments. The objective of this work was to evaluate whether the local breed ancestry percentages in the BoLA locus of a Brangus population diverged from the estimated genome-wide proportions and to identify signatures of positive selection in this genomic region. For this, 167 animals (100 Brangus, 45 Angus and 22 Brahman) were genotyped using a high-density single nucleotide polymorphism array. The local ancestry analysis showed that more than half of the haplotypes (55.0%) shared a Brahman origin. This value was significantly different from the global genome-wide proportion estimated by cluster analysis (34.7% Brahman), and the proportion expected by pedigree (37.5% Brahman). The analysis of selection signatures by genetic differentiation (F_st) and extended haplotype homozygosity-based methods (iHS and Rsb) revealed 10 and seven candidate regions, respectively. The analysis of the genes located within these candidate regions showed mainly genes involved in immune response-related pathway, while other genes and pathways were also observed (cell surface signalling pathways, membrane proteins and ion-binding proteins). Our results suggest that the BoLA region of Brangus cattle may have been enriched with Brahman haplotypes as a consequence of selection processes to promote adaptation to subtropical environments.     

Fine‐mapping the POLL locus in Brahman cattle yields the diagnostic marker CSAFG29

The POLL locus has been mapped to the centromeric region of bovine chromosome 1 (BTA1) in both taurine breeds and taurine–indicine crosses in an interval of approximately 1 Mb. It has not yet been mapped in pure‐bred zebu cattle. Despite several efforts, neither causative mutations in candidate genes nor a singular diagnostic DNA marker has been identified. In this study, we genotyped a total of 68 Brahman cattle and 20 Hereford cattle informative for the POLL locus for 33 DNA microsatellites, 16 of which we identified de novo from the bovine genome sequence, mapping the POLL locus to the region of the genes IFNAR2 and SYNJ1. The 303‐bp allele of the new microsatellite, CSAFG29, showed strong association with the POLL allele. We then genotyped 855 Brahman cattle for CSAFG29 and confirmed the association between the 303‐bp allele and POLL. To determine whether the same association was found in taurine breeds, we genotyped 334 animals of the Angus, Hereford and Limousin breeds and 376 animals of the Brangus, Droughtmaster and Santa Gertrudis composite taurine–zebu breeds. The association between the 303‐bp allele and POLL was confirmed in these breeds; however, an additional allele (305 bp) was also associated but not fully predictive of POLL. Across the data, CSAFG29 was in sufficient linkage disequilibrium to the POLL allele in Australian Brahman cattle that it could potentially be used as a diagnostic marker in that breed, but this may not be the case in other breeds. Further, we provide confirmatory evidence that the scur phenotype generally occurs in animals that are heterozygous for the POLL allele.     

Ultra-low-density genotype panels for breed assignment of Angus and Hereford cattle

Angus and Hereford beef is marketed internationally for apparent superior meat quality attributes; DNA-based breed authenticity could be a useful instrument to ensure consumer confidence on premium meat products. The objective of this study was to develop an ultra-low-density genotype panel to accurately quantify the Angus and Hereford breed proportion in biological samples. Medium-density genotypes (13 306 single nucleotide polymorphisms (SNPs)) were available on 54 703 commercial and 4042 purebred animals. The breed proportion of the commercial animals was generated from the medium-density genotypes and this estimate was regarded as the gold-standard breed composition. Ten genotype panels (100 to 1000 SNPs) were developed from the medium-density genotypes; five methods were used to identify the most informative SNPs and these included the Delta statistic, the fixation (F_st) statistic and an index of both. Breed assignment analyses were undertaken for each breed, panel density and SNP selection method separately with a programme to infer population structure using the entire 13 306 SNP panel (representing the gold-standard measure). Breed assignment was undertaken for all commercial animals (n=54 703), animals deemed to contain some proportion of Angus based on pedigree (n=5740) and animals deemed to contain some proportion of Hereford based on pedigree (n=5187). The predicted breed proportion of all animals from the lower density panels was then compared with the gold-standard breed prediction. Panel density, SNP selection method and breed all had a significant effect on the correlation of predicted and actual breed proportion. Regardless of breed, the Index method of SNP selection numerically (but not significantly) outperformed all other selection methods in accuracy (i.e. correlation and root mean square of prediction) when panel density was ⩾300 SNPs. The correlation between actual and predicted breed proportion increased as panel density increased. Using 300 SNPs (selected using the global index method), the correlation between predicted and actual breed proportion was 0.993 and 0.995 in the Angus and Hereford validation populations, respectively. When SNP panels optimised for breed prediction in one population were used to predict the breed proportion of a separate population, the correlation between predicted and actual breed proportion was 0.034 and 0.044 weaker in the Hereford and Angus populations, respectively (using the 300 SNP panel). It is necessary to include at least 300 to 400 SNPs (per breed) on genotype panels to accurately predict breed proportion from biological samples.     

Response to GnRH on day 6 of the estrous cycle is diminished as the percentage of Bos indicus breeding increases in Angus, Brangus, and Brahman x Angus heifers

Angus (n=6), Brangus (5/8 Angus x 3/8 Brahman, n=6), and Brahman x Angus (3/8 Angus x 5/8 Brahman, n=6) heifers exhibiting estrous cycles at regular intervals were used to determine if the percentage of Bos indicus breeding influenced the secretory patterns of LH in response to a GnRH treatment on Day 6 of the estrous cycle. Heifers were pre-synchronized with a two-injection PGF(2 alpha) protocol (25 mg i.m. Day -14 and 12.5 mg i.m. Day -3 and -2 of experiment). Heifers received 100 microg GnRH i.m. on Day 6 of the subsequent estrous cycle. Blood samples were collected at -60, -30, and -1 min before GnRH and 15, 30, 60, 90, 120, 150, 180, 240, 300, 360, 420, and 480 min after GnRH to determine concentrations of serum LH. Estradiol concentrations were determined at -60, -30, and -1 min before GnRH. On Day 6 and 8, ovaries were examined by ultrasonography to determine if ovulation occurred. On Day 13, heifers received 25 mg PGF(2 alpha) i.m. and blood samples were collected daily until either the expression of estrus or Day 20 for heifers not exhibiting estrus to determine progesterone concentrations. There was no effect (P>0.10) of breed on ovulation rate to GnRH as well as size of the largest follicle, mean estradiol, and mean corpus luteum volume at GnRH. Mean LH was greater (P<0.05) for Angus (7.0+/-0.8 ng/mL) compared to Brangus (4.6+/-0.8 ng/mL) and Brahman x Angus (2.9+/-0.8 ng/mL), which were similar (P>0.10). Mean LH peak-height was similar (P>0.10) for Brangus (13.9+/-3.4 ng/mL) compared to Angus (21.9+/-3.4 ng/mL) and Brahman x Angus (8.0+/-3.4 ng/mL), but was greater (P<0.05) for Angus compared to Brahman x Angus. Interval from GnRH to LH peak was similar (P>0.10) between breeds. As the percentage of Bos indicus breeding increased the amount of LH released in response to GnRH on Day 6 of the estrous cycle decreased.     

Genotype imputation in the domestic dog

Application of imputation methods to accurately predict a dense array of SNP genotypes in the dog could provide an important supplement to current analyses of array-based genotyping data. Here, we developed a reference panel of 4,885,283 SNPs in 83 dogs across 15 breeds using whole genome sequencing. We used this panel to predict the genotypes of 268 dogs across three breeds with 84,193 SNP array-derived genotypes as inputs. We then (1) performed breed clustering of the actual and imputed data; (2) evaluated several reference panel breed combinations to determine an optimal reference panel composition; and (3) compared the accuracy of two commonly used software algorithms (Beagle and IMPUTE2). Breed clustering was well preserved in the imputation process across eigenvalues representing 75 % of the variation in the imputed data. Using Beagle with a target panel from a single breed, genotype concordance was highest using a multi-breed reference panel (92.4 %) compared to a breed-specific reference panel (87.0 %) or a reference panel containing no breeds overlapping with the target panel (74.9 %). This finding was confirmed using target panels derived from two other breeds. Additionally, using the multi-breed reference panel, genotype concordance was slightly higher with IMPUTE2 (94.1 %) compared to Beagle; Pearson correlation coefficients were slightly higher for both software packages (0.946 for Beagle, 0.961 for IMPUTE2). Our findings demonstrate that genotype imputation from SNP array-derived data to whole genome-level genotypes is both feasible and accurate in the dog with appropriate breed overlap between the target and reference panels.     

Variation in the XKR4 gene was significantly associated with subcutaneous rump fat thickness in indicine and composite cattle

Variation in the XK, Kell blood group complex subunit–related family, member 4 (XKR4) gene on BTA14 was associated with rump fat thickness in a recent genome‐wide association study. This region is also of interest because it is known to show evidence of a signature of population genetic selection. In this study, additional variation in this gene was genotyped in a sample of a total of 1283 animals of the Belmont Red (BEL) and Santa Gertrudis (SGT) breeds. The SNP rs41724387 was significantly (P < 0.001) associated with rump fat thickness and explained 5.9% of the genetic variance for the trait in this sample. Using the 4466 genotypes for the SNP rs42646708 from several data sets to estimate effects in seven breeds, this relatively large quantitative trait locus effect appears to be a result of the variation in indicine and taurine–indicine composite cattle. However, the only DNA variant found in Brahman cattle that altered the predicted amino acid sequence of XKR4 was not associated with rump fat thickness. This suggests that causative mutations lie outside the coding sequence of this gene.     

Characterization of experimental cattle populations from Argentina with a low-density SNP genotyping panel

The objectives of the present experiment were to evaluate a low-density SNP array designed for the molecular characterisation of gene banks and to assess the genetic diversity and population structure of beef cattle herds from an Argentinean research station. Forty-nine animals from three breeds (Angus, Hereford, and Argentinean Creole) were genotyped using the multi-species IMAGE001 60-K SNP array (10 K for cattle). Genotypes of other 19 cattle populations from Argentina, other American countries, and Europe were included in the study. Of special interest was the characterization of the Argentinean Creole, the only autochthonous cattle breed in the country. Due to the merging of different datasets, approximately 5 K SNPs were effectively used. Genetic differentiation (F_ST), principal component analysis, neighbour-joining tree of Reynolds distances and ancestry analysis showed that autochthonous American breeds are clearly differentiated, but all have genetic influences of Iberian cattle. The analysed herds of Argentinean Creole showed no evidence of recent admixture and represent a unique genetic pool within local American breeds. An experimental herd and the local commercial Hereford population have also diverged, probably due to the influence of current selection objectives in the breed. Our results illustrate the utility of using low-cost, low density SNP arrays in the evaluation of animal genetic resources. This type of panels could become a very useful resource in developing countries, where most endangered cattle breeds are located. The results also reinforce the importance of experimental herds as reservoir of genetic diversity, particularly in the case of local breeds under-represented in traditional production systems.     

Growth and reproductive development from weaning through 20 months of age among breeds of bulls in subtropical Florida

To determine the effect of breed on growth and reproductive development, weaned bulls in each of 2 yr were managed as a single group for approximately a year. In Year 1, the study group consisted of 24 Angus, 24 Brahman, 20 Hereford and 14 Senepol bulls, while in Year 2, it contained 25 Angus, 17 Brahman. 13 Romosinuano and 9 Nellore x Brahman bulls. Body and testicular growth measurements were recorded at 6-wk intervals. At approximately 1 yr of age and quarterly thereafter (4 periods), bulls were evaluated for libido, pubertal status, and GnRH-induced LH and testosterone secretion. Significant breed-by-age interactions occurred for most growth measurements. Brahman bulls (Bos indicus ) were (P < 0.05) older and heavier at puberty than Angus, Hereford, Senepol and Romosinuano bulls (Bos taurus ). Libido scores were lowest for Brahman and Nell ore x Brahman bulls (Bos indicus ). highest for Angus and Hereford bulls (temperate Bos taurus breeds) and intermediate for Senepol and Romosinuano bulls (tropical Bos taurus breeds; P < 0.05). Differences were not consistent among breeds or between years for GnRH-induced LH secretion. In both years, basal testosterone concentrations and areas under the GnRH-induced testosterone curve were higher (P < 0.05) for Angus and Hereford bulls (temperate breeds) than for Brahman, Senepol, Romosinuano and Nellore x Brahman bulls (tropical breeds). In conclusion, reproductive development of Senepol and Romosinuano bulls (tropical Bos taurus breeds) was more similar to Angus and Hereford bulls (temperate Bos taurus breeds) than to Brahman and Nellore x Brahman bulls (Bos indicus ).     

Bull selection and use in northern Australia. 4. Calf output and predictors of fertility of bulls in multiple-sire herds

On 10 northern Australian properties, the number of calves sired (calf output) by individual bulls in multiple-sire matings was measured by DNA typing for paternity. There were 235 bulls (92 Santa Gertrudis, 25 5/8 Brahman and 119 Brahman) from 37 multiple-sire mating groups. Number of bulls in groups ranged from 2 to 25 and ages of bulls ranged from 2 to 5 years. Mating periods were for 3-12 months and bull mating percentages were 2.5-6%. In all, there were 4251 calves tested and the resolution of paternity ranged from 92.5 to 100% and averaged 97.7% across all sites. This included 9.9% of calves with no potential sires in any of the mating groups. Of the 235 bulls mated, 58% sired 10% or less calves in each of their respective mating groups with 6% not siring any calves. In contrast, 14% sired over 30% of the calves in each of the respective mating groups. When bulls were mated in groups of 8-24, the maximum percent of calves sired by individual bulls was 26+/-7% (mean+/-S.D.) with a range 11-36%. However, when bulls were mated in groups of 2-7, the maximum percent of calves sired by individual bulls was 59+/-19% with a range 24-94%. Calf output of bulls was moderately repeatable across years at four of five sites. Multiple regression models relating pre-mating measures of physical, seminal and behavioural traits to calf output were developed for the three breed groups. In all, only 138 of the 235 bulls were included in the models (40 Santa Gertrudis, 24 5/8 Brahman and 74 Brahman). Sheath and testicular traits, such as scrotal circumference and testicular tone, were generally not related to calf output, the exceptions being sheath depth in Brahman bulls which was negatively related (P<0.05) and scrotal circumference in 5/8 Brahmans which was positively related (P=0.08) to calf output. Dominance was only included in the 5/8 Brahman model but there was no significant relationship between dominance hierarchy and calf output. Semen motility was only related (P<0.05) to calf output in 5/8 Brahmans. However, measures of semen quality based on spermatozoa morphology were important contributors to calf output in the Santa Gertrudis and Brahman models where percent normal spermatozoa was positively related (P<0.01) to calf output. In Santa Gertrudis and Brahman bulls, measures of sexual behaviour in the serving capacity test were related to calf output. In Santa Gertrudis, these were for the number of displays of sexual interest (P<0.05), and mounts (P<0.01), but not number of serves, whilst in Brahman bulls, libido score was positively related to calf output (P<0.05). The models only explained 35-57% of the variation in calf output.     

Bull selection and use in northern Australia. 1. Physical traits

A systematic breeding soundness examination, including assessment of serving capacity was conducted on 5/8 Brahman, Brahman and Belmont Red bulls, most aged 2-4 years, at 12 different cattle properties across northern Australia. A subset of bulls (n=235) were subsequently mated in various groups, as multiple-sires, to cows and heifers at bull:female ratios of 2.5-6%. The number of calves sired by individual bulls (calf output) was determined by DNA typing for paternity. Overall, the incidence of physical abnormalities that were judged likely to result in reduced reproductive performance was low (5.6-12% of bulls). Measures of key physical traits scrotal circumference, sheath depth (vertical distance from ventral abdominal wall to preputial orifice), and sheath score were moderately to highly repeatable. Except in Belmont Red bulls, liveweight was positively correlated with scrotal circumference (r=0.36-0.78; P<0.01), and both traits increased with age. In 2-year-old Brahman bulls only, percent normal spermatozoa was correlated with scrotal circumference (r=0.34; P<0.05). However, 12 and 15% of Santa Gertrudis and Brahman bulls, respectively, with greater than or equal to the recommended threshold values for scrotal circumferences, had less than 50% morphologically normal spermatozoa. Size and conformation of the umbilicus was associated with conformation of the sheath, and influenced mating ability in 2-year-old Brahman and 3-year-old Santa Gertrudis bulls. In 2-year-old Brahman bulls umbilical cord thickness was positively related (r=0.36; P<0.05) to sheath depth, and negatively related (r=-0.65; P<0.05) to number of mounts and serves in a serving capacity test. In 3-year-old Santa Gertrudis bulls navel (external umbilical scar) score was negatively related to the number of serves (r=-0.53; P<0.01) in a serving capacity test. None of the physical traits measured were consistently included in the final regression model for calf output for each breed.     

Aspects of the sexual development of Brahman versus Angus bulls in Florida

Postweaning growth and reproductive traits were studied in 10 Brahman and 12 Angus bulls from 8 through 20 months of age. Brahman bulls reached puberty at 15.9 +/- .4 months of age, weighed 432 +/- 16 kg, had a scrotal circumference (SC) of 33.4 +/- 1.2 cm, and plasma testosterone of 3.96 +/- 1.03 ng/ml. Breed differences in SC averaged over the entire study were not significant. However, the breed x day interaction (BxD) (P<.01) showed that, initially, the Brahman SC was smaller than the Angus SC; however, by the end of the study, the Brahman SC was larger than the Angus. When SC was adjusted for body weight, breed differences (P<.01) and BxD (P <.01) for SC/body weight (BW) reflected the later age and heavier weight at which the Brahman bull reached puberty. Plasma testosterone differed between breeds (Angus > Brahman, P< .01) and increased at a linear (P< .01) rate with age. There was no BxD in plasma testosterone. No breed differences in sperm concentration were observed. However, other semen traits were different (P< .01), i.e., rate of forward movement, sperm motility, total abnormalities and semen volume. A BxD (P< .01) was also evident for breed differences in these semen traits. Sexual development of the Brahman bull occurred at a later chronological age and in a nonparallel pattern to that of the Angus. Between animal variation in SC within the Brahmans and differences between this study and other reports suggest that differences in SC exist for various populations of Brahman bulls and should provide opportunities for progress in selection for this trait.     

Dynamics of genomic architecture during composite breed development in cattle

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.     

Bull selection and use in northern Australia. 3. Serving capacity tests

Serving capacity studies were conducted on Santa Gertrudis, 5/8 Brahman, Brahman and Belmont Red bulls at 10 sites. Modifications, such as providing prior sexual experience, using females in oestrus and comparing restrained with unrestrained females, were made to the standard Bos taurus serving capacity test in a series of experiments with Santa Gertrudis bulls. Providing sexual experience to virgin bulls prior to testing improved the number of serves and libido scores in a serving capacity test. In serving capacity tests that used restrained females, expressions of sexual behaviour of bulls were not improved by having females in oestrus compared with non-oestrus females. There were more sexual behaviours, particularly mounts, in serving capacity tests involving restrained rather than unrestrained females. Summary statistics were collated from over 1100 serving capacity tests, in the 2 months prior to mating, across the 10 sites. In serving capacity tests using restrained females, in bulls of the same age, there were more expressions of sexual behaviours (except interest) in Belmont Red than in Brahman or 5/8 Brahman bulls with Santa Gertrudis bulls intermediate. The percentage of bulls displaying serves tended to increase with age; 82, 83 and 86% of 2-, 3- and > or =4-year-old Belmont Red bulls, respectively, completed 1 or more serves. This compared with 50, 50 and 66% of 2-, 3- and > or =4-year-old Santa Gertrudis bulls and 31% of 2-year-old Brahman bulls. Between 2 and 9% of bulls did not display any sexual behaviour in serving capacity tests with restrained females. Across genotypes, there were few and inconsistent relationships between either mounts or serves and liveweight, condition score, scrotal circumference, testicular tone and sheath measurements. The exceptions were in 2-year-old Brahmans where mounts were negatively correlated (r=-0.61; P<0.05) with umbilicus thickness and in 3-year-old Santa Gertrudis bulls, where serves (r=-0.53; P<0.01) and mounts+serves (r=-0.49; P<0.01) were negatively correlated with navel score. Measures of sexual behaviour were unable to consistently predict calf output of bulls in multiple-sire mating groups in extensive pastoral conditions. The main value of the serving capacity test in these genotypes is to identify whether a bull is capable of natural service.     

Genome‐wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds

Uganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium‐density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (H_O) and expected (H_E) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub‐structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f₄ statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within‐breed diversity and heterozygote advantage in crossbreeding schemes.     

A comparison between cervical dimensions of pregnant and nonpregnant Santa Gertrudis and Bos taurus cows

Cervical measurements were recorded on palpation per rectum for 603 Bos taurus , Santa Gertrudis, and crossbred cattle. Postmortem cervical measurements were recorded for 77 of these cows. Antemortem measurements were positively correlated with those taken at postmortem examination (cervical diameter: r=0.75, P<0.01 and cervical length: r=0.29, P<0.06). Cervical diameters from antemortem and postmortem specimens were greater in Santa Gertrudis cattle than in Bos taurus and crossbred cattle (P<0.05). Cervical shape was different between Santa Gertrudis cattle and Bos taurus or crossbred cattle (P<0.01), with Santa Gertrudis cows having a higher incidence of conical shaped cervices than other breeds. Within Santa Gertrudis cattle, the incidence of conical cervical shape increased with age (P<0.05). Conical cervical shape was associated with a lower pregnancy rate in all breeds (P<0.01). Within Santa Gertrudis cattle, conical cervical shape was associated with a lower pregnancy rate in 1- to 4-year-old cows (P<0.05) but not in cows more than 4 years of age (P>0.05).     

Accuracy of genotype imputation in sheep breeds

Although genomic selection offers the prospect of improving the rate of genetic gain in meat, wool and dairy sheep breeding programs, the key constraint is likely to be the cost of genotyping. Potentially, this constraint can be overcome by genotyping selection candidates for a low density (low cost) panel of SNPs with sparse genotype coverage, imputing a much higher density of SNP genotypes using a densely genotyped reference population. These imputed genotypes would then be used with a prediction equation to produce genomic estimated breeding values. In the future, it may also be desirable to impute very dense marker genotypes or even whole genome re‐sequence data from moderate density SNP panels. Such a strategy could lead to an accurate prediction of genomic estimated breeding values across breeds, for example. We used genotypes from 48 640 (50K) SNPs genotyped in four sheep breeds to investigate both the accuracy of imputation of the 50K SNPs from low density SNP panels, as well as prospects for imputing very dense or whole genome re‐sequence data from the 50K SNPs (by leaving out a small number of the 50K SNPs at random). Accuracy of imputation was low if the sparse panel had less than 5000 (5K) markers. Across breeds, it was clear that the accuracy of imputing from sparse marker panels to 50K was higher if the genetic diversity within a breed was lower, such that relationships among animals in that breed were higher. The accuracy of imputation from sparse genotypes to 50K genotypes was higher when the imputation was performed within breed rather than when pooling all the data, despite the fact that the pooled reference set was much larger. For Border Leicesters, Poll Dorsets and White Suffolks, 5K sparse genotypes were sufficient to impute 50K with 80% accuracy. For Merinos, the accuracy of imputing 50K from 5K was lower at 71%, despite a large number of animals with full genotypes (2215) being used as a reference. For all breeds, the relationship of individuals to the reference explained up to 64% of the variation in accuracy of imputation, demonstrating that accuracy of imputation can be increased if sires and other ancestors of the individuals to be imputed are included in the reference population. The accuracy of imputation could also be increased if pedigree information was available and was used in tracking inheritance of large chromosome segments within families. In our study, we only considered methods of imputation based on population‐wide linkage disequilibrium (largely because the pedigree for some of the populations was incomplete). Finally, in the scenarios designed to mimic imputation of high density or whole genome re‐sequence data from the 50K panel, the accuracy of imputation was much higher (86–96%). This is promising, suggesting that in silico genome re‐sequencing is possible in sheep if a suitable pool of key ancestors is sequenced for each breed.     

A genomewide single-nucleotide-polymorphism panel with high ancestry information for African American admixture mapping

Admixture mapping requires a genomewide panel of relatively evenly spaced markers that can distinguish the ancestral origins of chromosomal segments in admixed individuals. Through use of the results of the International HapMap Project and specific selection criteria, the current study has examined the ability of selected single-nucleotide polymorphisms (SNPs) to extract continental ancestry information in African American subjects and to explore parameters for admixture mapping. Genotyping of two linguistically diverse West African populations (Bini and Kanuri Nigerians, who are Niger-Congo [Bantu] and Nilo-Saharan speakers, respectively), European Americans, and African Americans validated a genomewide set of >4,000 SNP ancestry-informative markers with mean and median F(ST) values >0.59 and mean and median Fisher's information content >2.5. This set of SNPs extracted a larger amount of ancestry information in African Americans than previously reported SNP panels and provides nearly uniform coverage of the genome. Moreover, in the current study, simulations show that this more informative panel improves power for admixture mapping in African Americans when ethnicity risk ratios are modest. This is particularly important in the application of admixture mapping in complex genetic diseases for which only modest ethnicity risk ratios of relevant susceptibility genes are expected.     

Selection and use of SNP markers for animal identification and paternity analysis in U.S. beef cattle

DNA marker technology represents a promising means for determining the genetic identity and kinship of an animal. Compared with other types of DNA markers, single nucleotide polymorphisms (SNPs) are attractive because they are abundant, genetically stable, and amenable to high-throughput automated analysis. In cattle, the challenge has been to identify a minimal set of SNPs with sufficient power for use in a variety of popular breeds and crossbred populations. This report describes a set of 32 highly informative SNP markers distributed among 18 autosomes and both sex chromosomes. Informativity of these SNPs in U.S. beef cattle populations was estimated from the distribution of allele and genotype frequencies in two panels: one consisting of 96 purebred sires representing 17 popular breeds, and another with 154 purebred American Angus from six herds in four Midwestern states. Based on frequency data from these panels, the estimated probability that two randomly selected, unrelated individuals will possess identical genotypes for all 32 loci was 2.0 × 10⁻¹³ for multi-breed composite populations and 1.9 × 10⁻¹⁰ for purebred Angus populations. The probability that a randomly chosen candidate sire will be excluded from paternity was estimated to be 99.9% and 99.4% for the same respective populations. The DNA immediately surrounding the 32 target SNPs was sequenced in the 96 sires of the multi-breed panel and found to contain an additional 183 polymorphic sites. Knowledge of these additional sites, together with the 32 target SNPs, allows the design of robust, accurate genotype assays on a variety of high-throughput SNP genotyping platforms.