Single-marker and haplotype-based genome-wide association studies for the number of teats in two heavy pig breeds.

The number of teats is a reproductive-related trait of great economic relevance as it affects the mothering ability of the sows and thus the number of properly weaned piglets. Moreover, genetic improvement of this trait is fundamental to parallelly help the selection for increased litter size. We present the results of single-marker and haplotypes-based genome-wide association studies for the number of teats in two large cohorts of heavy pig breeds (Italian Large White and Italian Landrace) including 3990 animals genotyped with the 70K GGP Porcine BeadChip and other 1927 animals genotyped with the Illumina PorcineSNP60 BeadChip. In the Italian Large White population, genome scans identified three genome regions (SSC7, SSC10, and SSC12) that confirmed the involvement of the VRTN gene (as we previously reported) and highlighted additional loci known to affect teat counts, including the FRMD4A and HOXB1 gene regions. A different picture emerged in the Italian Landrace population, with a total of 12 genome regions in eight chromosomes (SSC3, SSC6, SSC8, SSC11, SSC13, SSC14, SSC15, and SSC16) mainly detected via the haplotype-based genome scan. The most relevant QTL was close to the ARL4C gene on SSC15. Markers in the VRTN gene region were not significant in the Italian Landrace breed. The use of both single-marker and haplotype-based genome-wide association analyses can be helpful to exploit and dissect the genome of the pigs of different populations. Overall, the obtained results supported the polygenic nature of the investigated trait and better elucidated its genetic architecture in Italian heavy pigs.

Exploiting within-breed variability in the autochthonous Reggiana breed identified several candidate genes affecting pigmentation-related traits, stature and udder defects in cattle.

Autochthonous cattle breeds constitute important reservoirs of genetic diversity. Reggiana is an Italian local cattle breed reared in the north of Italy for the production of a mono-breed Parmigiano-Reggiano cheese. Reggiana cattle usually have a classical solid red coat colour and pale muzzle. As part of the strategies designed for the sustainable conservation of this genetic resource, we investigated at the genome-wise level the within-breed detected variability of three pigmentation-related traits (intensity of red coat colour, based on three classes - light/diluted, normal and dark; spotted patterns/piebaldism that sometime emerge in the breed; muzzle colour - pink/pale, grey and black), stature, presence/absence and number of supernumerary teats and teat length. A total of 1776 Reggiana cattle (about two-thirds of the extant breed population) were genotyped with the GeneSeek GGP Bovine 150k SNP array and single-marker and haplotype-based GWASs were carried out. The results indicated that two main groups of genetic factors affect the intensity of red coat colour: darkening genes (including EDN3 and a few other genes) and diluting genes (including PMEL and a few other genes). Muzzle colour was mainly determined by MC1R gene markers. Piebaldism was mainly associated with KIT gene markers. Stature was associated with BTA6 markers upstream of the NCAPG-LCORL genes. Teat defects were associated with TBX3/TBX5, MCC and LGR5 genes. Overall, the identified genomic regions not only can be directly used in selection plans in the Reggiana breed, but also contribute to clarifying the genetic mechanisms involved in determining exterior traits in cattle.

Runs of homozygosity provide a genome landscape picture of inbreeding and genetic history of European autochthonous and commercial pig breeds.

ROHs are long stretches of DNA homozygous at each polymorphic position. The proportion of genome covered by ROHs and their length are indicators of the level and origin of inbreeding. Frequent common ROHs within the same population define ROH islands and indicate hotspots of selection. In this work, we investigated ROHs in a total of 1131 pigs from 20 European local pig breeds and in three cosmopolitan breeds, genotyped with the GGP Porcine HD Genomic Profiler. plink software was used to identify ROHs. Size classes and genomic inbreeding parameters were evaluated. ROH islands were defined by evaluating different thresholds of homozygous SNP frequency. A functional overview of breed-specific ROH islands was obtained via over-representation analyses of GO biological processes. Mora Romagnola and Turopolje breeds had the largest proportions of genome covered with ROH (~1003 and ~955 Mb respectively), whereas Nero Siciliano and Sarda breeds had the lowest proportions (~207 and 247 Mb respectively). The highest proportion of long ROH (>16 Mb) was in Apulo-Calabrese, Mora Romagnola and Casertana. The largest number of ROH islands was identified in the Italian Landrace (n = 32), Cinta Senese (n = 26) and Lithuanian White Old Type (n = 22) breeds. Several ROH islands were in regions encompassing genes known to affect morphological traits. Comparative ROH structure analysis among breeds indicated the similar genetic structure of local breeds across Europe. This study contributed to understanding of the genetic history of the investigated pig breeds and provided information to manage these pig genetic resources.

A comparative whole genome sequencing analysis identified a candidate locus for lack of operculum in cultivated gilthead seabream (Sparus aurata).

The gilthead seabream (Sparus aurata, Sparidae family) is commonly used for aquaculture. Despite its great economic value, several problems in its cultivation remain. One of the major concerns is the high frequency of morphological abnormalities occurring during the early developmental stages. Partial and/or total lack of operculum is the most frequent anomaly affecting the fish cranial region. The existence of genetic factors that can at least partially determine this defect has been hypothesized. In this work, two DNA pools of highly related fry, one composed of normal-looking (control) fish and the other lacking an operculum (case), were constructed and whole-genome resequencing data produced from the two were compared. The analysis revealed a 1 Mb region on chromosome 2 with higher heterozygosity in the lack of operculum DNA pool than in the control DNA pool, consistent with the enrichment, in the first DNA pool, of one or more haplotypes causing or predisposing to the defect together with other normal haplotypes. A window-based FST analysis between the two DNA pools indicated that the same region had the highest divergence score. This region contained 2921 SNVs, 10 of which, with predicted high impacts (three splice donor and seven stop-gained variants), were detected in novel genes that are homologous to calcium-sensing receptor-like genes, probably involved in bone development. Other studies are needed to clarify the genetic mechanisms involved in predisposing fry to this deformity and then to identify associated markers that could be used in breeding programs to reduce the frequency of this defect in the broodstock.

A genome-wide association study for the number of teats in European rabbits (Oryctolagus cuniculus) identifies several candidate genes affecting this trait.

In the European rabbit (Oryctolagus cuniculus), a polytocous livestock species, the number of teats indirectly impacts the doe reproduction efficiency and, in turn, the sustainable production of rabbit meat. In this study, we carried out a genome-wide association study (GWAS) for the total number of teats in 247 Italian White does included in the Italian White rabbit breed selection program, by applying a selective genotyping approach. Does had either 8 (n = 121) or 10 teats (n = 126). All rabbits were genotyped with the Affymetrix Axiom OrcunSNP Array. Genomic data from the two extreme groups of rabbits were also analysed with the single-marker fixation index statistic and combined with the GWAS results. The GWAS identified 50 significant SNPs and the fixation index analysis identified a total of 20 SNPs that trespassed the 99.98th percentile threshold, 19 of which confirmed the GWAS results. The most significant SNP (P = 4.31 × 10-11 ) was located on OCU1, close to the NUDT2 gene, a breast carcinoma cells proliferation promoter. Another significant SNP identified as candidate gene NR6A1, which is well known to play an important role in affecting the correlated number of vertebrae in pigs. Other significant markers were close to candidate genes involved in determining body length in mice. Markers associated with increased number of teats could be included in selection programmes to speed up the improvement for this trait in rabbit lines that need to increase maternal performances.

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states.

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30-35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.

Genome-wide association studies for the number of teats and teat asymmetry patterns in Large White pigs.

The number of teats is a morphological trait that influences the mothering ability of the sows and thus their reproduction performances. In this study, we carried out GWASs for the total number of teats and other 12 related parameters in 821 Italian Large White heavy pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip array. For four investigated parameters (total number of teats, the number of teats of the left line, the number of teats of the right line and the maximum number of teats comparing the two sides), significant markers were identified on SSC7, in the region of the vertnin (VRTN) gene. Significant markers for the numbers of posterior teats and the absolute difference between anterior and posterior teat numbers were consistently identified on SSC6. The most significant SNP for these parameters was an intron variant in the TOX high mobility group box family member 3 (TOX3) gene. For the other four parameters (absolute difference between the two sides; anterior teats; the ratio between the posterior and the anterior number of teats; and the absence or the presence of extra teats) only suggestively significant markers were identified on several other chromosomes. This study further supported the role of the VRTN gene region in affecting the recorded variability of the number of teats in the Italian Large White pig population and identified a genomic region potentially affecting the biological mechanisms controlling the developmental programme of morphological features in pigs.

Identification of a major locus determining a pigmentation defect in cultivated gilthead seabream (Sparus aurata).

The gilthead seabream (Sparus aurata) is an important cultivated species in the Mediterranean area. A major problem for the gilthead seabream aquaculture sector derives from the high frequency of phenotypic abnormalities, including discolorations. In this study, we applied a whole-genome resequencing approach to identify a genomic region affecting a pigmentation defect that occurred in a cultivated S. aurata population. Two equimolar DNA pools were constructed using DNA extracted from 30 normally coloured and 21 non-pigmented fish collected among the offspring of the same broodstock nucleus. Whole-genome resequencing reads from the two DNA pools were aligned to the S. aurata draft genome and variant calling was performed. A whole-genome heterozygosity scan from single pool sequencing data highlighted a peak of reduced heterozygosity of approximately 5 Mbp on chromosome 6 in the non-pigmented pool that was not present in the normally coloured pool. The comparison of the non-pigmented with the normally coloured fish using a whole-genome FST analysis detected three main regions within the coordinates previously detected with the heterozygosity analysis. The results support the presence of a major locus affecting this discoloration defect in this fish population. The results of this study have practical applications, including the possibility of eliminating this defect from the breeding stock, with direct economic advantages derived from the reduction of discarded fry. Other studies are needed to identify the candidate gene and the causative mutation, which could add information to understand the complex biology of fish pigmentation.

Genomic diversity, linkage disequilibrium and selection signatures in European local pig breeds assessed with a high density SNP chip.

Genetic characterization of local breeds is essential to preserve their genomic variability, to advance conservation policies and to contribute to their promotion and sustainability. Genomic diversity of twenty European local pig breeds and a small sample of Spanish wild pigs was assessed using high density SNP chips. A total of 992 DNA samples were analyzed with the GeneSeek Genomic Profiler (GGP) 70 K HD porcine genotyping chip. Genotype data was employed to compute genetic diversity, population differentiation and structure, genetic distances, linkage disequilibrium and effective population size. Our results point out several breeds, such as Turopolje, Apulo Calabrese, Casertana, Mora Romagnola and Lithuanian indigenous wattle, having the lowest genetic diversity, supported by low heterozygosity and very small effective population size, demonstrating the need of enhanced conservation strategies. Principal components analysis showed the clustering of the individuals of the same breed, with few breeds being clearly isolated from the rest. Several breeds were partially overlapped, suggesting genetic closeness, which was particularly marked in the case of Iberian and Alentejana breeds. Spanish wild boar was also narrowly related to other western populations, in agreement with recurrent admixture between wild and domestic animals. We also searched across the genome for loci under diversifying selection based on FST outlier tests. Candidate genes that may underlie differences in adaptation to specific environments and productive systems and phenotypic traits were detected in potentially selected genomic regions.

Signatures of de-domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution.

Autochthonous pig breeds are usually reared in extensive or semi-extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south-eastern European countries (Krskopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro-geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild-type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of 'de-domestication' process, and wild resources are challenged by a 'domestication' drift. Both need to be further investigated and managed.

Taking advantage from phenotype variability in a local animal genetic resource: identification of genomic regions associated with the hairless phenotype in Casertana pigs.

Casertana is an endangered autochthonous pig breed (raised in south-central Italy) that is considered to be the descendant of the influential Neapolitan pig population that was used to improve British breeds in the 19th century. Casertana pigs are characterized by a typical, almost complete, hairless phenotype, even though a few Casertana pigs are normal haired. In this work, using Illumina PorcineSNP60 BeadChip data, we carried out a genome-wide association study and an FST analysis with this breed by comparing animals showing the classical hairless phenotype (n = 81) versus pigs classified as haired (n = 15). Combining the results obtained with the two approaches, we identified two significant regions: one on porcine chromosome (SSC) 7 and one on SSC15. The SSC7 region contains the forkhead box N3 (FOXN3) gene, the most plausible candidate gene of this region, considering that mutations in another gene of the same family (forkhead box N1; Foxn1 or FOXN1) are responsible for the nude locus in rodents and alopecia in humans. Another potential candidate gene, rho guanine nucleotide exchange factor 10 (ARHGEF10), is located in the SSC15 region. FOXN3 and ARHGEF10 have been detected as differentially expressed in androgenetic and senescent alopecia respectively. This study on an autochthonous pig breed contributes to shed some light on novel genes potentially involved in hair development and growth and demonstrates that local animal breeds can be valuable genetic resources for disclosing genetic factors affecting unique traits, taking advantage of phenotype variability segregating in small populations.

Few mitochondrial DNA sequences are inserted into the turkey (Meleagris gallopavo) nuclear genome: evolutionary analyses and informativity in the domestic lineage.

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many eukaryotes. These sequences are pseudogenes originated by horizontal transfer of mtDNA fragments into the nuclear genome, producing nuclear DNA sequences of mitochondrial origin (numt). In this study we determined the frequency and distribution of mtDNA-originated pseudogenes in the turkey (Meleagris gallopavo) nuclear genome. The turkey reference genome (Turkey_2.01) was aligned with the reference linearized mtDNA sequence using last. A total of 32 numt sequences (corresponding to 18 numt regions derived by unique insertional events) were identified in the turkey nuclear genome (size ranging from 66 to 1415 bp; identity against the modern turkey mtDNA corresponding region ranging from 62% to 100%). Numts were distributed in nine chromosomes and in one scaffold. They derived from parts of 10 mtDNA protein-coding genes, ribosomal genes, the control region and 10 tRNA genes. Seven numt regions reported in the turkey genome were identified in orthologues positions in the Gallus gallus genome and therefore were present in the ancestral genome that in the Cretaceous originated the lineages of the modern crown Galliformes. Five recently integrated turkey numts were validated by PCR in 168 turkeys of six different domestic populations. None of the analysed numts were polymorphic (i.e. absence of the inserted sequence, as reported in numts of recent integration in other species), suggesting that the reticulate speciation model is not useful for explaining the origin of the domesticated turkey lineage.

Next generation semiconductor based sequencing of bitter taste receptor genes in different pig populations and association analysis using a selective DNA pool-seq approach.

Taste perception in animals affects feed intake and may influence production traits. In particular, bitter is sensed by receptors encoded by the family of TAS2R genes. In this research, using a DNA pool-seq approach coupled with next generation semiconductor based target resequencing, we analysed nine porcine TAS2R genes (TAS2R1, TAS2R3, TAS2R4, TAS2R7, TAS2R9, TAS2R10, TAS2R16, TAS2R38 and TAS2R39) to identify variability and, at the same time, estimate single nucleotide polymorphism (SNP) allele frequencies in several populations and testing differences in an association analysis. Equimolar DNA pools were prepared for five pig breeds (Italian Duroc, Italian Landrace, Pietrain, Meishan and Casertana) and wild boars (5-10 individuals each) and for two groups of Italian Large White pigs with extreme and divergent back fat thickness (50 + 50 pigs). About 1.8 million reads were obtained by sequencing amplicons generated from these pools. A total of 125 SNPs were identified, of which 37 were missense mutations. Three of them (p.Ile53Phe and p.Trp85Leu in TAS2R4; p.Leu37Ser in TAS2R39) could have important effects on the function of these bitter taste receptors, based on in silico predictions. Variability in wild boars seems lower than that in domestic breeds potentially as a result of selective pressure in the wild towards defensive bitter taste perception. Three SNPs in TAS2R38 and TAS2R39 were significantly associated with back fat thickness. These results may be important to understand the complexity of taste perception and their associated effects that could be useful to develop nutrigenetic approaches in pig breeding and nutrition.

The albinism of the feral Asinara white donkeys (Equus asinus) is determined by a missense mutation in a highly conserved position of the tyrosinase (TYR) gene deduced protein.

A feral donkey population (Equus asinus), living in the Asinara National Park (an island north-west of Sardinia, Italy), includes a unique white albino donkey subpopulation or colour morph that is a major attraction of this park. Disrupting mutations in the tyrosinase (TYR) gene are known to cause recessive albinisms in humans (oculocutaneous albinism Type 1; OCA1) and other species. In this study, we analysed the donkey TYR gene as a strong candidate to identify the causative mutation of the albinism of these donkeys. The TYR gene was sequenced from 13 donkeys (seven Asinara white albino and six coloured animals). Seven single nucleotide polymorphisms were identified. A missense mutation (c.604C>G; p.His202Asp) in a highly conserved amino acid position (even across kingdoms), which disrupts the first copper-binding site (CuA) of functional protein, was identified in the homozygous condition (G/G or D/D) in all Asinara white albino donkeys and in the albino son of a trio (the grey parents had genotype C/G or H/D), supporting the recessive mode of inheritance of this mutation. Genotyping 82 donkeys confirmed that Asinara albino donkeys had genotype G/G whereas all other coloured donkeys had genotype C/C or C/G. Across-population association between the c.604C>G genotypes and the albino coat colour was highly significant (P = 6.17E-18). The identification of the causative mutation of the albinism in the Asinara white donkeys might open new perspectives to study the dynamics of this putative deleterious allele in a feral population and to manage this interesting animal genetic resource.

Differentiation of meat from European wild boars and domestic pigs using polymorphisms in the MC1R and NR6A1 genes.

Wild boar meat cannot be easily distinguished from domestic pig meat, especially in processed products, thus it can be fraudulently substituted with cheaper domestic pork. In this study we genotyped polymorphisms in two genes (MC1R, affecting coat color and NR6A1, associated with number of vertebrae) in 293 domestic pigs of five commercial breeds, 111 wild boars sampled in Italy, and 90 in Slovenia and other Western Balkan regions. Allele and genotype frequency data were used to set up a DNA-based method to distinguish meat of wild boars and domestic pigs. Genotyping results indicated that domesticated genes were introgressed into wild boar populations. This complicated the determination of the origin of the meat and would cause a high error rate if markers of only one gene were used. The combined use of polymorphisms in the two analyzed genes substantially reduced false negative results.

A premature stop codon in the TYRP1 gene is associated with brown coat colour in the European rabbit (Oryctolagus cuniculus).

Classical genetic studies in European rabbits (Oryctolagus cuniculus) suggested the presence of two alleles at the brown coat colour locus: a wild-type B allele that gives dense black pigment throughout the coat and a recessive b allele that in the homozygous condition (b/b genotype) produces brown rabbits that are unable to develop black pigmentation. In several other species, this locus is determined by mutations in the tyrosinase-related protein 1 (TYRP1) gene, encoding a melanocyte enzyme needed for the production of dark eumelanin. In this study, we investigated the rabbit TYRP1 gene as a strong candidate for the rabbit brown coat colour locus. A total of 3846 bp of the TYRP1 gene were sequenced in eight rabbits of different breeds and identified 23 single nucleotide polymorphisms (SNPs; 12 in intronic regions, five in exons and six in the 3'-untranslated region) and an insertion/deletion of 13 bp, in the 3'-untranslated region, organised in a few haplotypes. A mutation in exon 2 (g.41360196G>A) leads to a premature stop codon at position 190 of the deduced amino acid sequence (p.Trp190ter). Therefore, translation predicts a truncated TYRP1 protein lacking almost completely the tyrosinase domain. Genotyping 203 rabbits of 32 different breeds identified this mutation only in brown Havana rabbits. Its potential functional relevance in disrupting the TYRP1 protein and its presence only in brown animals strongly argue for this non-sense mutation being a causative mutation for the recessive b allele at the brown locus in Oryctolagus cuniculus.

High-throughput SNP discovery in the rabbit (Oryctolagus cuniculus) genome by next-generation semiconductor-based sequencing.

The European rabbit (Oryctolagus cuniculus) is a domesticated species with one of the broadest ranges of economic and scientific applications and fields of investigation. Rabbit genome information and assembly are available (oryCun2.0), but so far few studies have investigated its variability, and massive discovery of polymorphisms has not been published yet for this species. Here, we sequenced two reduced representation libraries (RRLs) to identify single nucleotide polymorphisms (SNPs) in the rabbit genome. Genomic DNA of 10 rabbits belonging to different breeds was pooled and digested with two restriction enzymes (HaeIII and RsaI) to create two RRLs which were sequenced using the Ion Torrent Personal Genome Machine. The two RRLs produced 2 917 879 and 4 046 871 reads, for a total of 280.51 Mb (248.49 Mb with quality >20) and 417.28 Mb (360.89 Mb with quality >20) respectively of sequenced DNA. About 90% and 91% respectively of the obtained reads were mapped on the rabbit genome, covering a total of 15.82% of the oryCun2.0 genome version. The mapping and ad hoc filtering procedures allowed to reliably call 62 491 SNPs. SNPs in a few genomic regions were validated by Sanger sequencing. The Variant Effect Predictor Web tool was used to map SNPs on the current version of the rabbit genome. The obtained results will be useful for many applied and basic research programs for this species and will contribute to the development of cost-effective solutions for high-throughput SNP genotyping in the rabbit.