Francesca Beretti

An initial comparative map of copy number variations in the goat (Capra hircus) genome

BackgroundThe goat (Capra hircus) represents one of the most important farm animal species. It is reared in all continents with an estimated world population of about 800 million of animals. Despite its importance, studies on the goat genome are still in their infancy compared to those in other farm animal species. Comparative mapping between cattle and goat showed only a few rearrangements in agreement with the similarity of chromosome banding. We carried out a cross species cattle-goat array comparative genome hybridization (aCGH) experiment in order to identify copy number variations (CNVs) in the goat genome analysing animals of different breeds (Saanen, Camosciata delle Alpi, Girgentana, and Murciano-Granadina) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome.ResultsWe identified a total of 161 CNVs (an average of 17.9 CNVs per goat), with the largest number in the Saanen breed and the lowest in the Camosciata delle Alpi goat. By aggregating overlapping CNVs identified in different animals we determined CNV regions (CNVRs): on the whole, we identified 127 CNVRs covering about 11.47 Mb of the virtual goat genome referred to the bovine genome (0.435% of the latter genome). These 127 CNVRs included 86 loss and 41 gain and ranged from about 24 kb to about 1.07 Mb with a mean and median equal to 90,292 bp and 49,530 bp, respectively. To evaluate whether the identified goat CNVRs overlap with those reported in the cattle genome, we compared our results with those obtained in four independent cattle experiments. Overlapping between goat and cattle CNVRs was highly significant (P < 0.0001) suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Genes with environmental functions were over-represented in goat CNVRs as reported in other mammals.ConclusionsWe describe a first map of goat CNVRs. This provides information on a comparative basis with the cattle genome by identifying putative recurrent interspecies CNVs between these two ruminant species. Several goat CNVs affect genes with important biological functions. Further studies are needed to evaluate the functional relevance of these CNVs and their effects on behavior, production, and disease resistance traits in goats.

Copy Number Variation and Missense Mutations of the Agouti Signaling Protein (ASIP) Gene in Goat Breeds with Different Coat Colors

In goats, classical genetic studies reported a large number of alleles at the Agouti locus with effects on coat color and pattern distribution. From these early studies, the dominant AWt (white/tan) allele was suggested to cause the white color of the Saanen breed. Here, we sequenced the coding region of the goat ASIP gene in 6 goat breeds (Girgentana, Maltese, Derivata di Siria, Murciano-Granadina, Camosciata delle Alpi, and Saanen), with different coat colors and patterns. Five single nucleotide polymorphisms (SNPs) were identified, 3 of which caused missense mutations in conserved positions of the cysteine-rich carboxy-terminal domain of the protein (p.Ala96Gly, p.Cys126Gly, and p.Val128Gly). Allele and genotype frequencies suggested that these mutations are not associated or not completely associated with coat color in the investigated goat breeds. Moreover, genotyping and sequencing results, deviation from Hardy-Weinberg equilibrium, as well as allele copy number evaluation from semiquantitative fluorescent multiplex PCR, indicated the presence of copy number variation (CNV) in all investigated breeds. To confirm the presence of CNV and evaluate its extension, we applied a bovine-goat cross-species array comparative genome hybridization (aCGH) experiment using a custom tiling array based on bovine chromosome 13. aCGH results obtained for 8 goat DNA samples confirmed the presence of CNV affecting a region of less that 100 kb including the ASIP and AHCY genes. In Girgentana and Saanen breeds, this CNV might cause the AWt allele, as already suggested for a similar structural mutation in sheep affecting the ASIP and AHCY genes, providing evidence for a recurrent interspecies CNV. However, other mechanisms may also be involved in determining coat color in these 2 breeds.

A first comparative map of copy number variations in the sheep genome

We carried out a cross species cattle-sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P<0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.

Missense and nonsense mutations in melanocortin 1 receptor (MC1R) gene of different goat breeds: association with red and black coat colour phenotypes but with unexpected evidences

BackgroundAgouti and Extension loci control the relative amount of eumelanin and pheomelanin production in melanocytes that, in turn, affects pigmentation of skin and hair. The Extension locus encodes the melanocortin 1 receptor (MC1R) whose permanent activation, caused by functional mutations, results in black coat colour, whereas other inactivating mutations cause red coat colour in different mammals.ResultsThe whole coding region of the MC1R gene was sequenced in goats of six different breeds showing different coat colours (Girgentana, white cream with usually small red spots in the face; Maltese, white with black cheeks and ears; Derivata di Siria, solid red; Murciano-Granadina, solid black or solid brown; Camosciata delle Alpi, brown with black stripes; Saanen, white; F1 goats and the parental animals). Five single nucleotide polymorphisms (SNPs) were identified: one nonsense mutation (p.Q225X), three missense mutations (p.A81V, p.F250V, and p.C267W), and one silent mutation. The stop codon at position 225 should cause the production of a shorter MC1R protein whose functionality may be altered. These SNPs were investigated in a larger sample of animals belonging to the six breeds. The Girgentana breed was almost fixed for the p.225X allele. However, there was not complete association between the presence of red spots in the face and the presence of this allele in homozygous condition. The same allele was identified in the Derivata di Siria breed. However, its frequency was only 33%, despite the fact that these animals are completely red. The p.267W allele was present in all Murciano-Granadina black goats, whereas it was never identified in the brown ones. Moreover, the same substitution was present in almost all Maltese goats providing evidence of association between this mutation and black coat colour.ConclusionAccording to the results obtained in the investigated goat breeds, MC1R mutations may determine eumelanic and pheomelanic phenotypes. However, they are probably not the only factors. In particular, the surprising not complete association of the nonsense mutation (p.Q225X) with red coat colour raises a few hypotheses on the determination of pheomelanic phenotypes in goats that should be further investigated.

Investigation of candidate genes for glycolytic potential of porcine skeletal muscle: Association with meat quality and production traits in Italian Large White pigs

The objective of this study was to investigate the association of DNA markers in candidate genes for glycolytic potential on meat quality parameters (pH(1), pH(u), glycogen and lactate content and glycolytic potential of semimembranosus muscle) and estimated breeding values (EBVs) for average daily gain, lean cuts, back fat thickness, ham weight, and feed:gain ratio in 272 Italian Large White pigs. Three mutations in the PRKAG3 gene (T30N, G52S and I199V) were investigated as well as single nucleotide polymorphisms in two other skeletal muscle genes (PGAM2 and PKM2) involved in the glycolytic pathway. Association analysis with the PRKAG3 markers showed significant results (P<0.05) only for pH(1) (I199V, with significant additive effect) and lactate content (T30N), confirming, at least in part, the effects of this gene on meat quality traits. Significant association (P<0.05) was also observed for PGAM2 and ham weight EBV with significant additive and dominance effects. PKM2 was associated with average daily gain, lean cuts (P<0.001), back fat thickness and feed:gain ratio (P<0.05), with significant additive and/or dominance effects on these traits. PKM2 encodes for a key enzyme of the muscle glycolytic pathway and maps on porcine chromosome 7 where other studies have reported important QTL for the same traits. These data might suggest an important function of this gene in the mechanisms that produce the observed effects. The results will be important to evaluate the inclusion of some of these DNA polymorphisms in marker assisted selection programs.

Single nucleotide polymorphisms in several porcine cathepsin genes are associated with growth, carcass, and production traits in Italian Large White pigs

To identify DNA markers associated with performance, carcass, and meat production traits including muscle postmortem cathepsin activity, several porcine genes encoding for lysosomal proteinases (cathepsin B, CTSB; cathepsin D, CTSD; cathepsin F, CTSF; cathepsin H, CTSH; cathepsin L, CTSL; and cathepsin Z, CTSZ) and for a cathepsin inhibitor (cystatin B) were investigated. Single nucleotide polymorphisms were identified in CTSD, CTSH, CTSL, and CTSZ genes with a combination of in silico expressed sequence tag database mining and single-strand conformation polymorphism analysis. Sequencing and PCR-RFLP protocols were used to validate the identified polymorphisms. Allele frequencies at these loci were investigated in Italian Large White, Landrace, Duroc, Piétrain, Belgian Landrace, Hampshire, and Meishan breeds. Genotyping CTSD and CTSH markers made it possible to genetically map these genes to SSC 2 and 7, respectively. Markers in CTSD, CTSH, CTSL, and CTSZ genes, together with mutations we previously reported in cystatin B, CTSB, and CTSF genes, were genotyped in an Italian Large White sib-tested population (272 or 482 animals). For these animals, meat quality traits (cathepsin B activity, pH measured at 2 h postmortem, pH measured at 24 h postmortem, glycogen, lactate, and glycolytic potential of semimembranosus muscle) and EBV for ADG, lean cuts (LC), backfat thickness (BFT), ham weight (HW), and feed:gain ratio (FGR) were determined. Analyzed markers did not show any association with muscle cathepsin B activity. Thus, it could be possible that different genes, other than these investigated candidates, affect this trait, which is correlated with the excessive softness defect of dry-cured hams. The results of association analysis confirmed the effects we already reported in another study for CTSF on ADG (P = 0.008), LC (P = 0.001), and BFT (P = 0.02). Moreover, CTSD was associated with ADG, LC (P < 0.0001), BFT, HW, and FGR (P < 0.001); CTSH was associated with FGR (P = 0.026); and CTSZ was associated with ADG (P = 0.006), LC (P = 0.01), HW (P = 0.024), and FGR (P = 0.029). The biochemical and physiological functions of the lysosomal proteinases, together with the results obtained in our investigation, suggest that the cathepsin gene family might play important roles affecting economic traits in pigs.

Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes

Massese is an Italian dairy sheep breed characterized by animals with black skin and horns and black or apparent grey hairs. Owing to the presence of these two coat colour types, this breed can be considered an interesting model to evaluate the effects of coat colour gene polymorphisms on this phenotypic trait. Two main loci have been already shown to affect coat colour in sheep: Agouti and Extension coding for the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes, respectively. The Agouti locus is affected by a large duplication including the ASIP gene that may determine the Agouti white and tan allele (A(Wt)). Other disrupting or partially inactivating mutations have been identified in exon 2 (a deletion of 5 bp, D(5); and a deletion of 9 bp, D(9)) and in exon 4 (g.5172T>A, p.C126S) of the ASIP gene. Three missense mutations in the sheep MC1R gene cause the dominant black E(D) allele (p.M73K and p.D121N) and the putative recessive e allele (p.R67C). Here, we analysed these ASIP and MC1R mutations in 161 Massese sheep collected from four flocks. The presence of one duplicated copy allele including the ASIP gene was associated with grey coat colour (P = 9.4E-30). Almost all animals with a duplicated copy allele (37 out of 41) showed uniform apparent grey hair and almost all animals without a duplicated allele (117 out of 120) were completely black. Different forms of duplicated alleles were identified in Massese sheep including, in almost all cases, copies with exon 2 disrupting or partially inactivating mutations making these alleles different from the A(Wt) allele. A few exceptions were observed in the association between ASIP polymorphisms and coat colour: three grey sheep did not carry any duplicated copy allele and four black animals carried a duplicated copy allele. Of the latter four sheep, two carried the E(D) allele of the MC1R gene that may be the cause of their black coat colour. The coat colour of all other black animals may be determined by non-functional ASIP alleles (non-agouti alleles, A(a)) and in a few cases by the E(D) Extension allele. At least three frequent ASIP haplotypes ([D(5):g.5172T], [N:g.5172A] and [D(5):g.5172A]) were detected (organized into six different diplotypes). In conclusion, the results indicated that coat colours in the Massese sheep breed are mainly derived by combining ASIP and MC1R mutations.

Evaluation of post mortem stability of porcine skeletal muscle RNA

The objective of this study was to evaluate the effect of postmortem times on the quality of porcine skeletal muscle total RNA in order to consider the possibility to use postmortem material for gene expression analysis. Samples of Musculus semimembranosus were collected at 20min, 2h, 6h, 24h and 48h postmortem from the left legs of four commercial heavy pigs. Total RNA was analysed by agarose gel electrophoresis stained with ethidium bromide and by microfluidic capillary electrophoresis on an Agilent 2100 Bioanalyzer instrument obtaining 28S:18S rRNA peak ratios and RIN values. The average RIN values of the analysed samples were 7.45±0.13, 7.43±0.15, 7.45±0.10, 7.33±0.15 and 3.95±0.58 for the same postmortem times, respectively, indicating that RNA degradation was present at 48h postmortem. In a similar experiment, carried out by other authors on beef cattle muscle total RNA extracted at different postmortem times, RNA was stable up to 8days after death as indicated by intact 28S and 18S rRNA bands. Thus, differences among species or other environmental factors might affect the level of RNA degradation. In the porcine postmortem samples, qualitative assessment of GAPDH transcripts by PCR amplification of different cDNA fragments indicated that postmortem stages did not affect the possibility of analysing this housekeeping gene. Thus, postmortem porcine skeletal muscle can be an useful tissue to obtain gene expression based information.

A melanocortin 1 receptor ( MC1R ) gene polymorphism is useful for authentication of Massese sheep dairy products

Massese is an Italian sheep breed, with black or grey coat colour, mainly reared in the Tuscany and Emilia Romagna regions. Recently, the emerging interests in this breed have resulted in the production of Pecorino cheese obtained with only Massese milk. In order to be profitable, this marketing link between Massese breed and its products should be defended against fraudsters who could include milk of other sheep breeds or cow milk in Massese labelled productions. To identify the genetic factors affecting coat colour in sheep, we have recently analysed the melanocortin 1 receptor (MC1R) gene and identified several single nucleotide polymorphisms (SNPs). In this work, as a first step to set up a DNA based protocol for authentication of Massese dairy products, we further investigated the presence and distribution of one of these SNPs (c.-31G>A) in 143 Massese sheep and in another 13 sheep breeds (for a total of 351 animals). The Massese breed was fixed for allele c.-31A, whereas in all other breeds allele c.-31 G was the most frequent or with frequency of 0·50. At the same nucleotide position the cattle MC1R gene carries the G nucleotide. Using these data we developed a method to detect adulterating milk (from other sheep breeds or from cow) in Massese dairy products based on the analysis of the c.-31G>A SNP. We first tested the sensitivity of the protocol and then applied it to analyse DNA extracted from ricotta and Pecorino cheese obtained with only Massese milk or obtained with unrestricted sheep and cattle milk. To our knowledge, this system represents the first one that can be used for breed authentication of a sheep production and that, at the same time, can reveal frauds derived from the admixture of milk of an unreported species.

Investigation of allele frequencies of the growth hormone receptor (GHR) F279Y mutation in dairy and dual purpose cattle breeds

Abstract A major QTL for milk production traits was reported in the middle of bovine chromosome 20 and, for its map position, the growth hormone receptor (GHR) gene was considered a strong positional and functional candidate gene. A missense mutation in exon VIII (F279Y amino acid substitution) showed highly significant effects mainly on milk composition (protein percentage and fat percentage) as well as on milk yield in several dairy cattle populations. As no information about the frequency of these two GHR alleles is available in any population, we studied their distribution in dairy and dual purpose cattle breeds reared in Italy. A total of 679 animals belonging to seven cattle breeds (Italian Holstein-Friesian, n=108; Italian Brown, n=104; Italian Simmental, n=104; Jersey, n=104; Reggiana, n=108; Modenese, n=66; Rendena, n=85) were sampled. A new PCR-RFLP protocol was designed to analyse this mutation inserting an artificial restriction site for enzyme SspI. In all investigated breeds, allele F was the most frequent and ranged from 0.947 (Italian Brown and Jersey) to 0.727 (Italian Holstein-Friesian). In Rendena, Italian Simmental, Reggiana and Modenese it was 0.824, 0.909, 0.921 and 0.924, respectively. For all breeds no significant deviation from the Hardy-Weinberg equilibrium was observed. Differences in allele frequencies were statistically significant between the Italian Holstein Friesian and Rendena against all other breeds. Due to the high frequency of the putative positive allele for milk protein percentage the use of the GHR F279Y marker in marker assisted selection plans should not have a great impact on this trait in the studied breeds.