Application of Microsatellites to Trace the Dairy Products Back to the Farm of Origin

The increasing number of food frauds, mainly targeting high quality products, is a rising concern among producers and authorities appointed to food controls. Therefore, the development or implementation of methods to reveal frauds is desired. The genetic traceability of traditional or high-quality dairy products (i.e., products of protected designation of origin, PDO) represents a challenging issue due to the technical problems that arise. The aim of the study was to set up a genetic tool for the origin traceability of dairy products. We investigated the use of Short Tandem Repeats (STRs) to assign milk and cheese to the corresponding producer. Two farms were included in the study, and the blood of the cows, bulk milk, and derived cheese were sampled monthly for one year. Twenty STRs were selected and Polymerase Chain Reactions for each locus were carried out. The results showed that bulk milk and derived cheese express an STR profile composed of a subset of STRs of the lactating animals. A bioinformatics tool was used for the exclusion analysis. The study allowed the identification of a panel of 20 markers useful for the traceability of milk and cheeses, and its effectiveness in the traceability of dairy products obtained from small producers was demonstrated.

Genetic and Phenotypic Characteristics of Belted Pig Breeds: A Review

Belted pig breeds have unique, distinguishing phenotypic characteristics. This review summarises the current knowledge on pig breeds displaying a belted coat pattern. Belts of different widths and positions around the animal's trunk characterise specific pig breeds from all around the world. All the breeds included in the present paper have been searched through the FAO domestic animal diversity information system (DAD-IS), Every country was checked to identify all breeds described as having black or red piebald coat pattern variations. Advances in genomic technologies have made it possible to identify the specific genes and genetic markers associated with the belted phenotype and explore the genetic relationships between different local breeds. Thus, the origin, history, and production traits of these breeds, together with all the genomic information related to the mechanism of skin pigmentation, are discussed. By increasing our understanding of these breeds, we can appreciate the richness of our biological and cultural heritage and work to preserve the biodiversity of the world's animals.

Genomic regions underlying positive selection in local, Alpine cattle breeds

We used genome-wide SNP data from 18 local cattle breeds from six countries of the Alpine region to characterize population structure and identify genomic regions underlying positive selection. The geographically close breeds Evolèner, Eringer, Valdostana Pezzata Nera, and Valdostana Castana were found to differ from all other Alpine breeds. In addition, three breeds, Simmental, and Original Braunvieh from Switzerland and Pinzgauer from Austria built three separate clusters. Of the 18 breeds studied, the intra-alpine Swiss breed Evolèner had the highest average inbreeding based on runs of homozygosity (F_ROH ) and the highest average genomic relationship within the breed. In contrast, Slovenian Cika cattle had the lowest average genomic inbreeding and the lowest average genomic relationship within the breed. We found selection signatures on chromosome 6 near known genes such as KIT and LCORL explaining variation in coat color and body size in cattle. The most prominent selection signatures were similar regardless of marker density and the breeds in the data set. In addition, using available high-density SNP data from 14 of the breeds we identified 47 genome regions as ROH islands. The proportion of homozygous animals was higher in all studied animals of local breeds than in Holstein and Brown Swiss cattle, the two most important commercial breeds in the Alpine region. We report ROH islands near genes related to thermoregulation, coat color, production, and stature. The results of this study serve as a basis for the search for causal variants underlying adaptation to the alpine environment and other specific characteristics selected during the evolution of local Alpine cattle breeds.

MicroRNA-based signatures impacting clinical course and biology of ovarian cancer: a miRNOmics study

Background

In Western countries, ovarian cancer (OC) still represents the leading cause of gynecological cancer-related deaths, despite the remarkable gains in therapeutical options. Novel biomarkers of early diagnosis, prognosis definition and prediction of treatment outcomes are of pivotal importance. Prior studies have shown the potentials of micro-ribonucleic acids (miRNAs) as biomarkers for OC and other cancers.

Methods

We focused on the prognostic and/or predictive potential of miRNAs in OC by conducting a comprehensive array profiling of miRNA expression levels in ovarian tissue samples from 17 non-neoplastic controls, and 60 tumor samples from OC patients treated at the Regina Elena National Cancer Institute (IRE). A set of 54 miRNAs with differential expression in tumor versus normal samples (T/N-deregulated) was identified in the IRE cohort and validated against data from the Cancer Genoma Atlas (TCGA) related to 563 OC patients and 8 non-neoplastic controls. The prognostic/predictive role of the selected 54 biomarkers was tested in reference to survival endpoints and platinum resistance (P-res).

Results

In the IRE cohort, downregulation of the 2 miRNA-signature including miR-99a-5p and miR-320a held a negative prognostic relevance, while upregulation of miR-224-5p was predictive of less favorable event free survival (EFS) and P-res. Data from the TCGA showed that downregulation of 5 miRNAs, i.e., miR-150, miR-30d, miR-342, miR-424, and miR-502, was associated with more favorable EFS and overall survival outcomes, while miR-200a upregulation was predictive of P-res. The 9 miRNAs globally identified were all included into a single biologic signature, which was tested in enrichment analysis using predicted/validated miRNA target genes, followed by network representation of the miRNA-mRNA interactions.

Conclusions

Specific dysregulated microRNA sets in tumor tissue showed predictive/prognostic value in OC, and resulted in a promising biological signature for this disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40364-021-00289-6.

Genomic Analyses of Unveil Helmeted Guinea Fowl (Numida meleagris) Domestication in West Africa

Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300-5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.

A genome-wide scan of copy number variants in three Iranian indigenous river buffaloes

BACKGROUND

In Iran, river buffalo is of great importance. It plays an important role in the economy of the Country, because its adaptation to harsh climate conditions and long productive lifespan permitting its farming across the Country and to convert low-quality feed into valuable milk. The genetic variability in Iranian buffalo breeds have been recently studied using SNPs genotyping data, but a whole genome Copy Number Variants (CNVs) mapping was not available. The aim of this study was to perform a genome wide CNV scan in 361 buffaloes of the three Iranian river breeds (Azeri, Khuzestani and Mazandarani) through the analysis of data obtained using the Axiom® Buffalo Genotyping Array 90 K.

RESULTS

CNVs detection resulted in a total of 9550 CNVs and 302 CNVRs identified in at least 5% of samples within breed, covering around 1.97% of the buffalo genome. and A total of 22 CNVRs were identified in all breeds and a different proportion of regions were in common among the three populations. Within the more represented CNVRs (n = 302) mapped a total of 409 buffalo genes, some of which resulted associated with morphological, healthy, milk, meat and reproductive traits, according to Animal Genome Cattle database.

CONCLUSIONS

This work provides a step forward in the interpretation of genomic variation within and among the buffalo populations, releasing a first map of CNVs and providing insights about their recent selection and adaptation to environment. The presence of the set of genes and QTL traits harbored in the CNVRs could be possibly linked with the buffalo's natural adaptive history together to a recent selection for milk used as primary food source from this species.

The Genomic Variation in the Aosta Cattle Breeds Raised in an Extensive Alpine Farming System

The Aosta Red Pied (Valdostana Pezzata Rossa (VRP)), the Aosta Black Pied (Valdostana Pezzata Nera (VBP)) and the Aosta Chestnut (Valdostana Castana (CAS)) are dual-purpose cattle breeds (meat and milk), very well adapted to the harsh environmental conditions of alpine territories: their farming is in fact characterized by summer pasture at very high altitude. A total of 728 individuals were genotyped with the GeenSeek Genomic Profiler® (GGP) Bovine 150K Illumina SNP chip as a part of the DUALBREEDING-PSRN Italian-funded research project. The genetic diversity among populations showed that the three breeds are distinct populations based on the F_ST values, ADMIXTURE and Principal Component Analysis (PCA) results. Runs of Homozygosity (ROH) were obtained for the three populations to disclose recent autozygosity. The genomic inbreeding based on the ROH was calculated and coupled with information derived from the F (inbreeding coefficient) and F_ST parameters. The mean F_ROH values were low: CAS = 0.06, VBP = 0.05 and VRP = 0.07, while the average F values were -0.003, -0.01 and -0.003, respectively. The annotation and enrichment analysis, performed in the identified most frequent ROH (TOP_ROH), showed genes that can be linked to the resilience capacity of these populations to harsh environmental farming conditions, and to the peculiar characteristics searched for by farmers in each breed.

The German Shorthair Pointer Dog Breed (Canis lupus familiaris): Genomic Inbreeding and Variability

The German Shorthaired Pointer (GSHP) is a breed worldwide known for its hunting versatility. Dogs of this breed are appreciated as valuable companions, effective trackers, field trailers and obedience athletes. The aim of the present work is to describe the genomic architecture of the GSHP breed and to analyze inbreeding levels under a genomic and a genealogic perspective. A total of 34 samples were collected (24 Italian, 10 USA), and the genomic and pedigree coefficients of inbreeding have been calculated. A total of 3183 runs of homozygosity (ROH) across all 34 dogs have been identified. The minimum and maximum number of Single Nucleotide Polymorphisms (SNPs) defining all ROH are 40 and 3060. The mean number of ROH for the sample was 93.6. ROH were found on all chromosomes. A total of 854 SNPs (TOP_SNPs) defined 11 ROH island regions (TOP_ROH), in which some gene already associated with behavioral and morphological canine traits was annotated. The proportion of averaged observed homozygotes estimated on total number of SNPs was 0.70. The genomic inbreeding coefficient based on ROH was 0.17. The mean inbreeding based on genealogical information resulted 0.023. The results describe a low inbred population with quite a good level of genetic variability.

Canine fertility: The consequences of selection for special traits

Pedigree dogs and cats are bred aiming to conform breed standards with very poor consideration for breeding stock fertility. At the same time, the genetic asset underlining reproductive traits could be effectively analysed like in other species under selection. The definition of selection targets is very important in breeding protocols determination. The aim of the present work is to present an overview of the different correlations between reproduction and genetics, starting from selection procedure and inbreeding coefficient moving to genomic and the application of SNPs and GWAS on population study and identification of genes involved in phenotypical variation of reproductive traits in dogs. Particular relevance has been given to the concept of inbreeding which effects on canine reproduction have been presented. The use of genomic information in inbreeding coefficient calculation can be considered an improved effective procedure in the evaluation of the genetic variability loss in canine population and its negative effects on reproductive traits.

Copy Number Variation Mapping and Genomic Variation of Autochthonous and Commercial Turkey Populations

This study aims at investigating genomic diversity of several turkey populations using Copy Number Variants (CNVs). A total of 115 individuals from six Italian breeds (Colle Euganei, Bronzato Comune Italiano, Parma e Piacenza, Brianzolo, Nero d'Italia, and Ermellinato di Rovigo), seven Narragansett, 38 commercial hybrids, and 30 Mexican turkeys, were genotyped with the Affymetrix 600K single nucleotide polymorphism (SNP) turkey array. The CNV calling was performed with the Hidden Markov Model of PennCNV software and with the Copy Number Analysis Module of SVS 8.4 by Golden Helix^®. CNV were summarized into CNV regions (CNVRs) at population level using BEDTools. Variability among populations has been addressed by hierarchical clustering (pvclust R package) and by principal component analysis (PCA). A total of 2,987 CNVs were identified covering 4.65% of the autosomes of the Turkey_5.0/melGal5 assembly. The CNVRs identified in at least two individuals were 362-189 gains, 116 losses, and 57 complexes. Among these regions the 51% contain annotated genes. This study is the first CNV mapping of turkey population using 600K chip. CNVs clustered the individuals according to population and their geographical origin. CNVs are known to be indicators also of adaptation, as some researches in different species are suggesting.

Genome-Wide Association Study in Mexican Holstein Cattle Reveals Novel Quantitative Trait Loci Regions and Confirms Mapped Loci for Resistance to Bovine Tuberculosis

Simple Summary

Bovine tuberculosis is an infectious disease of cattle caused by Mycobacterium bovis characterized by the formation of tubercles in any organ or tissue. Bovine tuberculosis represents a significant veterinary and public health problem in many parts of the world. It is zoonotic, transmitted to humans through consumption of infected milk and other cattle products. Although many factors influence infection and progression of the disease, there must be an important host genetic component that explains why some animals get sick and others remain healty. We present evidence of genetic variants associated with resistance to tuberculosis in Mexican Holstein dairy cattle using a case-control approach with a selective DNA pooling. Here, we identified novel quantitative trait loci regions harboring genes involved in Mycobacterium spp. immune response. This is a first screening about resistance to tuberculosis infection on Mexican Holstein cattle based on a dense single nucleotide polymorphism chip. The identified genes belong to both, the already known, and the undisclosed quantitative trait loci regions.

Abstract

Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

Genome-wide association study between CNVs and milk production traits in Valle del Belice sheep

Copy number variation (CNV) is a major source of genomic structural variation. The aim of this study was to detect genomic CNV regions (CNVR) in Valle del Belice dairy sheep population and to identify those affecting milk production traits. The GO analysis identified possible candidate genes and pathways related to the selected traits. We identified CNVs in 416 individuals genotyped using the Illumina OvineSNP50 BeadChip array. The CNV association using a correlation-trend test model was examined with the Golden Helix SVS 8.7.0 tool. Significant CNVs were detected when their adjusted p-value was <0.01 after false discovery rate (FDR) correction. We identified 7,208 CNVs, which gave 365 CNVRs after aggregating overlapping CNVs. Thirty-one CNVRs were significantly associated with one or more traits included in the analysis. All CNVRs, except those on OAR19, overlapped with quantitative trait loci (QTL), even if they were not directly related to the traits of interest. A total of 222 genes were annotated within the significantly associated CNVRs, most of which played important roles in biological processes related to milk production and health-related traits. Identification of the genes in the CNVRs associated with the studied traits will provide the basis for further investigation of their role in the metabolic pathways related to milk production and health traits.

First genome-wide CNV mapping in FELIS CATUS using next generation sequencing data

Background

Copy Number Variations (CNVs) have becoming very significant variants, representing a major source of genomic variation. CNVs involvement in phenotypic expression and different diseases has been widely demonstrated in humans as well as in many domestic animals. However, genome wide investigation on these structural variations is still missing in Felis catus. The present work is the first CNV mapping from a large data set of Next Generation Sequencing (NGS) data in the domestic cat, performed within the 99 Lives Consortium.

Results

Reads have been mapped on the reference assembly_6.2 by Maverix Biomics. CNV detection with cn.MOPS and CNVnator detected 592 CNVs. These CNVs were used to obtain 154 CNV Regions (CNVRs) with BedTools, including 62 singletons. CNVRs covered 0.26% of the total cat genome with 129 losses, 19 gains and 6 complexes. Cluster Analysis and Principal Component Analysis of the detected CNVRs showed that breeds tend to cluster together as well as cats sharing the same geographical origins. The 46 genes identified within the CNVRs were annotated.

Conclusion

This study has improved the genomic characterization of 14 cat breeds and has provided CNVs information that can be used for studies of traits in cats. It can be considered a sound starting point for genomic CNVs identification in this species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5297-2) contains supplementary material, which is available to authorized users.

A copy number variant scan in the autochthonous Valdostana Red Pied cattle breed and comparison with specialized dairy populations

Copy number variants (CNVs) are an important source of genomic structural variation, recognized to influence phenotypic variation in many species. Many studies have focused on identifying CNVs within and between human and livestock populations alike, but only few have explored population-genetic properties in cattle based on CNVs derived from a high-density SNP array. We report a high-resolution CNV scan using Illumina’s 777k BovineHD Beadchip for Valdostana Red Pied (VRP), an autochthonous Italian dual-purpose cattle population reared in the Alps that did not undergo strong selection for production traits. After stringent quality control and filtering, CNVs were called across 108 bulls using the PennCNV software. A total of 6,784 CNVs were identified, summarized to 1,723 CNV regions (CNVRs) on 29 autosomes covering a total of ~59 Mb of the UMD3.1 assembly. Among the mapped CNVRs, there were 812 losses, 832 gains and 79 complexes. We subsequently performed a comparison of CNVs detected in the VRP and those available from published studies in the Italian Brown Swiss (IBS) and Mexican Holstein (HOL). A total of 171 CNVRs were common to all three breeds. Between VRP and IBS, 474 regions overlapped, while only 313 overlapped between VRP and HOL, indicating a more similar genetic background among populations with common origins, i.e. the Alps. The principal component, clustering and admixture analyses showed a clear separation of the three breeds into three distinct clusters. In order to describe the distribution of CNVs within and among breeds we used the pair V_ST statistic, considering only the CNVRs shared to more than 5 individuals (within breed). We identified unique and highly differentiated CNVs (n = 33), some of which could be due to specific breed selection and adaptation. Genes and QTL within these regions were characterized.

Genetic variability of Akhal-Teke horses bred in Italy

Background

The Akhal-Teke horse (AKH) is native of the modern Turkmenistan area. It was introduced in Italy from 1991 to 2000 mainly as an endurance horse. This paper characterizes the genetic variability of the whole Italian AKH horse population and evaluates their inbreeding level by analyzing microsatellite markers and mitochondrial D-Loop sequences.

Methods

Seventeen microsatellite marker loci were genotyped on 95 DNA samples from almost all the AKH horses bred in Italy in the last 20 years. Standard genetic variability measures (H_o, H_e, F_IS) were compared against the same variables published on other eight AKH populations. In addition, 397 bp of mtDNA D-loop region were sequenced on a sub-group of 22 unrelated AKH out of the 95 sampled ones, and on 11 unrelated Arab horses. The haplotypes identified in the Italian population were aligned to sequences of AKH (56), Arab (five), Caspian Pony (13), Przewalskii (two) and Barb (15) horses available in GenBank. The Median Joining Network (MJN), Principal Component Analysis (PCA) and Neighbor-joining (NJ) tree were calculated on the total 126 sequences.

Results

Nucleic markers showed a high degree of polymorphism (H_o = 0.642; H_e = 0.649) and a low inbreeding level (F_IS = 0.016) in Italian horses, compared to other AKH populations (ranged from −0.103 AKH from Estonia to 0.114 AKH from Czech Republic). High variability was also recorded in the D-Loop region. 11 haplotypes were identified with haplotype diversity (hd), nucleotide diversity (π) and average number of nucleotide differences (k) of 0.938, 0.021 and 6.448, respectively. When all the 126 D-Loop sequences were compared, 51 haplotypes were found, and four were here found only in the Italian AKH horses. The 51 haplotypes were conformed to eight recognized mtDNA haplogroups (A, C, F, G, L, M, P and Q) and confirmed by MJN analysis, Italian horses being assigned to five haplogroups (A, C, G, L and M). Using a PCA approach to the same data, the total haplotypes were grouped into two clusters including A+C+M+P and G+F haplogroups, while L and Q haplogroups remained ungrouped. Finally, the NJ algorithm effectively discretizes only the L haplogroup. All the above data univocally indicate good genetic variability and accurate management of the Akhal-Teke population in Italy.

Looking at genetic structure and selection signatures of the Mexican chicken population using single nucleotide polymorphism markers

Genetic variation enables both adaptive evolutionary changes and artificial selection. Genetic makeup of populations is the result of a long-term process of selection and adaptation to specific environments and ecosystems. The aim of this study was to characterize the genetic variability of México's chicken population to reveal any underlying population structure. A total of 213 chickens were sampled in different rural production units located in 25 states of México. Genotypes were obtained using the Affymetrix Axiom® 600 K Chicken Genotyping Array. The Identity by Descent (IBD) and the principal components analysis (PCA) were performed by SVS software on pruned single nucleotide polymorphisms (SNPs).ADMIXTURE analyses identified 3 ancestors and the proportion of the genetic contribution of each of them has been determined in each individual. The results of the Neighbor-Joining (NJ) analysis resulted consistent with those obtained by the PCA. All methods utilized in this study did not allow a classification of Mexican chicken in distinct clusters or groups. A total of 3,059 run of homozygosity (ROH) were identified and, being mainly short in length (<4 Mb), these regions are indicative of a low inbreeding level in the population. Finally, findings from the ROH analysis indicated the presence of natural selective pressure in the population of Mexican chicken.The study indicates that the Mexican chicken clearly appear to be a unique creole chicken population that was not subjected to a specific artificial selection. Results provide a genetic knowledge that can be used as a basis for the genetic management of a unique and very large creole population, especially in the view of using it in production of hybrids to increase the productivity and economic revenue of family farming agriculture, which is widely present in México.

Conservation status and historical relatedness of Italian cattle breeds

Background

In the last 50 years, the diversity of cattle breeds has experienced a severe contraction. However, in spite of the growing diffusion of cosmopolite specialized breeds, several local cattle breeds are still farmed in Italy. Genetic characterization of breeds represents an essential step to guide decisions in the management of farm animal genetic resources. The aim of this work was to provide a high-resolution representation of the genome-wide diversity and population structure of Italian local cattle breeds using a medium-density single nucleotide polymorphism (SNP) array.

Results

After quality control filtering, the dataset included 31,013 SNPs for 800 samples from 32 breeds. Our results on the genetic diversity of these breeds agree largely with their recorded history. We observed a low level of genetic diversity, which together with the small size of the effective populations, confirmed that several breeds are threatened with extinction. According to the analysis of runs of homozygosity, evidence of recent inbreeding was strong in some local breeds, such as Garfagnina, Mucca Pisana and Pontremolese. Patterns of genetic differentiation, shared ancestry, admixture events, and the phylogenetic tree, all suggest the presence of gene flow, in particular among breeds that originate from the same geographical area, such as the Sicilian breeds. In spite of the complex admixture events that most Italian cattle breeds have experienced, they have preserved distinctive characteristics and can be clearly discriminated, which is probably due to differences in genetic origin, environment, genetic isolation and inbreeding.

Conclusions

This study is the first exhaustive genome-wide analysis of the diversity of Italian cattle breeds. The results are of significant importance because they will help design and implement conservation strategies. Indeed, efforts to maintain genetic diversity in these breeds are needed. Improvement of systems to record and monitor inbreeding in these breeds may contribute to their in situ conservation and, in view of this, the availability of genomic data is a fundamental resource.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0406-x) contains supplementary material, which is available to authorized users.

Evaluación morfológica de gallinas de traspatio mexicanas (Gallus gallus domesticus)

El objetivo fue evaluar algunas características morfológicas de gallinas de traspatio (n= 255) provenientes de 65 unidades rurales de producción, localizadas en 52 municipios de 18 estados de la República Mexicana. El modelo estadístico incluyó los efectos de sexo, estado y municipio anidado en estado. El efecto de sexo fue importante para todas las variables analizadas (P<0.0001), excepto para robustez (P=0.33). El efecto de estado fue importante (P<0.0001) para todas las variables analizadas, mientras que municipio solo influyó significativamente a envergadura (P<0.0001), peso corporal (P=0.0271) y solidez (P=0.0267). Los machos tuvieron mayor (P<0.0001) longitud corporal (4.85 cm), envergadura (5.66 cm), circunferencia de la pechuga (2.94 cm), longitud del tarso (1.73 cm), peso corporal (0.65 kg) y solidez (0.99 puntos porcentuales). El peso corporal aumentó 126 g (P<0.01) en machos y 61 g en hembras (P<0.01) por cada centímetro que aumentó la circunferencia de la pechuga. El peso corporal mostró estar moderadamente correlacionado con longitud del cuerpo (r= 0.65), envergadura (r= 0.49) y longitud del tarso (r= 0.67), y altamente correlacionada con circunferencia de la pechuga (r= 0.76) en machos. La gallina de traspatio en México presentó una diferenciación morfológica importante entre machos y hembras para las características estudiadas, excepto para robustez, en la cual no se detectaron diferencias.

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.

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals

Stature is affected by many polymorphisms of small effect in humans ¹ . In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes^2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10^-8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP-seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals.

Short communication: Genomic prediction using imputed whole-genome sequence variants in Brown Swiss Cattle

The accuracy of genomic prediction determines response to selection. It has been hypothesized that accuracy of genomic breeding values can be increased by a higher density of variants. We used imputed whole-genome sequence data and various single nucleotide polymorphism (SNP) selection criteria to estimate genomic breeding values in Brown Swiss cattle. The extreme scenarios were 50K SNP chip data and whole-genome sequence data with intermediate scenarios using linkage disequilibrium-pruned whole-genome sequence variants, only variants predicted to be missense, or the top 50K variants from genome-wide association studies. We estimated genomic breeding values for 3 traits (somatic cell score, nonreturn rate in heifers, and stature) and found differences in accuracy levels between traits. However, among different SNP sets, accuracy was very similar. In our analyses, sequence data led to a marginal increase in accuracy for 1 trait and was lower than 50K for the other traits. We concluded that the inclusion of imputed whole-genome sequence data does not lead to increased accuracy of genomic prediction with the methods.

Genomic variability in Mexican chicken population using copy number variants

BACKGROUND

Copy number variations are genome polymorphism that influence phenotypic variation and are an important source of genetic variation in populations. The aim of this study was to investigate genetic variability in the Mexican Creole chicken population using CNVs.

RESULTS

The Hidden Markov Model of the PennCNV software detected a total of 1924 CNVs in the genome of the 256 samples processed with Axiom® Genome-Wide Chicken Genotyping Array (Affymetrix). The mapped CNVs comprised 1538 gains and 386 losses, resulting at population level in 1216 CNV regions (CNVRs), of which 959 gains, 226 losses and 31 complex (i.e. containing both losses and gains). The CNVRs covered a total of 47 Mb of the whole genome sequence length, corresponding to 5.12% of the chicken galGal4 autosome assembly.

CONCLUSIONS

This study allowed a deep insight into the structural variation in the genome of unselected Mexican chicken population, which up to now has not been genetically characterized. The genomic study disclosed that the population, even if presenting extreme morphological variation, cannot be organized in differentiated genetic subpopulations. Finally this study provides a chicken CNV map based on the 600 K SNP chip array jointly with a genome-wide gene copy number estimates in a native unselected for more than 500 years chicken population.

Mitochondrial DNA genetic diversity in six Italian donkey breeds (Equus asinus)

Donkeys have played an important role in agricultural land practices and in human historical periods of recent past and, still today, are used as a working power in several world areas. The objective of this study was to identify genetic variability in six Italian donkey breeds using mtDNA D-loop. Fifteen haplotypes, grouped in three haplogroups, were identified. The genetic indices were informative and showed a high population genetic variability. The results of AMOVA analyses based on geographic structuring of Italian populations highlighted that the majority of the observed variance is due to differences among samples within breeds. Comparison among Italian haplotypes and mtDNA D-loop sequences belonging to European domestic and Ethiopian donkeys and wild asses, clearly define two clades referred to Nubian lineage. The results can be useful to complement safeguard planes for donkey breeds that are considered to extinction endangered.

Microvesicles secreted from equine amniotic-derived cells and their potential role in reducing inflammation in endometrial cells in an in-vitro model

Background

It is known that a paracrine mechanism exists between mesenchymal stem cells and target cells. This process may involve microvesicles (MVs) as an integral component of cell-to-cell communication.

Methods

In this context, this study aims to understand the efficacy of MVs in in-vitro endometrial stressed cells in view of potential healing in in-vivo studies. For this purpose, the presence and type of MVs secreted by amniotic mesenchymal stem cells (AMCs) were investigated and the response of endometrial cells to MVs was studied using a dose-response curve at different concentrations and times. Moreover, the ability of MVs to counteract the in vitro stress in endometrial cells induced by lipopolysaccharide was studied by measuring the rate of apoptosis and cell proliferation, the expression of some pro-inflammatory genes such as tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), interleukin 1β (IL-1β), and metalloproteinases (MMP) 1 and 13, and the release of some pro- or anti-inflammatory cytokines.

Results

MVs secreted by the AMCs ranged in size from 100 to 200 nm. The incorporation of MVs was gradual over time and peaked at 72 h. MVs reduced the apoptosis rate, increased cell proliferation values, downregulated pro-inflammatory gene expression, and decreased the secretion of pro-inflammatory cytokines.

Conclusion

Our data suggest that some microRNAs could contribute to counteracting in-vivo inflammation of endometrial tissue.

Genome-wide association study for milk somatic cell score in holstein cattle using copy number variation as markers

Mastitis, the most common and expensive disease in dairy cows, implies significant losses in the dairy industry worldwide. Many efforts have been made to improve genetic mastitis resistance in dairy populations, but low heritability of this trait made this process not as effective as desired. The purpose of this study was to identify genomic regions explaining genetic variation of somatic cell count using copy number variations (CNVs) as markers in the Holstein population, genotyped with the Illumina BovineHD BeadChip. We found 24 and 47 copy number variation regions significantly associated with estimated breeding values for somatic cell score (SCS_EBVs) using SVS 8.3.1 and PennCNV-CNVRuler software, respectively. The association analysis performed with these two software allowed the identification of 18 candidate genes (TERT, NOTCH1, SLC6A3, CLPTM1L, PPARα, BCL-2, ABO, VAV2, CACNA1S, TRAF2, RELA, ELF3, DBH, CDK5, NF2, FASN, EWSR1 and MAP3K11) that result classified in the same functional cluster. These genes are also part of two gene networks, whose genes share the 'stress', 'cell death', 'inflammation' and 'immune response' GO terms. Combining CNV detection/association analysis based on two different algorithms helps towards a more complete identification of genes linked to phenotypic variation of the somatic cell count.

Análisis de variables morfológicas de pavos de traspatio mexicanos (Meleagris gallopavo gallopavo)

El objetivo fue evaluar algunas características morfológicas de pavos de traspatio (n=248) provenientes de 126 unidades rurales de producción localizadas en 75 municipios de 24 estados de la República Mexicana. El modelo estadístico incluyó sexo, estado y municipio anidado en estado. Las tres variables explicativas afectaron a todas las variables de respuesta (P<0.01), excepto municipio, que no afectó circunferencia de la pechuga (P>0.05). Los machos tuvieron mayor (P<0.001) longitud corporal (10.4 cm más), envergadura (11.4 cm más), circunferencia de la pechuga (13.8 cm más), longitud de tarso (2.5 cm más), peso corporal (2.5 kg más), robustez (9.0 puntos porcentuales más) y solidez (2.8 puntos porcentuales más) que las hembras. Peso corporal mostró estar altamente correlacionado fenotípicamente (P<0.01) con circunferencia de la pechuga tanto en machos (r=0.74) como en hembras (r=0.71). En machos, longitud corporal mostró una correlación baja con longitud del tarso (r=0.25; P<0.01), pero en hembras no estuvo correlacionada (r=0.05; P>0.05). El peso corporal aumentó 143 g (P<0.01) en machos y 113 g en hembras (P<0.01) por cada centímetro que aumentó la circunferencia de la pechuga. Los colores predominantes en el plumaje, la piel y el tarso fueron negro, blanco y café, respectivamente. El pavo de traspatio mexicano presentó dimorfismo sexual significativo y alta correlación fenotípica entre circunferencia de la pechuga y peso corporal.

The Use of Kosher Phenotyping for Mapping QTL Affecting Susceptibility to Bovine Respiratory Disease

Bovine respiratory disease (BRD) is the leading cause of morbidity and mortality in feedlot cattle, caused by multiple pathogens that become more virulent in response to stress. As clinical signs often go undetected and various preventive strategies failed, identification of genes affecting BRD is essential for selection for resistance. Selective DNA pooling (SDP) was applied in a genome wide association study (GWAS) to map BRD QTLs in Israeli Holstein male calves. Kosher scoring of lung adhesions was used to allocate 122 and 62 animals to High (Glatt Kosher) and Low (Non-Kosher) resistant groups, respectively. Genotyping was performed using the Illumina BovineHD BeadChip according to the Infinium protocol. Moving average of -logP was used to map QTLs and Log drop was used to define their boundaries (QTLRs). The combined procedure was efficient for high resolution mapping. Nineteen QTLRs distributed over 13 autosomes were found, some overlapping previous studies. The QTLRs contain polymorphic functional and expression candidate genes to affect kosher status, with putative immunological and wound healing activities. Kosher phenotyping was shown to be a reliable means to map QTLs affecting BRD morbidity.

191 SEARCHING FOR THE IN VIVO TRANSCRIPTOME BLUEPRINT OF COMPETENT BOVINE OOCYTES

Gene expression in early stage embryos relies mostly on post-transcriptional control of maternal transcripts accumulated during oocyte maturation. However, while the building process to obtain a competent oocyte is now better understood, it is still not clear what transcriptome blueprint composes a competent oocyte. The aim of the study was to compare the mRNA expression pattern between oocytes collected from fertile heifers and repeat breeders by using RNAseq. Oocytes were collected by ovum pickup from 3 heifers that were 11–15 months of age and became pregnant at the following oestrus and from 4 adult cows with an age of 4 to 7 years, classified as repeat breeders after they failed to become pregnant for a minimum of 3 consecutive AI. To obtain oocytes from follicles with the same degree of development, at time 0 all follicles visible through transrectal ultrasound examination were removed by transvaginal aspiration. Five days later oocytes were collected by ovum pick up from the newly formed follicles with diameters >5 mm. Oocytes from each animal were pooled and sequenced as a single sample. Total RNA was extracted by RNeasy Micro Kit (Qiagen, Valencia, CA, USA). Amplified cDNA, was prepared starting from total RNA using the Ovation RNA-Seq System V2 (Nugen Technologies, San Carlos, CA, USA). After library preparation with TruSeq DNA Sample Prep kit (Illumina, Madison, WI, USA), sequencing was performed on an Illumina HiSEqn 2000. Galaxy and Chipster open web-based platforms were used to analyse the data. We identified 49 differentially expressed genes. Heifers’ oocytes mRNA pattern indicated greater potential to sustain cell division. In particular, oocytes expressed more Keratin 14 (a gene involved in cell proliferation) and kinesin family member 20B (a protein involved in cytokinesis). More competent oocytes also have a greater ability to repair single-strand breaks due to the high levels of endo/exonuclease (5′-3′), endonuclease G-like. This may reflect greater capacity to neutralise DNA damage and, therefore, greater ability to preserve and transmit high-quality DNA. Repeat breeders portray a different landscape; their greater expression of Jun oncogene, Heat shock protein 1, Stimulated by retinoic acid gene6, arylhydrocarbon receptor nuclear translocator, fibromodulin, and aryl-hydrocarbon receptor repressor suggest that these oocytes have been subjected to environmental stress during oocyte maturation. Their greater expressions of inhibin α, stearoyl-CoA desaturase, junctional adhesion molecule 2 have been previously shown to correlate with a reduced oocyte developmental potential. Furthermore, the Cannabinoid receptor protein 1 expression suggests a compromised ion function that can lead to a failed activation of the development program. Finally, the greater expression of Ubiquilin3 and Heat shock protein 1 led to high protein and mRNA degradation, respectively, suggesting that these oocytes are deprived of essential components to sustain embryo growth. In conclusion our data provide the first detailed snapshot of the mRNA pattern defining the differences between a competent oocyte and an incompetent oocyte in vivo. Study supported by PRIN 2008, 2009 and EU-Quantomics.

Endothelin A receptor drives invadopodia function and cell motility through the β-arrestin/PDZ-RhoGEF pathway in ovarian carcinoma

The endothelin-1 (ET-1)/ET A receptor (ETAR) signalling pathway is a well-established driver of epithelial ovarian cancer (EOC) progression. One key process promoted by ET-1 is tumor cell invasion, which requires the scaffolding functions of β-arrestin-1 (β-arr1) downstream of the receptor; however, the potential role of ET-1 in inducing invadopodia, which are crucial for cellular invasion and tumor metastasis, is completely unknown. We describe here that ET-1/ETAR, through β-arr1, activates RhoA and RhoC GTPase and downstream ROCK (Rho-associated coiled coil-forming kinase) kinase activity, promoting actin-based dynamic remodelling and enhanced cell invasion. This is accomplished by the direct interaction of β-arr1 with PDZ-RhoGEF (postsynaptic density protein 95/disc-large/zonula occludens-RhoGEF). Interestingly, ETAR-mediated invasive properties are related to the regulation of invadopodia, as evaluated by colocalization of actin with cortactin, as well as with TKS5 and MT1-MMP (membrane type 1-matrix metalloproteinase) with areas of matrix degradation, and activation of cofilin pathway, which is crucial for regulating invadopodia activity. Depletion of PDZ-RhoGEF, or β-arr1, or RhoC, as well as the treatment with the dual ET-1 receptor antagonist macitentan, significantly impairs invadopodia function, MMP activity and invasion, demonstrating that β-arr1/PDZ-RhoGEF interaction mediates ETAR-driven ROCK-LIMK-cofilin pathway through the control of RhoC activity. In vivo, macitentan is able to inhibit metastatic dissemination and cofilin phosphorylation. Collectively, our data unveil a noncanonical activation of the RhoC/ROCK pathway through the β-arr1/PDZ-RhoGEF complex as a regulator of ETAR-induced motility and metastasis, establishing ET-1 axis as a novel regulator of invadopodia protrusions through the RhoC/ROCK/LIMK/cofilin pathway during the initial steps of EOC invasion.

Genome-wide association study for somatic cell score in Valdostana Red Pied cattle breed using pooled DNA

BACKGROUND

Mastitis is a major disease of dairy cattle occurring in response to environmental exposure to infective agents with a great economic impact on dairy industry. Somatic cell count (SCC) and its log transformation in somatic cell score (SCS) are traits that have been used as indirect measures of resistance to mastitis for decades in selective breeding. A selective DNA pooling (SDP) approach was applied to identify Quantitative Trait Loci (QTL) for SCS in Valdostana Red Pied cattle using the Illumina Bovine HD BeadChip.

RESULTS

A total of 171 SNPs reached the genome-wide significance for association with SCS. Fifty-two SNPs were annotated within genes, some of those involved in the immune response to mastitis. On BTAs 1, 2, 3, 4, 9, 13, 15, 17, 21 and 22 the largest number of markers in association to the trait was found. These regions identified novel genomic regions related to mastitis (1-Mb SNP windows) and confirmed those already mapped. The largest number of significant SNPs exceeding the threshold for genome-wide significant signal was found on BTA 15, located at 50.43-51.63 Mb.

CONCLUSIONS

The genomic regions identified in this study contribute to a better understanding of the genetic control of the mastitis immune response in cattle and may allow the inclusion of more detailed QTL information in selection programs.

Genomic prediction based on runs of homozygosity

Background

Genomic prediction is based on the accurate estimation of the genomic relationships among and between training animals and selection candidates in order to obtain accurate estimates of the genomic estimated breeding values (GEBV). Various methods have been used to predict GEBV based on population-wide linkage disequilibrium relationships (G_IBS) or sometimes on linkage analysis relationships (G_LA). Here, we propose a novel method to predict GEBV based on a genomic relationship matrix using runs of homozygosity (G_ROH). Runs of homozygosity were used to derive probabilities of multi-locus identity by descent chromosome segments. The accuracy and bias of the prediction of GEBV using G_ROH were compared to those using G_IBS and G_LA. Comparisons were performed using simulated datasets derived from a random pedigree and a real pedigree of Italian Brown Swiss bulls. The comparison of accuracies of GEBV was also performed on data from 1086 Italian Brown Swiss dairy cattle.

Results

Simulations with various thresholds of minor allele frequency for markers and quantitative trait loci showed that G_ROH achieved consistently more accurate GEBV (0 to 4% points higher) than G_IBS and G_LA. The bias of GEBV prediction for simulated data was higher based on the real pedigree than based on a random pedigree. In the analyses with real data, G_ROH and G_LA had similar accuracies. However, G_LA achieved a higher accuracy when the prediction was done on the youngest animals. The G_IBS matrices calculated with and without standardized marker genotypes resulted in similar accuracies.

Conclusions

The present study proposes G_ROH as a novel method to estimate genomic relationship matrices and predict GEBV based on runs of homozygosity and shows that it can result in higher or similar accuracies of GEBV prediction than G_LA, except for the real data analysis with validation of young animals. Compared to G_IBS, G_ROH resulted in more accurate GEBV predictions.