Search

16 Pág. Departamento de ​Mejora Genética Animal, 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., This work has been funded by University of Bologna RFO 2016-2019 programs and by the European Union’s Horizon 2020 research and innovation programme (TREASURE project – grant agreement no. 634476). The content of this work reflects only the authors’ view and the European Union Agency is not responsible for any use that may be made of the information it contains.

European sea bass (Dicentrarchus labrax) is an important marine species for commercial and sport fisheries and aquaculture production. Recently, the European sea bass genome has been sequenced and assembled. This resource can open new opportunities to evaluate and monitor variability and identify variants that could contribute to the adaptation to farming conditions. In this work, two DNA pools constructed from cultivated European sea bass were sequenced using a next generation semiconductor sequencing approach based on Ion Proton sequencer. Using the first draft version of the D. labrax genome as reference, sequenced reads obtained a total of about 1.6 million of single nucleotide polymorphisms (SNPs), spread all over the chromosomes. Transition/transversion (Ti/Tv) was equal to 1.28, comparable to what was already reported in Salmon species. A pilot homozygosity analysis across the D. labrax genome using DNA pool sequence datasets indicated that this approach can identify chromosome regions with putative signatures of selection, including genes involved in ion transport and chloride channel functions, amino acid metabolism and circadian clock and related neurological systems. This is the first study that reported genome wide polymorphisms in a fish species obtained with the Ion Proton sequencer. Moreover, this study provided a methodological approach for selective sweep analysis in this species.

The gilthead seabream (Sparus aurata) is a species of relevance for the Mediterranean aquaculture industry. Despite the advancement of genetic tools for the species, breeding programs still do not often include genomics. In this study, we designed a genomic strategy to identify signatures of selection and genomic regions of high differentiation among populations of farmed fish stocks. A comparative DNA pooling sequencing approach was applied to identify signatures of selection in gilthead seabream from the same hatchery and from different nuclei that had not been subjected to genetic selection. Identified genomic regions were further investigated to detect SNPs with predicted high impact. The analyses underlined major genomic differences in the proportion of fixed alleles among the investigated nuclei. Some of these differences highlighted genomic regions, including genes involved in general metabolism and development already detected in QTL for growth, size, skeletal deformity, and adaptation to variation of oxygen levels in other teleosts. The obtained results pointed out the need to control the genetic effect of breeding programs in this species to avoid the reduction of genetic variability within populations and the increase in inbreeding level that, in turn, might lead to an increased frequency of alleles with deleterious effects.

Summary 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.

Coat colour in livestock speciesis one of the most distinctive traits that characterise many different breeds. Pigmentation does not affect only skin and hairs but also the eyes. Few genomic studies in cattle and horses also investigated iris pigmentation, but no similar studies was conducted in pigs thus far. In this study we analysed eye colour diversity in a Large White pig population (n=897) and reported the results of genome-wide association studies based on several comparisonsincluding pigs having four main eye colour categories(three with both pigmented eyes of different brown grades: pale, 17.9%; medium, 14.8%; and dark, 54.3%; another one with both eyes completely depigmented, 3.8%) and heterochromia patterns (heterochromia iridis, i.e. depigmented iris sectors in pigmented irises, 3.2%; heterochromia iridium, i.e. a whole eye iris of depigmented phenotype and the other eye with the iris completely pigmented, 5.9%). Pigs were genotyped with the Illumina PorcineSNP60 BeadChip and GEMMA was used for the association analyses. The results indicated that the eye pigmented patterns (different grades of brown pigmentation), the total absence of pigmentation in the both eyes, and heterochromia iridis defect were under controlled of SLC45A2 (on chromosome 16, SSC16), EDNRB (SSC11) and KITLG (SSC5), respectively. In addition, to new candidate associations for the eye depigmented patterns were also identified for SNPs on two SSC4 regions (including two candidate genes: NOTCH2 and PREX2) and on SSC6, SSC8 and SSC14 (including COL17A1 as candidate gene). This study provided useful information to understand the genetic mechanisms affecting eye pigmentation in pigs.

Whole genome sequencing (WGS) datasets, usually generated for the investigation of the individual animal genome, can be used for additional mining of the fraction of sequencing reads that remains unmapped to the respective reference genome. A significant proportion of these reads contains viral DNA derived from viruses that infected the sequenced animals. In this study, we mined more than 480 billion sequencing reads derived from 1471 WGS datasets produced from cattle, pigs, chickens and rabbits. We identified 367 different viruses among which 14, 11, 12 and 1 might specifically infect the cattle, pig, chicken and rabbit, respectively. Some of them are ubiquitous, avirulent, highly or potentially damaging for both livestock and humans. Retrieved viral DNA information provided a first unconventional and opportunistic landscape of the livestock viromes that could be useful to understand the distribution of some viruses with potential deleterious impacts on the animal food production systems.

Reggiana and Modenese are autochthonous cattle breeds, reared in the North of Italy, that can be mainly distinguished for their standard coat color (Reggiana is red, whereas Modenese is white with some pale gray shades). Almost all milk produced by these breeds is transformed into 2 mono-breed branded Parmigiano-Reggiano cheeses, from which farmers receive the economic incomes needed for the sustainable conservation of these animal genetic resources. After the setting up of their herd books in 1960s, these breeds experienced a strong reduction in the population size that was subsequently reverted starting in the 1990s (Reggiana) or more recently (Modenese) reaching at present a total of about 2,800 and 500 registered cows, respectively. Due to the small population size of these breeds, inbreeding is a very important cause of concern for their conservation programs. Inbreeding is traditionally estimated using pedigree data, which are summarized in an inbreeding coefficient calculated at the individual level (FPED). However, incompleteness of pedigree information and registration errors can affect the effectiveness of conservation strategies. High-throughput SNP genotyping platforms allow investigation of inbreeding using genome information that can overcome the limits of pedigree data. Several approaches have been proposed to estimate genomic inbreeding, with the use of runs of homozygosity (ROH) considered to be the more appropriate. In this study, several pedigree and genomic inbreeding parameters, calculated using the whole herd book populations or considering genotyping information (GeneSeek GGP Bovine 150K) from 1,684 Reggiana cattle and 323 Modenese cattle, were compared. Average inbreeding values per year were used to calculate effective population size. Reggiana breed had generally lower genomic inbreeding values than Modenese breed. The low correlation between pedigree-based and genomic-based parameters (ranging from 0.187 to 0.195 and 0.319 to 0.323 in the Reggiana and Modenese breeds, respectively) reflected the common problems of local populations in which pedigree records are not complete. The high proportion of short ROH over the total number of ROH indicates no major recent inbreeding events in both breeds. ROH islands spread over the genome of the 2 breeds (15 in Reggiana and 14 in Modenese) identified several signatures of selection. Some of these included genes affecting milk production traits, stature, body conformation traits (with a main ROH island in both breeds on BTA6 containing the ABCG2, NCAPG, and LCORL genes) and coat color (on BTA13 in Modenese containing the ASIP gene). In conclusion, this work provides an extensive comparative analysis of pedigree and genomic inbreeding parameters and relevant genomic information that will be useful in the conservation strategies of these 2 iconic local cattle breeds.

14 Pág. Departamento de ​Mejora Genética Animal, Coronaviruses silently circulate in human and animal populations, causing mild to severe diseases. Therefore, livestock are important components of a "One Health" perspective aimed to control these viral infections. However, at present there is no example that considers pig genetic resources in this context. In this study, we investigated the variability of four genes (ACE2, ANPEP and DPP4 encoding for host receptors of the viral spike proteins and TMPRSS2 encoding for a host proteinase) in 23 European (19 autochthonous and three commercial breeds and one wild boar population) and two Asian Sus scrofa populations. A total of 2229 variants were identified in the four candidate genes: 26% of them were not previously described; 29 variants affected the protein sequence and might potentially interact with the infection mechanisms. The results coming from this work are a first step towards a "One Health" perspective that should consider conservation programs of pig genetic resources with twofold objectives: (i) genetic resources could be reservoirs of host gene variability useful to design selection programs to increase resistance to coronaviruses; (ii) the described variability in genes involved in coronavirus infections across many different pig populations might be part of a risk assessment including pig genetic resources., This work has received funding from the University of Bologna RFO 2016-2019 programmes, the Italian MIUR 2017 PigPhenomics project, from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 634476 for the project with acronym TREASURE, from the European Open Science Cloud (EOSC) Secretariat, project “Application of animal genomics and data mining to predict and monitor novel coronavirus potential infections (VirAnimalOne)”, the EGI call for COVID-19 research projects (AnGen1H project) and from the Por Fesr Emilia-Romagna 2014-2020 (actions 1.1.4 and 1.2.2—Bando per sostenere progetti di ricerca ed innovazione per lo sviluppo di soluzioni finalizzate al contrasto dell’epidemia da COVID-19—Project LIVESTOCK-STOP-COVI). The content of this article reflects only the authors' view and the European Union Agency is not responsible for any use that may be made of the information it contains.

Simple Summary Honey contains traces of the DNA of the honey bees that produced it. This environmental DNA can therefore be used to investigate the genome of the honey bees. In this study, we used a next generation sequencing technology to analyze the variability of a key gene of Apis mellifera L., the complementary sex determiner (csd) gene, using honey environmental DNA as a source of honey bee DNA. This gene determines the sex of the bees. Two different alleles at this locus are needed to produce females whereas males have only one copy of this gene as they are haploid. In case two identical alleles are present in a diploid individual, the larvae are not vital and are discarded by the workers. Therefore, there is an advantage in maintaining a large csd diversity in honey bee populations. In light of the recent decline in honey bee populations, it is important to monitor the allele variability at this gene. The applied methodology provided a new strategy to disclose the genetic diversity at the csd gene at the population-wide level and identify most, if not all, csd alleles present in the colonies in a single analysis. Abstract The complementary sex determiner (csd) gene plays an essential role in the sex determination of Apis mellifera L. Females develop only if fertilized eggs have functional heterozygous genotypes at this gene whereas males, being haploids, are hemizygous. Two identical csd alleles produce non vital males. In light of the recent decline in honey bee populations, it is therefore important to monitor the allele variability at this gene. In this study, we tested the application of next generation semiconductor-based sequencing technology (Ion Torrent) coupled with environmental honey DNA as a source of honey bee genome information to retrieve massive sequencing data for the analysis of variability at the hypervariable region (HVR) of the csd gene. DNA was extracted from 12 honey samples collected from honeycombs directly retrieved from 12 different colonies. A specifically designed bioinformatic pipeline, applied to analyze a total of about 1.5 million reads, identified a total of 160 different csd alleles, 55% of which were novel. The average number of alleles per sample was compatible with the number of expected patrilines per colony, according to the mating behavior of the queens. Allele diversity at the csd could also provide information useful to reconstruct the history of the honey.

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.

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 1Mb 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.

Summary 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.

Disrupting mutations affecting the tyrosinase (TYR) gene cause different forms of albinism in mice, humans and several other mammals. Classical genetic studies have already reported five alleles at the European rabbit Albino locus, indicated to be part of the C series, each of them with different actions on pheomelanin and eumelanin production, as well as on the eye colour. A few of these alleles have been already characterised at the DNA level by sequencing the coding region of the TYR gene in few rabbit breeds or strains with specific alleles at this locus. In this study, we further characterised the TYR gene by sequencing all coding and flanking regions in a total of 25 rabbits from 11 domestic breeds (Belgian Hare, Burgundy Fawn, Californian, Champagne d’Argent, Giant Chinchilla, Giant Grey, Havana, Leprino di Viterbo, New Zealand White, Silver and White Vienna) and in 11 wild rabbits from Sardinia. Sequencing data identified a total of 15 polymorphisms. We confirmed five missense mutations already detected by other authors, three of which associated with different coat colour phenotypes: p.T373K determining the albino allele; p.E294G causing the Himalayan and the chinchilla alleles; p.T358I observed only in Chinchilla rabbits. In addition to seven other synonymous mutations and one polymorphism in the 3’-untranslated region, two novel missense mutations were identified (p.T144S and p.K224T, the latest was detected only in wild rabbits). This study further contributed to disclose variability in the rabbit TYR gene and confirmed the effects on coat colours of missense mutations.Highlights Mutations in the tyrosinase (TYR) gene cause several coat colour alleles at the Albino locus in mammals. The TYR gene was sequenced in 11 domestic rabbit breeds and in a wild rabbit population. Fifteen single nucleotide polymorphisms (SNPs) were identified: seven were missense mutations. In silico analyses and information already available predicted the effects of some SNPs on coat colour. This study further contributed to disclose variability in the rabbit TYR gene.

Honey contains DNA from many different organisms that are part of hive micro-environmental niches and honey bee pathospheres. In this study, we recovered and sequenced mite mitochondrial DNA (mtDNA) from honey from different locations around the world (Europe, Asia, Africa, North and South America). DNA extracted from 17 honey samples was amplified with eight primer pairs targeting three mite mtDNA genes, obtaining 88 amplicons that were sequenced with an Ion Torrent sequencing platform. A bioinformatic pipeline compared produced reads with Varroa spp. mtDNA sequence entries available in GenBank and assigned them to different mitotypes. In all honey samples, the highest percentage of reads was attributed to the K1 lineage, including a few variants derived from it, in addition to J1 reads observed in the two South American samples and C1-1 reads obtained from the Chinese honey. This study opens new possibilities to analyse mite lineages and variants and monitor their geographical and temporal distribution, simplifying surveillance against this damaging honey bee parasite.

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.

Mitochondrial DNA (mtDNA) insertions have been detected in the nuclear genome of many mammals. These mtDNA related regions (called nuclear DNA sequences of mitochondrial origin or NUMTs) are pseudogenes derived from mtDNA fragments that have been integrated into the nuclear genome. Unintentionally amplified NUMTs can create biases in detecting true heteroplasmy and in phylogenetic analyses based on mtDNA as also demonstrated by the large fraction of mtDNA sequences available in GenBank/EMBL databases that are estimated to contain errors. In this study we obtained a detailed genome map of NUMTs in the Bos taurus genome and compared their distribution between the latest assembled versions, UMD3.1 and ARS-UCD1.2. These two genome versions were aligned with the reference linearized mtDNA sequence using LAST software. The highest proportion of NUMTs had an identity of about 70-80% with the modern mtDNA sequence, suggesting that most integration events occurred about 10 to 30 million of years before present in ancient genomes that evolved in the current Bos taurus species. A total of 422 and 441 NUMTs (for a total of 242,409 and 238,345 bp, which covered the whole mtDNA sequence) were identified in the two versions.These NUMT sequences could be grouped into 171 and 178 NUMT regions that represent putative unique insertional events. Two gff files were produced to support the annotation of the two cattle genome versions. There data will be helpful to plan studies based on mtDNA sequences in cattle and to interpret previous mtDNA sequence results obtained in this species.

Commercial single nucleotide polymorphism (SNP) arrays have been recently developed for several species and can be used to identify informative markers to differentiate breeds or populations for several downstream applications. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this work, we compared several methods of SNPs preselection (Delta, F st and principal component analyses (PCA)) in addition to Random Forest classifications to analyse SNP data from six dairy cattle breeds, including cosmopolitan (Holstein, Brown and Simmental) and autochthonous Italian breeds raised in two different regions and subjected to limited or no breeding programmes (Cinisara, Modicana, raised only in Sicily and Reggiana, raised only in Emilia Romagna). From these classifications, two panels of 96 and 48 SNPs that contain the most discriminant SNPs were created for each preselection method. These panels were evaluated in terms of the ability to discriminate as a whole and breed-by-breed, as well as linkage disequilibrium within each panel. The obtained results showed that for the 48-SNP panel, the error rate increased mainly for autochthonous breeds, probably as a consequence of their admixed origin lower selection pressure and by ascertaining bias in the construction of the SNP chip. The 96-SNP panels were generally more able to discriminate all breeds. The panel derived by PCA-chrom (obtained by a preselection chromosome by chromosome) could identify informative SNPs that were particularly useful for the assignment of minor breeds that reached the lowest value of Out Of Bag error even in the Cinisara, whose value was quite high in all other panels. Moreover, this panel contained also the lowest number of SNPs in linkage disequilibrium. Several selected SNPs are located nearby genes affecting breed-specific phenotypic traits (coat colour and stature) or associated with production traits. In general, our results demonstrated the usefulness of Random Forest in combination to other reduction techniques to identify population informative SNPs.

Background Natural and artificial directional selection in cosmopolitan and autochthonous pig breeds and wild boars have shaped their genomes and resulted in a reservoir of animal genetic diversity. Signatures of selection are the result of these selection events that have contributed to the adaptation of breeds to different environments and production systems. In this study, we analysed the genome variability of 19 European autochthonous pig breeds (Alentejana, Bísara, Majorcan Black, Basque, Gascon, Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda, Krškopolje pig, Black Slavonian, Turopolje, Moravka, Swallow-Bellied Mangalitsa, Schwäbisch-Hällisches Schwein, Lithuanian indigenous wattle and Lithuanian White old type) from nine countries, three European commercial breeds (Italian Large White, Italian Landrace and Italian Duroc), and European wild boars, by mining whole-genome sequencing data obtained by using a DNA-pool sequencing approach. Signatures of selection were identified by using a single-breed approach with two statistics [within-breed pooled heterozygosity (HP) and fixation index (FST)] and group-based FST approaches, which compare groups of breeds defined according to external traits and use/specialization/type. Results We detected more than 22 million single nucleotide polymorphisms (SNPs) across the 23 compared populations and identified 359 chromosome regions showing signatures of selection. These regions harbour genes that are already known or new genes that are under selection and relevant for the domestication process in this species, and that affect several morphological and physiological traits (e.g. coat colours and patterns, body size, number of vertebrae and teats, ear size and conformation, reproductive traits, growth and fat deposition traits). Wild boar related signatures of selection were detected across all the genome of several autochthonous breeds, which suggests that crossbreeding (accidental or deliberate) occurred with wild boars. Conclusions Our findings provide a catalogue of genetic variants of many European pig populations and identify genome regions that can explain, at least in part, the phenotypic diversity of these genetic resources.

We report haplotype-based GWASs for 33 blood parameters measured in 843 Italian Large White pigs. In the single-trait analysis, a total of 30 QTL for number of basophils, six erythrocyte traits (haemoglobin, haematocrit, mean corpuscular haemoglobin, mean corpuscular haemoglobin concentration, mean corpuscular volume and red blood cell count) and two clinical-biochemical traits (alkaline phosphatase and Ca2+ contents) were identified. In the multiple-trait analysis, a total of five QTL affected three different clusters of traits. Only four of these QTL were already reported in the single-marker and multi-marker GWASs we previously carried out on the same pig population. QTL on SSC11 and SSC17 showed effects on multiple traits. These results further dissected the genetic architecture of parameters that could be used as proxies in breeding programmes for more complex traits. In addition, these results might help to better define the pig as an animal model for several blood-related biological functions.

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. info:eu-repo/semantics/acceptedVersion

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 5Mbp 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.

Italian heavy pig selection programs are based on production (average daily gain, ADG; feed gain ratio, FGR), carcass (weight of lean cuts, LC; ham weight, HW; backfat thickness, BFT) and meat quality traits of the green hams (ham weight loss at first salting, HWLFS; and visible intermuscular fat, VIF). In this study, we evaluated whether a selective genotyping approach previously designed for BFT could identify markers associated with other traits (ADG, FGR, HW, HWLFS, LC and VIF), most of which correlated with BFT. Italian Large White gilts were chosen according to their extreme estimated breeding values for BFT (275 with the most negative and 278 with the most positive EBVs) within a population of about 12,000 sib-tested pigs. A total of 736 single nucleotide polymorphisms (SNPs) in candidate genes for fat deposition and/or related traits were genotyped. After filtering, a total of 279 SNPs in 234 genes remained for association analyses that were carried out using random residuals (RRs) for the six mentioned traits. Four SNPs in three candidate genes were associated with some of the considered traits. A SNP in the exosome component 1 (EXOSC1) gene and a SNP in the glycogen phosphorylase L (PYGL) gene were associated with HWLFS RRs. Two SNPs within the TBC1 (tre-2/USP6, BUB2, cdc16) domain family, member 1 (TBC1D1) gene were associated with LC RRs and suggestively associated with other performance, carcass and meat traits. Association analysis using six TBC1D1 haplotypes, inferred from four genotyped SNPs within this gene, confirmed the results of the single-marker association study. This study contributes to dissect the biological mechanisms underlying the genetic correlations among BFT, HWLFS and LC, which are critical for the sustainable implementation of tailored selection programs in heavy pigs.

Haematological and clinical-biochemical parameters are considered indicators of the physiological/health status of animals and might serve as intermediate phenotypes to link physiological aspects to production and disease resistance traits. The dissection of the genetic variability affecting these phenotypes might be useful to describe the resilience of the animals and to support the usefulness of the pig as animal model. Here, we analysed 15 haematological and 15 clinical-biochemical traits in 843 Italian Large White pigs, via three genome-wide association scan approaches (single-trait, multi-trait and Bayesian). We identified 52 quantitative trait loci (QTLs) associated with 29 out of 30 analysed blood parameters, with the most significant QTL identified on porcine chromosome 14 for basophil count. Some QTL regions harbour genes that may be the obvious candidates: QTLs for cholesterol parameters identified genes (ADCY8, APOB, ATG5, CDKAL1, PCSK5, PRL and SOX6) that are directly involved in cholesterol metabolism; other QTLs highlighted genes encoding the enzymes being measured [ALT (known also as GPT) and AST (known also as GOT)]. Moreover, the multivariate approach strengthened the association results for several candidate genes. The obtained results can contribute to define new measurable phenotypes that could be applied in breeding programs as proxies for more complex traits.

Intermuscular fat content in protected designations of origin dry-cured hams is a very important meat quality trait that affects the acceptability of the product by the consumers. An excess in intermuscular fat (defined as the level of fat deposition between leg muscles) is a defect that depreciates the final product. In this study we carried out a genome-wide association study for visible intermuscular fat (VIF) of hams in the Italian Large White pig breed. This trait was evaluated on the exposed muscles of green legs in 1122 performance-tested gilts by trained personnel, according to a classification scale useful for routine and cheap evaluation. All animals were genotyped with the Illumina PorcineSNP60 BeadChip. The genome-wide association study identified three QTL regions on porcine chromosome 1 (SSC1; accounting for ~79% of the SNPs below the 5.0E−04 threshold) and SSC2, two on SSC7 and one each on SSC3, SSC6, SSC9, SSC11, SSC13, SSC15, SSC16 and SSC17. The most significant SNP (ALGA0004143 on SSC1 at 77.3 Mb; PFDRÂ

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.

Summary Calcium, magnesium and phosphorus are essential electrolytes involved in a large number of biological processes. Imbalance of these minerals in blood may indicate clinically relevant conditions and are important in inferring acute or chronic pathologies in humans and animals. In this work, we carried out a genome-wide association study (GWAS) for the level of these three electrolytes in the serum of 843 performance-tested Italian Large White pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip, and GWAS was carried out using genome-wide efficient mixed-model association. For the level of Ca2+, eight single nucleotide polymorphisms (SNPs) were significant, considering a false discovery rate (FDR)

Knowledge of genes involved in the variability of performance of racehorses is, in general, limited. In this study, six single nucleotide polymorphisms (SNPs) of five genes were genotyped in 412 Italian Trotter horses (ITHs) and possible associations between five segregating SNPs and harness racing performance (best racing time in career, earnings, placings, starts, and wins) were investigated. These markers were in genes related to exercise physiology (creatine kinase muscle, CKM; and cytochrome c oxidase, subunit 4, isoform 2, COX412), gait patterns (DMRT3), carbohydrates metabolism (GCK) and morphological phenotypes (LCORL). Two SNPs of the GCK gene were used to predict haplotypes. The mutate A allele of the DMRT3 was almost fixed in the ITHs. Single marker association analyses showed that the COX412 SNP was significantly associated with variation of number of starts (P = 0.0003), and wins (P = 0.0045). This COX4I2 polymorphic site seems therefore to be a promising marker for racing longevity and durability. A haplotype of the GCK was instead positively associated with best racing time (P = 0.0118), earnings (P = 0.0370) and wins (P = 0.0061). This is the first study highlighting the effect of GCK on harness racing performance in horses. Confirmation of these results in independent populations would be important to evaluate the inclusion of the SNPs in marker-assisted selection programmes in ITHs.

Casertana is an autochthonous pig breed reared in the Central-South of Italy, mainly in extensive or semi-extensive production systems. A small number of animals are registered to the breed herd-book. The population of this breed is not completely fixed for several exterior traits, including ear size (small vs large), ear bearing (forward vs floppy), presence or absence of wattles and coat colour (black vs slate-grey) that we recorded in a total of 101 pigs from six different farms. All these pigs were also genotyped with the Illumina PorcineSNP60 BeadChip. Then, with the aim to identify DNA markers associated with these traits, we run four genome-wide association studies (GWAS) comparing pigs having different phenotypes for each of the four external traits. Genome-wide fixation index (FST) analyses were run to confirm the results. GEMMA software was used in the GWAS. Single nucleotide polymorphisms (SNPs) associated to the wattle trait were identified on porcine chromosome (SSC) 1 and SSC11. Markers associated to the coat colour trait were located on SSC6, SSC8, SSC14 and SSC15. These chromosome regions were confirmed by the independent FST analyses. FST did not confirm the signals detected for the two ear conformation traits. Candidate genes were identified according to the function of the annotated genes in the detected genome regions. Despite the small number of pigs that could be investigated, this study showed that autochthonous animal genetic resources might be useful to disclose genetic factors affecting traits that are not fixed in a population.

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.

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.

In this study, we investigated at the genome-wide level if 20 years of artificial directional selection based on boar genetic evaluation obtained with a classical BLUP animal model shaped the genome of the Italian Large White pig breed. The most influential boars of this breed (n = 192), born from 1992 (the beginning of the selection program of this breed) to 2012, with an estimated breeding value reliability of >0.85, were genotyped with the Illumina Porcine SNP60 BeadChip. After grouping the boars in eight classes according to their year of birth, filtered single nucleotide polymorphisms (SNPs) were used to evaluate the effects of time on genotype frequency changes using multinomial logistic regression models. Of these markers, 493 had a PBonferroni

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.

The acquisition of epithelial-neuroendocrine differentiation (ND) is a peculiarity of human advanced, androgen-independent, prostate cancers. The HOX genes are a network of transcription factors controlling embryonal development and playing an important role in crucial adult eukaryotic cell functions. The molecular organization of this 39-gene network is unique in the genome and probably acts by regulating phenotype cell identity. The expression patterns of the HOX gene network in human prostate cell phenotypes, representing different stages of prostate physiology and prostate cancer progression, make it possible to discriminate between different human prostate cell lines and to identify loci and paralogous groups harboring the HOX genes mostly involved in prostate organogenesis and cancerogenesis. Exposure of prostate epithelial phenotypes to cAMP alters the expression of lumbo-sacral HOX D genes located on the chromosomal region 2q31-33 where the cAMP effector genes CREB1, CREB2, and cAMP-GEFII are present. Interestingly, this same chromosomal area harbors: (i) a global cis-regulatory DNA control region able to coordinate the expression of HOX D and contiguous phylogenetically unrelated genes; (ii) a prostate specific ncRNA gene associated with high-risk prostate cancer (PCGEM1); (iii) a series of neurogenic-related genes involved with epithelial-neuronal cell conversion. We report the expression of neurexin 1, Neuro D1, dlx1, and dlx2 in untreated and cAMP treated epithelial prostate cells. The in vivo expression of Neuro D1 in human advanced prostate cancers correlate with the state of tumor differentiation as measured by Gleason score. Thus, we suggest that the chromosomal area 2q 31-33 might be involved in the epithelial-ND characteristic of human advanced prostate cancers.

Few studies investigated the donkey (Equus asinus) at the whole genome level so far. Here, we sequenced the genome of two male donkeys using a next generation semiconductor based sequencing platform (the Ion Proton sequencer) and compared obtained sequence information with the available donkey draft genome (and its Illumina reads from which it was originated) and with the EquCab2.0 assembly of the horse genome. Moreover, the Ion Torrent Personal Genome Analyzer was used to sequence reduced representation libraries (RRL) obtained from a DNA pool including donkeys of different breeds (Grigio Siciliano, Ragusano and Martina Franca). The number of next generation sequencing reads aligned with the EquCab2.0 horse genome was larger than those aligned with the draft donkey genome. This was due to the larger N50 for contigs and scaffolds of the horse genome. Nucleotide divergence between E. caballus and E. asinus was estimated to be similar to 0.52-0.57%. Regions with low nucleotide divergence were identified in several autosomal chromosomes and in the whole chromosome X. These regions might be evolutionally important in equids. Comparing Y-chromosome regions we identified variants that could be useful to track donkey paternal lineages. Moreover, about 4.8 million of single nucleotide polymorphisms (SNPs) in the donkey genome were identified and annotated combining sequencing data from Ion Proton (whole genome sequencing) and Ion Torrent (RRL) runs with Illumina reads. A higher density of SNPs was present in regions homologous to horse chromosome 12, in which several studies reported a high frequency of copy number variants. The SNPs we identified constitute a first resource useful to describe variability at the population genomic level in E. asinus and to establish monitoring systems for the conservation of donkey genetic resources.

Summary The genetic identification of the population of origin of individuals, including animals, has several practical applications in forensics, evolution, conservation genetics, breeding and authentication of animal products. Commercial high-density single nucleotide polymorphism (SNP) genotyping tools that have been recently developed in many species provide information from a large number of polymorphic sites that can be used to identify population-/breed-informative markers. In this study, starting from Illumina BovineSNP50 v1 BeadChip array genotyping data available from 3711 cattle of four breeds (2091 Italian Holstein, 738 Italian Brown, 475 Italian Simmental and 407 Marchigiana), principal component analysis (PCA) and random forests (RFs) were combined to identify informative SNP panels useful for cattle breed identification. From a PCA preselected list of 580 SNPs, RFs were computed using ranking methods (Mean Decrease in the Gini Index and Mean Accuracy Decrease) to identify the most informative 48 and 96 SNPs for breed assignment. The out-of-bag (OOB) error rate for both ranking methods and SNP densities ranged from 0.0 to 0.1% in the reference population. Application of this approach in a test population (10% of individuals pre-extracted from the whole data set) achieved 100% of correct assignment with both classifiers. Linkage disequilibrium between selected SNPs was relevant (r2 > 0.6) only in few pairs of markers indicating that most of the selected SNPs captured different fractions of variance. Several informative SNPs were in genes/QTL regions that affect or are associated with phenotypes or production traits that might differentiate the investigated breeds. The combination of PCA and RF to perform SNP selection and breed assignment can be easily implemented and is able to identify subsets of informative SNPs useful for population assignment starting from a large number of markers derived by high-throughput genotyping platforms.

SummaryIn this study, we investigated whether a selection programme based onboar genetic evaluation obtained with a classical BLUP animal model canchange allele frequencies in a pig population. All Italian Large White boarsborn from 1992 to 2012 with estimated breeding value reliability >0.85(n = 200) were selected among all boars of this breed. Boars were geno-typed with markers in major genes (IGF2 intron3-g.3072G>A, MC4Rp.D298N, VRTN PRE1 insertion, PRKAG3 p.I199V and FTO g.276T>G).Genotyping data were analysed grouping boars in eight classes accordingto their year of birth. To evaluate the influence of time on allele frequen-cies of the genotyped markers, multinomial logistic regression modelswere computed. Four of five polymorphic sites (IGF2, MC4R, VRTN andFTO) showed significant (p < 0.01) changes in allele frequencies over timedue to a progressive and continuous increase of one allele (associated withhigher lean meat content, higher average daily gain and favourable feed:gain ratio) and, consequently, decrease of the other one, following thedirectional selection of the selection programme of this pig breed. The ret-rospective analysis that was carried out in Italian Large White boars sug-gests that selection based on methodologies assuming the infinitesimalmodel is able to modify in a quite short period of time allele frequencies inmajor genes, increasing the frequency of alleles explaining a relevant(non-infinitesimal) fraction of the overall genetic variability for produc-tion traits.IntroductionSelection in livestock populations has been mainlybased on methodologies assuming the infinitesimalmodel that describes quantitative traits as determinedby an infinite number of genes, each of them with avery small effect. By this approach, genetic improve-ment is not expected to significantly change allelefrequencies at any particular locus in a population(Barton & Keightley 2002). However, experimentsdissecting genetic variability in livestock genomeshave substantially challenged the assumptions of theinfinitesimal model showing that some loci explain arelevant part of the genetic variance of quantitativetraits under selection. These loci have been termedquantitative trait loci (QTL) or major genes when thecausative genes and mutations are known and theireffects are particularly relevant. Therefore, directional

Background Carcass 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. Results We identified 4 single nucleotide polymorphisms with P≤5.0E-07 and additional 119 ones with 5.0E-07 Conclusions Further 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.

The expression of the HOX gene network in mid-stage human tooth development mostly concerns the epithelial tooth germ compartment and involves the C and D HOX loci. To further dissect the HOX gene implication with tooth epithelium differentiation we compared the expression of the whole HOX network in human ameloblastomas, as paradigm of epithelial odontogenic tumors, with tooth germs. We identified two ameloblastoma molecular types with respectively low and high number of active HOX C genes. The highly expressing HOX C gene ameloblastomas were characterized by a strong keratinized phenotype. Locus C HOX genes are located on chromosome 12q13–15 in physical contiguity with one of the two keratin gene clusters included in the human genome. The most posterior HOX C gene, HOX C13, is capable to interact with hair keratin genes located on the other keratin gene cluster in physical contiguity with the HOX B locus on chromosome 17q21– 22. Inside the HOX C locus, a 2.2 kb ncRNA (HOTAIR) able to repress transcription, in cis, along the entire HOX C locus and, in trans, at the posterior region of the HOX D locus has recently been identified. Interestingly both loci are deregulated in ameloblastomas. Our finding support an important role of the HOX network in characterizing the epithelial tooth compartment. Furthermore, the physical contiguity between locus C HOX and keratin genes in normal tooth epithelium and their deregulation in the neoplastic counterparts suggest they may act on the same mechanism potentially involved with epithelial tumorigenesis. J. Cell. Biochem. 112: 3206–3215, 2011.

Homeobox containing genes are a family of transcription factors regulating normal development and controlling primary cellular processes (cell identity, cell division and differentiation) recently enriched by the discovery of their interaction with miRNAs and ncRNAs. Class I human homeobox genes (HOX genes) are characterized by a unique genomic network organization: four compact chromosomal loci where 39 sequence corresponding genes can be aligned with each other in 13 antero-posterior paralogous groups. The cardiovascular system is the first mesoderm organ-system to be generated during embryonic development; subsequently it generates the blood and lymphatic vascular systems. Cardiovascular remodelling is involved through homeobox gene regulation and deregulation in adult physiology (menstrual cycle and wound healing) and pathology (atherosclerosis, arterial restenosis, tumour angiogenesis and lymphangiogenesis). Understanding the role played by homeobox genes in endothelial and smooth muscle cell phenotype determination will be crucial in identifying the molecular processes involved in vascular cell differentiation, as well as to support future therapeutic strategies. We report here on the current knowledge of the role played by homeobox genes in normal and abnormal vasculogenesis and postulate a common molecular mechanism accounting for the involvement of homeobox genes in the regulation of the nuclear export of specific transcripts potentially capable of generating endothelial phenotype modification involved in new vessel formation.

Homeobox-containing genes play a crucial role in odontogenesis. After the detection of Dlx and Msx genes in overlapping domains along maxillary and mandibular processes, a homeobox odontogenic code has been proposed to explain the interaction between different homeobox genes during dental lamina patterning. No role has so far been assigned to the Hox gene network in the homeobox odontogenic code due to studies on specific Hox genes and evolutionary considerations. Despite its involvement in early patterning during embryonal development, the HOX gene network, the most repeat-poor regions of the human genome, controls the phenotype identity of adult eukaryotic cells. Here, according to our results, the HOX gene network appears to be active in human tooth germs between 18 and 24 weeks of development. The immunohistochemical localization of specific HOX proteins mostly concerns the epithelial tooth germ compartment. Furthermore, only a few genes of the network are active in embryonal retromolar tissues, as well as in ectomesenchymal dental pulp cells (DPC) grown in vitro from adult human molar. Exposure of DPCs to cAMP induces the expression of from three to nine total HOX genes of the network in parallel with phenotype modifications with traits of neuronal differentiation. Our observations suggest that: (i) by combining its component genes, the HOX gene network determines the phenotype identity of epithelial and ectomesenchymal cells interacting in the generation of human tooth germ; (ii) cAMP treatment activates the HOX network and induces, in parallel, a neuronal-like phenotype in human primary ectomesenchymal dental pulp cells.

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 (Krškopolje 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.