Luca Fontanesi

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication

Rabbits softly swept to domestication When people domesticate animals, they select for tameness and tolerance of humans. What else do they look for? To identify the selective pressures that led to rabbit domestication, Carneiro et al. sequenced a domestic rabbit genome and compared it to that of its wild brethren (see the Perspective by Lohmueller). Domestication did not involve a single gene changing, but rather many gene alleles changing in frequency between tame and domestic rabbits, known as a soft selective sweep. Many of these alleles have changes that may affect brain development, supporting the idea that tameness involves changes at multiple loci. Science, this issue p. 1074; see also p. 1000 The domestication of rabbits primarily shifted the frequencies of alleles represented, rather than creating new genes. [Also see Perspective by Lohmueller] The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

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.

Identification and association analysis of several hundred single nucleotide polymorphisms within candidate genes for back fat thickness in Italian Large White pigs using a selective genotyping approach1

Combining different approaches (resequencing of portions of 54 obesity candidate genes, literature mining for pig markers associated with fat deposition or related traits in 77 genes, and in silico mining of porcine expressed sequence tags and other sequences available in databases), we identified and analyzed 736 SNP within candidate genes to identify markers associated with back fat thickness (BFT) in Italian Large White sows. Animals were chosen using a selective genotyping approach according to their EBV for BFT (276 with most negative and 279 with most positive EBV) within a population of ≈ 12,000 pigs. Association analysis between the SNP and BFT has been carried out using the MAX test proposed for case-control studies. The designed assays were successful for 656 SNP: 370 were excluded (low call rate or minor allele frequency <5%), whereas the remaining 286 in 212 genes were taken for subsequent analyses, among which 64 showed a P(nominal) value <0.1. To deal with the multiple testing problem in a candidate gene approach, we applied the proportion of false positives (PFP) method. Thirty-eight SNP were significant (P(PFP) < 0.20). The most significant SNP was the IGF2 intron3-g.3072G>A polymorphism (P(nominal) < 1.0E-50). The second most significant SNP was the MC4R c.1426A>G polymorphism (P(nominal) = 8.0E-05). The third top SNP (P(nominal) = 6.2E-04) was the intronic TBC1D1 g.219G>A polymorphic site, in agreement with our previous results obtained in an independent study. The list of significant markers also included SNP in additional genes (ABHD16A, ABHD5, ACP2, ALMS1, APOA2, ATP1A2, CALR, COL14A1, CTSF, DARS, DECR1, ENPP1, ESR1, GH1, GHRL, GNMT, IKBKB, JAK3, MTTP, NFKBIA, NT5E, PLAT, PPARG, PPP2R5D, PRLR, RRAGD, RFC2, SDHD, SERPINF1, UBE2H, VCAM1, and WAT). Functional relationships between genes were obtained using the Ingenuity Pathway Analysis (IPA) Knowledge Base. The top scoring pathway included 19 genes with a P(nominal) < 0.1, 2 of which (IKBKB and NFKBIA) are involved in the hypothalamic IKKβ/NFκB program that could represent a key axis to affect fat deposition traits in pigs. These results represent a starting point to plan marker-assisted selection in Italian Large White nuclei for BFT. Because of similarities between humans and pigs, this study might also provide useful clues to investigate genetic factors affecting human obesity.

A genome wide association study for backfat thickness in Italian Large White pigs highlights new regions affecting fat deposition including neuronal genes

BackgroundCarcass fatness is an important trait in most pig breeding programs. Following market requests, breeding plans for fresh pork consumption are usually designed to reduce carcass fat content and increase lean meat deposition. However, the Italian pig industry is mainly devoted to the production of Protected Designation of Origin dry cured hams: pigs are slaughtered at around 160 kg of live weight and the breeding goal aims at maintaining fat coverage, measured as backfat thickness to avoid excessive desiccation of the hams. This objective has shaped the genetic pool of Italian heavy pig breeds for a few decades. In this study we applied a selective genotyping approach within a population of ~ 12,000 performance tested Italian Large White pigs. Within this population, we selectively genotyped 304 pigs with extreme and divergent backfat thickness estimated breeding value by the Illumina PorcineSNP60 BeadChip and performed a genome wide association study to identify loci associated to this trait.ResultsWe identified 4 single nucleotide polymorphisms with P≤5.0E-07 and additional 119 ones with 5.0E-07<P≤5.0E-05. These markers were located throughout all chromosomes. The largest numbers were found on porcine chromosomes 6 and 9 (n=15), 4 (n=13), and 7 (n=12) while the most significant marker was located on chromosome 18. Twenty-two single nucleotide polymorphisms were in intronic regions of genes already recognized by the Pre-Ensembl Sscrofa10.2 assembly. Gene Ontology analysis indicated an enrichment of Gene Ontology terms associated with nervous system development and regulation in concordance with results of large genome wide association studies for human obesity.ConclusionsFurther investigations are needed to evaluate the effects of the identified single nucleotide polymorphisms associated with backfat thickness on other traits as a pre-requisite for practical applications in breeding programs. Reported results could improve our understanding of the biology of fat metabolism and deposition that could also be relevant for other mammalian species including humans, confirming the role of neuronal genes on obesity.

Genetic heterogeneity at the bovine KIT gene in cattle breeds carrying different putative alleles at the spotting locus

According to classical genetic studies, piebaldism in cattle is largely influenced by the allelic series at the spotting locus (S), which includes the S(H) (Hereford pattern), S(+) (non-spotted) and s (spotted) alleles. The S locus was mapped on bovine chromosome 6 in the region containing the KIT gene. We investigated the KIT gene, analysing its variability and haplotype distribution in cattle of three breeds (Angus, Hereford and Holstein) with different putative alleles (S(+), S(H) and s respectively) at the S locus. Resequencing of a whole of 0.485 Mb revealed 111 polymorphisms. The global nucleotide diversity was 0.087%. Tajimas D-values were negative for all breeds, indicating putative directional selection. Of the 28 inferred haplotypes, only five were observed in the Hereford breed, in which one was the most frequent. Coalescent simulation showed that it is highly unlikely (P < 10E-6) to obtain this low number of haplotypes conditionally on the observed number of segregating SNPs. Therefore, the neutral model could be rejected for the Hereford breed, suggesting that a selection sweep occurred at the KIT locus. Twelve haplotypes were inferred in Holstein and Angus. For these two breeds, the neutral model could not be rejected. High heterogeneity of the KIT gene was confirmed from a phylogenetic analysis. Our results suggest a role of the KIT gene in determining the S(H) allele(s) in the Hereford, but no evidence of selective sweep was obtained in Holstein, suggesting that complex mechanisms (or other genes) might be the cause of the spotted phenotype in this breed.

Analysis of Horse Myostatin Gene and Identification of Single Nucleotide Polymorphisms in Breeds of Different Morphological Types

Myostatin (MSTN) is a negative modulator of muscle mass. We characterized the horse (Equus caballus) MSTN gene and identified and analysed single nucleotide polymorphisms (SNPs) in breeds of different morphological types. Sequencing of coding, untranslated, intronic, and regulatory regions of MSTN gene in 12 horses from 10 breeds revealed seven SNPs: two in the promoter, four in intron 1, and one in intron 2. The SNPs of the promoter (GQ183900:g.26T>C and GQ183900:g.156T>C, the latter located within a conserved TATA-box like motif) were screened in 396 horses from 16 breeds. The g.26C and the g.156C alleles presented higher frequency in heavy (brachymorphic type) than in light breeds (dolichomorphic type such as Italian Trotter breed). The significant difference of allele frequencies for the SNPs at the promoter and analysis of molecular variance (AMOVA) on haplotypes indicates that these polymorphisms could be associated with variability of morphology traits in horse breeds.

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.

Genetic heterogeneity and selection signature at the KIT gene in pigs showing different coat colours and patterns

Mutations in the porcine KIT gene (Dominant white locus) have been shown to affect coat colours and colour distribution in pigs. We analysed this gene in several pig breeds and populations (Sicilian black, completely black or with white patches; Cinta Senese; grey local population; Large White; Duroc; Hampshire; Pietrain; wild boar; Meishan) with different coat colours and patterns, genotyping a few polymorphisms. The 21 exons and parts of the intronic regions were sequenced in these pigs and 69 polymorphisms were identified. The grey-roan coat colour observed in a local grey population was completely associated with a 4-bp deletion of intron 18 in a single copy KIT gene, providing evidence that this mutation characterizes the I(d) allele described in the early genetic literature. The white patches observed in black Sicilian pigs were not completely associated with the presence of a duplicated KIT allele (I(p) ), suggesting that genetic heterogeneity is a possible cause of different coat colours in this breed. Selection signature was evident at the KIT gene in two different belted pig breeds, Hampshire and Cinta Senese. The same mutation(s) may cause the belted phenotype in these breeds that originated in the 18th-19th centuries from English pigs (Hampshire) and in Tuscany (Italy) in the 14th century (Cinta Senese). Phylogenetic relationships of 28 inferred KIT haplotypes indicated two clades: one of Asian origin that included Meishan and a few Sicilian black haplotypes and another of European origin.