Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes.

BACKGROUND: Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected. RESULTS: A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3. CONCLUSIONS: The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Replacement of microsatellite markers by imputed medium-density SNP arrays for parentage control in German warmblood horses.

In horses, parentage control is currently performed based on an internationally standardized panel of 17 microsatellite (MS) markers comprising 12 mandatory and five optional markers. Unlike MS, single nucleotide polymorphism (SNP) profiles support a wider portfolio of genomic applications, including parentage control. A transition to SNP-based parentage control is favorable, but requires additional efforts for ensuring generation-overlapping availability of marker genotypes of the same type. To avoid double genotyping of either parents or offspring for changing to SNP technology and enable efficient transition, we tested whether MS genotypes used for parentage control could be reliably imputed from a medium-density SNP panel in German warmblood horses. Imputation accuracy was tested in a tenfold cross-validation with two approaches: within breed (option A) and across breeds (option B). Average imputation accuracies of 97.98% (A) and 96.17% (B) were achieved, respectively. Due to interbreed differences in genotyping rates, five MS markers of low genotyping rate (GTR; < 90%) could be imputed with higher accuracy within breed (98.18%) than across breeds (90.73%). MS markers with high GTR performed homogeneously well in option B (98.44%) and showed slightly lower accuracy in option A (97.90%). Among these markers, AHT5 proved to be problematic for imputation regardless of the approach, revealing accuracies of 86.40% (A) and 88.70% (B). Better results for MS markers with high GTR and savings in computational processing justified the choice of option B for routine implementation. To date, more than 9500 horses have undergone the new parentage control based on imputed MS genotypes.

Genetic and genomic characterization followed by single-step genomic evaluation of withers height in German Warmblood horses.

Reliability of genomic predictions is influenced by the size and genetic composition of the reference population. For German Warmblood horses, compilation of a reference population has been enabled through the cooperation of five German breeding associations. In this study, preliminary data from this joint reference population were used to genetically and genomically characterize withers height and to apply single-step methodology for estimating genomic breeding values for withers height. Using data on 2113 mares and their genomic information considering about 62,000 single nucleotide polymorphisms (SNPs), analysis of the genomic relationship revealed substructures reflecting breed origin and different breeding goals of the contributing breeding associations. A genome-wide association study confirmed a known quantitative trait locus (QTL) for withers height on equine chromosome (ECA) 3 close to LCORL and identified a further significant peak on ECA 1. Using a single-step approach with a combined relationship matrix, the estimated heritability for withers height was 0.31 (SE = 0.08) and the corresponding genomic breeding values ranged from - 2.94 to 2.96 cm. A mean reliability of 0.38 was realized for these breeding values. The analyses of withers height showed that compiling a reference population across breeds is a suitable strategy for German Warmblood horses. The single-step method is an appealing approach for practical genomic prediction in horses, because not many genotypes are available yet and animals without genotypes can by this way directly contribute to the estimation system.

Quantitative trait loci for meat quality traits in pigs considering imprinting and epistatic effects.

The aim of the research was to gain a better understanding of the genomic regulation of meat quality by investigating individual and epistatic QTL in a three-generation full-sib population (Pietrain x crossbred dam line). In total, 386 animals were genotyped for 96 markers. Analysed traits included pH, reflectance value, conductivity, and meat colour. Thirteen significant individual QTL were identified. The most significant QTL were detected on SSC1 and SSC9 for pH, on SSC4 for meat colour, and on SSC8 for conductivity, accounting for 3.4% to 4.7% of the phenotypic variance. Nine significant epistatic QTL pairs were detected accounting for between 5.7% and 10.9% of the phenotypic variance. Epistatic QTL pairs showing the largest effects were for reflectance value between two locations of SSC4, and for pH between SSC10 and SSC13, explaining 9.5% and 10.9% of the phenotypic variance, respectively. This study indicates that meat quality traits are influenced by numerous QTL as well as a complex network of interactions.

Genomic scan for quantitative trait loci of chemical and physical body composition and deposition on pig chromosome X including the pseudoautosomal region of males.

A QTL analysis of pig chromosome X (SSCX) was carried out using an approach that accurately takes into account the specific features of sex chromosomes i.e. their heterogeneity, the presence of a pseudoautosomal region and the dosage compensation phenomenon. A three-generation full-sib population of 386 animals was created by crossing Pietrain sires with a crossbred dam line. Phenotypic data on 72 traits were recorded for at least 292 and up to 315 F2 animals including chemical body composition measured on live animals at five target weights ranging from 30 to 140 kg, daily gain and feed intake measured throughout growth, and carcass characteristics obtained at slaughter weight (140 kg). Several significant and suggestive QTL were detected on pig chromosome X: (1) in the pseudoautosomal region of SSCX, a QTL for entire loin weight, which showed paternal imprinting, (2) closely linked to marker SW2456, a suggestive QTL for feed intake at which Pietrain alleles were found to be associated with higher feed intake, which is unexpected for a breed known for its low feed intake capacity, (3) at the telomeric end of the q arm of SSCX, QTL for jowl weight and lipid accretion and (4) suggestive QTL for chemical body composition at 30 kg. These results indicate that SSCX is important for physical and chemical body composition and accretion as well as feed intake regulation.

Antimicrobial activity of bovine psoriasin.

Human psoriasin (S100A7) has originally been described as a member of the family of S100 calcium-binding proteins which is overexpressed in patients suffering from psoriasis. The bovine homolog was first identified as a cow-derived respiratory allergen. As Escherichia coli mastitis is a common problem in dairy cattle, and human psoriasin was found to exhibit antimicrobial activity preferentially against E. coli, we examined whether the bovine mRNA is expressed in the mammary gland. To demonstrate the antimicrobial activity of bovine psoriasin, we isolated cDNA from the udder, cloned the bovine psoriasin gene in a bacterial expression vector, and the recombinant protein was expressed in BL21 cells. The in vitro antibacterial activity was tested by performing microdilution susceptibility tests and radial diffusion assays with eight different bacterial strains, thereof three different E. coli strains, and one yeast. The antimicrobial activity of the recombinant bovine psoriasin is comparable with human psoriasin and also limited to E. coli. Psoriasin appears to be a part of the local host defense mechanism in the udder, is a putative candidate for a cow-specific factor influencing mastitis susceptibility, and a possible alternative to conventional antibiotics.

A novel horse alpha-defensin: gene transcription, recombinant expression and characterization of the structure and function.

Defensins are a predominant class of antimicrobial peptides, which act as endogenous antibiotics. Defensins are classified into three distinct sub-families: theta-, beta-, and alpha-defensins. Synthesis of alpha-defensin has been confirmed only in primates and glires to date and is presumably unique for a few tissues, including neutrophils and Paneth cells of the small intestine. Antimicrobial activities of these peptides were shown against a wide variety of microbes including bacteria, fungi, viruses and protozoan parasites. In the present study, we report the characterization of the equine alpha-defensin DEFA (defensin alpha) 1. Transcription analysis revealed that the transcript of the gene is present in the small intestine only. An alignment with known alpha-defensins from primates and glires displayed a homology with Paneth-cell-specific alpha-defensins. DEFA1 was recombinantly expressed in Escherichia coli and subsequently analysed structurally by CD and molecular modelling. To examine the antimicrobial properties, a radial diffusion assay was performed with 12 different micro-organisms and the LD90 (lethal dose killing > or =90% of target organism) and MBC (minimal bactericidal concentration) values were examined. DEFA1 showed an antimicrobial activity against different Gram-positive and Gram-negative bacteria and against the yeast Candida albicans. Using viable bacteria in combination with a membrane-impermeable fluorescent dye, as well as depolarization of liposomes as a minimalistic system, it became evident that membrane permeabilization is at least an essential part of the peptide's mode of action.

Estimation of breed contributions to present and future genetic diversity of 44 North Eurasian cattle breeds using core set diversity measures.

Extinction of breeds threatens genetic diversity of livestock species. The need to conserve genetic diversity is widely accepted but involves in general two questions: (i) is the expected loss of diversity in a set of breeds within a defined future time horizon large enough to establish a conservation plan, and if so (ii) which breeds should be prioritised for such a conservation plan? The present study uses a marker assisted methodology to address these questions. The methodology combines core set diversity measures with a stochastic method for the estimation of expected future diversity and breed marginal diversities. The latter is defined as the change in the total diversity of all breeds caused by a one unit decrease in extinction probability of a particular breed. The stochastic method was validated by means of simulations. A large field data set consisting of 44 North Eurasian cattle breeds was analysed using simplified determined extinction probabilities. The results show that the expected loss of diversity in this set within the next 20 to 50 years is between 1 and 3% of the actual diversity, provided that the extinction probabilities which were used are approximately valid. If this loss is to be reduced, it is sufficient to include those three to five breeds with the highest marginal diversity in a conservation scheme.

Bovine beta-defensins: identification and characterization of novel bovine beta-defensin genes and their expression in mammary gland tissue.

beta-Defensin genes code for multifunctional peptides with a broad-range antimicrobial activity. In this project we hypothesized that beta-defensin genes may be candidate genes for resistance to mastitis. In this article we describe the identification and genomic characterization of eight bovine beta-defensin genes, including six novel defensin genes and two pseudogenes. Expression in the bovine mammary gland of one of the novel genes, DEFB401, has been demonstrated, as well as the expression of LAP, TAP, DEFB1, BNBD3, BNBD9, and BNBD12. For genomic characterization, 20 BACs from two different bovine BAC libraries (RZPD numbers 750 and 754) were isolated by PCR screening with beta-defensin consensus primers derived from published sequences. PCR products from BACs generated with consensus primers have been subcloned and sequenced, revealing a total of 16 genes and two pseudogenes. Six novel beta-defensin genes share the typical exon-intron structure and are highly homologous to published bovine beta-defensin genes. They are named DEFB401- DEFB405 and LAP-like, and two novel pseudogenes are named EBD-P and EBD-P2. Analysis of mammary gland tissue-derived cDNA from nine cows with different clinical findings demonstrated the expression of several beta-defensin genes mentioned above. First results indicate that the lactational status of the cow presumably has no influence on gene expression. Competent knowledge of antimicrobial activity of beta-defensins from literature, the abundance of beta-defensin mRNA in the bovine mammary gland, and the inducibility of some genes give first evidence that beta-defensins may play a role in local host defense during udder infections.

Multiple quantitative trait loci mapping with cofactors and application of alternative variants of the false discovery rate in an enlarged granddaughter design.

The experimental power of a granddaughter design to detect quantitative trait loci (QTL) in dairy cattle is often limited by the availability of progeny-tested sires, by the ignoring of already identified QTL in the statistical analysis, and by the application of stringent experimentwise significance levels. This study describes an experiment that addressed these points. A large granddaughter design was set up that included sires from two countries (Germany and France), resulting in almost 2000 sires. The animals were genotyped for markers on nine different chromosomes. The QTL analysis was done for six traits separately using a multimarker regression that included putative QTL on other chromosomes as cofactors in the model. Different variants of the false discovery rate (FDR) were applied. Two of them accounted for the proportion of truly null hypotheses, which were estimated to be 0.28 and 0.3, respectively, and were therefore tailored to the experiment. A total of 25 QTL could be mapped when cofactors were included in the model-7 more than without cofactors. Controlling the FDR at 0.05 revealed 31 QTL for the two FDR methods that accounted for the proportion of truly null hypotheses. The relatively high power of this study can be attributed to the size of the experiment, to the QTL analysis with cofactors, and to the application of an appropriate FDR.

[Amplification of pig microsatellite markers using multiplex PCR].

In order to rapidly amplify pig microsatellite markers and save materials,multiplex PCR was used and its reaction condition was optimized. Forty-six combinations of multiplex PCR with good effects were obtained. Thirty of them are duplex-PCRs, sixteen are triplex-PCRs. The results of multiplexes showed that the concentration of primers varied among 0.06 approximately 0.3 micromol/L, the Mg(2+) concentration among 1.5 approximately 3.0 mmol/L; 0.2 approximately 0.4 U of Taq polymerase and 1.0-, 1.2-, 1.4-, 1.6-fold buffer were used, the annealing temperature and the cycle number varied among 52 approximately 60 degrees and 32 approximately 50 degrees, respectively. All multiplexes were further combined into 17 sets for the electrophoresis on ABI 377 sequencer.

[Comparison between the female- and male-linkage maps in pigs].

The difference between the length of female- and male-linkage map, which was created with a reference pedigree based on a commercial porcine population and using 163 microsatellite markers as well as 3 type-I markers (RYR1, PRKAG3, PIT1), was statistic analyzed. The results showed that the total length of female linkage map of autosomes is 2625.9 cm and the total length of the male linkage map is 2259.7 cm; the ratio between the total length of the female- and male-linkage maps is 1.16 : 1; except for the chromosomes 1 and 14, the female linkage maps of the other chromosomes are longer than the male linkage maps. The difference between the length of female- and male-linkage maps of chromosomes 1, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 16, 17 and 18 is very significant (P<0.01) and the difference of chromosome 9 is significant (P<0.05); but there is no significance on chromosomes 2, 4, 12 and 15.

[A porcine linkage map of microsatellite markers using a commercial population as a reference].

A commercial population created with 19 hybrid boars [Piétrain x (Piétrain x Hampshire)], 52 hybrid sows [Leicoma x (Large White x Landrace)] and their 332 offspring was used to construct a whole genome porcine linkage map. The genetic markers used in this study include 172 microsatellite markers and 3 type- I markers (RYR1, PIT1 and PRKAG3). All microsatellite markers were genotyped using multiplex PCR reactions and visualized on ABI 377 sequencer. Three type I markers were assayed using PCR-RFLP technique. CRIMAP (2.4) analysis revealed a total length of the sex-averaged map for SSC1-SSC18 by 2449.2 cM and the length of SSCX by 143.1 cM. The average interval distance accounts to 16.3 cM. The heterozygosity of the parents at the microsatellite loci averaged 0.70. This map will play an important role in screening of QTL for growth, carcass and meat quality and reproduction in commercial populations.

Combined analysis of data from two granddaughter designs: A simple strategy for QTL confirmation and increasing experimental power in dairy cattle.

A joint analysis of five paternal half-sib Holstein families that were part of two different granddaughter designs (ADR- or Inra-design) was carried out for five milk production traits and somatic cell score in order to conduct a QTL confirmation study and to increase the experimental power. Data were exchanged in a coded and standardised form. The combined data set (JOINT-design) consisted of on average 231 sires per grandsire. Genetic maps were calculated for 133 markers distributed over nine chromosomes. QTL analyses were performed separately for each design and each trait. The results revealed QTL for milk production on chromosome 14, for milk yield on chromosome 5, and for fat content on chromosome 19 in both the ADR- and the Inra-design (confirmed within this study). Some QTL could only be mapped in either the ADR- or in the Inra-design (not confirmed within this study). Additional QTL previously undetected in the single designs were mapped in the JOINT-design for fat yield (chromosome 19 and 26), protein yield (chromosome 26), protein content (chromosome 5), and somatic cell score (chromosome 2 and 19) with genomewide significance. This study demonstrated the potential benefits of a combined analysis of data from different granddaughter designs.

Improved confidence intervals in quantitative trait loci mapping by permutation bootstrapping.

The nonparametric bootstrap approach is known to be suitable for calculating central confidence intervals for the locations of quantitative trait loci (QTL). However, the distribution of the bootstrap QTL position estimates along the chromosome is peaked at the positions of the markers and is not tailed equally. This results in conservativeness and large width of the confidence intervals. In this study three modified methods are proposed to calculate nonparametric bootstrap confidence intervals for QTL locations, which compute noncentral confidence intervals (uncorrected method I), correct for the impact of the markers (weighted method I), or both (weighted method II). Noncentral confidence intervals were computed with an analog of the highest posterior density method. The correction for the markers is based on the distribution of QTL estimates along the chromosome when the QTL is not linked with any marker, and it can be obtained with a permutation approach. In a simulation study the three methods were compared with the original bootstrap method. The results showed that it is useful, first, to compute noncentral confidence intervals and, second, to correct the bootstrap distribution of the QTL estimates for the impact of the markers. The weighted method II, combining these two properties, produced the shortest and less biased confidence intervals in a large number of simulated configurations.