Stefano Frattini

A method for single nucleotide polymorphism selection for parentage assessment in goats

Accurate pedigrees are essential to optimize genetic improvement and conservation of animal genetic resources. In goats, the use of mating groups and kidding management procedures hamper the identification of parentage. Small panels of single nucleotide polymorphisms (SNP) have been proposed in other species to substitute microsatellites for parentage assessment. Using data from the current GoatSNP50 chip, we developed a new 3-step procedure to identify a low-density SNP panel for highly accurate parentage assessment. Methodologies for SNP selection used in other species are less suitable in the goat because of uncertainties in the genome assembly. The procedure developed in this study is based on parent-offspring identification and on estimation of Mendelian errors, followed by canonical discriminant analysis identification and stepwise regression reduction. Starting from a reference sample of 109 Alpine goats with known pedigree relationships, we first identified a panel of 200 SNP that was further reduced to 2 final panels of 130 and 114 SNP with random coincidental match inclusion of 1.51×10(-57) and 2.94×10(-34), respectively. In our reference data set, all panels correctly identified all parent-offspring combinations, revealing a 40% pedigree error rate in the information provided by breeders. All reference trios were confirmed by official tests based on microsatellites. Panels were also tested on Saanen and Teramana breeds. Although the testing on a larger set of breeds in the reference population is still needed to validate these results, our findings suggest that our procedure could identify SNP panels for accurate parentage assessment in goats or in other species with unreliable marker positioning.

Early career researchers want Open Science

Open Science is encouraged by the European Union and many other political and scientific institutions. However, scientific practice is proving slow to change. We propose, as early career researchers, that it is our task to change scientific research into open scientific research and commit to Open Science principles.

Genomic analysis suggests KITLG is responsible for a roan pattern in two Pakistani goat breeds

The roan coat color pattern is described as the presence of white hairs intermixed with pigmented hairs. This kind of pigmentation pattern has been observed in many domestic species, including the goat. The molecular mechanisms and inheritance that underlie this pattern are known for some species and the KITLG gene has been shown associated with this phenotype. To date, no research effort has been carried out to find the gene(s) that control(s) roan coat color pattern in goats. In the present study, after genotyping with the GoatSNP50 BeadChip, 35 goats that showed a roan pattern and that belonged to two Pakistan breeds (Group A) were analyzed and then compared to 740 goats of 39 Italian and Pakistan goat breeds that did not have the same coat color pattern (Group B). Runs of homozygosity-based and XP-EHH analyses were used to identify unique genomic regions potentially associated with the roan pattern. A total of 3 regions on chromosomes 5, 6, and 12 were considered unique among the group A versus group B comparisons. The A region > 1.7 Mb on chromosome 5 was the most divergent between the two groups. This region contains six genes, including the KITLG gene. Our findings support the hypothesis that the KITLG gene may be associated with the roan phenotype in goats.

Polymorphism of the STAT5A, MTNR1A and TNFα genes and their effect on dairy production in Bubalus bubalis

Abstract The water buffalo is a fundamental resource, especially in developing countries, however, differently from other species, its genetic potential is still poorly investigated. In this work, we performed a candidate gene association study for milk composition in 491 female buffaloes. Animals were from four farms located in Southern Italy, where the Out-of-Breeding-Season-Mating technique is usually performed. We analysed three genes: (1) the signal transducer and activator of transcription 5A (STAT5A), (2) the tumour necrosis factor alpha (TNFα) and (3) the melatonin receptor 1A (MTNR1A). We confirmed the mutation at the MTNR1A gene and we found five novel single nucleotide polymorphisms (SNPs): one in the TNFα and four in the STAT5A. No associations were found for the SNPs in the MTNR1A and TNFα genes, while we identified a marked association with milk protein % for a C > T substitution at the STAT5A gene. At this locus, the TT buffaloes showed significantly higher protein percentage in milk. Conversely, this genotype class was the less frequent in the population. Moreover, an A > G substitution at the STAT5A showed an influence on reproductive seasonality, with the advantageous allele most frequent in the population, suggesting a possible effect of selection for this trait. The C > T substitution on STAT5A detected in present study could be used in marker assisted selection of Mediterranean Italian buffalo, and should be monitored to understand the reasons behind the low frequency of the favourable genotype at this locus and to stop this unfavourable trend in the population.

Genome-wide analysis of DNA methylation in hypothalamus and ovary of Capra hircus

BackgroundDNA methylation is a frequently studied epigenetic modification due to its role in regulating gene expression and hence in biological processes and in determining phenotypic plasticity in organisms. Rudimentary DNA methylation patterns for some livestock species are publically available: among these, goat methylome deserves to be further explored.ResultsGenome-wide DNA methylation maps of the hypothalamus and ovary from Saanen goats were generated using Methyl-CpG binding domain protein sequencing (MBD-seq). Analysis of DNA methylation patterns indicate that the majority of methylation peaks found within genes are located gene body regions, for both organs. Analysis of the distribution of methylated sites per chromosome showed that chromosome X had the lowest number of methylation peaks. The X chromosome has one of the highest percentages of methylated CpG islands in both organs, and approximately 50% of the CpG islands in the goat epigenome are methylated in hypothalamus and ovary. Organ-specific Differentially Methylated Genes (DMGs) were correlated with the expression levels.ConclusionsThe comparison between transcriptome and methylome in hypothalamus and ovary showed that a higher level of methylation is not accompanied by a higher gene suppression. The genome-wide DNA methylation map for two goat organs produced here is a valuable starting point for studying the involvement of epigenetic modifications in regulating goat reproduction performance.

Initial genomic characterization of Italian, Egyptian and Pakistani goat breeds.

Selection and breeding practices in goats have differed greatly among countries and populations. These processes, together with natural selection and regional drift, have shaped the phenotypic variability of goat breeds (Kim et al., 2015). The availability of improved genomic analysis tools for this species may provide useful information on the history of selection, adaptation and differentiation of goats from different areas of the world, that can be evaluated by the study of gene frequencies and length of the Runs of Homozigosity (contiguous length of homozygous genotypes, ROH; Purfield et al., 2012). In current study, we examined using a goat medium density SNP chip animals from three different countries: Egypt (with lack of selection scheme), Italy (with several standardized breeds; Nicoloso et al., 2015) and Pakistan (with several breeds showing peculiar phenotypes) to produce a genomic landscape of goats breeds in these countries. A total of 1,123 animals of 39 different populations, and 48,895 SNP markers were analyzed. Genotypes were imputed on a country-based approach, and markers without known position in the genome were excluded from the analysis. MDS and ADMIXTURE plots confirmed the good differentiation among populations from the three countries. Runs of Homozygosity (ROH) were performed for each country and population allowed the detection of genomic regions with high homozygosity levels, common in at least two out of three sampling areas. These results provide new insights into goat genome structure within and among breeds and countries. The detection of conserved regions with different lengths may explain recent selection strategies or adaptation to different, extreme environmental conditions. The research was funded by INNOVAGEN project. Support by Iowa State University and the Ensminger funds for AE and AT as well as support by the Fulbright Foundation for AE are gratefully acknowledged. Sampling from Pakistan was funded by PAK-USAID project.

Studies of modern Italian dog populations reveal multiple patterns for domestic breed evolution

Abstract Through thousands of years of breeding and strong human selection, the dog (Canis lupus familiaris) exists today within hundreds of closed populations throughout the world, each with defined phenotypes. A singular geographic region with broad diversity in dog breeds presents an interesting opportunity to observe potential mechanisms of breed formation. Italy claims 14 internationally recognized dog breeds, with numerous additional local varieties. To determine the relationship among Italian dog populations, we integrated genetic data from 263 dogs representing 23 closed dog populations from Italy, seven Apennine gray wolves, and an established dataset of 161 globally recognized dog breeds, applying multiple genetic methods to characterize the modes by which breeds are formed within a single geographic region. Our consideration of each of five genetic analyses reveals a series of development events that mirror historical modes of breed formation, but with variations unique to the codevelopment of early dog and human populations. Using 142,840 genome‐wide SNPs and a dataset of 1,609 canines, representing 182 breeds and 16 wild canids, we identified breed development routes for the Italian breeds that included divergence from common populations for a specific purpose, admixture of regional stock with that from other regions, and isolated selection of local stock with specific attributes.

Distribution of ncRNAs expression across hypothalamic-pituitary-gonadal axis in Capra hircus

BackgroundMolecular regulation of the hypothalamic-pituitary-gonadal (HPG) axis plays an essential role in the fine tuning of seasonal estrus in Capra hircus. Noncoding RNAs (ncRNAs) are emerging as key regulators in sexual development and mammalian reproduction. In order to identify ncRNAs and to assess their expression patterns, along the HPG axis, we sequenced ncRNA libraries from hypothalamus, pituitary and ovary of three goats.ResultsAmong the medium length noncoding RNAs (mncRNAs) identified, small nucleolar RNAs (snoRNAs) and transfer RNAs (tRNAs) were found to be more abundant in ovary and hypothalamus, respectively. The observed GC content was representative for different classes of ncRNAs, allowing the identification of a tRNA-derived RNA fragments (tRFs) subclass, which had a peak distribution around 32–38% GC content in the hypothalamus. Differences observed among organs confirmed the specificity of microRNA (miRNA) profiles for each organ system.ConclusionsData on ncRNAs in organs constituting the HPG axis will contribute to understanding their role in the physiological regulation of reproduction in goats.

Investigating the population structure and genetic differentiation of livestock guard dog breeds

Livestock guarding dogs are a valuable adjunct to the pastoral community. Having been traditionally selected for their working ability, they fulfil their function with minimal interaction or command from their human owners. In this study, the population structure and the genetic differentiation of three Italian livestock guardian breeds (Silas Dog, Maremma and Abruzzese Sheepdog and Mannaras Dog) and three functionally and physically similar breeds (Cane Corso, Central Asian Shepherd Dog and Caucasian Shepherd Dog), totalling 179 dogs unrelated at the second generation, were investigated with 18 autosomal microsatellite markers. Values for the number of alleles per locus, observed and expected heterozygosity, Hardy-Weinberg Equilibrium, F stats, Neis and Reynolds genetic distances, clustering and sub-population formation abilities and individual genetic structures were calculated. Our results show clear breed differentiation, whereby all the considered breeds show reasonable genetic variability despite small population sizes and variable selection schemes. These results provide meaningful data to stakeholders in specific breed and environmental conservation programmes.

The methylome of the hypothalamus of prepubertal and pubertal goats

Puberty is the fulfillment of fertility, a process involving physiological and morphological development. It is well known that the increased hypothalamic secretion of the gonadotropin-releasing hormone (GnRH) is essential for the activation of this process, even if the elements coordinating the timing of puberty have not been fully identified1,2. Recent studies provide proof that there is an epigenetic regulation of female puberty, and DNA methylation, the most studied epigenetic modification, plays a major role in it3. We analyzed DNA methylation patterns of 5 Alpine goats at their prepubertal stage and 5 that reached puberty in order to highlight differences in their methylome. Detection of methylated regions across the goat genome involved a Methyl Binding Domain (MBD) enrichment followed by deep sequencing (Hiseq2000 Illumina). The software ChIPseeqer4 permitted the identification of peaks corresponding to hyper-methylated regions. We have observed a higher methylation level in prepubertal goats. The distribution of the methylation peaks across the genome and within CpG islands per chromosome per group of animals has been analyzed. Furthermore, we have investigated differential methylation in genes associated with puberty. Specifically, Cbx7, coding for a core component of the Polycomb group silencing complex, and GnRHR, the gene coding for GnRH receptor, showed a higher number of peaks into two intragenic fragments within prepubertal goats. These results, accompanied by transcriptome analysis, provide a foundation for elucidating the role of DNA methylation in the complex mechanisms that drive puberty in goat species.