A second KRT71 allele in curly coated dogs

Major characteristics of coat variation in dogs can be explained by variants in only a few genes. Until now, only one missense variant in the KRT71 gene, p.Arg151Trp, has been reported to cause curly hair in dogs. However, this variant does not explain the curly coat in all breeds as the mutant ¹⁵¹Trp allele, for example, is absent in Curly Coated Retrievers. We sequenced the genome of a Curly Coated Retriever at 22× coverage and searched for variants in the KRT71 gene. Only one protein‐changing variant was present in a homozygous state in the Curly Coated Retriever and absent or present in a heterozygous state in 221 control dogs from different dog breeds. This variant, NM_001197029.1:c.1266_1273delinsACA, was an indel variant in exon 7 that caused a frameshift and an altered and probably extended C‐terminus of the KRT71 protein NP_001183958.1:p.(Ser422ArgfsTer?). Using Sanger sequencing, we found that the variant was fixed in a cohort of 125 Curly Coated Retrievers and segregating in five of 14 additionally tested breeds with a curly or wavy coat. KRT71 variants cause curly hair in humans, mice, rats, cats and dogs. Specific KRT71 variants were further shown to cause alopecia. Based on this knowledge from other species and the predicted molecular consequence of the newly identified canine KRT71 variant, it is a compelling candidate causing a second curly hair allele in dogs. It might cause a slightly different coat phenotype than the previously published p.Arg151Trp variant and could potentially be associated with follicular dysplasia in dogs.     

ARSB gene variants causing Mucopolysaccharidosis VI in Miniature Pinscher and Miniature Schnauzer dogs

Mucopolysaccharidosis (MPS) VI is a lysosomal storage disease caused by a deficiency of N-acetylgalactosamine-4-sulfatase, also called arylsulfatase B (ARSB, EC 3.1.6.12). Dogs with MPS VI show progressive predominantly oculoskeletal signs homologous to those in human and feline patients. We report herein two pathogenic ARSB gene variants in Miniature Pinscher and Miniature Schnauzer dogs with MPS VI and a genotyping survey in these breeds. All exons and adjacent regions of the ARSB gene were sequenced from three affected Miniature Pinschers and three affected Miniature Schnauzers. Allelic discrimination assays were used for genotyping. A missense variant (NM_001048133.1:c.910G>A) was found in exon 5 of MPS VI-affected Miniature Pinschers that is predicted to result in a deleterious amino acid substitution of a highly conserved glycine to arginine (NP_001041598.1:p.Gly304Arg). In MPS VI-affected Miniature Schnauzers, a 56 bp deletion (NM_001048133.1:c.-24_32del) was found at the junction of exon 1 and its upstream region, predicting no enzyme synthesis. All clinically affected Miniature Pinschers and Miniature Schnauzers were homozygous for the respective variants, and screened healthy dogs in each breed were either heterozygous or homozygous for the wt allele. Whereas the Miniature Pinscher variant seemed to occur commonly (0.133 allele frequency), the Miniature Schnauzer variant was presumed to be rare. In conclusion, two breed-specific pathogenic ARSB gene variants were identified in Miniature Pinscher and Miniature Schnauzer dogs with MPS VI, allowing for genotyping and informed breeding to prevent the production of affected offspring.     

Identification of compound heterozygous DNAH11 variants in a Han-Chinese family with primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) is a group of genetically and clinically heterogeneous disorders with motile cilia dysfunction. It is clinically characterized by oto-sino-pulmonary diseases and subfertility, and half of the patients have situs inversus (Kartagener syndrome). To identify the genetic cause in a Han-Chinese pedigree, whole-exome sequencing was conducted in the 37-year-old proband, and then, Sanger sequencing was performed on available family members. Minigene splicing assay was applied to verify the impact of the splice-site variant. Compound heterozygous variants including a splice-site variant (c.1974-1G>C, rs1359107415) and a missense variant (c.7787G>A, p.(Arg2596Gln), rs780492669), in the dynein axonemal heavy chain 11 gene (DNAH11) were identified and confirmed as the disease-associated variants of this lineage. The minigene expression in vitro revealed that the c.1974-1G>C variant could cause skipping over exon 12, predicted to result in a truncated protein. This discovery may enlarge the DNAH11 variant spectrum of PCD, promote accurate genetic counselling and contribute to PCD diagnosis.     

A novel KRT71 variant in curly‐coated dogs

Curly fur is a common phenotype in many dog breeds, known to result from a missense variant (c.451C>T) in exon 2 of the keratin 71 (KRT71) gene. During screening for this variant across various breeds, we found that Curly Coated Retrievers (CCRs) fixed with the trait did not carry the known variant. By analysis of whole‐genome sequencing data of one CCR we identified a novel genetic cause for curly fur. We found a novel structural variant in exon 7 of the KRT71 gene (c.1266_1273delinsACA) that was predicted to result in a frameshift and stop loss, therefore significantly affecting the structure of the protein, if translated. The variant was also found at lower frequencies in five other breeds, including Lagotto Romagnolo, Bichon Frise, Spanish Water Dog, Chesapeake Bay Retriever and Irish Terrier. One curly‐coated Lagotto carried neither of the two KRT71 variants. These results identify a second variant for curly coat in KRT71 and suggest the existence of additional alleles. This study enables the development of an additional KRT71 gene test for breeders to understand and manage coat types.     

Utilizing the Dog Genome in the Search for Novel Candidate Genes Involved in Glioma Development—Genome Wide Association Mapping followed by Targeted Massive Parallel Sequencing Identifies a Strongly Associated Locus

Gliomas are the most common form of malignant primary brain tumors in humans and second most common in dogs, occurring with similar frequencies in both species. Dogs are valuable spontaneous models of human complex diseases including cancers and may provide insight into disease susceptibility and oncogenesis. Several brachycephalic breeds such as Boxer, Bulldog and Boston Terrier have an elevated risk of developing glioma, but others, including Pug and Pekingese, are not at higher risk. To identify glioma-associated genetic susceptibility factors, an across-breed genome-wide association study (GWAS) was performed on 39 dog glioma cases and 141 controls from 25 dog breeds, identifying a genome-wide significant locus on canine chromosome (CFA) 26 (p = 2.8 x 10⁻⁸). Targeted re-sequencing of the 3.4 Mb candidate region was performed, followed by genotyping of the 56 SNVs that best fit the association pattern between the re-sequenced cases and controls. We identified three candidate genes that were highly associated with glioma susceptibility: CAMKK2, P2RX7 and DENR. CAMKK2 showed reduced expression in both canine and human brain tumors, and a non-synonymous variant in P2RX7, previously demonstrated to have a 50% decrease in receptor function, was also associated with disease. Thus, one or more of these genes appear to affect glioma susceptibility.     

A missense variant in the NSDHL gene in a Chihuahua with a congenital cornification disorder resembling inflammatory linear verrucous epidermal nevi

Congenital hemidysplasia with ichthyosiform nevus and limb defects syndrome in humans is a genodermatosis characterized by inflammatory linear verrucous epidermal nevi (ILVEN), often showing a striking lateralization pattern. It is caused by variants in the NSDHL gene encoding a 3β‐hydroxysteroid dehydrogenase involved in the cholesterol biosynthesis pathway. In the present study, we investigated a female Chihuahua, which showed clinical and histological signs of ILVEN. We performed a candidate gene analysis in the affected animal. This analysis revealed a single missense variant in the NSDHL gene in the affected dog (XM_014111859.2:c.700G>A). The variant is predicted to cause a non‐conservative amino acid change from glycine to arginine, XP_013967334.1:p.(Gly234Arg). The mutant allele was absent from WGS data of 594 genetically diverse dogs and eight wolves. Sanger sequencing confirmed that the variant was heterozygous in the affected dog and absent from 22 control Chihuahuas. Based on the knowledge about the functional impact of NSDHL variants in dogs and other species, c.700G>A is probably pathogenic and a convincing candidate causative variant for the observed skin lesions in the affected Chihuahua.     

Two Independent Mutations in ADAMTS17 Are Associated with Primary Open Angle Glaucoma in the Basset Hound and Basset Fauve de Bretagne Breeds of Dog

PurposeMutations in ADAMTS10 (CFA20) have previously been associated with primary open angle glaucoma (POAG) in the Beagle and Norwegian Elkhound. The closely related gene, ADAMTS17, has also been associated with several different ocular phenotypes in multiple breeds of dog, including primary lens luxation and POAG. We investigated ADAMTS17 as a candidate gene for POAG in the Basset Hound and Basset Fauve de Bretagne dog breeds.MethodsWe performed ADAMTS17 exon resequencing in three Basset Hounds and three Basset Fauve de Bretagne dogs with POAG. Identified variants were genotyped in additional sample cohorts of both breeds and dogs of other breeds to confirm their association with disease.ResultsAll affected Basset Hounds were homozygous for a 19 bp deletion in exon 2 that alters the reading frame and is predicted to lead to a truncated protein. Fifty clinically unaffected Basset Hounds were genotyped for this mutation and all were either heterozygous or homozygous for the wild type allele. Genotyping of 223 Basset Hounds recruited for a different study revealed a mutation frequency of 0.081 and predicted frequency of affected dogs in the population to be 0.007. Based on the entire genotyping dataset the association statistic for the POAG-associated deletion was p = 1.26 x 10⁻¹⁰. All affected Basset Fauve de Bretagne dogs were homozygous for a missense mutation in exon 11 causing a glycine to serine amino acid substitution (G519S) in the disintegrin-like domain of ADAMTS17 which is predicted to alter protein function. Unaffected Basset Fauve de Bretagne dogs were either heterozygous for the mutation (5/24) or homozygous for the wild type allele (19/24). Based on the entire genotyping dataset the association statistic for the POAG-associated deletion was p = 2.80 x 10⁻⁷. Genotyping of 85 dogs of unrelated breeds and 90 dogs of related breeds for this variant was negative.ConclusionThis report documents strong associations between two independent ADAMTS17 mutations and POAG in two different dog breeds.     

Co-Expression of Wild-Type P2X7R with Gln460Arg Variant Alters Receptor Function

The P2X7 receptor is a member of the P2X family of ligand-gated ion channels. A single-nucleotide polymorphism leading to a glutamine (Gln) by arginine (Arg) substitution at codon 460 of the purinergic P2X7 receptor (P2X7R) has been associated with mood disorders. No change in function (loss or gain) has been described for this SNP so far. Here we show that although the P2X7R-Gln460Arg variant per se is not compromised in its function, co-expression of wild-type P2X7R with P2X7R-Gln460Arg impairs receptor function with respect to calcium influx, channel currents and intracellular signaling in vitro. Moreover, co-immunoprecipitation and FRET studies show that the P2X7R-Gln460Arg variant physically interacts with P2X7R-WT. Specific silencing of either the normal or polymorphic variant rescues the heterozygous loss of function phenotype and restores normal function. The described loss of function due to co-expression, unique for mutations in the P2RX7 gene so far, explains the mechanism by which the P2X7R-Gln460Arg variant affects the normal function of the channel and may represent a mechanism of action for other mutations.     

DRD4 and TH gene polymorphisms are associated with activity,impulsivity and inattention in Siberian Husky dogs

Both dopamine receptor D4 (DRD4) exon 3 and tyrosine hydroxylase (TH) intron 4 repeat polymorphisms have been linked to activity and impulsivity in German Shepherd dogs (GSDs). However, the results in GSDs may not be generalisable to other breeds, as allelic frequencies vary markedly among breeds. We selected the Siberian Husky for further study, because it is highly divergent from most dog breeds, including the GSD. The study sample consisted of 145 racing Siberian Huskies from Europe and North America. We found that this breed possesses seven DRD4 length variants, two to five more variants than found in other breeds. Among them was the longest known allele, previously described only in wolves. Short alleles of the DRD4 and TH repeat polymorphisms were associated with higher levels of activity, impulsivity and inattention. Siberian Huskies possessing at least one short allele of the DRD4 polymorphism displayed greater activity in a behavioural test battery than did those with two long alleles. However, the behavioural test was brief and may not have registered variation in behaviour across time and situations. Owners also completed the Dog‐Attention Deficit Hyperactivity Disorder Rating Scale (Dog‐ADHD RS), a more general measure of activity and attention. Siberian Huskies from Europe with two short alleles of the TH polymorphism received higher ratings of inattention on the Dog‐ADHD RS than did those with the long allele. Investigation of the joint effect of DRD4 and TH showed that dogs possessing long alleles at both sites were scored as less active–impulsive than were others. Our results are aligned with previous studies showing that DRD4 and TH polymorphisms are associated with activity–impulsivity related traits in dogs. However, the prevalence of variants of these genes differs across breeds, and the functional role of specific variants is unclear.     

Missense mutations in the BMP15 gene are associated with ovarian failure

Premature ovarian failure (POF) is an unexplained amenorrhoea (>6 months) with raised levels of gonadotropins (FSH>40 U/L) occurring before the age of 40 years. Recent studies have elucidated the role of oocyte derived growth factors (BMP15 and GDF9) in maintenance of folliculogenesis, granulosa cell (GC) proliferation and overall fertility. Our recently published work showed presence of two rare missense variants in the GDF9 gene associated with ovarian failure (Dixit et al. 2005, Menopause 12:749–754). The present case-control study has been structured to establish the role of BMP15 germline status associated with ovarian failure. Sequence analysis of the coding region of the BMP15 gene was carried out in a cohort of women with POF (n=133), primary amenorrhoea (n=60), and secondary amenorrhoea (n=9) compared with control females (n=197). This study revealed a total of 18 germline variants in the coding region of BMP15 gene, including 16 novel variants. These novel variants include one intronic variant, one 3’ flanking variant, one silent variant, and 13 missense variants. Eleven missense variants were present only in cases with complete absence in the control females. The remaining two missense variants viz. c.308A>G (p.Asn103Ser) and c.788_789insTCT (p.Leu263_Arg264insLeu) were present both in the cases and in the controls. The c.788_789insTCT variant was significantly higher in primary amenorrhoea cases than in the controls (Fisher’s exact test, P=0.034). Three frequent variants c.-9C>G, c.308A>G, and c.852C>T were chosen for haplotyping. The haplotype G-G-C was found to be significantly associated with ovarian failure (P=0.0075). In a nutshell, the BMP15 gene is highly associated with etiology of ovarian failure.     

Genetic analysis of <Emphasis Type="Italic">COL11A2</Emphasis> in Korean patients with autosomal dominant non-syndromic hearing loss

The collagen type XI alpha 2 gene (COL11A2) is associated with autosomal dominant non-syndromic hearing loss (ADNSHL), and all mutations of this gene in ADNSHL are missense mutations. To evaluate its potential as a major causative gene of ADNSHL in the Korean population, we performed genetic analysis of COL11A2 in 75 unrelated Korean patients with ADNSHL. Consequently, 5 non-synonymous variants, 7 synonymous variants, and 6 intronic variants were identified in COL11A2. Among them, a novel variant, p.G829R (c.2485G>C) was found in a patient as a heterozygote. However, pedigree analysis showed this variation was not co-segregated with hearing loss. Previously reported variants p.G230W (c.688G>T) and p.P1422L (c.4265C>T) were discovered in Korean patients. However, these variants were also detected in normal individuals. These results suggest that COL11A2 is not a major causative gene of ADNSHL in the Korean population.     

Genetic Variants on Chromosome 1p13.3 Are Associated with Non-ST Elevation Myocardial Infarction and the Expression of DRAM2 in the Finnish Population

Myocardial infarction (MI) is divided into either ST elevation MI (STEMI) or non-ST elevation MI (NSTEMI), differing in a number of clinical characteristics. We sought to identify genetic variants conferring risk to NSTEMI or STEMI by conducting a genome-wide association study (GWAS) of MI stratified into NSTEMI and STEMI in a consecutive sample of 1,579 acute MI cases with 1,576 controls. Subsequently, we followed the results in an independent population-based sample of 562 cases and 566 controls, a partially independent prospective cohort (N = 16,627 with 163 incident NSTEMI cases), and examined the effect of disease-associated variants on gene expression in 513 healthy participants. Genetic variants on chromosome 1p13.3 near the damage-regulated autophagy modulator 2 gene DRAM2 associated with NSTEMI (rs656843; odds ratio 1.57, P = 3.11 × 10⁻¹⁰) in the case-control analysis with a consistent but not statistically significant effect in the prospective cohort (rs656843; hazard ratio 1.13, P = 0.43). These variants were not associated with STEMI (rs656843; odds ratio, 1.11, P = 0.20; hazard ratio 0.97, P = 0.87), appearing to have a pronounced effect on NSTEMI risk. A majority of the variants at 1p13.3 associated with NSTEMI were also associated with the expression level of DRAM2 in blood leukocytes of healthy controls (top-ranked variant rs325927, P = 1.50 × 10⁻¹²). The results suggest that genetic factors may in part influence whether coronary artery disease results in NSTEMI rather than STEMI.     

P2RX7: Expression Responds to Sleep Deprivation and Associates with Rapid Cycling in Bipolar Disorder Type 1

Context

Rapid cycling is a severe form of bipolar disorder with an increased rate of episodes that is particularly treatment-responsive to chronotherapy and stable sleep-wake cycles. We hypothesized that the P2RX7 gene would be affected by sleep deprivation and be implicated in rapid cycling.

Objectives

To assess whether P2RX7 expression is affected by total sleep deprivation and if variation in P2RX7 is associated with rapid cycling in bipolar patients.

Design

Gene expression analysis in peripheral blood mononuclear cells (PBMCs) from healthy volunteers and case-case and case-control SNP/haplotype association analyses in patients.

Participants

Healthy volunteers at the sleep research center, University of California, Irvine Medical Center (UCIMC), USA (n = 8) and Swedish outpatients recruited from specialized psychiatric clinics for bipolar disorder, diagnosed with bipolar disorder type 1 (n = 569; rapid cycling: n = 121) and anonymous blood donor controls (n = 1,044).

Results

P2RX7 RNA levels were significantly increased during sleep deprivation in PBMCs from healthy volunteers (p = 2.3*10⁻⁹). The P2RX7 rs2230912 _A allele was more common (OR = 2.2, p = 0.002) and the ACGTTT haplotype in P2RX7 (rs1718119 to rs1621388) containing the protective rs2230912_G allele (OR = 0.45–0.49, p = 0.003–0.005) was less common, among rapid cycling cases compared to non-rapid cycling bipolar patients and blood donor controls.

Conclusions

Sleep deprivation increased P2RX7 expression in healthy persons and the putatively low-activity P2RX7 rs2230912 allele A variant was associated with rapid cycling in bipolar disorder. This supports earlier findings of P2RX7 associations to affective disorder and is in agreement with that particularly rapid cycling patients have a more vulnerable diurnal system.     

Exome sequencing reveals novel rare variants in Iranian familial multiple sclerosis: The importance of POLD2 in the disease pathogenesis

The prevalence of familial multiple sclerosis (FMS) is increasing worldwide which endorses the heritability of the disease. Given that many genome variations are ethnicity-specific and consanguineous marriage could affect genetic diseases, hereditary disease gene analysis among FMS patients from Iran, a country with high rates of parental consanguinity, could be highly effective in finding mutations underlying disease pathogenesis. To examine rare genetic mutations, we selected three Iranian FMS cases with ≥3 MS patients in more than one generation and performed whole exome sequencing. We identified a homozygous rare missense variant in POLD2 (p. Arg141Cys; rs372336011). Molecular dynamics analysis showed reduced polar dehydration energy and conformational changes in POLD2 mutant. Further, we found a heterozygote rare missense variant in NBFP1 (p. Gly487Asp; rs778806175). Our study revealed the possible role of novel rare variants in FMS. Molecular dynamic simulation provided the initial evidence of the structural changes behind POLD2 mutant.     

A CDH23 missense variant in Beauceron dogs with non-syndromic deafness

Congenital coat-colour-related deafness is common among certain canine breeds especially those exhibiting extreme white spotting or merle patterning. We identified a non-syndromic deafness in Beauceron dogs characterised by a bilateral hearing loss in puppies that is not linked to coat colour. Pedigree analysis suggested an autosomal recessive transmission. By combining homozygosity mapping with whole genome sequencing and variant filtering in affected dogs we identified a CDH23:c.700C>T variant. The variant, located in the CHD23 (cadherin related 23) gene, was predicted to induce a CDH23:p.(Pro234Ser) change in the protein. Proline-234 of CDH23 protein is highly conserved across different vertebrate species. In silico tools predicted the CDH23:p.(Pro234Ser) change to be deleterious. CDH23 encodes a calcium-dependent transmembrane glycoprotein localised near the tips of hair-cell stereocilia in the mammalian inner ear. Intact function of these cilia is mandatory for the transformation of the acoustical wave into a neurological signal, leading to sensorineural deafness when impaired. By genotyping a cohort of 90 control Beauceron dogs sampled in France, we found a 3.3% carrier frequency. The CDH23:c.[700C>T] allele is easily detectable with a genetic test to avoid at-risk matings.     

Variants in Nebulin (NEB) Are Linked to the Development of Familial Primary Angle Closure Glaucoma in Basset Hounds

Several dog breeds are susceptible to developing primary angle closure glaucoma (PACG), which suggests a genetic basis for the disease. We have identified a four-generation Basset Hound pedigree with characteristic autosomal recessive PACG that closely recapitulates PACG in humans. Our aim is to utilize gene mapping and whole exome sequencing approaches to identify PACG-causing sequence variants in the Basset. Extensive clinical phenotyping of all pedigree members was conducted. SNP-chip genotyping was carried out in 9 affected and 15 unaffected pedigree members. Two-point and multipoint linkage analyses of genome-wide SNP data were performed using Superlink-Online SNP-1.1 and a locus was mapped to chromosome 19q with a maximum LOD score of 3.24. The locus contains 12 Ensemble predicted canine genes and is syntenic to a region on chromosome 2 in the human genome. Using exome-sequencing analysis, a possibly damaging, non-synonymous variant in the gene Nebulin (NEB) was found to segregate with PACG which alters a phylogenetically conserved Lysine residue. The association of this variants with PACG was confirmed in a secondary cohort of unrelated Basset Hounds (p = 3.4 × 10^-4, OR = 15.3 for homozygosity). Nebulin, a protein that promotes the contractile function of sarcomeres, was found to be prominently expressed in the ciliary muscles of the anterior segment. Our findings may provide insight into the molecular mechanisms that underlie PACG. The phenotypic similarities of disease presentation in dogs and humans may enable the translation of findings made in this study to patients with PACG.     

The genetic association study of TP53 polymorphisms in Saudi obese patients

Obesity is a multifactorial metabolic disorder characterized by low grade chronic inflammation. Rare and novel mutations in genes which are vital in several key pathways have been reported to alter the energy expenditure which regulates body weight. The TP53 or p53 gene plays a prominent role in regulating various metabolic activities such as glycolysis, lipolysis, and glycogen synthesis. Recent genome-wide association studies reported that tumor suppressor gene p53 variants play a critical role in the predisposition of type 2 diabetes and obesity. Till date, no reports are available from the Arabian population; hence the present study was intended to assess the association between p53 variants with risk of obesity development in the Saudi population. We have selected three p53 polymorphisms, rs1642785 (C > G), and rs9894946 (A > G), and rs1042522 (Pro72Arg; C > G) and assessed their association with obesity risk in the Saudi population. Phenotypic and biochemical parameters were also evaluated to check their association with p53 genotypes and obesity. Genotyping was carried out on 136 obese and 122 normal samples. We observed that there is significantly increased prevalence p52 Pro72Arg (rs1042522) polymorphism in obese persons when compared to controls at GG genotype in overall comparison (OR: 2.169, 95% CI: 1.086-4.334, p = 0.02716). Male obese subjects showed three-fold higher risk at GG genotype (OR: 3.275, 95% CI: 1.230-8.716, p = 0.01560) and two-fold risk at G allele (OR: 1.827, 95% CI: 1.128-2.958, p = 0.01388) of p53 variant Pro72Arg respectively. This variant has also shown significant influence on cholesterol, LDL level, and random insulin levels in obese subjects (p ≤ 0.05). In conclusion, p53 Pro72Arg variant is highly prevalent among obese individuals and may act as a genetic modifier for obesity development among Saudis.     

Association of the melanocortin 4 receptor with feed intake and daily gain in F2 Mangalitsa x Piétrain pigs

The melanocortin 4 receptor (MC4R) is a key factor in the regulation of energy balance and body weight. Hence it is a candidate for feed intake and energy homeostasis-related traits. Studies in humans and swine have revealed several sequence variants in the gene that are associated with some of these traits. In pigs the coding non-synonymous missense variant Asp298Asn in MC4R has been associated with feed intake, fatness and growth. Here we confirm the association of this Piétrain-derived polymorphism with feed intake and daily gain in the F2 generation of a Mangalitsa x Piétrain cross. In one Piétrain founder animal, we detected an additional non-synonymous missense variant Arg236His. Thus, the MC4R gene could be a useful marker for increased growth in the relatively slow-growing Piétrain breed.