Molecular characterization of the Montecristo feral goats in the Mediterranean context

The Montecristo wild goat is an endangered feral population occurring on the homonymous island in the Tuscan Archipelago since a long time. The origins of Montecristo goats are still debated, with authors dating their introduction either back to Neolithic times or between the 6th and 13th century of the Common Era. To investigate the evolutionary history and relationships of this population we assembled a 50K SNP dataset including 55 Mediterranean breeds and two nuclei of Montecristo goats sampled on the island and from an ex situ conservation project on the Italian mainland, respectively. Diversity levels, gene flow, population structure and relationships were assessed through multiple approaches. The insular population scored the lowest values of both observed and expected heterozygosity, highlighting reduced genetic variation, while the ex situ nucleus highlighted a less severe reduction. Multivariate statistics, Neighbour-network and population structure analyses clearly separated the insular nucleus from all other breeds, but also the two Montecristo populations from each other. Treemix software analysis pinpointed possible genetic inputs received by the two Montecristo goat nuclei from different sources, while Runs Of Homozygosity (ROHs) indicated an ancient bottleneck/founder effect in the insular population and recent inbreeding in the ex situ one. Overall, our results suggest that Montecristo goats experienced several demographic fluctuations combined with admixture events over time, and highlighted a noticeable differentiation between the two nuclei. This evidence can serve as a starting point to implement marker-assisted conservation plans for the endangered Montecristo feral goat

The SNP-Based Profiling of Montecristo Feral Goat Populations Reveals a History of Isolation, Bottlenecks, and the Effects of Management

The Montecristo wild goat is an endangered feral population that has been on the homony-mous island in the Tuscan Archipelago since ancient times. The origins of Montecristo goats are still debated, with authors dating their introduction either back to Neolithic times or between the 6th and 13th century of the Common Era. To investigate the evolutionary history and relationships of this population we assembled a 50K SNP dataset including 55 Mediterranean breeds and two nuclei of Montecristo goats sampled on the island and from an ex situ conservation project. Diversity levels, gene flow, population structure, and genetic relationships were assessed through multiple approaches. The insular population scored the lowest values of both observed and expected heterozygosity, high-lighting reduced genetic variation, while the ex situ nucleus highlighted a less severe reduction. Multivariate statistics, network, and population structure analyses clearly separated the insular nucleus from all other breeds, including the population of Montecristo goats from the mainland. Moreover, admixture and gene flow analyses pinpointed possible genetic inputs received by the two Montecristo goat nuclei from different sources, while Runs of Homozygosity (ROHs) indicated an ancient bottleneck/founder effect in the insular population and recent extensive inbreeding in the ex situ one. Overall, our results suggest that Montecristo goats experienced several demographic fluctuations combined with admixture events over time and highlighted a noticeable differentiation between the two nuclei

Local adaptations of Mediterranean sheep and goats through an integrative approach

Small ruminants are suited to a wide variety of habitats and thus represent promising study models for identifying genes underlying adaptations. Here, we considered local Mediterranean breeds of goats (n = 17) and sheep (n = 25) from Italy, France and Spain. Based on historical archives, we selected the breeds potentially most linked to a territory and defined their original cradle (i.e., the geographical area in which the breed has emerged), including transhumant pastoral areas. We then used the programs PCAdapt and LFMM to identify signatures of artificial and environmental selection. Considering cradles instead of current GPS coordinates resulted in a greater number of signatures identified by the LFMM analysis. The results, combined with a systematic literature review, revealed a set of genes with potentially key adaptive roles in relation to the gradient of aridity and altitude. Some of these genes have been previously implicated in lipid metabolism (SUCLG2, BMP2), hypoxia stress/lung function (BMPR2), seasonal patterns (SOX2, DPH6) or neuronal function (TRPC4, TRPC6). Selection signatures involving the PCDH9 and KLH1 genes, as well as NBEA/NBEAL1, were identified in both species and thus could play an important adaptive role

Identification of differential gene expression in in vitro FSH treated pig granulosa cells using suppression subtractive hybridization

FSH, which binds to specific receptors on granulosa cells in mammals, plays a key role in folliculogenesis. Its biological activity involves stimulation of intercellular communication and upregulation of steroidogenesis, but the entire spectrum of the genes regulated by FSH has yet to be fully characterized. In order to find new regulated transcripts, however rare, we have used a Suppression Subtractive Hybridization approach (SSH) on pig granulosa cells in primary culture treated or not with FSH. Two SSH libraries were generated and 76 clones were sequenced after selection by differential screening. Sixty four different sequences were identified, including 3 novel sequences. Experiments demonstrated the presence of 25 regulated transcripts. A gene ontology analysis of these 25 genes revealed (1) catalytic; (2) transport; (3) signal transducer; (4) binding; (5) anti-oxidant and (6) structural activities. These findings may deepen our understanding of FSH's effects. Particularly, they suggest that FSH is involved in the modulation of peroxidase activity and remodelling of chromatin

Methods for interpreting lists of affected genes obstained in a DNA microarray experiment

Background - The aim of this paper was to describe and compare the methods used and the results obtained by the participants in a joint EADGENE (European Animal Disease Genomic Network of Excellence) and SABRE (Cutting Edge Genomics for Sustainable Animal Breeding) workshop focusing on post analysis of microarray data. The participating groups were provided with identical lists of microarray probes, including test statistics for three different contrasts, and the normalised log-ratios for each array, to be used as the starting point for interpreting the affected probes. The data originated from a microarray experiment conducted to study the host reactions in broilers occurring shortly after a secondary challenge with either a homologous or heterologous species of Eimeria. Results - Several conceptually different analytical approaches, using both commercial and public available software, were applied by the participating groups. The following tools were used: Ingenuity Pathway Analysis, MAPPFinder, LIMMA, GOstats, GOEAST, GOTM, Globaltest, TopGO, ArrayUnlock, Pathway Studio, GIST and AnnotationDbi. The main focus of the approaches was to utilise the relation between probes/genes and their gene ontology and pathways to interpret the affected probes/genes. The lack of a well-annotated chicken genome did though limit the possibilities to fully explore the tools. The main results from these analyses showed that the biological interpretation is highly dependent on the statistical method used but that some common biological conclusions could be reached. Conclusion - It is highly recommended to test different analytical methods on the same data set and compare the results to obtain a reliable biological interpretation of the affected genes in a DNA microarray experimen

Signatures of selection and environmental adaptation across the goat genome post-domestication

Background: Since goat was domesticated 10,000 years ago, many factors have contributed to the differentiation of goat breeds and these are classified mainly into two types: (i) adaptation to different breeding systems and/or purposes and (ii) adaptation to different environments. As a result, approximately 600 goat breeds have developed worldwide; they differ considerably from one another in terms of phenotypic characteristics and are adapted to a wide range of climatic conditions. In this work, we analyzed the AdaptMap goat dataset, which is composed of data from more than 3000 animals collected worldwide and genotyped with the CaprineSNP50 BeadChip. These animals were partitioned into groups based on geographical area, production uses, available records on solid coat color and environmental variables including the sampling geographical coordinates, to investigate the role of natural and/or artificial selection in shaping the genome of goat breeds. Results: Several signatures of selection on different chromosomal regions were detected across the different breeds, sub-geographical clusters, phenotypic and climatic groups. These regions contain genes that are involved in important biological processes, such as milk-, meat- or fiber-related production, coat color, glucose pathway, oxidative stress response, size, and circadian clock differences. Our results confirm previous findings in other species on adaptation to extreme environments and human purposes and provide new genes that could explain some of the differences between goat breeds according to their geographical distribution and adaptation to different environments. Conclusions: These analyses of signatures of selection provide a comprehensive first picture of the global domestication process and adaptation of goat breeds and highlight possible genes that may have contributed to the differentiation of this species worldwide

VarGoats project: a dataset of 1159 whole-genome sequences to dissect Capra hircus global diversity

Background: Since their domestication 10,500 years ago, goat populations with distinctive genetic backgrounds have adapted to a broad variety of environments and breeding conditions. The VarGoats project is an international 1000-genome resequencing program designed to understand the consequences of domestication and breeding on the genetic diversity of domestic goats and to elucidate how speciation and hybridization have modeled the genomes of a set of species representative of the genus Capra. Findings: A dataset comprising 652 sequenced goats and 507 public goat sequences, including 35 animals representing eight wild species, has been collected worldwide. We identified 74,274,427 single nucleotide polymorphisms (SNPs) and 13,607,850 insertion-deletions (InDels) by aligning these sequences to the latest version of the goat reference genome (ARS1). A Neighbor-joining tree based on Reynolds genetic distances showed that goats from Africa, Asia and Europe tend to group into independent clusters. Because goat breeds from Oceania and Caribbean (Creole) all derive from imported animals, they are distributed along the tree according to their ancestral geographic origin. Conclusions: We report on an unprecedented international effort to characterize the genome-wide diversity of domestic goats. This large range of sequenced individuals represents a unique opportunity to ascertain how the demographic and selection processes associated with post-domestication history have shaped the diversity of this species. Data generated for the project will also be extremely useful to identify deleterious mutations and polymorphisms with causal effects on complex traits, and thus will contribute to new knowledge that could be used in genomic prediction and genome-wide association studies

SNPchiMp v.3: integrating and standardizing single nucleotide polymorphism data for livestock species

Published: 10 April 2015In recent years, the use of genomic information in livestock species for genetic improvement, association studies and many other fields has become routine. In order to accommodate different market requirements in terms of genotyping cost, manufacturers of single nucleotide polymorphism (SNP) arrays, private companies and international consortia have developed a large number of arrays with different content and different SNP density. The number of currently available SNP arrays differs among species: ranging from one for goats to more than ten for cattle, and the number of arrays available is increasing rapidly. However, there is limited or no effort to standardize and integrate array- specific (e.g. SNP IDs, allele coding) and species-specific (i.e. past and current assemblies) SNP information.Here we present SNPchiMp v.3, a solution to these issues for the six major livestock species (cow, pig, horse, sheep, goat and chicken). Original data was collected directly from SNP array producers and specific international genome consortia, and stored in a MySQL database. The database was then linked to an open-access web tool and to public databases. SNPchiMp v.3 ensures fast access to the database (retrieving within/across SNP array data) and the possibility of annotating SNP array data in a user-friendly fashion.This platform allows easy integration and standardization, and it is aimed at both industry and research. It also enables users to easily link the information available from the array producer with data in public databases, without the need of additional bioinformatics tools or pipelines. In recognition of the open-access use of Ensembl resources, SNPchiMp v.3 was officially credited as an Ensembl E!mpowered tool. Availability at http://bioinformatics.tecnoparco.org/SNPchimp.Ezequiel L Nicolazzi, Andrea Caprera, Nelson Nazzicari, Paolo Cozzi, Francesco Strozzi, Cindy Lawley, Ali Pirani, Chandrasen Soans, Fiona Brew, Hossein Jorjani, Gary Evans, Barry Simpson, Gwenola Tosser-Klopp, Rudiger Brauning, John L Williams and Alessandra Stell

On the Effect of Thermodynamic Equilibrium on the Assembly Efficiency of Complex Multi-Layered Virus-Like Particles (VLP): the Case of Rotavirus VLP

Previous studies have reported the production of malformed virus-like-particles (VLP) in recombinant host systems. Here we computationally investigate the case of a large triple-layered rotavirus VLP (RLP). In vitro assembly, disassembly and reassembly data provides strong evidence of microscopic reversibility of RLP assembly. Light scattering experimental data also evidences a slow and reversible assembly untypical of kinetic traps, thus further strengthening the fidelity of a thermodynamically controlled assembly. In silico analysis further reveals that under favourable conditions particles distribution is dominated by structural subunits and completely built icosahedra, while other intermediates are present only at residual concentrations. Except for harshly unfavourable conditions, assembly yield is maximised when proteins are provided in the same VLP protein mass composition. The assembly yield decreases abruptly due to thermodynamic equilibrium when the VLP protein mass composition is not obeyed. The latter effect is more pronounced the higher the Gibbs free energy of subunit association is and the more complex the particle is. Overall this study shows that the correct formation of complex multi-layered VLPs is restricted to a narrow range of association energies and protein concentrations, thus the choice of the host system is critical for successful assembly. Likewise, the dynamic control of intracellular protein expression rates becomes very important to minimize wasted proteins