Epistatic interactions on chromosome 14 influencing stillbirth in Fleckvieh cattle

Single nucleotide polymorphism (SNP) data of 7384 Fleckvieh bulls was analyzed to identify epistatic interactions influencing stillbirth. Deregressed breeding values were used as phenotypes. The epistatic effects were identified as significant interaction terms from pairwise linear regressions performed for each SNP after accounting for multiple testing. Majority of the detected epistatic effects were located in the 9-31Mb region of chromosome 14, corresponding to the most significant region from the genome wide association. Additional epistatic SNPs at 50.5Mb and 80.5 Mb at the same chromosome were detected. The region around 25 Mb contained genes connected to height and body size such as PLAG1, CHCHD7, LYN, RDHE2 (SDR16C5) and PENK. The other interesting region at 50.5Mb contained the TRPS1 gene influencing bone malformations. Both regions have been identified as candidates influencing stillbirth

Comparison of gluteus medius muscle activity in Haflinger and Noriker horses with polysaccharide storage myopathy

Type 1 polysaccharide storage myopathy caused by genetic mutation in the glycogen synthase 1 gene is present in many breeds including the Noriker and Haflinger horses. In humans, EMG has already been used to document changes in the muscle activity patterns of patients affected by human glycogen storage disorders. Therefore, the aim of the present study was to describe gluteus muscle activity with surface electromyography (sEMG) in Haflinger and Noriker horses with known GYS1 mutation status during walk and trot. Thirty‐two horses (11 Haflinger and 21 Noriker horses) with homozygous non‐affected (GG), heterozygous affected (GA) and homozygous affected (AA) status of GYS1 mutation without overt clinical signs of any myopathy were selected for the current study. Using surface electromyography gluteus medius muscle activity at walk and at trot was measured, and muscle activity was described in relation to the maximum observed value at the same sensor and the same gait. In order to further describe the signals in detail comprising both frequencies and amplitudes, the crossings through the baseline and the 25, 50 and 75 percentile lines were determined. The result of the relative muscle activity did not show a consistent difference between affected and non‐affected horses. Genetically affected (GA and AA) horses showed significantly less density of muscle activity for both gaits and horse breeds except for the crossings per second at the baseline and 75 percentile at walk in the Haflinger horses and 75 percentile at trot in the Noriker horses. The medians of all calculated density values were significantly lower in the GA Haflingers compared to the GG Haflingers (p = 0.012) and also in the AA Norikers compared to the GG Norikers (p = 0.011). Results indicate that the GYS1 mutation reduces the number of functional muscle fibres detected by sEMG measurements even in the absence of overt clinical signs

Recovery of Native Genetic Background in Admixed Populations Using Haplotypes, Phenotypes, and Pedigree Information - Using Cika Cattle as a Case Breed

The aim of this study was to obtain unbiased estimates of the diversity parameters, the population history, and the degree of admixture in Cika cattle which represents the local admixed breeds at risk of extinction undergoing challenging conservation programs. Genetic analyses were performed on the genome-wide Single Nucleotide Polymorphism (SNP) Illumina Bovine SNP50 array data of 76 Cika animals and 531 animals from 14 reference populations. To obtain unbiased estimates we used short haplotypes spanning four markers instead of single SNPs to avoid an ascertainment bias of the BovineSNP50 array. Genome-wide haplotypes combined with partial pedigree and type trait classification show the potential to improve identification of purebred animals with a low degree of admixture. Phylogenetic analyses demonstrated unique genetic identity of Cika animals. Genetic distance matrix presented by rooted Neighbour-Net suggested long and broad phylogenetic connection between Cika and Pinzgauer. Unsupervised clustering performed by the admixture analysis and two-dimensional presentation of the genetic distances between individuals also suggest Cika is a distinct breed despite being similar in appearance to Pinzgauer. Animals identified as the most purebred could be used as a nucleus for a recovery of the native genetic background in the current admixed population. The results show that local well-adapted strains, which have never been intensively managed and differentiated into specific breeds, exhibit large haplotype diversity. They suggest a conservation and recovery approach that does not rely exclusively on the search for the original native genetic background but rather on the identification and removal of common introgressed haplotypes would be more powerful. Successful implementation of such an approach should be based on combining phenotype, pedigree, and genome-wide haplotype data of the breed of interest and a spectrum of reference breeds which potentially have had direct or indirect historical contribution to the genetic makeup of the breed of interest

Population Structure and Genetic Diversity of Sheep Breeds in the Kyrgyzstan

Sheep are a main livestock species of Kyrgyzstan, a Central Asian country with predominating mountain terrain. The current gene pool of local sheep resources has been forming under diverse climate conditions from the era of the trading caravans of the Great Silk Road, through the Soviet period of large-scale livestock improvements, which was followed by the deep crisis at the end of the 20th century, up to now. However, not much is known about the genetic background and variability of the local sheep populations. Therefore, our aims were to provide a characterization of the population structure and genetic relations within the Kyrgyz sheep breeds and to study their genetic connections with the global sheep breeds using SNP analysis. Samples of the Alai (n = 31), Gissar (n = 30), Kyrgyz coarse wool (n = 13), Aykol (n = 31), and Tien-Shan (n = 24) breeds were genotyped with the OvineSNP50 BeadChip or the Ovine Infinium HD BeadChip (Illumina Inc., USA). The measure of inbreeding based on runs of homozygosity showed a minimum value in the Aykol breed (FROH = 0.034), while the maximum was found in the Alai breed (FROH = 0.071). Short ROH segments (ROH ≤ 4 Mb) were predominant in all breeds. Long ROH segments (ROH > 16 Mb) were absent in the Gissar breed. The Gissar and Aykol breeds had the highest values of the effective population sizes estimated for five generations ago (Ne5 = 660 and 563), whereas the Alai and Kyrgyz coarse wool displayed lower values (Ne5 = 176 and 128, respectively). The synthetic origin of the Aykol breed was clearly evidenced by all analyses applied. Based on the network and admixture analyses of the Kyrgyz and global sheep breeds, the Tien-Shan and the Russian semi-fine wool breeds demonstrated a common ancestry that most likely is due to a contribution of the Lincoln breed. The Gissar, Aykol, and Kyrgyz coarse wool breeds showed a genetic background predominating in sheep populations from Iran and China whereas the Alai demonstrated the different ancestry type. The revealed admixture patterns probably resulted from the exchange and trade during the era of the Great Silk Road, which partly overlapped with historical and archeological findings

Using the community-based breeding program (CBBP) model as a collaborative platform to develop the African Goat Improvement Network—Image collection protocol (AGIN-ICP) with mobile technology for data collection and management of livestock phenotypes

Introduction: The African Goat Improvement Network Image Collection Protocol (AGIN-ICP) is an accessible, easy to use, low-cost procedure to collect phenotypic data via digital images. The AGIN-ICP collects images to extract several phenotype measures including health status indicators (anemia status, age, and weight), body measurements, shapes, and coat color and pattern, from digital images taken with standard digital cameras or mobile devices. This strategy is to quickly survey, record, assess, analyze, and store these data for use in a wide variety of production and sampling conditions.Methods: The work was accomplished as part of the multinational African Goat Improvement Network (AGIN) collaborative and is presented here as a case study in the AGIN collaboration model and working directly with community-based breeding programs (CBBP). It was iteratively developed and tested over 3 years, in 12 countries with over 12,000 images taken.Results and discussion: The AGIN-ICP development is described, and field implementation and the quality of the resulting images for use in image analysis and phenotypic data extraction are iteratively assessed. Digital body measures were validated using the PreciseEdge Image Segmentation Algorithm (PE-ISA) and software showing strong manual to digital body measure Pearson correlation coefficients of height, length, and girth measures (0.931, 0.943, 0.893) respectively. It is critical to note that while none of the very detailed tasks in the AGIN-ICP described here is difficult, every single one of them is even easier to accidentally omit, and the impact of such a mistake could render a sample image, a sampling day’s images, or even an entire sampling trip’s images difficult or unusable for extracting digital phenotypes. Coupled with tissue sampling and genomic testing, it may be useful in the effort to identify and conserve important animal genetic resources and in CBBP genetic improvement programs by providing reliably measured phenotypes with modest cost. Potential users include farmers, animal husbandry officials, veterinarians, regional government or other public health officials, researchers, and others. Based on these results, a final AGIN-ICP is presented, optimizing the costs, ease, and speed of field implementation of the collection method without compromising the quality of the image data collection

Genetic association between somatic cell score and milk lactose in early- to mid-lactation of first calving Fleckvieh cows

The present study aimed to assess genetic correlations of milk lactose (percentage and yield), lactose to fat ratio and lactose to protein ratio with somatic cell score (SCS) in the first half of lactation, the period with the highest risk of mastitis in cows. More than 35,000 first lactation records from Austrian Fleckvieh breed were available for this purpose. A linear animal model was adopted for the genetic analysis, with inclusion of age at calving and year-month of calving as fixed effects, and additive genetic animal, herd-year of calving and residual as random effects. The average SCS in the first 150 days in milk was low (1.53 \ub11.29) compared with values reported in literature for other breeds in first lactation. Heritabilities of lactose percentage, lactose to fat ratio and lactose to protein ratio were high, ranging between 0.65 and 0.71. According to lactose yield and SCS, their heritabilities were equal to 0.3 and 0.2, respectively. Since the focus was on the first lactation, where the heritability of lactose is reported to be the highest, these results agreed with expectations. The strongest genetic correlation was found between SCS and lactose percentage, and the weakest between SCS and lactose yield, suggesting that lactose percentage and yield have different variability at animal level and different potential roles at genetic level. Both lactose to fat ratio and lactose to protein ratio showed weak and negative genetic correlation with SCS. Further research is required to assess relations between lactose and recorded mastitis data, with the final goal of validating lactose level as udder inflammation indicator and explore its potential role in breeding programs to reduce cow\u2019s susceptibility to mastitis

Copy number expansion of the STX17 duplication in melanoma tissue from Grey horses

Background: Greying with age in horses is an autosomal dominant trait, associated with loss of hair pigmentation, melanoma and vitiligo-like depigmentation. We recently identified a 4.6 kb duplication in STX17 to be associated with the phenotype. The aims of this study were to investigate if the duplication in Grey horses shows copy number variation and to exclude that any other polymorphism is uniquely associated with the Grey mutation. Results: We found little evidence for copy number expansion of the duplicated sequence in blood DNA from Grey horses. In contrast, clear evidence for copy number expansions was indicated in five out of eight tested melanoma tissues or melanoma cell lines. A tendency of a higher copy number in aggressive tumours was also found. Massively parallel resequencing of the similar to 350 kb Grey haplotype did not reveal any additional mutations perfectly associated with the phenotype, confirming the duplication as the true causative mutation. We identified three SNP alleles that were present in a subset of Grey haplotypes within the 350 kb region that shows complete linkage disequilibrium with the causative mutation. Thus, these three nucleotide substitutions must have occurred subsequent to the duplication, consistent with our interpretation that the Grey mutation arose more than 2,000 years before present. Conclusions: These results suggest that the mutation acts as a melanoma-driving regulatory element. The elucidation of the mechanistic features of the duplication will be of considerable interest for the characterization of these horse melanomas as well as for the field of human melanoma research

Data from: Conservation of a domestic metapopulation structured into related and partly admixed strains

Preservation of genetic diversity is one of the most pressing challenges in the planetary boundaries concept. Within this context, we focused on genetic diversity in a native, unselected and highly admixed domesticated metapopulation. A set of 1828 individuals from 60 different cattle breeds was analysed using a medium density SNP chip. Among these breeds, 14 Buša strains formed a metapopulation represented by 350 individuals, while the remaining 46 breeds represented the global cattle population. Genetic analyses showed that the scarcely selected and less differentiated Buša metapopulation contributed a substantial proportion (52.6%) of the neutral allelic diversity to this global taurine population. Consequently, there is an urgent need for synchronised maintenance of this highly fragmented domestic metapopulation, which is distributed over several countries without sophisticated infrastructure and highly endangered by continuous replacement crossing as part of the global genetic homogenisation process. This study collected and evaluated samples, data and genome-wide information and developed genome-assisted cross-border conservation concepts. To detect and maintain genetic integrity of the metapopulation strains, we designed and applied a composite test that combines six metrics based on additive genetic relationships, a nearest neighbour graph and the distribution of semi-private alleles. Each metric provides distinct information components about past admixture events and offers an objective and powerful tool for the detection of admixed outliers. The here developed conservation methods and presented experiences could easily be adapted to comparable conservation programmes of domesticated or other metapopulations bred and kept in captivity or under some other sort of human control

Accuracy of genotype imputation in Nelore cattle

Background: Genotype imputation from low-density (LD) to high-density single nucleotide polymorphism (SNP) chips is an important step before applying genomic selection, since denser chips tend to provide more reliable genomic predictions. Imputation methods rely partially on linkage disequilibrium between markers to infer unobserved genotypes. Bos indicus cattle (e.g. Nelore breed) are characterized, in general, by lower levels of linkage disequilibrium between genetic markers at short distances, compared to taurine breeds. Thus, it is important to evaluate the accuracy of imputation to better define which imputation method and chip are most appropriate for genomic applications in indicine breeds.Methods: Accuracy of genotype imputation in Nelore cattle was evaluated using different LD chips, imputation software and sets of animals. Twelve commercial and customized LD chips with densities ranging from 7 K to 75 K were tested. Customized LD chips were virtually designed taking into account minor allele frequency, linkage disequilibrium and distance between markers. Software programs Flmpute and BEAGLE were applied to impute genotypes. From 995 bulls and 1247 cows that were genotyped with the Illumina (R) BovineHD chip (HD), 793 sires composed the reference set, and the remaining 202 younger sires and all the cows composed two separate validation sets for which genotypes were masked except for the SNPs of the LD chip that were to be tested.Results: Imputation accuracy increased with the SNP density of the LD chip. However, the gain in accuracy with LD chips with more than 15 K SNPs was relatively small because accuracy was already high at this density. Commercial and customized LD chips with equivalent densities presented similar results. Flmpute outperformed BEAGLE for all LD chips and validation sets. Regardless of the imputation software used, accuracy tended to increase as the relatedness between imputed and reference animals increased, especially for the 7 K chip.Conclusions: If the Illumina (R) BovineHD is considered as the target chip for genomic applications in the Nelore breed, cost-effectiveness can be improved by genotyping part of the animals with a chip containing around 15 K useful SNPs and imputing their high-density missing genotypes with Flmpute

The genotypes from 55,266 SNP scored in 1,210 cattle

This file contains the genotypes from 55,266 SNP markers scored in 1,210 cattle. Samples were genotyped with the Illumina BovineSNP50 BeadChip version 1 and version 2. Genotypes were called using GenomeStudio