Transcription profiling provides insights into gene pathways involved in horn and scurs development in cattle

IL-1β promotes A7r5 and HASMC migration and invasion via the p38-MAPK/Angpt-2 pathway

The migration, proliferation, and inflammatory factor secretion of vascular smooth muscle cells (VSMCs) are involved in the important pathological processes of several vascular occlusive diseases, including coronary atherosclerosis (CAS). Interleukin 1β(IL-1β), as a bioactive mediator of VSMC synthesis and secretion, can promote the pathological progress of CAS. In this study, we further explored the underlying molecular mechanisms by which IL-1β regulates VSMC migration, invasion. We pretreated A7r5 and HASMC with IL-1β for 24 h, and measured the expression of IL-1β, proliferating cell nuclear antigen (PCNA), cyclin D1, matrix metalloproteinase 2 (MMP2) and matrix metalloproteinase 2 (MMP9) in the cells by Western blotting. Cell migration and invasion ability were measured by Transwell and wound healing assays. Cell viability was measured by an MTT assay. We found that IL-1β upregulated the expression of proliferation-related proteins (PCNA and Cyclin D1) in A7r5 and HASMC, and induces the secretion of MMP2 and MMP9, promotes cell invasion and migration. In addition, in A7r5 and HASMCs treated with IL-1β, the expression of Angiopoietin-2 (Angpt-2) increased in a time-dependent manner, transfection with si-Angpt-2 suppressed cell migration and invasion, with downregulated MMP2 and MMP9 expression. Parallelly, we further found that the p38-MAPK pathway is activated in cells induced by IL-1β, p38-MAPK inhibitors can down-regulate the expression of Angpt-2. Collectively, these data demonstrated that IL-1β promotes A7r5 and HASMC migration and invasion via the p38-MAPK/Angpt-2 pathway.

The Complex and Diverse Genetic Architecture of the Absence of Horns (Polledness) in Domestic Ruminants, including Goats and Sheep

Horns are the most obvious common feature of Bovidae. The naturally occurring absence of horns in these species, also known as polledness, is of surprisingly heterogeneous nature, although they are Mendelian traits. This review compares in detail the molecular differences among the causes of inherited polledness in the domestic ruminant species of cattle, yak, sheep, and goat based on the causal gene variants that have been discovered in recent years. The genetic causes for the lack of horns in small ruminants seem not only to be more complex, e.g., in sheep, breed-specific characteristics are still unexplained, but in goats, there is also the associated disorder of intersexuality—polled intersex syndrome (PIS). In connection with animal welfare and the associated discussion about a legal ban on the dehorning of all farm animals, naturally hornless animals and the causal genetic variants are of increasing research interest in the age of genome editing. However, the low acceptance of genetic engineering in livestock, especially in European societies, limits its use in food-producing animals. Therefore, genotype-based targeted selection of naturally occurring variants is still a widely used method for spreading this desired trait within and across populations, at least in cattle and sheep.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Migration of BEAS-2B cells enhanced by H1299 cell derived-exosomes

Previous studies reported that exosomes (Exos) secreted by tumor cells could affect the tumor cells themselves and normal cells. However, the effects of exosomes derived from tumor cells on normal cells' migration and mechanical characteristics are rarely reported. This work explores the effects of H1299 cell-derived exosomes (H1299-Exos) on the migration of BEAS-2B cells, and analyzes possible mechanical mechanisms. In the experiments, exosomes were isolated from the culture supernatants of H1299 cells by ultracentrifugation. The H1299-Exos were confirmed by scanning electron microscope (SEM) and western blotting (WB). The BEAS-2B cell migration was assessed using scratch assays. Cytoskeletal structure changes were detected by immunofluorescence. Surface morphology and mechanical properties were measured by atomic force microscopy (AFM). After incubation with H1299-Exos for 48 h, BEAS-2B cells enhanced migration ability, with increased filopodia and cytoskeletal rearrangements. The changes in the morphology and mechanical properties of the cells caused by H1299-Exos were detected using AFM, including the increase in cell length and the decrease in cell height, Young's modulus and adhesion. In short, H1299-Exos promoted the BEAS-2B cell migrations. It indicates that the morphological and mechanical properties can be used as a means to assess normal cell alterations induced by tumor cell derived-exosomes. This provides a method for studying the effects of exosomes secreted by tumor cells on normal cells and the changes in their physical properties.

Heterozygosity of the Celtic polled locus in Canadian scurred beef cattle

Polled cattle are preferable to horned or scurred animals because they are safer for handling and cause less bruising. Although DNA testing can determine horned/polled genotype, scurs may appear in polled animals. The inheritance of scurs is complex because it is a sex-influenced trait that interacts with the polled locus. We demonstrate that in 685 purebred and crossbred Canadian beef cattle, all 153 scurred animals were heterozygous polled at the Celtic variant. In addition, male obligate carriers of scurs were smooth polled when homozygous for the polled mutation. Scurred and non-scurred males were sequenced for five genes (CTDNEP1, SHBG, SOX15, FGF11, and DHRS7C) within the scur candidate region on BTA19 that are functionally related to bone development and hormone regulation. Multipoint linkage analysis was conducted using 18 microsatellite markers and two informative variants (DHRS7C g.29594018G>C and CTDNEP1 c.462G>A) in the scurred families and further supported mapping on BTA19 between BMS2142 (logarithm of the odds (LOD) = 5.42) and IDVGA46 (LOD = 3.47). These data indicate epistatic interactions between the scurred and polled loci and emphasise the necessity for a scurred DNA test to assist purebred beef producers in eradicating the scur trait.

Are scurs in heterozygous polled (Pp) cattle a complex quantitative trait?

Background

Breeding genetically hornless, i.e. polled, cattle provides an animal welfare-friendly and non-invasive alternative to the dehorning of calves. However, the molecular regulation of the development of horns in cattle is still poorly understood. Studying genetic characters such as polledness and scurs, can provide valuable insights into this process. Scurs are hornlike formations that occur occasionally in a wide variety of sizes and forms as an unexpected phenotype when breeding polled cattle.

Methods

We present a unique dataset of 885 Holstein–Friesian cattle with polled parentage. The horn phenotype was carefully examined, and the phenotypic heterogeneity of the trait is described. Using a direct gene test for polledness, the polled genotype of the animals was determined. Subsequently, the existence of a putative scurs locus was investigated using high-density genotype data of a selected subset of 232 animals and two mapping approaches: mixed linear model-based association analyses and combined linkage disequilibrium and linkage analysis.

Results

The results of an exploratory data analysis indicated that the expression of scurs depends on age at phenotyping, sex and polled genotype. Scurs were more prevalent in males than in females. Moreover, homozygous polled animals did not express any pronounced scurs and we found that the Friesian polled allele suppresses the development of scurs more efficiently than the Celtic polled allele. Combined linkage and linkage disequilibrium mapping revealed four genome-wide significant loci that affect the development of scurs, one on BTA5 and three on BTA12. Moreover, suggestive associations were detected on BTA16, 18 and 23. The mixed linear model-based association analysis supports the results of the combined linkage and linkage disequilibrium analysis. None of the mapping approaches provided convincing evidence for a monogenic inheritance of scurs.

Conclusions

Our results contradict the initial and still broadly accepted model for the inheritance of horns and scurs. We hypothesise an oligogenetic model to explain the development of scurs and polledness.

Optimized Genetic Testing for Polledness in Multiple Breeds of Cattle

Many breeds of modern cattle are naturally horned, and for sound husbandry management reasons the calves frequently undergo procedures to physically remove the horns by disbudding or dehorning. These procedures are however a welfare concern. Selective breeding for polledness - absence of horns - has been effective in some cattle breeds but not in others (Bos indicus genotypes) due in part to the complex genetics of horn phenotype. To address this problem different approaches to genetic testing which provide accurate early-in-life prediction of horn phenotype have been evaluated, initially using microsatellites (MSAT) and more recently single nucleotide polymorphism (SNP). A direct gene test is not effective given the genetic heterogeneity and large-sized sequence variants associated with polledness in different breeds. The current study investigated 39,943 animals of multiple breeds to assess the accuracy of available poll testing assays. While the standard SNP-based test was an improvement on the earlier MSAT haplotyping method, 1999 (9.69%) out of 20,636 animals tested with this SNP-based assay did not predict a genotype, most commonly associated with the Indicus-influenced breeds. The current study has developed an optimized poll gene test that resolved the vast majority of these 1999 unresolved animals, while the predicted genotypes of those previously resolved remained unchanged. Hence the optimized poll test successfully predicted a genotype in 99.96% of samples assessed. We demonstrated that a robust set of 5 SNPs can effectively determine P_C and P_F alleles and eliminate the ambiguous and undetermined results of poll gene testing previously identified as an issue in cattle.

Understanding the effects of the bovine POLLED variants

Horns are paired appendages on the head of bovine species, comprising an inner bony core and outer keratin sheath. The horn bud forms during early fetal development but ossification of the developing horn does not occur until approximately 1 month after birth. Little is known about the genetic pathways that lead to horn growth. Hornless, or polled, animals are found in all domestic bovids. Histological studies of bovine fetuses have shown that the horn bud does not form in polled individuals. There are currently four known genetic variants for polledness in cattle on BTA1. All of the variants are intergenic, but probably affect regulation of nearby genes or long non-coding RNAs. Transcriptomic studies suggest that the expression of two nearby long non-coding RNAs are affected by the Celtic POLLED variant, but further studies are required to confirm these data. Candidate genes located elsewhere in the genome are involved in regulating bone formation and epithelial-to-mesenchymal transition. Expression of one of these candidate genes, RXFP2, appears to be reduced in the fetal horn bud of polled animals carrying the Celtic variant compared with horned individuals. Investigating horn ontogenesis and the genetic pathway by which the POLLED variants prevent horn development has implications for cattle breeding. If the genetic basis of horn bud formation and polledness is better understood, then new targets may be identified for precision genome editing to create polled individuals.

The effectiveness of clove oil and two different cautery disbudding methods on preventing horn growth in dairy goat kids

The effectiveness of clove oil and cautery disbudding on horn growth was evaluated in goat kids. The study used 243 Saanen doe kids (4±1 days old; mean±SD) on two goat farms that were disbudded with either (i) clove oil injection (CLOVE), (ii) a cautery iron and bud removed (BUDOFF), or (iii) a cautery iron with bud left intact (BUDON). Each kid received a different treatment per bud, which were balanced between buds (left/right) and randomly allocated. A trained observer monitored bud growth following treatment for 3 months recording either: N: no growth, H: normal horn, S: abnormal horn (scur), or SC: soft, fibrous lump (scorn). After the final observation, buds were assessed for the probability of detecting (i) success (no growth), (ii) scurs, (iii) horns or (iv) scorns [with 95% CI]. The probability of success for BUDOFF (0.77 [0.63, 0.87]) was higher than for BUDON (0.20 [0.11, 0.34]) and CLOVE (0.09 [0.04, 0.18]; P ≤ 0.05). Furthermore, the probability of success for BUDON was higher than for CLOVE (P ≤ 0.05). The probability of scurs was higher for CLOVE (0.72 [0.63, 0.80]) than BUDOFF (0.25 [0.17, 0.34]) and BUDON (0.30 [0.21, 0.39]; P ≤ 0.05). There was no difference in the probability of scurs for BUDOFF and BUDON (P > 0.05). The probability of horns was higher for CLOVE (0.21 [0.15, 0.29]) than BUDON (0.02 [0.01, 0.06]; P ≤ 0.05); horns were not observed for BUDOFF. The probability of scorns for BUDON, the only treatment that led to scorns, was 0.41 (0.25, 0.60). These results suggest that BUDOFF was more effective at preventing growth than CLOVE and BUDON and appears the most effective method, of the methods tested, for disbudding kids. Future research should explore other alternatives to cautery disbudding that may be both efficacious and cause less pain.

Comparative iTRAQ proteomics revealed proteins associated with horn development in yak

Background

The practice of dehorning yak raises animal safety concerns, which have been addressed by selective breeding to obtain genetically hornless yak. The POLLED locus in yak has been studied extensively; however, little is known regarding the proteins that regulate horn bud development.

Methods

A differential proteomic analysis was performed to compare the skin from the horn bud region of polled yak fetuses and the horn bud tissue of horned yak fetuses using isobaric tags for relative and absolute quantitation (iTRAQ) technology coupled with 2D LC-MS/MS.

Results

One hundred differentially abundant proteins (DAPs) were identified. Of these, 29 were up-regulated and 71 were down-regulated in skin from the horn bud region of polled fetuses when compared to the horn bud tissue of horned fetuses. Bioinformatics analyses showed that the up-regulated DAPs were mainly associated with metabolic activities, while the down-regulated DAPs were significantly enriched in cell adhesion and cell movement activities.

Conclusions

We concluded that some important proteins were associated with cell adhesion, cell motility, keratinocyte differentiation, cytoskeleton organization, osteoblast differentiation, and fatty acid metabolism during horn bud development. These results advance our understanding of the molecular mechanisms underlying horn development.

Simultaneous introgression of three POLLED mutations into a synthetic breed of Chinese cattle

The polled phenotype of cattle is increasingly becoming favourable mainly because of the enhanced emphasis on animal welfare, for which the causative mutations have been reported during the past years. The Shuxuan cattle are a new synthetic breed by crossing the indigenous cattle with both Simmental and Holstein semen in Sichuan of Southwest China, in which about 15% of polled individuals have newly emerged. Because official record about POLLED genotypes for the historically imported sires is unavailable, we therefore genotyped the proposed POLLED variants of P_202ID, P_80kbID and P_219ID among 48 polled and 16 horned Shuxuan cattle. It was first revealed that all three candidate mutations have been simultaneously introgressed into Shuxuan cattle, whereas the P_202ID mutation is dominant. Furthermore, one polled animal still remains to carry none of the three candidate mutations, which suggests that further mutation(s) would also exist. Additionally, we sequenced mitochondrial DNA and found that Shuxuan cattle are composed of two matrilineal origins of Bos taurus (65.6%) and B. indicus (34.4%); and there is no origin-biased distribution of polled phenotype. In conclusion, our study first supports the recently reported novel candidate mutation of P_219ID and detects simultaneous presences of all three known POLLED mutations within a cattle breed.

Genome-wide Association Study Identifies Loci for the Polled Phenotype in Yak

The absence of horns, known as the polled phenotype, is an economically important trait in modern yak husbandry, but the genomic structure and genetic basis of this phenotype have yet to be discovered. Here, we conducted a genome-wide association study with a panel of 10 horned and 10 polled yaks using whole genome sequencing. We mapped the POLLED locus to a 200-kb interval, which comprises three protein-coding genes. Further characterization of the candidate region showed recent artificial selection signals resulting from the breeding process. We suggest that expressional variations rather than structural variations in protein probably contribute to the polled phenotype. Our results not only represent the first and important step in establishing the genomic structure of the polled region in yak, but also add to our understanding of the polled trait in bovid species.

Evaluation of breeding strategies for polledness in dairy cattle using a newly developed simulation framework for quantitative and Mendelian traits

Background

Intensified selection of polled individuals has recently gained importance in predominantly horned dairy cattle breeds as an alternative to routine dehorning. The status quo of the current polled breeding pool of genetically-closely related artificial insemination sires with lower breeding values for performance traits raises questions regarding the effects of intensified selection based on this founder pool.

Methods

We developed a stochastic simulation framework that combines the stochastic simulation software QMSim and a self-designed R program named QUALsim that acts as an external extension. Two traits were simulated in a dairy cattle population for 25 generations: one quantitative (QMSim) and one qualitative trait with Mendelian inheritance (i.e. polledness, QUALsim). The assignment scheme for qualitative trait genotypes initiated realistic initial breeding situations regarding allele frequencies, true breeding values for the quantitative trait and genetic relatedness. Intensified selection for polled cattle was achieved using an approach that weights estimated breeding values in the animal best linear unbiased prediction model for the quantitative trait depending on genotypes or phenotypes for the polled trait with a user-defined weighting factor.

Results

Selection response for the polled trait was highest in the selection scheme based on genotypes. Selection based on phenotypes led to significantly lower allele frequencies for polled. The male selection path played a significantly greater role for a fast dissemination of polled alleles compared to female selection strategies. Fixation of the polled allele implies selection based on polled genotypes among males. In comparison to a base breeding scenario that does not take polledness into account, intensive selection for polled substantially reduced genetic gain for this quantitative trait after 25 generations. Reducing selection intensity for polled males while maintaining strong selection intensity among females, simultaneously decreased losses in genetic gain and achieved a final allele frequency of 0.93 for polled.

Conclusions

A fast transition to a completely polled population through intensified selection for polled was in contradiction to the preservation of high genetic gain for the quantitative trait. Selection on male polled genotypes with moderate weighting, and selection on female polled phenotypes with high weighting, could be a suitable compromise regarding all important breeding aspects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0228-7) contains supplementary material, which is available to authorized users.

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.

Associations of single nucleotide polymorphisms in candidate genes with the polled trait in Datong domestic yaks

The domestic yak (Bos grunniens) is an iconic symbol of animal husbandry at high altitudes. Yaks exhibit unique external characteristics including long hair and large horns. However, hornless yaks can be found in different breeds and different populations. The hornless trait is also known as polled, and the POLL locus has been fine-mapped to chromosome 1 in cattle (Bos taurus), although the underlying genetic basis of the polled trait is still unclear in the yak. Thus, we performed an association study to identify the genetic polymorphisms responsible for the polled trait in the yak. Fifty polled Datong domestic yaks and 51 horned individuals were selected randomly from a huge herd and were used as the case and control groups respectively for the association analysis. Twelve genes located in the candidate region of the POLL locus in cattle were used as references to detect DNA polymorphisms related to yak polledness, which were analyzed by sequencing and a high-resolution melting test. We applied Fisher's exact test and haplotype analysis to show that a 147-kb segment that included three protein-coding genes C1H21orf62, GCFC1 and SYNJ1 was the most likely location of the POLL mutation in domestic yaks.

A Single Nucleotide Polymorphism within the Interferon Gamma Receptor 2 Gene Perfectly Coincides with Polledness in Holstein Cattle

Polledness is a high impact trait in modern milk and beef production to meet the demands of animal welfare and work safety. Previous studies have mapped the polled-locus to the proximal region of the bovine chromosome 1 (BTA1) and narrowed it down to approximately 1 Mb. Sequencing of the positional candidate genes within the 1 Mb polled region and whole genome sequencing of Holsteins revealed a single nucleotide polymorphism (SNP) AC000158: g.1390292G>A within intron 3 of the interferon gamma receptor 2 gene (IFNGR2) in perfect co-segregation with polledness in Holsteins. This complete association was validated in 443 animals of the same breed. This SNP allows reliable genotyping of horned, heterozygous and homozygous polled Holsteins, even in animals that could not be resolved using the previously published haplotype for Holstein.

Novel insights into the bovine polled phenotype and horn ontogenesis in Bovidae

Despite massive research efforts, the molecular etiology of bovine polledness and the developmental pathways involved in horn ontogenesis are still poorly understood. In a recent article, we provided evidence for the existence of at least two different alleles at the Polled locus and identified candidate mutations for each of them. None of these mutations was located in known coding or regulatory regions, thus adding to the complexity of understanding the molecular basis of polledness. We confirm previous results here and exhaustively identify the causative mutation for the Celtic allele (P_C) and four candidate mutations for the Friesian allele (P_F). We describe a previously unreported eyelash-and-eyelid phenotype associated with regular polledness, and present unique histological and gene expression data on bovine horn bud differentiation in fetuses affected by three different horn defect syndromes, as well as in wild-type controls. We propose the ectopic expression of a lincRNA in P_C/p horn buds as a probable cause of horn bud agenesis. In addition, we provide evidence for an involvement of OLIG2, FOXL2 and RXFP2 in horn bud differentiation, and draw a first link between bovine, ovine and caprine Polled loci. Our results represent a first and important step in understanding the genetic pathways and key process involved in horn bud differentiation in Bovidae.

Inferring the in vivo cellular program of developing bovine skeletal muscle from expression data

We outline an in vivo cellular program of bovine longissimus muscle development inferred from expression data from 60 days post conception to 3months postnatal. Analytic challenges included changes in cellular composition, ambiguous 'diagnostic' markers of cell type and contrasts between cattle human and mouse myogenesis. Nevertheless, the expression profiles of the myosin isoforms support slow and fast muscle fibres emanating from primary and secondary myogenesis respectively, while expression of the prenatal myosin subunits is down regulated prior to birth. Of the canonical pro-myogenic transcription factors (TF), MYF6 and MYF5 are negatively co-expressed, with MYF6 displaying higher expression in the post-natal samples and MYF5, MYOG, HES6 and PAX7 displaying higher expression in early development. A set of TFs (SIX1, EYA2 and DACH2) considered important in undifferentiated murine cells were equally abundant in differentiated bovine cells. An examination of mammalian regulators of fibre composition, muscle mass and muscle metabolism, underscored the roles of PPARGC1A, TGFβ signalling and the NHR4 Nuclear Hormone Receptors on bovine muscle development. Enriched among the most variably expressed genes from the entire data set were molecules regulating mitochondrial metabolism of carbohydrate (PDK4), fat (UCP3), protein (AGXT2L1) and high energy phosphate (CKMT2). The dramatic increase in the expression of these transcripts, which may enable the peri-natal transition to metabolic independence critical for new-born herbivores, provides surprising evidence for substantial developmental remodelling of muscle mitochondria and reflects changes in nutrient availability. Overall, despite differences in size, metabolism and physiology, the muscle structural subunit expression program appears very similar in ruminants, rodents and humans.

Gene expression phenotypes for lipid metabolism and intramuscular fat in skeletal muscle of cattle

Gene expression phenotypes were evaluated for intramuscular fat (IMF) in bovine skeletal muscle as an alternative to traditional estimates of IMF%. Gene expression data from a time course of LM development in high- and low-marbling Bos taurus cattle crosses were compared to identify genes involved in intramuscular adipocyte lipid metabolism with developmentally similar gene expression profiles. Three sets of genes were identified: triacylglyceride (TAG) synthesis and storage, fatty acid (FA) synthesis, and PPARγ-related genes. In an independent analysis in the LM of 48 Bos indicus cattle, TAG and FA gene sets were enriched in the top 100 genes of which expression was most correlated with IMF% (P = 1.2 × 10(-24) and 3.5 × 10(-9), respectively). In general, genes encoding enzymes involved in the synthesis of FA and TAG in the intramuscular adipocytes were present in the top 100 genes. In B. indicus, effects of a steroid hormone growth promotant (HGP), 2 experimental sites [New South Wales (NSW) and Western Australia (WA)], and 3 tenderness genotypes on the expression levels of genes in the TAG gene set and the correlation of gene expression with IMF% were investigated. Although correlation between expression of 12 individual TAG genes and IMF% was observed in HGP-treated animals in both experimental sites (mean r = 0.43), correlation was not observed for untreated animals at the NSW site (mean r = -0.07, P < 3 × 10(-6)). However, TAG genes showed an average 1.6-fold (P < 0.0004) reduction in expression in the LM of HGP-treated cattle relative to untreated cattle, an effect consistent across both experimental sites. Cattle possessing the favored tenderness calpain 1 and 3 and calpastatin alleles exhibited a greater (P = 0.008) reduction in expression in NSW (1.8-fold reduction, P = 0.0002) compared with WA (1.2-fold reduction, P = 0.03). Tenderness genotype had no impact (P > 0.05) on the correlation of TAG genes with IMF%. In general, the interactions among genotype, treatment and location, and TAG gene set gene expression were consistent with the interactions among the same factors and IMF% detected using 255 animals, of which the 48 in this study were a subset. Thus, the TAG gene set constitutes a gene expression phenotype able to predict effects of different genotypes and treatments on IMF% using much smaller groups than current approaches, even in animals with very low IMF%.

An Always Correlated gene expression landscape for ovine skeletal muscle, lessons learnt from comparison with an "equivalent" bovine landscape

Background

We have recently described a method for the construction of an informative gene expression correlation landscape for a single tissue, longissimus muscle (LM) of cattle, using a small number (less than a hundred) of diverse samples. Does this approach facilitate interspecies comparison of networks?

Findings

Using gene expression datasets from LM samples from a single postnatal time point for high and low muscling sheep, and from a developmental time course (prenatal to postnatal) for normal sheep and sheep exhibiting the Callipyge muscling phenotype gene expression correlations were calculated across subsets of the data comparable to the bovine analysis. An “Always Correlated” gene expression landscape was constructed by integrating the correlations from the subsets of data and was compared to the equivalent landscape for bovine LM muscle. Whilst at the high level apparently equivalent modules were identified in the two species, at the detailed level overlap between genes in the equivalent modules was limited and generally not significant. Indeed, only 395 genes and 18 edges were in common between the two landscapes.

Conclusions

Since it is unlikely that the equivalent muscles of two closely related species are as different as this analysis suggests, within tissue gene expression correlations appear to be very sensitive to the samples chosen for their construction, compounded by the different platforms used. Thus users need to be very cautious in interpretation of the differences. In future experiments, attention will be required to ensure equivalent experimental designs and use cross-species gene expression platform to enable the identification of true differences between different species.

A 3.7 Mb deletion encompassing ZEB2 causes a novel polled and multisystemic syndrome in the progeny of a somatic mosaic bull

Polled and Multisystemic Syndrome (PMS) is a novel developmental disorder occurring in the progeny of a single bull. Its clinical spectrum includes polledness (complete agenesis of horns), facial dysmorphism, growth delay, chronic diarrhea, premature ovarian failure, and variable neurological and cardiac anomalies. PMS is also characterized by a deviation of the sex-ratio, suggesting male lethality during pregnancy. Using Mendelian error mapping and whole-genome sequencing, we identified a 3.7 Mb deletion on the paternal bovine chromosome 2 encompassing ARHGAP15, GTDC1 and ZEB2 genes. We then produced control and affected 90-day old fetuses to characterize this syndrome by histological and expression analyses. Compared to wild type individuals, affected animals showed a decreased expression of the three deleted genes. Based on a comparison with human Mowat-Wilson syndrome, we suggest that deletion of ZEB2, is responsible for most of the effects of the mutation. Finally sperm-FISH, embryo genotyping and analysis of reproduction records confirmed somatic mosaicism in the founder bull and male-specific lethality during the first third of gestation. In conclusion, we identified a novel locus involved in bovid horn ontogenesis and suggest that epithelial-to-mesenchymal transition plays a critical role in horn bud differentiation. We also provide new insights into the pathogenicity of ZEB2 loss of heterozygosity in bovine and humans and describe the first case of male-specific lethality associated with an autosomal locus in a non-murine mammalian species. This result sets PMS as a unique model to study sex-specific gene expression/regulation.

Bovine polledness--an autosomal dominant trait with allelic heterogeneity

The persistent horns are an important trait of speciation for the family Bovidae with complex morphogenesis taking place briefly after birth. The polledness is highly favourable in modern cattle breeding systems but serious animal welfare issues urge for a solution in the production of hornless cattle other than dehorning. Although the dominant inhibition of horn morphogenesis was discovered more than 70 years ago, and the causative mutation was mapped almost 20 years ago, its molecular nature remained unknown. Here, we report allelic heterogeneity of the POLLED locus. First, we mapped the POLLED locus to a ∼381-kb interval in a multi-breed case-control design. Targeted re-sequencing of an enlarged candidate interval (547 kb) in 16 sires with known POLLED genotype did not detect a common allele associated with polled status. In eight sires of Alpine and Scottish origin (four polled versus four horned), we identified a single candidate mutation, a complex 202 bp insertion-deletion event that showed perfect association to the polled phenotype in various European cattle breeds, except Holstein-Friesian. The analysis of the same candidate interval in eight Holsteins identified five candidate variants which segregate as a 260 kb haplotype also perfectly associated with the POLLED gene without recombination or interference with the 202 bp insertion-deletion. We further identified bulls which are progeny tested as homozygous polled but bearing both, 202 bp insertion-deletion and Friesian haplotype. The distribution of genotypes of the two putative POLLED alleles in large semi-random sample (1,261 animals) supports the hypothesis of two independent mutations.

SNP-based association mapping of the polled gene in divergent cattle breeds

Naturally, hornless cattle are called polled. Although the POLL locus could be assigned to a c. 1.36-Mb interval in the centromeric region of BTA1, the underlying genetic basis for the polled trait is still unknown. Here, an association mapping design was set up to refine the candidate region of the polled trait for subsequent high-throughput sequencing. The case group comprised 101 homozygous polled animals from nine divergent cattle breeds, the majority represented by Galloway, Angus, Fleckvieh and Holstein Friesian. Additionally, this group included some polled individuals of Blonde d'Aquitaine, Charolais, Hereford, Jersey and Limousin breeds. The control group comprised horned Belgian Blue, Fleckvieh, Holstein Friesian and Illyrian Buša cattle. A genome-wide scan using 49,163 SNPs was performed, which revealed one shared homozygous haplotype block consisting of nine neighbouring SNPs in all polled animals. This segment defines a 381-kb interval on BTA1 that we consider to be the most likely location of the POLL mutation. Our results further demonstrate that the polled-associated haplotype is also frequent in horned animals included in this study, and thus the haplotype as such cannot be used for population-wide genetic testing. The actual trait-associated haplotype may be revealed by using higher-density SNP arrays. For the final identification of the causal mutation, we suggest high-throughput sequencing of the entire candidate region, because the identification of functional candidate genes is difficult owing to the lack of a comparable model.

A newly described bovine type 2 scurs syndrome segregates with a frame-shift mutation in TWIST1

The developmental pathways involved in horn development are complex and still poorly understood. Here we report the description of a new dominant inherited syndrome in the bovine Charolais breed that we have named type 2 scurs. Clinical examination revealed that, despite a strong phenotypic variability, all affected individuals show both horn abnormalities similar to classical scurs phenotype and skull interfrontal suture synostosis. Based on a genome-wide linkage analysis using Illumina BovineSNP50 BeadChip genotyping data from 57 half-sib and full-sib progeny, this locus was mapped to a 1.7 Mb interval on bovine chromosome 4. Within this region, the TWIST1 gene encoding a transcription factor was considered as a strong candidate gene since its haploinsufficiency is responsible for the human Saethre-Chotzen syndrome, characterized by skull coronal suture synostosis. Sequencing of the TWIST1 gene identified a c.148_157dup (p.A56RfsX87) frame-shift mutation predicted to completely inactivate this gene. Genotyping 17 scurred and 20 horned founders of our pedigree as well as 48 unrelated horned controls revealed a perfect association between this mutation and the type 2 scurs phenotype. Subsequent genotyping of 32 individuals born from heterozygous parents showed that homozygous mutated progeny are completely absent, which is consistent with the embryonic lethality reported in Drosophila and mouse suffering from TWIST1 complete insufficiency. Finally, data from previous studies on model species and a fine description of type 2 scurs symptoms allowed us to propose different mechanisms to explain the features of this syndrome. In conclusion, this first report on the identification of a potential causal mutation affecting horn development in cattle offers a unique opportunity to better understand horn ontogenesis.