A nonsense mutation in PLD4 is associated with a zinc deficiency-like syndrome in Fleckvieh cattle
- Simone Jung†1,
- Hubert Pausch†1,
- Martin C Langenmayer3,
- Hermann Schwarzenbacher4,
- Monir Majzoub-Altweck3,
- Nicole S Gollnick2 and
- Ruedi Fries1Email author
© Jung et al.; licensee BioMed Central Ltd. 2014
Received: 28 February 2014
Accepted: 28 May 2014
Published: 22 July 2014
Bovine hereditary zinc deficiency (BHZD) is an autosomal recessive disorder of cattle, first described in Holstein-Friesian animals. Affected calves suffer from severe skin lesions and show a poor general health status. Recently, eight calves with the phenotypic appearance of BHZD have been reported in the Fleckvieh cattle population.
In spite of the similar disease phenotypes, SLC39A4, the gene responsible for BHZD in Holstein-Friesian was excluded as underlying gene for the disorder in the affected Fleckvieh calves. In order to identify the disease-associated region, genotypes of eight affected calves obtained with the Illumina BovineHD BeadChip comprising 777,962 SNPs were contrasted with the genotypes of 1,339 unaffected animals. A strong association signal was observed on chromosome 21 (P = 5.87 × 10-89). Autozygosity mapping in the eight affected animals revealed a common segment of extended homozygosity encompassing 1,023 kb (BTA 21: 70,550,045 - 71,573,501). This region contains 17 genes/transcripts, among them two genes encoding gastro-intestinal zinc transporters (CRIP1, CRIP2). However, no mutation that was compatible with recessive inheritance could be detected in these candidate genes. One of the affected calves was re-sequenced together with 42 unaffected Fleckvieh animals. Analysis of the sequencing data revealed a nonsense mutation (p.W215X) in a phospholipase encoding gene (PLD4) as candidate causal polymorphism. To confirm the causality, genotypes of the p.W215X-mutation were obtained from 3,650 animals representing three different breeds. None of the unaffected animals was homozygous for the defect allele, while all eight affected calves were homozygous. The deleterious effect of the mutation is manifested in a significantly lower survival rate of descendants from risk matings when compared with the survival rate of descendants from non-risk matings. The deleterious allele has an estimated frequency of 1.1% in the Fleckvieh population.
Our results provide strong evidence that a newly identified recessive disorder in the Fleckvieh population is caused by a nonsense mutation in PLD4, most likely resulting in an impaired function of the encoded protein. Although the phenotype of affected calves strongly resembles BHZD, a zinc deficiency resulting from malabsorption is unlikely to be responsible for the diseased Fleckvieh calves.
KeywordsAutozygosity mapping Next-generation sequencing Recessive inheritance Zinc deficiency-like Fleckvieh Nonsense mutation PLD4
So far, 415 congenital defects have been identified in cattle, 109 of them being inherited in an autosomal recessive mode . The widespread use of elite sires by means of artificial insemination makes cattle populations highly susceptible to the propagation of recessively inherited disorders . However, the availability of dense SNP genotyping arrays facilitates the rapid identification of underlying genomic regions .
This paper details a disorder that resembles Bovine Hereditary Zinc Deficiency (BHZD, OMIA 000593-9913), an autosomal recessive condition, primarily seen in Holstein-Friesian calves . This disease is characterized by an impaired function of the immune system, growth retardation and skin alterations as a result of a deficient gastrointestinal zinc absorption . While affected animals are born without apparent clinical symptoms, first skin lesions emerge between the age of four and eight weeks . Impairment of immune functions makes affected animals more susceptible to infectious pathogens and leads to an increased incidence of common calf diseases, e.g., enteritis and pneumonia . Highly-dosed oral zinc supplementation ameliorates clinical symptoms in affected Holstein-Friesian animals, however, if untreated, BHZD is lethal . Inherited zinc absorption disorders, caused by mutations in the zinc transporter encoding gene SLC39A4, have been reported for several mammalian species, including mouse and human [7–9]. In Holstein-Friesian cattle, BHZD is caused by a splice-site variant in SLC39A4. Hereditary zinc absorption disorders have not yet been described in the Fleckvieh population.
Here we show that a phenotype of Fleckvieh calves that resembles BHZD is not caused by a defective zinc transporter gene. Genome-wide association analysis, autozygosity mapping and analysis of whole-genome sequencing data enabled to identify a nonsense mutation in the phospholipase D family member 4 encoding gene (PLD4) which is likely causal for the observed phenotype.
Phenotypic manifestation of the defect
In addition to the seven cases described above, another Fleckvieh calf (eleven weeks old) with severe alterations of the skin and diarrhoea was reported by a veterinarian. The described lesions were very similar to the phenotype of the previously reported seven calves. The calf received dietary zinc (750 mg/day) for a period of two weeks. However, the overall condition of the animal did not improve and the calf died at the age of 13 weeks.
Inspection of the pedigrees of the affected calves revealed a common ancestor suggesting a genetic background. On the basis of the patients’ history, the clinical and pathological findings, the calves were tentatively diagnosed to suffer from BHZD.
Analysis of SLC39A4 –the gene causing BHZD in Holstein-Friesian
Mutations in SLC39A4 are known to cause defects resembling the phenotypic appearance of the eight affected Fleckvieh calves in various species including cattle . SLC39A4 is located at the proximal region of bovine chromosome 14 (BTA 14: 1,719,732 bp – 1,724,221 bp). The gene was re-sequenced in a case–control panel consisting of all affected animals, all available dams and sires and randomly selected, unaffected control animals. Totally ~7 kb of genomic sequence was screened, resulting in the detection of ten SNPs (Additional file 1). The mutation causing BHZD in Holstein-Friesian was not present in the diseased animals and none of the detected polymorphisms was associated with the disease phenotype, nor was any of the polymorphic sites compatible with the supposed pattern of recessive inheritance.
Identification of the disease-associated region
Autozygosity mapping within the distal region of BTA 21 revealed a common 1,023 kb segment of extended homozygosity in the eight affected animals (70,550,045 bp – 71,573,501 bp) (Figure 3B). However, one out of 1,339 animals of the control group was homozygous for the disease-associated region as well. Signal intensities obtained from high-density genotyping revealed no indication for the presence of large structural variants (deletions, copy number variations) within the associated region (Additional file 2).
The frequency of the associated haplotype was estimated in a sample of 10,355 unaffected Fleckvieh animals. Among them, 380 were heterozygous carriers and three were homozygous, yielding a frequency for the associated haplotype of 1.86%. The haplotype distribution shows no deviation from the Hardy-Weinberg equilibrium (P = 0.75).
Analysis of CRIP1 and CRIP2 – genes with similar function as SLC39A4
The segment of extended homozygosity contains 17 transcripts and protein encoding genes, respectively (UMD3.1 annotation , Figure 3C). Two of these genes encode cysteine-rich proteins (CRIP1 and CRIP2). CRIPs are highly expressed in the intestine  and are essential for zinc absorption as well as for immune system related cytokine regulation [12, 13]. Thus, they represent excellent candidate genes for the observed disease phenotype.
CRIP1 and CRIP2 (71,390,596 bp – 71,392,110 bp and 71,365,858 bp – 71,382,666 bp, respectively) were re-sequenced in the case–control panel. We screened ~6 kb of the genomic sequence of CRIP1, resulting in the detection of 18 SNPs and two insertion/deletion (InDel) polymorphisms (Additional file 1). Approximately 8 kb of the genomic sequence of CRIP2 were screened, resulting in the identification of 30 SNPs (Additional file 1). None of the variants in CRIP1 and CRIP2 was compatible with the disease phenotype. Thus, variation in the two positional and functional candidate genes is likely not causal for the observed disease.
Identification of the underlying mutation by exploiting whole-genome sequencing data
Genotypes of four variants compatible with recessive inheritance
NCBI reference ID
Validation of the p.W215X-mutation in PLD4
Genotypes for the p.W215X-mutation in PLD4 for 3,650 animals
Genotypes for p.W215X
Dams of the calves
Unaffected animals being homozygous for the disease-associated haplotype
Animals carrying the disease-associated haplotype
Animals not carrying the disease-associated haplotype
Animals with unknown haplotype status
Randomly selected animals
Randomly selected animals
Survival rate of progeny descending from bulls carrying the p.W215X-mutation
Survival rate of calves descending from different mating types
Carrier-state for the p.W215X-mutation
Number of calves
Survival rate at different days
Eight calves with severe skin lesions and poor general health status were recently identified in the Fleckvieh cattle population. Clinical and pathological examinations of affected animals showed striking similarities to findings described for individuals suffering from zinc deficiency in cattle [4, 6, 18] and humans . Inherited zinc absorption disorders, resulting from mutations in SLC39A4, can be ameliorated by highly-dosed oral zinc supplementation . However the Fleckvieh calves did not respond to oral zinc supplementation. Variants in SLC39A4 could be excluded from being causal for the phenotype of the affected Fleckvieh animals, suggesting a different aetiological basis of the syndrome observed in the Fleckvieh calves.
Genome-wide association analysis followed by autozygosity mapping identified a common 1,023 kb segment of extended homozygosity on BTA 21 in the affected calves. Variants in two genes (CRIP1 and CRIP2) with functional similarity to SLC39A4[20, 21] were excluded to be associated with the phenotype of the affected calves. Finally, genome-wide re-sequencing of one of the affected calves revealed a putatively causal loss-of-function mutation in the phospholipase D family member 4 encoding gene PLD4. Given the relatively low fold sequence coverage, some variants might have been wrongly called and lost, respectively. However, multi-sample variant calling followed by population-based genotype imputation as applied in the present study provides high quality genotypes even for lowly covered sites . The p.W125X-mutation was perfectly associated in a panel of 3,650 animals, further substantiating our suspicion of causality.
Phospholipase D4 is a member of the family of phospholipid signalling enzymes  and is mainly expressed in spleen and early microglia, suggesting a role in immunological pathways [22, 23]. Although PLD4 seems to be involved in phagocytosis, little is known about its enzymatic function in cells . Recent studies revealed an association between variants in PLD4 and two autoimmune diseases in human, namely systemic sclerosis  and rheumatoid arthritis , disorders with inflammatory skin lesions  and interstitial lung diseases [26, 27]. Furthermore, knock-out mice with a nonsense mutation in PLD4 manifest a phenotype with thin hair and growth retardation . In addition, Arun et al. demonstrated the importance of phospholipase family D members in mediating the repair of plasma membrane disruptions in mice keratinocytes emphasising the role of these phospholipid signalling enzymes in membrane function and wound healing. A mutation in a phospholipase-domain containing protein (PNPLA1) results in a severe cornification disorder in dogs, providing evidence for a key role of lipases in the keratinisation process and the metabolism of the epidermal barrier in general . Taken together, it seems very likely that an impaired function of PLD4 is causal for the severe skin lesions and the poor general health status of the affected calves.
Compared to the wild-type protein, the PLD4 protein is predicted to be shortened by 273 amino acids in the affected calves. The truncated protein may be retained with an impaired function and/or the transcript may be degraded via nonsense-mediated mRNA decay . If the truncated protein is retained, its function is likely to be severely compromised as it lacks essential domains for enzymatic activity [16, 17]. However, the actual effect of the mutation on the expression of PLD4 needs to be unravelled in subsequent studies.
Although there is a striking similarity in the phenotypic appearance of the Fleckvieh calves and Holstein-Friesian animals suffering from BHZD , there are no clues for a connection between PLD4 and zinc metabolism. Furthermore, zinc supplementation did not ameliorate any of the symptoms of the affected calves, corroborating that impaired zinc metabolism might not be causal for the phenotype. Based on our findings the tentative diagnosis of zinc deficiency can no longer be maintained.
Haplotype analysis revealed that three unaffected animals are homozygous for the associated haplotype. However, none of them was homozygous for the p.W215X-mutation. It seems likely that the mutation might have occurred in the germline of a recent founder animal, resulting in two identical haplotypes differing for the deleterious allele only. A similar situation has been observed for the arachnomelia syndrome (OMIA 000059-9913) in Brown Swiss cattle . Using array-derived genotypes does not allow distinguishing between animals carrying the affected and unaffected haplotype version. Thus, the haplotype frequency derived from array-based genotypes (1.86%) is higher than the frequency of the p.W125X-mutation. The frequency of the p.W215X-mutation was estimated to 1.1% only in the current Fleckvieh population. However, the widespread use of a single bull carrying the defect allele could lead to a rapid increase of the frequency within few generations. Since array-based genotypes are routinely obtained for all candidate bulls, animals carrying the disease-associated haplotype can now be easily identified. Such animals should be directly genotyped for the p.W215X-mutation and animals carrying the mutation should be excluded from artificial insemination. This cost-effective approach will prevent unneeded animal suffering and economic losses by avoiding inadvertent carrier x carrier matings.
Assuming a frequency of 1.1% of the deleterious allele in the Fleckvieh population, equal use of all bulls and ~1,000,000 annual births, one would expect 121 affected calves per year. The actual number of expected cases is probably somewhat lower, as most carriers are related and farmers avoid close inbreeding. Furthermore, most of these cases will probably not be reported since the affected calves die due to rather unspecific diseases (e.g., enteritis or diarrhoea) before the manifestation of characteristic symptoms (i.e., hyperkeratotic and fissured skin). The significantly reduced survival rate of descendants from risk matings supports this assumption. However, incomplete penetrance could also reduce the observed incidences. Therefore it will be important to study animals being homozygous for the identified mutation from birth in a controlled environment for a better characterisation of the disorder.
A recessively inherited condition in Fleckvieh cattle resembles the phenotypic appearance of bovine hereditary zinc deficiency. Our results strongly support that a nonsense mutation (p.W215X) in the PLD4 encoding gene is causative for this disease in the Fleckvieh breed and that zinc deficiency is not involved in the aetiology of the disease. Further studies are necessary to unravel the detailed genotype-phenotype relationship. However, the identification of the causal variant allows for broad testing in the Fleckvieh population. Hence, an efficient management of this new defect is now possible.
Animal ethics statement
Semen samples were collected by approved commercial artificial insemination stations as part of their regular breeding and reproduction measures in cattle industry. The collection of blood samples was carried out by trained veterinarians during treatment of affected animals following standard veterinary protocols in Germany. No ethical approval was required for this study.
Animals and DNA extraction
In order to confirm a genetic predisposition to the observed phenotype, a case–control panel consisting of the eight affected calves, their dams and sires and twelve unaffected control animals was set up. Seven affected animals were identified and examined by veterinarians of the Clinic for Ruminants. One affected calf was reported by a herd veterinarian. Blood samples of the affected calves and their dams were collected by trained veterinarians following standard procedures and DNA was extracted using proteinase K digestion and salt-out extraction. For the control group, genomic DNA was prepared from semen straws following standard protocols using proteinase K digestion and phenol-chloroform extraction.
Annotation and polymorphism screening of SLC39A4, CRIP1 and CRIP2
The GENOMETHREADER software tool  was used to predict the genomic structure and localization of SLC39A4, CRIP1 and CRIP2 based on the University of Maryland UMD3.1 assembly of the bovine genome sequence  and the Dana–Farber Cancer Institute bovine gene index release 12.0  together with the annotated RNA sequences of the UMD3.1 assembly . The GENOMETHREADER output was viewed and edited using the Apollo sequence annotation editor . The genes were PCR amplified (the primers are listed in Additional file 5), including exons, introns and the flanking 3'- and 5'-regions. Sequencing reactions were done using BigDye® Terminator v1.1 Cycle Sequencing Kit (Applied Biosystems, Foster, CA., USA; Life Technologies Corporation, USA). Electrophoresis of purified sequencing reactions was performed on the ABI 3130x1 Genetic Analyzer (Applied Biosystems, Foster, CA., USA; Life Technologies Corporation, USA). The Phred/Phrap/Polyphred software suite [36–38] was used for base calling, sequence alignment and polymorphism detection. Sequences were viewed with consed.
High- density genotyping and quality control
The eight affected calves and 1,339 unaffected Fleckvieh bulls were genotyped with the Illumina BovineHD BeadChip comprising 777,962 SNPs. Genotype calling was performed using default parameters of Illumina's BeadStudio. The chromosomal position of the SNPs was determined based on the UMD3.1 assembly of the bovine genome . 343 mitochondrial, 1,224 Y-chromosomal and 1,735 SNPs with unknown chromosome position were not considered for subsequent analyses. Quality control was carried out with PLINK v1.07 . 19,734 SNPs for which genotyping failed in more than 5% of the individuals and 110,746 SNPs with minor allele frequency < 0.5% were omitted. Sporadically missing genotypes were imputed using Beagle genetic analysis software . The final dataset comprised 1,347 animals (eight affected calves and 1,339 controls) and 644,450 SNPs.
Genome-wide association study
To account for population stratification and the resulting inflation of false positive associations, a mixed model based association analysis was performed. We used GEMMA to fit the linear mixed model y = μ + xb + Zu + e, where y denotes the affection status (coded as 1 and 2 for affected and unaffected animals, respectively), μ is the intercept, x is a vector of marker genotypes, b is the SNP effect, Z is an incidence matrix, u is a vector of random polygenic effects ~ N(0, G), where is the additive genetic variance and G is the genomic relationship matrix (GRM) among the 1,347 animals built based upon 627,627 autosomal SNPs following VanRaden's approach  and e is a vector of error terms.
Haplotypes of 1,347 animals of the initial genome-wide scan were inferred using default parameters of Beagle genetic analysis software (see above). To assess the frequency of the associated haplotype in a larger sample of the Fleckvieh population, array-derived genotypes of another 9,016 unaffected animals were analysed. Of the 9,016 animals, 7,000 were genotyped with the BovineSNP50 BeadChip (50 K) and 2,016 were genotyped with the Illumina BovineHD BeadChip (777 K). High-density genotypes and haplotypes of the 50 K data set were inferred based on haplotypes of the 777 K data set using a combination of Beagle (see above) and Minimac, which yields high imputation accuracy in cattle .
Exploiting whole-genome re-sequencing data for mutation screening
The genomes of 43 animals of the Fleckvieh population were sequenced to an average coverage of 7.46x, among them an affected calf and a healthy animal being homozygous for the disease-associated haplotype with coverages of 7.8x and 6.2x, respectively. Sequencing on Illumina GA IIX and Hiseq 2000 instruments was performed as detailed by Jansen et al.. Paired-end reads were obtained and mapped to the bovine reference sequence (see above) using the Burrows-Wheeler Aligner (BWA) . PICARD (http://picard.sourceforge.net) was used to mark PCR-duplicates. Subsequent multi-sample variant calling with mpileup yielded genotypes at 17.17 million sites. Beagle phasing and imputation was applied to improve the primary genotype calls. A detailed overview of the entire variant calling pipeline and all obtained variants is presented in Jansen et al.. Of 17.17 million sites, 7,086 SNPs and 574 InDels were located within the 1,023 kb segment (70,550,045 bp – 71,573,501 bp) of extended homozygosity on BTA 21. To account for inaccurately genotyped variants due to the low-coverage sequence data (e.g., mis-calling of heterozygous genotypes for rare variants [49, 50]), we filtered for variants that were segregating (heterozygous or homozygous for the non-reference allele) in the affected animal, heterozygous or homozygous for the reference allele in the healthy animal being homozygous for the disease-associated haplotype and homozygous for the reference allele in 41 healthy animals.
5'-exonuclease diagnostic assay of the c.G645A-PLD4mutation
A 5'-exonuclease assay was developed to obtain genotypes for the nonsense mutation (c.G645A, p.W215X, Chr21:71,001,232 bp, rs378824791), using 5'-GGG CGG CAC ATC TAC GT-3' and 5'-CCA GGG CGG ACG AAC TC-3' as PCR primers, and 5'-CAT GGA CTG GCG GTC C-3' (wild type G allele) and 5'-CAT GGA CTG ACG GTC C-3' (mutant A allele) as probes (TaqMan®, Life Technologies Corporation, USA). Reactions were carried out on an ABI7500 Real-Time PCR system (Life Technologies Corporation, USA) using standard procedures and analysed using the allelic discrimination endpoint analysis mode of the 7500 software package v2.0.5. Genotypes for the polymorphism were obtained in 3,650 animals representing three different breeds (Fleckvieh, Holstein-Friesian, Braunvieh).
The survival rate of calves from two different mating types was estimated using a Kaplan-Meier estimator as implemented in the R package ('survival') . For the control group the survival rate of 2,552 calves descending from matings where the sire does not carry the p.W215X-mutation whereas the maternal grandsire is a carrier S5(non-risk matings) was estimated. The case group comprised 1,213 calves descending from matings where both sire and maternal grandsire carry the p.W215X-mutation (risk-matings).
Availability of supporting data
All relevant SNPs (with accession numbers of dbSNP) supporting the results of this article are included within the article and its additional files.
We thank the Arbeitsgemeinschaft Süddeutscher Rinderzüchter und Besamungsorganisation e.V (ASR) and the Förderverein Biotechnologieforschung e.V. (FBF) for sharing genotype data and the Bavarian State Research Center for Agriculture for providing additional genotypes and DNA samples. The study was supported by the German Federal Ministry of Education and Research (BMBF) within the AgroClustEr “Synbreed – Synergistic plant and animal breeding” (FKZ: 0315528A). We further acknowledge financial support from the Förderverein Biotechnologieforschung e.V. (FBF) and the Universität Bayern e.V.
- Online Mendelian Inheritance in Animals, OMIA, Faculty of Veterinary Science, University of Sydney (January 2014). http://omia.angis.org.au/,
- Van Vleck LD: Effect of artificial insemination on frequency of undesirable recessive genes. J Dairy Sci. 1967, 50: 201-204. 10.3168/jds.S0022-0302(67)87388-0.PubMedView ArticleGoogle Scholar
- Charlier C, Coppieters W, Rollin F, Desmecht D, Agerholm JS, Cambisano N, Carta E, Dardano S, Dive M, Fasquelle C, Frennet J-C, Hanset R, Hubin X, Jorgensen C, Karim L, Kent M, Harvey K, Pearce BR, Simon P, Tama N, Nie H, Vandeputte S, Lien S, Longeri M, Fredholm M, Harvey RJ, Georges M: Highly effective SNP-based association mapping and management of recessive defects in livestock. Nat Genet. 2008, 40: 449-454. 10.1038/ng.96.PubMedView ArticleGoogle Scholar
- Machen M, Montgomery T, Holland R, Braselton E, Dunstan R, Brewer G, Yuzbasiyan-Gurkan V: Bovine hereditary zinc deficiency: lethal trait A 46. J Vet Diagn Invest. 1996, 8: 219-227. 10.1177/104063879600800212.PubMedView ArticleGoogle Scholar
- Brummerstedt E: Animal model of human disease: acrodermatitis enteropathica, zinc malabsorption. Am J Pathol. 1977, 87: 725-728.PubMed CentralPubMedGoogle Scholar
- Stöber M: Parakeratosis in black pied calves: 1: clinical findings and etiology. DTW Dtsch Tierarztl Wochenschr. 1971, 78: 257-265.PubMedGoogle Scholar
- Wang K, Zhou B, Kuo Y-M, Zemansky J, Gitschier J: A novel member of a zinc transporter family is defective in acrodermatitis enteropathica. Am J Hum Genet. 2002, 71: 66-73. 10.1086/341125.PubMed CentralPubMedView ArticleGoogle Scholar
- Yuzbasiyan-Gurkan V, Bartlett E: Identification of a unique splice site variant in SLC39A4 in bovine hereditary zinc deficiency, lethal trait A46: an animal model of acrodermatitis enteropathica. Genomics. 2006, 88: 521-526. 10.1016/j.ygeno.2006.03.018.PubMedView ArticleGoogle Scholar
- Dufner-Beattie J, Weaver BP, Geiser J, Bilgen M, Larson M, Xu W, Andrews GK: The mouse acrodermatitis enteropathica gene Slc39a4 (Zip4) is essential for early development and heterozygosity causes hypersensitivity to zinc deficiency. Hum Mol Genet. 2007, 16: 1391-1399. 10.1093/hmg/ddm088.PubMedView ArticleGoogle Scholar
- Florea L, Souvorov A, Kalbfleisch TS, Salzberg SL: Genome assembly has a major impact on gene content: a comparison of annotation in two bos taurus assemblies. PLoS One. 2011, 6: e21400-10.1371/journal.pone.0021400.PubMed CentralPubMedView ArticleGoogle Scholar
- Hallquist NA, Khoo C, Cousins RJ: Lipopolysaccharide regulates cysteine-rich intestinal protein, a zinc-finger protein, in immune cells and plasma. J Leukoc Biol. 1996, 59: 172-177.PubMedGoogle Scholar
- Hempe JM, Cousins RJ: Cysteine-rich intestinal protein binds zinc during transmucosal zinc transport. Proc Natl Acad Sci USA. 1991, 88: 9671-9674. 10.1073/pnas.88.21.9671.PubMed CentralPubMedView ArticleGoogle Scholar
- Cousins RJ, Lanningham-Foster L: Regulation of cysteine-rich intestinal protein, a zinc finger protein, by mediators of the immune response. J Infect Dis. 2000, 182 (Suppl 1): S81-84.PubMedView ArticleGoogle Scholar
- Jansen S, Aigner B, Pausch H, Wysocki M, Eck S, Benet-Pagès A, Graf E, Wieland T, Strom TM, Meitinger T, Fries R: Assessment of the genomic variation in a cattle population by re-sequencing of key animals at low to medium coverage. BMC Genomics. 2013, 14: 446-10.1186/1471-2164-14-446.PubMed CentralPubMedView ArticleGoogle Scholar
- Daetwyler HD, Capitan A, Pausch H, Stothard P, van Binsbergen R, Brøndum RF, Liao X, Djari A, Rodriguez SC, Grohs C, Esquerré D, Bouchez O, Rossignol M-N, Klopp C, Rocha D, Fritz S, Eggen A, Bowman PJ, Coote D, Chamberlain AJ, Anderson C, VanTassell CP, Hulsegge I, Goddard ME, Guldbrandtsen B, Lund MS, Veerkamp RF, Boichard DA, Fries R, Hayes BJ: Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. Nat Genet. 2014, advance online publication: doi:10.1038/ng.3034Google Scholar
- Xie Z, Ho WT, Exton JH: Association of the N- and C-terminal domains of phospholipase D: contribution of the conserved HKD motifs to the interaction and the requirement of the association for Ser/Thr phosphorylation of the enzyme. J Biol Chem. 2000, 275: 24962-24969. 10.1074/jbc.M909745199.PubMedView ArticleGoogle Scholar
- Yoshikawa F, Banno Y, Otani Y, Yamaguchi Y, Nagakura-Takagi Y, Morita N, Sato Y, Saruta C, Nishibe H, Sadakata T, Shinoda Y, Hayashi K, Mishima Y, Baba H, Furuichi T: Phospholipase D family member 4, a transmembrane glycoprotein with no phospholipase D activity, expression in spleen and early postnatal microglia. PLoS One. 2010, 5: e13932-10.1371/journal.pone.0013932.PubMed CentralPubMedView ArticleGoogle Scholar
- Trautwein G: Parakeratosis in black pied calves: 2: pathological findings. DTW Dtsch Tierarztl Wochenschr. 1971, 78: 265-270.PubMedGoogle Scholar
- Maverakis E, Lynch PJ, Fazel N: Acrodermatitis enteropathica. Dermatol Online J. 2007, 13: 11-PubMedGoogle Scholar
- O’Dell BL: Cysteine-rich intestinal protein (CRIP): a new intestinal zinc transport protein. Nutr Rev. 1992, 50: 232-233.PubMedView ArticleGoogle Scholar
- Cousins RJ: Gastrointestinal factors influencing zinc absorption and homeostasis. Int J Vitam Nutr Res. 2010, 80: 243-248. 10.1024/0300-9831/a000030.PubMed CentralPubMedView ArticleGoogle Scholar
- Otani Y, Yamaguchi Y, Sato Y, Furuichi T, Ikenaka K, Kitani H, Baba H: PLD4 is involved in phagocytosis of microglia: expression and localization changes of PLD4 are correlated with activation state of microglia. PLoS One. 2011, 6: e27544-10.1371/journal.pone.0027544.PubMed CentralPubMedView ArticleGoogle Scholar
- Terao C, Ohmura K, Kawaguchi Y, Nishimoto T, Kawasaki A, Takehara K, Furukawa H, Kochi Y, Ota Y, Ikari K, Sato S, Tohma S, Yamada R, Yamamoto K, Kubo M, Yamanaka H, Kuwana M, Tsuchiya N, Matsuda F, Mimori T: PLD4 as a novel susceptibility gene for systemic sclerosis in a Japanese population. Arthritis & Rheumatism. 2013, 65: 472-480. 10.1002/art.37777.View ArticleGoogle Scholar
- Okada Y, Terao C, Ikari K, Kochi Y, Ohmura K, Suzuki A, Kawaguchi T, Stahl EA, Kurreeman FAS, Nishida N, Ohmiya H, Myouzen K, Takahashi M, Sawada T, Nishioka Y, Yukioka M, Matsubara T, Wakitani S, Teshima R, Tohma S, Takasugi K, Shimada K, Murasawa A, Honjo S, Matsuo K, Tanaka H, Tajima K, Suzuki T, Iwamoto T, Kawamura Y, et al: Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population. Nat Genet. 2012, 44: 511-516. 10.1038/ng.2231.PubMedView ArticleGoogle Scholar
- Krieg T, Takehara K: Skin disease: a cardinal feature of systemic sclerosis. Rheumatology (Oxford). 2009, 3: iii14-18.Google Scholar
- Picchianti Diamanti A, Germano V, Bizzi E, Laganà B, Migliore A: Interstitial lung disease in rheumatoid arthritis in the era of biologics. Pulm Med. 2011, 2011: 931342-PubMed CentralPubMedView ArticleGoogle Scholar
- Bussone G, Mouthon L: Interstitial lung disease in systemic sclerosis. Autoimmun Rev. 2011, 10: 248-255. 10.1016/j.autrev.2010.09.012.PubMedView ArticleGoogle Scholar
- Harris B, Ward-Bailey PF, Bergstrom DE, Bronson RT, Donahue LR: Thin hair with small size (thss) is a new recessive hair mutation on Chromosom 12. Bar Harbor, Maine: Mouse Mutant Resource Web Site, The Jackson Laboratory, (http://mousemutant.jax.org/). [September 17, 2012]
- Arun SN, Xie D, Howard AC, Zhong Q, Zhong X, McNeil PL, Bollag WB: Cell wounding activates phospholipase D in primary mouse keratinocytes. J Lipid Res. 2013, 54: 581-591. 10.1194/jlr.M027060.PubMed CentralPubMedView ArticleGoogle Scholar
- Grall A, Guaguère E, Planchais S, Grond S, Bourrat E, Hausser I, Hitte C, Le Gallo M, Derbois C, Kim G-J, Lagoutte L, Degorce-Rubiales F, Radner FPW, Thomas A, Küry S, Bensignor E, Fontaine J, Pin D, Zimmermann R, Zechner R, Lathrop M, Galibert F, André C, Fischer J: PNPLA1 mutations cause autosomal recessive congenital ichthyosis in golden retriever dogs and humans. Nat Genet. 2012, 44: 140-147. 10.1038/ng.1056.PubMedView ArticleGoogle Scholar
- Frischmeyer PA, Dietz HC: Nonsense-mediated mRNA decay in health and disease. Hum Mol Genet. 1999, 8: 1893-1900. 10.1093/hmg/8.10.1893.PubMedView ArticleGoogle Scholar
- Drögemüller C, Tetens J, Sigurdsson S, Gentile A, Testoni S, Lindblad-Toh K, Leeb T: Identification of the Bovine Arachnomelia Mutation by Massively Parallel Sequencing Implicates Sulfite Oxidase (SUOX) in bone development. PLoS Genet. 2010, 6: e1001079-10.1371/journal.pgen.1001079.PubMed CentralPubMedView ArticleGoogle Scholar
- Gremme G, Brendel V, Sparks ME, Kurtz S: Engineering a software tool for gene structure prediction in higher organisms. Inform Software Tech. 2005, 47: 965-978. 10.1016/j.infsof.2005.09.005.View ArticleGoogle Scholar
- Quackenbush J, Cho J, Lee D, Liang F, Holt I, Karamycheva S, Parvizi B, Pertea G, Sultana R, White J: The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species. Nucl Acids Res. 2001, 29: 159-164. 10.1093/nar/29.1.159.PubMed CentralPubMedView ArticleGoogle Scholar
- Lewis SE, Searle SMJ, Harris N, Gibson M, Iyer V, Richter J, Wiel C, Bayraktaroglu L, Birney E, Crosby MA, Kaminker JS, Matthews BB, Prochnik SE, Smith CD, Tupy JL, Rubin GM, Misra S, Mungall CJ, Clamp ME: Apollo: a sequence annotation editor. Genome Biol. 2002, 3: research0082-1 – 0082.14PubMed CentralPubMedView ArticleGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred: II: error probabilities. Genome Res. 1998, 8: 186-194.PubMedView ArticleGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred: I: accuracy assessment. Genome Res. 1998, 8: 175-185. 10.1101/gr.8.3.175.PubMedView ArticleGoogle Scholar
- Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res. 1997, 25: 2745-2751. 10.1093/nar/25.14.2745.PubMed CentralPubMedView ArticleGoogle Scholar
- Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202. 10.1101/gr.8.3.195.PubMedView ArticleGoogle Scholar
- Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D, Hanrahan F, Pertea G, Van Tassell CP, Sonstegard TS, Marçais G, Roberts M, Subramanian P, Yorke JA, Salzberg SL: A whole-genome assembly of the domestic cow: bos taurus. Genome Biol. 2009, 10: R42-10.1186/gb-2009-10-4-r42.PubMed CentralPubMedView ArticleGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.PubMed CentralPubMedView ArticleGoogle Scholar
- Browning BL, Browning SR: A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009, 84: 210-223. 10.1016/j.ajhg.2009.01.005.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhou X, Stephens M: Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012, 44: 821-824. 10.1038/ng.2310.PubMed CentralPubMedView ArticleGoogle Scholar
- VanRaden PM: Efficient methods to compute genomic predictions. J Dairy Sci. 2008, 91: 4414-4423. 10.3168/jds.2007-0980.PubMedView ArticleGoogle Scholar
- Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR: Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet. 2012, 44: 955-959. 10.1038/ng.2354.PubMed CentralPubMedView ArticleGoogle Scholar
- Pausch H, Aigner B, Emmerling R, Edel C, Götz K-U, Fries R: Imputation of high-density genotypes in the Fleckvieh cattle population. Genet Sel Evol. 2013, 45: 3-10.1186/1297-9686-45-3.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralPubMedView ArticleGoogle Scholar
- Li Y, Sidore C, Kang HM, Boehnke M, Abecasis GR: Low-coverage sequencing: implications for design of complex trait association studies. Genome Res. 2011, 21: 940-951. 10.1101/gr.117259.110.PubMed CentralPubMedView ArticleGoogle Scholar
- DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, Philippakis AA, del Angel G, Rivas MA, Hanna M, McKenna A, Fennell TJ, Kernytsky AM, Sivachenko AY, Cibulskis K, Gabriel SB, Altshuler D, Daly MJ: A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011, 43: 491-498. 10.1038/ng.806.PubMed CentralPubMedView ArticleGoogle Scholar
- Therneau TM, Grambsch PM: Modeling Survival Data: Extending the Cox Model. 2000, New York: SpringerView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.