Jens Tetens

The pattern of linkage disequilibrium in German Holstein cattle

This study presents a second generation of linkage disequilibrium (LD) map statistics for the whole genome of the Holstein-Friesian population, which has a four times higher resolution compared with that of the maps available so far. We used DNA samples of 810 German Holstein-Friesian cattle genotyped by the Illumina Bovine SNP50K BeadChip to analyse LD structure. A panel of 40 854 (75.6%) markers was included in the final analysis. The pairwise r(2) statistic of SNPs up to 5 Mb apart across the genome was estimated. A mean value of r(2) = 0.30 +/- 0.32 was observed in pairwise distances of <25 kb and it dropped to 0.20 +/- 0.24 at 50-75 kb, which is nearly the average inter-marker space in this study. The proportion of SNPs in useful LD (r(2) > or = 0.25) was 26% for the distance of 50 and 75 kb between SNPs. We found a lower level of LD for SNP pairs at the distance < or =100 kb than previously thought. Analysis revealed 712 haplo-blocks spanning 4.7% of the genome and containing 8.0% of all SNPs. Mean and median block length were estimated as 164 +/- 117 kb and 144 kb respectively. Allele frequencies of the SNPs have a considerable and systematic impact on the estimate of r(2). It is shown that minimizing the allele frequency difference between SNPs reduces the influence of frequency on r(2) estimates. Analysis of past effective population size based on the direct estimates of recombination rates from SNP data showed a decline in effective population size to N(e) = 103 up to approximately 4 generations ago. Systematic effects of marker density and effective population size on observed LD and haplotype structure are discussed.

A genome‐wide scan for signatures of recent selection in Holstein cattle

The data from the newly available 50 K SNP chip was used for tagging the genome-wide footprints of positive selection in Holstein-Friesian cattle. For this purpose, we employed the recently described Extended Haplotype Homozygosity test, which detects selection by measuring the characteristics of haplotypes within a single population. To assess formally the significance of these results, we compared the combination of frequency and the Relative Extended Haplotype Homozygosity value of each core haplotype with equally frequent haplotypes across the genome. A subset of the putative regions showing the highest significance in the genome-wide EHH tests was mapped. We annotated genes to identify possible influence they have in beneficial traits by using the Gene Ontology database. A panel of genes, including FABP3, CLPN3, SPERT, HTR2A5, ABCE1, BMP4 and PTGER2, was detected, which overlapped with the most extreme P-values. This panel comprises some interesting candidate genes and QTL, representing a broad range of economically important traits such as milk yield and composition, as well as reproductive and behavioural traits. We also report high values of linkage disequilibrium and a slower decay of haplotype homozygosity for some candidate regions harbouring major genes related to dairy quality. The results of this study provide a genome-wide map of selection footprints in the Holstein genome, and can be used to better understand the mechanisms of selection in dairy cattle breeding.

Identification and Dissection of Four Major QTL Affecting Milk Fat Content in the German Holstein-Friesian Population

Milk composition traits exhibit a complex genetic architecture with a small number of major quantitative trait loci (QTL) explaining a large fraction of the genetic variation and numerous QTL with minor effects. In order to identify QTL for milk fat percentage (FP) in the German Holstein-Friesian (HF) population, a genome-wide association study (GWAS) was performed. The study population consisted of 2327 progeny-tested bulls. Genotypes were available for 44,280 SNPs. Phenotypes in the form of estimated breeding values (EBVs) for FP were used as highly heritable traits. A variance components-based approach was used to account for population stratification. The GWAS identified four major QTL regions explaining 46.18% of the FP EBV variance. Besides two previously known FP QTL on BTA14 (P = 8.91×10−198) and BTA20 (P = 7.03×10−12) within DGAT1 and GHR, respectively, we uncovered two additional QTL regions on BTA5 (P = 2.00×10−13) and BTA27 (P = 9.83×10−5) encompassing EPS8 and GPAT4, respectively. EPS8 and GPAT4 are involved in lipid metabolism in mammals. Re-sequencing of EPS8 and GPAT4 revealed 50 polymorphisms. Genotypes for five of them were inferred for the entire study population. Two polymorphisms affecting potential transcription factor binding sites of EPS8 (P = 1.40×10−12) and GPAT4 (P = 5.18×10−5), respectively, were highly significantly associated with the FP EBV. Our results provide evidence that alteration of regulatory sites is an important aspect of genetic variation of complex traits in cattle.

A genome-wide association study indicates LCORL/NCAPG as a candidate locus for withers height in German Warmblood horses

A genome-wide association scan for loci affecting withers height was conducted in 782 German Warmblood stallions, which were genotyped using the Illumina EquineSNP50 Bead Chip. A principal components approach was applied to correct for population structure. The analysis revealed a single major QTL on ECA3 explaining ~18 per cent of the phenotypic variance, which is in concordance with recent reports from other horse populations. The LCORL/NCAPG locus represents a strong candidate gene for this QTL. This locus is among a small number that have consistently been identified to influence human height in several large meta-analyses. Furthermore, a mutation within the NCAPG gene was found to affect growth and body frame size in cattle. Together with the results of this study in German Warmbloods, these findings strongly indicate LCORL/NCAPG as a candidate locus for withers height in horses. Further studies are, however, needed to confirm this.

Exploration of relationships between production and fertility traits in dairy cattle via association studies of SNPs within candidate genes derived by expression profiling

The objective of this work was to integrate findings from functional genomics studies with genome-wide association studies for fertility and production traits in dairy cattle. Association analyses of production and fertility traits with SNPs located within or close to 170 candidate genes derived from two gene expression studies and from the literature were performed. Data from 2294 Holstein bulls genotyped for 39557 SNPs were used. A total of 111 SNPs were located on chromosomal segments covered by a candidate gene. Allele substitution effects for each SNP were estimated using a mixed model with a fixed effect of marker and a random polygenic effect. Assumed covariance was derived either from marker or from pedigree information. Results from the analysis with the kinship matrix built from marker genotypes were more conservative than from the analysis with the pedigree-derived relationship matrix. From sixteen SNPs with significant effects on both classes of traits, ten provided evidence of an antagonistic relationship between productivity and fertility. However, we found four SNPs with favourable effects on fertility and on yield traits, one SNP with favourable effects on fertility and percentage traits, and one SNP with antagonistic effects on two fertility traits. While most quantitative genetic studies have proven genetic antagonisms between yield and functional traits, improvements in both production and functionality may be possible when focusing on a few relevant SNPs. Investigations combining input from quantitative genetics and functional genomics with association analysis may be applied for the identification of such SNPs.

Microphthalmia in Texel sheep is associated with a missense mutation in the paired-like homeodomain 3 (PITX3) gene.

Microphthalmia in sheep is an autosomal recessive inherited congenital anomaly found within the Texel breed. It is characterized by extremely small or absent eyes and affected lambs are absolutely blind. For the first time, we use a genome-wide ovine SNP array for positional cloning of a Mendelian trait in sheep. Genotyping 23 cases and 23 controls using Illuminas OvineSNP50 BeadChip allowed us to localize the causative mutation for microphthalmia to a 2.4 Mb interval on sheep chromosome 22 by association and homozygosity mapping. The PITX3 gene is located within this interval and encodes a homeodomain-containing transcription factor involved in vertebrate lens formation. An abnormal development of the lens vesicle was shown to be the primary event in ovine microphthalmia. Therefore, we considered PITX3 a positional and functional candidate gene. An ovine BAC clone was sequenced, and after full-length cDNA cloning the PITX3 gene was annotated. Here we show that the ovine microphthalmia phenotype is perfectly associated with a missense mutation (c.338G>C, p.R113P) in the evolutionary conserved homeodomain of PITX3. Selection against this candidate causative mutation can now be used to eliminate microphthalmia from Texel sheep in production systems. Furthermore, the identification of a naturally occurring PITX3 mutation offers the opportunity to use the Texel as a genetically characterized large animal model for human microphthalmia.

Identification of the Bovine Arachnomelia Mutation by Massively Parallel Sequencing Implicates Sulfite Oxidase (SUOX) in Bone Development

Arachnomelia is a monogenic recessive defect of skeletal development in cattle. The causative mutation was previously mapped to a ∼7 Mb interval on chromosome 5. Here we show that array-based sequence capture and massively parallel sequencing technology, combined with the typical family structure in livestock populations, facilitates the identification of the causative mutation. We re-sequenced the entire critical interval in a healthy partially inbred cow carrying one copy of the critical chromosome segment in its ancestral state and one copy of the same segment with the arachnomelia mutation, and we detected a single heterozygous position. The genetic makeup of several partially inbred cattle provides extremely strong support for the causality of this mutation. The mutation represents a single base insertion leading to a premature stop codon in the coding sequence of the SUOX gene and is perfectly associated with the arachnomelia phenotype. Our findings suggest an important role for sulfite oxidase in bone development.

Independent Polled Mutations Leading to Complex Gene Expression Differences in Cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

A deletion in the N-myc downstream regulated gene 1 (NDRG1) gene in Greyhounds with polyneuropathy.

The polyneuropathy of juvenile Greyhound show dogs shows clinical similarities to the genetically heterogeneous Charcot-Marie-Tooth (CMT) disease in humans. The pedigrees containing affected dogs suggest monogenic autosomal recessive inheritance and all affected dogs trace back to a single male. Here, we studied the neuropathology of this disease and identified a candidate causative mutation. Peripheral nerve biopsies from affected dogs were examined using semi-thin histology, nerve fibre teasing and electron microscopy. A severe chronic progressive mixed polyneuropathy was observed. Seven affected and 17 related control dogs were genotyped on the 50k canine SNP chip. This allowed us to localize the causative mutation to a 19.5 Mb interval on chromosome 13 by homozygosity mapping. The NDRG1 gene is located within this interval and NDRG1 mutations have been shown to cause hereditary motor and sensory neuropathy-Lom in humans (CMT4D). Therefore, we considered NDRG1 a positional and functional candidate gene and performed mutation analysis in affected and control Greyhounds. A 10 bp deletion in canine NDRG1 exon 15 (c.1080_1089delTCGCCTGGAC) was perfectly associated with the polyneuropathy phenotype of Greyhound show dogs. The deletion causes a frame shift (p.Arg361SerfsX60) which alters several amino acids before a stop codon is encountered. A reduced level of NDRG1 transcript could be detected by RT-PCR. Western blot analysis demonstrated an absence of NDRG1 protein in peripheral nerve biopsy of an affected Greyhound. We thus have identified a candidate causative mutation for polyneuropathy in Greyhounds and identified the first genetically characterized canine CMT model which offers an opportunity to gain further insights into the pathobiology and therapy of human NDRG1 associated CMT disease. Selection against this mutation can now be used to eliminate polyneuropathy from Greyhound show dogs.

The spatial expression pattern of antimicrobial peptides across the healthy bovine udder.

Antimicrobial peptides are key molecules in local host defense. With the aim to better understand the possible involvement of these peptides in the prevention of bovine mastitis, we determined, for the first time to our knowledge, the spatial pattern of constitutive expression of 5 bovine beta-defensins and bovine psoriasin (S100A7) across 5 localizations of the bovine mammary gland applying a quantitative real-time PCR approach. We observed 3 different expression patterns in the healthy udder: 1) constitutive expression of the lingual and tracheal antimicrobial peptides (LAP and TAP), as well that of bovine neutrophil beta-defensins 4 and 10 (BNBD4 and BNBD10), is essentially restricted to the mammary lymph node; 2) bovine beta-defensin 1 (DEFB1) is mainly expressed in the cisternal epithelium and the Rosette of Fürstenberg; 3) strong constitutive mRNA expression of the calcium-binding protein S100A7, which is also known as psoriasin and which has been reported to be highly active against Escherichia coli, was detected in the streak canal. These results indicate a crucial role of S100A7 in the early-stage prevention of coliform mastitis, and the analyzed beta-defensins might be regarded as inducible weapons against already invading pathogens.