A genomics-informed, SNP association study reveals FBLN1 and FABP4 as contributing to resistance to fleece rot in Australian Merino sheep

A Novel Missense SNP in the Fatty Acid-Binding Protein 4 (FABP4) Gene is Associated with Growth Traits in Karakul and Awassi Sheep

The fatty acid-binding protein 4 (FABP4) plays a crucial role in the transportation and metabolism of fatty acids. It binds to long-chain fatty acids and facilitates their transport within cells. FABP4 is involved in lipid metabolism, insulin sensitivity, inflammation, and energy homeostasis. This study was conducted to assess the association between the FABP4 gene and growth traits in Karakul and Awassi sheep. A PCR-single strand conformation polymorphism (SSCP) protocol was utilized to assess the polymorphism of FABP4 PCR products with growth traits measured at birth, three, six, nine, and twelve-month intervals. One non-synonymous SNP was identified in the second exon, in which threonine was converted to aspartate in the 61st position in FABP4 (p.61Thr > Asp). This novel SNP showed significant associations with all growth traits measured at all age intervals. The results showed that lambs with the TT genotype exhibited higher growth traits than those with the GT and GG genotypes, respectively. The conducted prediction showed a clearly deleterious effect of p.61Thr > Asp on FABP4 structure, which was accompanied by reduced fatty acid binding efficiency. Owing to the predicted damaging effects caused by p.61Thr > Asp on FABP, lower levels of lipid transport and its consequent increased weight gain and other growth trait indices are expected. Therefore, a putative mechanism through which lambs with these genotypes exhibit higher growth traits is proposed. The FABP4 gene is suggested as a promising marker to improve growth traits in Karakul and Awassi sheep. However, more research is required to validate this mechanism.

Traumatic sheep myiasis: A review of the current understanding

Myiasis, or the infestation of live humans and vertebrate animals by dipterous larvae, is a health issue worldwide. The economic impact and potential threat to animal health and wellbeing of this disease under the animal husbandry sector is considerable. Sheep are a highly vulnerable livestock category exposed to myiasis (sheep strike), due to several unique predisposing factors that attract flies. The successful mitigation of this disease relies on a thorough understanding of fly population dynamics associated with the change in weather patterns and the evaluation of this disease through different branches of science such as chemistry, molecular biology, and microbiology. The present review provides a summary of the existing knowledge of strike in sheep, discussed in relation to the application of volatile organic compounds, metagenomics, and molecular biology, and their use regarding implementing fly control strategies such as traps, and to increase the resilience of sheep to this disease through improving their health and wellbeing.

Two cSNPs sites in the fatty acid-binding protein 4 (FABP4) gene and their association analysis with body measurement data in five Chinese cattle breeds

Single nucleotide polymorphisms (SNPs) include the transition and transversion of a single base. In this study, we found two cSNPs (rs110757796 and rs110652478) in the FABP4 (fatty acid-binding protein 4) gene and tested their population genetic parameters in five Chinese cattle breeds. FABP4 is a key marker molecule for lipid production. It plays a crucial role in the growth and development of animals. Thus, we also analyze the association between the two cSNPs of FABP4 gene and body measurement data of Chinese cattle. Our results were indicated that rs110757796 was significantly associated with the chest width in Chinese cattle groups (p < 0.05). In addition, we tested the spatiotemporal expression profile of the bovine FABP4 gene and effect of genetic variation on its expression. It was found that bovine FABP4 gene has tissue-differential expression. Then, the genetic variation located on the exon has a certain effect on the expression of FABP4 gene in bovine muscle. Overall, our results provide that FABP4 could as a candidate gene to improve the progress of cattle molecular breeding.

Ovine Toll-like Receptor 9 (TLR9) Gene Variation and Its Association with Flystrike Susceptibility

Simple Summary

Flystrike is a major ectoparasitic disease of sheep and it creates both an economic impact and welfare issue for the sheep industry. Several factors control the responses of sheep to flystrike, and among these, immune response is regarded as an important factor. Toll-like receptors (TLRs) plays a crucial role in the innate immune system by recognizing pathogen-associated molecular patterns derived from various microbes. Of these receptors, TLR9 recognises unmethylated cytosine-phosphate-guanine (CpG) motifs that are known to be prevalent in bacterial genomes and are also reported in Dipteran insects, including Lucilia cuprina, one of the main species associated with flystrike in sheep. In this study, we used a polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) technique to investigate TLR9 variation in sheep with and without flystrike, and found that variation in a gene region containing the coding sequence of the putative CpG-DNA binding domain was associated with the likelihood of flystrike occurrence. This suggests that variation in ovine TLR9 may affect a sheep’s response to flystrike.

Abstract

Toll-like receptors (TLRs) are a family of proteins that play a role in innate immune responses by recognising pathogen-associated molecular patterns derived from various microbes. Of these receptors, TLR9 recognises bacterial and viral DNA containing unmethylated cytosine-phosphate-guanine (CpG) motifs, and variation in TLR9 has been associated with resistance to various infectious diseases. Flystrike is a problem affecting the sheep industry globally and the immune response of the sheep has been suggested as one factor that influences the response to the disease. In this study, variation in ovine TLR9 from 178 sheep with flystrike and 134 sheep without flystrike was investigated using a polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) approach. These sheep were collected from both commercial and stud farms throughout New Zealand and they were of 13 different breeds, cross-breds and composites. Four alleles of TLR9 were detected, including three previously identified alleles (*01, *02 and *03) and a new allele (*04). In total six single nucleotide polymorphisms (SNPs) were found. Of the three common alleles in the sheep studied, the presence of *03 was found to be associated with a reduced likelihood of flystrike being present (OR = 0.499, p = 0.024). This suggests that variation in ovine TLR9 may affect a sheep’s response to flystrike, and thus the gene may have value as a genetic marker for improving resistance to the disease.

Ovine FABP4 Variation and Its Association With Flystrike Susceptibility

Flystrike is a major cost and a welfare issue for the New Zealand sheep industry. There are several factors that can predispose sheep to flystrike, such as having fleecerot, a urine-stained breech, and “dags” (an accumulation of fecal matter in the wool of the breech). The FABP4 gene (FABP4) has been associated with variation in ovine fleecerot resistance, with a strong genetic correlation existing between fleecerot and flystrike occurrence. In this study, blood samples were collected from sheep with and without flystrike for DNA typing. PCR-SSCP analyses were used to genotype two regions of ovine FABP4. Sheep with the A₁ variant of FABP4 were found to be less likely (odds ratio 0.689, P = 0.014) to have flystrike than those without A₁. The likelihood of flystrike occurrence decreased as copy number of A₁ increased (odds ratio 0.695, P = 0.006). This suggests that FABP4 might be a candidate gene for flystrike resilience in sheep, although further research is required to verify this association.

Identification of Regulatory Functions of LncRNAs Associated With T. circumcincta Infection in Adult Sheep

Several recent studies have demonstrated the role of long non-coding RNAs (lncRNAs) in regulating the defense mechanism against parasite infections, but no studies are available that investigated their relevance for immune response to nematode infection in sheep. Thus, the aim of the current study was to (i) detect putative lncRNAs that are expressed in the abomasal lymph node of adult sheep after an experimental infection with the gastrointestinal nematode (GIN) Teladorsagia circumcincta and (ii) to elucidate their potential functional role associated with the differential host immune response. We hypothesized that putative lncRNAs differentially expressed (DE) between samples from animals that differ in resistance to infection may play a significant regulatory role in response to nematode infection in adult sheep. To obtain further support for our hypothesis, we performed co-expression and functional gene enrichment analyses with the differentially expressed lncRNAs (DE lncRNAs). In a conservative approach, we included for this predictive analysis only those lncRNAs that are confirmed and supported by documentation of expression in gastrointestinal tissues in the current sheep gene atlas. We identified 9,105 putative lncRNA transcripts corresponding to 7,124 gene loci. Of these, 457 were differentially expressed lncRNA loci (DELs) with 683 lncRNA transcripts. Based on a gene co-expression analysis via weighted gene co-expression network analysis, 12 gene network modules (GNMs) were found significantly correlated with at least one of 10 selected target DE lncRNAs. Based on the principle of “guilt-by-association,” the DE genes from each of the three most significantly correlated GNMs were subjected to a gene enrichment analysis. The significant pathways associated with DE lncRNAs included ERK5 Signaling, SAPK/JNK Signaling, RhoGDI Signaling, EIF2 Signaling, Regulation of eIF4 and p70S6K Signaling and Oxidative Phosphorylation pathways. They belong to signaling pathway categories like Cellular Growth, Proliferation and Development, Cellular Stress and Injury, Intracellular and Second Messenger Signaling and Apoptosis. Overall, this lncRNA study conducted in adult sheep after GIN infection provided first insights into the potential functional role of lncRNAs in the differential host response to nematode infection.

Discovery of genes and proteins possibly regulating mean wool fibre diameter using cDNA microarray and proteomic approaches

Wool fibre diameter (WFD) is one of the wool traits with higher economic impact. However, the main genes specifically regulating WFD remain unidentified. In this current work we have used Agilent Sheep Gene Expression Microarray and proteomic technology to investigate the gene expression patterns of body side skin, bearing more wool, in Aohan fine wool sheep, a Chinese indigenous breed, and compared them with that of small tail Han sheep, a sheep bread with coarse wool. Microarray analyses showed that most of the genes likely determining wool diameter could be classified into a few categories, including immune response, regulation of receptor binding and growth factor activity. Certain gene families might play a role in hair growth regulation. These include growth factors, immune cytokines, solute carrier families, cellular respiration and glucose transport amongst others. Proteomic analyses also identified scores of differentially expressed proteins.

Comparative proteomic analyses using iTRAQ-labeling provides insights into fiber diversity in sheep and goats

The structural component of wool and hair fibers, such as keratin-associated proteins (KAPs), has been well described, but the genetic determinants of fiber diameter are largely unknown. Here, we have used an iTRAQ-based proteomic approach to investigate differences in protein abundance among 18 samples from sheep and goats across a diverse range of fibers. We identified proteins with different abundance and are associated with variation in fiber features. Proteins with different abundance are mainly keratin or keratin-associated proteins (KRTAP11-1, KRT6A, KRT38), or are related to hair growth (DSC2, DSG3, EEF2, CALML5, TCHH, SELENBP1) and fatty acid synthesis (FABP4, FABP5). RNA-seq further confirmed the functional importance of the DSC2 gene in the determination of woolly phenotype in goat fibers. This comprehensive analysis of fibers from major fiber-producing animals is the first to provide a list of candidate proteins that are involved in fiber formation. This list will be valuable asset for future studies into the molecular mechanisms that underlie fiber diversity.

BIOLOGICAL SIGNIFICANCE

Proteins are the basis for animal-derived hair fibers, yet proteins conferring fiber structure and characteristics in sheep and goats are largely elusive. By examining 27 fibers samples representing 9 fiber types from sheep and goats through the iTRAQ approach, we show a list of differentially abundant proteins that are important to hair structural component, or genes related to hair growth and fatty acid synthesis. RNA-seq further validated the DSC2 gene is key to the woolly/straight hair phenotype in goats.

Single-nucleotide polymorphism-based genetic diversity analysis of the Kilakarsal and Vembur sheep breeds

AIM

The present study was thus undertaken to analyze the genetic diversity of Kilakarsal and Vembur sheep breeds using single-nucleotide polymorphism (SNP) markers within Toll-like receptor (TLR) 3, 5, 6, 9, and 10 genes.

MATERIALS AND METHODS

Competitive allele-specific polymerase chain reaction (PCR)-based end-point genotyping was performed using real-time PCR to type the SNPs. Allele discrimination module implemented in real-time PCR was utilized to call the genotypes based on fluorescence intensity recorded for each of the two alleles. Basic diversity indices, namely, gene frequencies, observed heterozygosity, expected heterozygosity, and inbreeding coefficient (F_IS), and testing for Hardy-Weinberg equilibrium (HWE) were estimated using package for elementary analysis of SNP data software program.

RESULTS

Of the 25 SNPs, 22 were found to be polymorphic, whereas two SNPs, namely, TLR3_1081_AC and TLR9_2036_CT, were monomorphic in both Kilakarsal and Vembur sheep populations. The SNP TLR10_1180_AG was monomorphic in Kilakarsal but polymorphic in Vembur sheep. The observed heterozygosities were estimated as 0.289 and 0.309 in Kilakarsal and Vembur sheep, respectively, whereas the expected heterozygosity values were 0.305 and 0.309 in the two breeds, respectively. The overall mean F_IS was 0.107 ranging from -0.005 to 0.241 in Kilakarsal sheep and -0.047 ranging from -0.005 to 0.255 in Vembur sheep. In Kilakarsal sheep, the test for HWE revealed TLR9_1308_GC SNP locus with significant deviation (p<0.05) due to heterozygosity deficit. In Vembur sheep, TLR10_82_CT and TLR10_292_CG loci showed significant deviation (p<0.05) due to heterozygosity excess. Other SNP loci did not deviate from HWE (p>0.05) revealing that the population was in HWE proportions.

CONCLUSIONS

The SNP markers within five TLR genes (TLR3, TLR5, TLR6, TLR9, and TLR10) utilized for genotyping in this study were highly polymorphic in Kilakarsal and Vembur breeds of sheep. This study on the genetic diversity analysis of the Kilakarsal and Vembur sheep breeds revealed considerable genetic variation within the breeds and it can be utilized to improve desirable traits.

The Battle Against Flystrike - Past Research and New Prospects Through Genomics

Flystrike, or cutaneous myiasis, is caused by blow fly larvae of the genus Lucilia. This disease is a major problem in countries with large sheep populations. In Australia, Lucilia cuprina (Wiedemann, 1830) is the principal fly involved in flystrike. While much research has been conducted on L. cuprina, including physical, chemical, immunological, genetic and biological investigations, the molecular biology of this fly is still poorly understood. The recent sequencing, assembly and annotation of the draft genome and analyses of selected transcriptomes of L. cuprina have given a first global glimpse of its molecular biology and insights into host-fly interactions, insecticide resistance genes and intervention targets. The present article introduces L. cuprina, flystrike and associated issues, details past control efforts and research foci, reviews salient aspects of the L. cuprina genome project and discusses how the new genomic and transcriptomic resources for this fly might accelerate fundamental molecular research of L. cuprina towards developing new methods for the treatment and control of flystrike.

Identification of genes and proteins associated with anagen wool growth

Identifying genes of major effect for wool growth would offer strategies for improving the quality and increasing the yield of fine wool. In this study, we employed the Agilent Sheep Gene Expression Microarray and proteomic technology to investigate the gene expression patterns of body side skin (more wool growing) in Aohan fine wool sheep (a Chinese indigenous breed) in comparison with groin skin (no wool growing) at the anagen stage of the wool follicle. A microarray study revealed that 4772 probes were differentially expressed, including 2071 upregulated and 2701 downregulated probes, in the comparisons of body side skin vs. groin skin (S/G). The microarray results were verified by means of quantitative PCR. A total of 1099 probes were assigned to unique genes/transcripts. The number of distinct genes/transcripts (annotated) was 926, of which 352 were upregulated and 574 were downregulated. In S/G, 13 genes were upregulated by more than 10 fold, whereas 60 genes were downregulated by more than 10 fold. Further analysis revealed that the majority of the genes possibly related to the wool growth could be assigned to categories including regulation of cell division, intermediate filament, cytoskeletal part and growth factor activity. Several potential gene families may participate in hair growth regulation, including fibroblast growth factors, transforming growth factor-β, WNTs, insulin-like growth factor, vascular endothelial growth factors and so on. Proteomic analysis also revealed 196 differentially expressed protein points, of which 121 were identified as single protein points.

Development of Vibrio spp. infection resistance related SNP markers using multiplex SNaPshot genotyping method in the clam Meretrix meretrix

The clam Meretrix meretrix is a commercially important mollusc species in the coastal areas of South and Southeast Asia. In the present study, large-scale SNPs were genotyped by the Multiplex SNaPshot genotyping method among the stocks of M. meretrix with different Vibrio spp. infection resistance profile. Firstly, the AUTOSNP software was applied to mine SNPs from M. meretrix transcriptome, and 323 SNP loci (including 120 indels) located on 64 contigs were selected based on Uniprot-GO associations. Then, 38 polymorphic SNP loci located on 15 contigs were genotyped successfully in the clam stocks with different resistance to Vibrio parahaemolyticus infection (11-R and 11-S groups). Pearson's Chi-square test was applied to compare the allele and genotype frequency distributions of the SNPs between the different stocks, and seven SNP markers located on three contigs were found to be associated with V. parahaemolyticus infection resistance trait. Haplotype-association analysis showed that six haplotypes had significantly different frequency distributions in 11-S and 11-R (P < 0.05). With selective genotyping between 09-R and 09-C populations, which had different resistance to Vibrio harveyi infection, four out of the seven selected SNPs had significantly different distributions (P < 0.05) and therefore they were considered to be associated with Vibrio spp. infection resistance. Sequence alignments and annotations indicated that the contigs containing the associated SNPs had high similarity to the immune related genes. All these results would be useful for the future marker-assisted selection of M. meretrix strains with high Vibrio spp. infection resistance.

Identification of skin-expressed genes possibly associated with wool growth regulation of Aohan fine wool sheep

Background

Sheep are valuable resources for the animal fibre industry. Therefore, identifying genes which regulate wool growth would offer strategies for improving the quality of fine wool. In this study, we employed Agilent sheep gene expression microarray and proteomic technology to compare the gene expression patterns of the body side (hair-rich) and groin (hairless) skins of Aohan fine wool sheep (a Chinese indigenous breed).

Results

Comparing the body side to the groin skins (S/G) of Aohan fine wool sheep, the microarray study revealed that 1494 probes were differentially expressed, including 602 more highly expressed and 892 less highly expressed probes. The microarray results were verified by means of quantitative PCR. Cluster analysis could distinguish the body side skin and the groin skin. Based on the Database for Annotation, Visualization and Integrated Discovery (DAVID), 38 of the differentially expressed genes were classified into four categories, namely regulation of receptor binding, multicellular organismal process, protein binding and macromolecular complex. Proteomic study revealed that 187 protein spots showed significant (p < 0.05) differences in their respective expression levels. Among them, 46 protein entries were further identified by MALDI-TOF/MS analyses.

Conclusions

Microarray analysis revealed thousands of differentially expressed genes, many of which were possibly associated with wool growth. Several potential gene families might participate in hair growth regulation. Proteomic analysis also indentified hundreds of differentially expressed proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0144-1) contains supplementary material, which is available to authorized users.

Differential expression of genes and proteins associated with wool follicle cycling

Sheep are valuable resources for the wool industry. Wool growth of Aohan fine wool sheep has cycled during different seasons in 1 year. Therefore, identifying genes that control wool growth cycling might lead to ways for improving the quality and yield of fine wool. In this study, we employed Agilent sheep gene expression microarray and proteomic technology to compare the gene expression patterns of the body side skins at August and December time points in Aohan fine wool sheep (a Chinese indigenous breed). Microarray study revealed that 2,223 transcripts were differentially expressed, including 1,162 up-regulated and 1,061 down-regulated transcripts, comparing body side skin at the August time point to the December one (A/D) in Aohan fine wool sheep. Then seven differentially expressed genes were selected to validated the reliability of the gene chip data. The majority of the genes possibly related to follicle development and wool growth could be assigned into the categories including regulation of receptor binding, extracellular region, protein binding and extracellular space. Proteomic study revealed that 84 protein spots showed significant differences in expression levels. Of the 84, 63 protein spots were upregulated and 21 were downregulated in A/D. Finally, 55 protein points were determined through MALDI-TOF/MS analyses. Furthermore, the regulation mechanism of hair follicle might resemble that of fetation.

Identification of more than two paternal haplotypes of the ovine fatty acid-binding protein 4 (FABP4) gene in half-sib families: evidence of intragenic meiotic recombination

The fatty acid binding protein 4 (FABP4) plays an important role in the regulation of lipid metabolism in mammals. In this study, two regions of ovine FABP4 spanning exon 2-intron 2 and exon 3-intron 3 were investigated in four hundred and twenty lambs derived from seven sires that were previously typed as having heterozygous genotypes in both these regions of the gene. These regions have been shown to be variable, with three SNPs plus one indel and four SNPs respectively constituting five and four allele variants in the two regions. Across these regions, fourteen haplotypes have been identified. The lambs were typed using a Polymerase Chain Reaction Single-Stranded Conformational Polymorphism (PCR-SSCP) method to identify the haplotypes inherited from the sires. Between three and four paternally-derived haplotypes were identified in the progeny of six of the seven sires, suggesting that meiotic recombination occurs within ovine FABP4. A number of sequence motifs associated with recombination “hotspots” were detected in the two regions of the gene that were analyzed and these may facilitate the recombination.

Allelic variation in ovine fatty acid-binding protein (FABP4) gene

The fatty acid-binding protein 4 (FABP4) plays a role in lipid metabolism and has been implicated in intra-cellular lipid transport. While FABP4 variation has been reported in some species, variation in the coding sequence has not been reported in sheep. In this study two regions of ovine FABP4 were analysed using PCR-SSCP and sequencing. Five different PCR-SSCP patterns, representing five specific sequences (A ( 1 )-E ( 1 )) were detected in region 1 (exon 2-intron 2) with sequence analysis revealing three nucleotide substitutions and one deletion in the intron. In region 2 (exon 3-intron 3), four different PCR-SSCP patterns (A ( 2 )-D ( 2 )) were observed and four nucleotide substitutions were revealed. In total, fourteen haplotypes through both regions were defined. There was a difference (P < 0.001) in allele frequencies between two selection lines of Coopworth sheep that have been bred over many generations to be lean or fat. In region 1, A ( 1 ) and B ( 1 ) were most common (at a frequency of 50 and 30 % respectively) in the fat line, whereas these two variants were absent or rare in the lean line in which C ( 1 ) predominated (89 %). In region 2, C ( 2 ) was the most common variant (59 %) in the lean line but was absent in the fat line, whereas B ( 2 ) was predominant (83 %) in the fat line but was rare (3 %) in the lean line. These results indicate that ovine FABP4 is polymorphic and suggest further analysis is required to see if the variation detected affects fat deposition in sheep.

miRNA-mediated relationships between Cis-SNP genotypes and transcript intensities in lymphocyte cell lines

In metazoans, miRNAs regulate gene expression primarily through binding to target sites in the 3' UTRs (untranslated regions) of messenger RNAs (mRNAs). Cis-acting variants within, or close to, a gene are crucial in explaining the variability of gene expression measures. Single nucleotide polymorphisms (SNPs) in the 3' UTRs of genes can affect the base-pairing between miRNAs and mRNAs, and hence disrupt existing target sites (in the reference sequence) or create novel target sites, suggesting a possible mechanism for cis regulation of gene expression. Moreover, because the alleles of different SNPs within a DNA sequence of limited length tend to be in strong linkage disequilibrium (LD), we hypothesize the variants of miRNA target sites caused by SNPs potentially function as bridges linking the documented cis-SNP markers to the expression of the associated genes. A large-scale analysis was herein performed to test this hypothesis. By systematically integrating multiple latest information sources, we found 21 significant gene-level SNP-involved miRNA-mediated post-transcriptional regulation modules (SNP-MPRMs) in the form of SNP-miRNA-mRNA triplets in lymphocyte cell lines for the CEU and YRI populations. Among the cognate genes, six including ALG8, DGKE, GNA12, KLF11, LRPAP1, and MMAB are related to multiple genetic diseases such as depressive disorder and Type-II diabetes. Furthermore, we found that ~35% of the documented transcript intensity-related cis-SNPs (~950) in a recent publication are identical to, or in significant linkage disequilibrium (LD) (p<0.01) with, one or multiple SNPs located in miRNA target sites. Based on these associations (or identities), 69 significant exon-level SNP-MPRMs and 12 disease genes were further determined for two populations. These results provide concrete in silico evidence for the proposed hypothesis. The discovered modules warrant additional follow-up in independent laboratory studies.

Using regulatory and epistatic networks to extend the findings of a genome scan: identifying the gene drivers of pigmentation in merino sheep

Extending genome wide association analysis by the inclusion of gene expression data may assist in the dissection of complex traits. We examined piebald, a pigmentation phenotype in both human and Merino sheep, by analysing multiple data types using a systems approach. First, a case control analysis of 49,034 ovine SNP was performed which confirmed a multigenic basis for the condition. We combined these results with gene expression data from five tissue types analysed with a skin-specific microarray. Promoter sequence analysis of differentially expressed genes allowed us to reverse-engineer a regulatory network. Likewise, by testing two-loci models derived from all pair-wise comparisons across piebald-associated SNP, we generated an epistatic network. At the intersection of both networks, we identified thirteen genes with insulin-like growth factor binding protein 7 (IGFBP7), platelet-derived growth factor alpha (PDGFRA) and the tetraspanin platelet activator CD9 at the kernel of the intersection. Further, we report a number of differentially expressed genes in regions containing highly associated SNP including ATRN, DOCK7, FGFR1OP, GLI3, SILV and TBX15. The application of network theory facilitated co-analysis of genetic variation with gene expression, recapitulated aspects of the known molecular biology of skin pigmentation and provided insights into the transcription regulation and epistatic interactions involved in piebald Merino sheep.