Prediction of detailed blood metabolic profile using milk infrared spectra and machine learning methods in dairy cattle

The adoption of preventive management decisions is crucial to dealing with metabolic impairments in dairy cattle. Various serum metabolites are known to be useful indicators of the health status of cows. In this study, we used milk Fourier-transform mid-infrared (FTIR) spectra and various machine learning (ML) algorithms to develop prediction equations for a panel of 29 blood metabolites, including those related to energy metabolism, liver function/hepatic damage, oxidative stress, inflammation/innate immunity, and minerals. For most traits, the data set comprised observations from 1,204 Holstein-Friesian dairy cows belonging to 5 herds. An exception was represented by β-hydroxybutyrate prediction, which contained observations from 2,701 multibreed cows pertaining to 33 herds. The best predictive model was developed using an automatic ML algorithm that tested various methods, including elastic net, distributed random forest, gradient boosting machine, artificial neural network, and stacking ensemble. These ML predictions were compared with partial least squares regression, the most commonly used method for FTIR prediction of blood traits. Performance of each model was evaluated using 2 cross-validation (CV) scenarios: 5-fold random (CV_r) and herd-out (CV_h). We also tested the best model's ability to classify values precisely in the 2 extreme tails, namely, the 25th (Q25) and 75th (Q75) percentiles (true-positive prediction scenario). Compared with partial least squares regression, ML algorithms achieved more accurate performance. Specifically, elastic net increased the R² value from 5% to 75% for CV_r and 2% to 139% for CV_h, whereas the stacking ensemble increased the R² value from 4% to 70% for CV_r and 4% to 150% for CV_h. Considering the best model, with the CV_r scenario, good prediction accuracies were obtained for glucose (R² = 0.81), urea (R² = 0.73), albumin (R² = 0.75), total reactive oxygen metabolites (R² = 0.79), total thiol groups (R² = 0.76), ceruloplasmin (R² = 0.74), total proteins (R² = 0.81), globulins (R² = 0.87), and Na (R² = 0.72). Good prediction accuracy in classifying extreme values was achieved for glucose (Q25 = 70.8%, Q75 = 69.9%), albumin (Q25 = 72.3%), total reactive oxygen metabolites (Q25 = 75.1%, Q75 = 74%), thiol groups (Q75 = 70.4%), total proteins (Q25 = 72.4%, Q75 = 77.2.%), globulins (Q25 = 74.8%, Q75 = 81.5%), and haptoglobin (Q75 = 74.4%). In conclusion, our study shows that FTIR spectra can be used to predict blood metabolites with relatively good accuracy, depending on trait, and are a promising tool for large-scale monitoring.

Short Communication: Maternal undernutrition during peri-conceptional period affects whole genome ovine muscle methylation in adult offspring

Experimental and epidemiological studies suggest that maternal nutritional status during early pregnancy, including the period around the time of conception, may induce long-lasting epigenetic changes in the offspring. However, this remains largely unexplored in livestock. Therefore, the objective of this study was to evaluate if modification of the maternal diet of sheep (CTR: control; UND: 50% undernutrition) during the peri-conceptional period (42 days in total: -14/+28 from mating), would impact CpG methylation in muscle tissue (Longissimus dorsi) of adult offspring (11.5 months old). Reduced Representation Bisulfite Sequencing (RRBS), identified 262 (Edge-R, FDR<0.05) and 686 (Logistic Regression, FDR <0.001) differentially methylated regions (DMRs) between the UND and CTR groups. Gene ontology (GO) analysis identified genes related to development, functions of the muscular system and steroid hormone receptor activity within the DMRs. The data reported here show that nutritional stress during early pregnancy leads to epigenetic modifications in the muscle of the resulting offspring, with possible implications for cardiac dysfunction, muscle physiology and meat production.

Associations between Milk Fatty Acid Profile and Body Condition Score, Ultrasound Hepatic Measurements and Blood Metabolites in Holstein Cows

Dairy cows have high incidences of metabolic disturbances, which often lead to disease, having a subsequent significant impact on productivity and reproductive performance. As the milk fatty acid (FA) profile represents a fingerprint of the cow’s nutritional and metabolic status, it could be a suitable indicator of metabolic status at the cow level. In this study, we obtained milk FA profile and a set of metabolic indicators (body condition score, ultrasound liver measurements, and 29 hematochemical parameters) from 297 Holstein–Friesian cows. First, we applied a multivariate factor analysis to detect latent structure among the milk FAs. We then explored the associations between these new synthetic variables and the morphometric, ultrasonographic and hematic indicators of immune and metabolic status. Significant associations were exhibited by the odd-chain FAs, which were inversely associated with β-hydroxybutyrate and ceruloplasmin, and positively associated with glucose, albumin, and γ-glutamyl transferase. Short-chain FAs were inversely related to predicted triacylglycerol liver content. Rumen biohydrogenation intermediates were associated with glucose, cholesterol, and albumin. These results offer new insights into the potential use of milk FAs as indicators of variations in energy and nutritional metabolism in early lactating dairy cows.

Runs of homozygosity in the Italian goat breeds: impact of management practices in low-input systems

Background

Climate and farming systems, several of which are considered as low-input agricultural systems, vary between goat populations from Northern and Southern Italy and have led to different management practices. These processes have impacted genome shaping in terms of inbreeding and regions under selection and resulted in differences between the northern and southern populations. Both inbreeding and signatures of selection can be pinpointed by the analysis of runs of homozygosity (ROH), which provides useful information to assist the management of this species in different rural areas.

Results

We analyzed the ROH distribution and inbreeding (F_ROH) in 902 goats from the Italian Goat Consortium2 dataset. We evaluated the differences in individual ROH number and length between goat breeds from Northern (NRD) and Central-southern (CSD) Italy. Then, we identified the signatures of selection that differentiate these two groups using three methods: ROH, ΔROH, and averaged F_ST. ROH analyses showed that some Italian goat breeds have a lower inbreeding coefficient, which is attributable to their management and history. ROH are longer in breeds that are undergoing non-optimal management or with small population size. In several small breeds, the ROH length classes are balanced, reflecting more accurate mating planning. The differences in climate and management between the NRD and CSD groups have resulted in different ROH lengths and numbers: the NRD populations bred in isolated valleys present more and shorter ROH segments, while the CSD populations have fewer and longer ROH, which is likely due to the fact that they have undergone more admixture events during the horizontal transhumance practice followed by a more recent standardization. We identified four genes within signatures of selection on chromosome 11 related to fertility in the NRD group, and 23 genes on chromosomes 5 and 6 related to growth in the CSD group. Finally, we identified 17 genes on chromosome 12 related to environmental adaptation and body size with high homozygosity in both groups.

Conclusions

These results show how different management practices have impacted the level of genomic inbreeding in two Italian goat groups and could be useful to assist management in a low-input system while safeguarding the diversity of small populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00685-4.

Associations between ultrasound measurements and hematochemical parameters for the assessment of liver metabolic status in Holstein-Friesian cows

Metabolic disorders, including hepatic lipidosis and ketosis, severely affect animal health status and welfare with a large economic burden in dairy herds. The gold standard for diagnosing hepatic lipidosis is the liver biopsy, which is impractical and invasive for the screening at farm level. Ultrasound (US) imaging is a promising technique for identifying liver dysfunction, but standardized specifications in physiological conditions are needed. Herein, we described the features of four US measurements, namely the liver predicted triacylglycerol (pTAG) content, liver depth (LD), and portal vein area (PVA) and depth (PVD) and we investigated their associations with a set of hematochemical (HC) indicators in 342 clinically healthy Holstein Friesian dairy cows. Liver pTAG content was negatively associated with hematocrit and positively with globulin, whereas PVA was negatively associated with thiol group levels, and LD positively with ceruloplasmin. We found significant interactions between some HC parameters and parity: in particular, creatinine, thiol groups and globulin for PVA, and aspartate aminotransferase, paraoxonase and ceruloplasmin for PVD. This study offers new insights on variations in liver function occurring after calving and pave the way for the potential use of minimally invasive techniques for prompt detection of metabolic disorders in dairy herds.

Maternal peri-conceptional undernourishment perturbs offspring sperm methylome

The genotype of an organism is stable throughout its life; however, its epigenome is dynamic and can be altered in response to environmental factors, such as diet. Inheritance of acquired epigenetic modifications by the next generation occurs through the germline, although the precise mechanisms remain to be elucidated. Here, we used a sheep model to evaluate if modification of the maternal diet (CTR; control, UND: undernutrition; FA: undernutrition and folic acid supplementation) during the peri-conceptional period affects the genome-wide methylation status of the gametes of male offspring. Sperm DNA methylation, measured by Reduced Representation Bisulfite Sequencing (RRBS), identified Differentially Methylated Regions (DMR) in offspring that experienced in utero undernutrition, both in UND (244) and FA (240), compared with CTR. Gene ontology (GO) analysis identified DMRs in categories related to sperm function, therefore we investigated whether the fertilizing capacity of the semen from the three groups differed in an in vitro fertilization assay. Spermatozoa from the undernourished groups showed lower motility and sperm chromatin structure abnormalities, represented by a higher percentage of DNA fragmentation and an increased number of immature cells, compared with CTR. While good quality blastocysts were obtained from all three groups, the proportion of embryos reaching the blastocyst stage was reduced in the UND vs CTR, an effect partially rescued by the FA treatment. The data reported here show that nutritional stress during early pregnancy leads to epigenetic modifications in the semen of the resulting offspring, the effects of which in next generation remain to be elucidated.

Genomic prediction for latent variables related to milk fatty acid composition in Holstein, Simmental and Brown Swiss dairy cattle breeds

Genomic selection (GS) reports on milk fatty acid (FA) profiles have been published quite recently and are still few despite this trait represents the most important aspect of milk nutritional and sensory quality. Reasons for this can be found in the high costs of phenotype recording but also in issues related to its nature of complex trait constituted by multiple genetically correlated variables with low heritabilities. One possible strategy to deal with such constraint is represented by the use of dimension reduction methods. We analysed 40 individual FAs from Italian Brown Swiss, Holstein and Simmental milk through multivariate factor analysis (MFA) to study the genetics of milk FA-related latent variables (factors) and assess their potential use in breeding. A total of nine factors were obtained, and their genetic parameters were inferred under a Bayesian framework using two statistical approaches: the classical pedigree best linear unbiased prediction (ABLUP) and the single-step genomic BLUP (ssGBLUP). The resulting factorial solutions were able to represent groups of FAs with common origin and function and can be considered concise pathway-level phenotypes. The heritability (h² ) values showed relevant variations across different factors in each breed (0.03 ≤ h²  ≤ 0.38). The accuracies of breeding values predicted were low to high, ranging from 0.13 to 0.72 and from 0.18 to 0.74 considering the pedigree and the genomic model, respectively. The gain in accuracy in genetic prediction due to the addition of genomic information was ~30% and ~5% in validation and training groups respectively, confirming the contribution of genomic information in yielding more accurate predictions compared to the traditional ABLUP methodology. Our results suggest that MFA in combination with GS can be a valuable tool in dairy cattle breeding and deserves to be further investigated for use in future breeding programs to improve cow's milk FA-related traits.

Genome-Wide DNA Methylation and Gene Expression Profiles in Cows Subjected to Different Stress Level as Assessed by Cortisol in Milk

Dairy cattle health, wellbeing and productivity are deeply affected by stress. Its influence on metabolism and immune response is well known, but the underlying epigenetic mechanisms require further investigation. In this study, we compared DNA methylation and gene expression signatures between two dairy cattle populations falling in the high- and low-variant tails of the distribution of milk cortisol concentration (MC), a neuroendocrine marker of stress in dairy cows. Reduced Representation Bisulfite Sequencing was used to obtain a methylation map from blood samples of these animals. The high and low groups exhibited similar amounts of methylated CpGs, while we found differences among non-CpG sites. Significant methylation changes were detected in 248 genes. We also identified significant fold differences in the expression of 324 genes. KEGG and Gene Ontology (GO) analysis showed that genes of both groups act together in several pathways, such as nervous system activity, immune regulatory functions and glucocorticoid metabolism. These preliminary results suggest that, in livestock, cortisol secretion could act as a trigger for epigenetic regulation and that peripheral changes in methylation can provide an insight into central nervous system functions.

Use of multivariate factor analysis of detailed milk fatty acid profile to perform a genome-wide association study in Italian Simmental and Italian Holstein

Milk fatty acid (FA) profile is a clear example of complex and multiple correlated traits whose genetic basis is difficult to assess. Although genome-wide association (GWA) studies have been successful in the identification of significant genetic variants for complex traits, when correlated phenotypes are analysed separately, the outcomes are difficult to compare and interpret in a metabolic context. Here, we performed a multivariate factor analysis (MFA) on Italian Simmental and Italian Holstein milk fat profiles to extract latent unobserved factors able to explain correlation structure and common metabolic function among different FAs. Individual factor scores obtained by MFA were used to perform a single-SNP based GWA. In both breeds, MFA was able to extract ten latent factors with specific biological meaning, notably: de novo synthesis, desaturation activity and biohydrogenation. The GWA result confirmed the increased power of joint association analysis on multiple correlated traits and allowed us to identify major candidate genes with well-documented function consistent with the metabolic classification of factors obtained, such as DGAT1, FASN and SCD.

A genetically unique Chinese cattle population shows evidence of common ancestry with wild species when analysed with a reduced ascertainment bias SNP panel

In Hong Kong, there is a cattle population of ~1,200 individuals of uncertain origin and genetic diversity. This population shows heterogeneous morphology, both in body type and pigmentation. Once used as draught animals by the local farmers, they were abandoned around the 1970s due to changes in the economy, and since then have lived as feral populations. To explore the origins of these cattle, we analysed ~50k genotype data of 21 Hong Kong feral cattle, along with data from 703 individuals of 36 cattle populations of European, African taurine, and Asian origin, the wild x domestic hybrid gayal, plus two wild bovine species, gaur and banteng. To reduce the effect of ascertainment bias ~4k loci that are polymorphic in the two wild species were selected for further analysis. The stringent SNP selection we applied resulted in increased heterozygosity across all populations studies, compared with the full panel of SNP, thus reducing the impact of ascertainment bias and facilitating the comparison of divergent breeds of cattle. Our results showed that Hong Kong feral cattle have relatively high levels of genetic distinctiveness, possibly due to the low level of artificial selection, and a likely common ancestry with wild species. We found signs of a putative taurine introgression, probably dating to the import of north European breeds during the British colonialism of Hong Kong. We showed that Hong Kong feral cattle, are distinct from Bos taurus and Bos indicus breeds. Our results highlight the distinctiveness of Hong Kong feral cattle and stress the conservation value of this indigenous breed that is likely to harbour adaptive genetic variation, which is a fundamental livestock resource in the face of climate change and diversifying market demands.

Occurrence, biological properties and potential effects on human health of β-casomorphin 7: Current knowledge and concerns

The genetic variant A1 of bovine β-casein (β-Cn) presents a His residue at a position 67 of the mature protein. This feature makes the Ile⁶⁶-His⁶⁷ bond more vulnerable to enzymatic cleavage, determining the release of the peptide β-Cn f(60-66), named β-casomorphin 7 (BCM7). BCM7 is an opioid-agonist for μ receptors, and it has been hypothesized to be involved in the development of different non-transmissible diseases in humans. In the last decade, studies have provided additional results on the potential health impact of β-Cn A1 and BCM7. These studies, here reviewed, highlighted a relation between the consumption of β-Cn A1 (and its derivative BCM7) and the increase of inflammatory response as well as discomfort at the gastrointestinal level. Conversely, the role of BCM7 and the effects of ingestion of β-Cn A1 on the onset or worsening of other non-transmissible diseases as caused or favored by still need proof of evidence. Overall, the reviewed literature demonstrates that the "β-Cn A1/BCM7 issue" remains an intriguing but not exhaustively explained topic in human nutrition. On this basis, policies in favor of breeding for β-Cn variants not releasing BCM7 and consumption of "A1-like" milk appear not yet sound for a healthier and safer nutrition.

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.

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

Background

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

Results

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

Conclusions

Data on ncRNAs in organs constituting the HPG axis will contribute to understanding their role in the physiological regulation of reproduction in goats.

Electronic supplementary material

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

Mitochondrial DNA variants of Podolian cattle breeds testify for a dual maternal origin

BACKGROUND

Over the past 15 years, 300 out of 6000 breeds of all farm animal species identified by the Food and Agriculture Organization of the United Nations (FAO) have gone extinct. Among cattle, many Podolian breeds are seriously endangered in various European areas. Podolian cattle include a group of very ancient European breeds, phenotypically close to the aurochs ancestors (Bos primigenius). The aim of the present study was to assess the genetic diversity of Podolian breeds and to reconstruct their origin.

METHODOLOGY

The mitochondrial DNA (mtDNA) control-regions of 18 Podolian breeds have been phylogenetically assessed. Nine non-Podolian breeds have been also included for comparison.

CONCLUSION

The overall analysis clearly highlights some peculiarities in the mtDNA gene pool of some Podolian breeds. In particular, a principal component analysis point to a genetic proximity between five breeds (Chianina, Marchigiana, Maremmana, Podolica Italiana and Romagnola) reared in Central Italy and the Turkish Grey. We here propose the suggestive hypothesis of a dual ancestral contribution to the present gene pool of Podolian breeds, one deriving from Eastern European cattle; the other arising from the arrival of Middle Eastern cattle into Central Italy through a different route, perhaps by sea, ferried by Etruscan boats. The historical migration of Podolian cattle from North Eastern Europe towards Italy has not cancelled the mtDNA footprints of this previous ancient migration.

Integrative analysis of methylomic and transcriptomic data in fetal sheep muscle tissues in response to maternal diet during pregnancy

BACKGROUND

Numerous studies have established a link between maternal diet and the physiological and metabolic phenotypes of their offspring. In previous studies in sheep, we demonstrated that different maternal diets altered the transcriptome of fetal tissues. However, the mechanisms underlying transcriptomic changes are poorly understood. DNA methylation is an epigenetic mark regulating transcription and is largely influenced by dietary components of the one-carbon cycle that generate the methyl group donor, SAM. Therefore, in the present study, we tested whether different maternal diets during pregnancy would alter the DNA methylation and gene expression patterns in fetal tissues.

RESULTS

Pregnant ewes were randomly divided into two groups which received either hay or corn diet from mid-gestation (day 67 ± 5) until day 131 ± 1 when fetuses were collected by necropsy. A total of 1516 fetal longissimus dorsi (LD) tissues were used for DNA methylation analysis and gene expression profiling. Whole genome DNA methylation using methyl-binding domain enrichment analysis revealed 60 differentially methylated regions (DMRs) between hay and corn fetuses with 39 DMRs more highly methylated in the hay fetuses vs. 21 DMRs more highly methylated in the corn fetuses. Three DMRs (LPAR3, PLIN5-PLIN4, and the differential methylation of a novel lincRNA) were validated using bisulfite sequencing. These DMRs were associated with differential gene expression. Additionally, significant DNA methylation differences were found at the single CpG level. Integrative methylome and transcriptome analysis revealed an association between gene expression and inter-/intragenic methylated regions. Furthermore, intragenic DMRs were found to be associated with expression of neighboring genes.

CONCLUSIONS

The findings of this study imply that maternal diet from mid- to late-gestation can shape the epigenome and the transcriptome of fetal tissues, and putatively affect phenotypes of the lambs.

Use of canonical discriminant analysis to study signatures of selection in cattle

BACKGROUND

Cattle include a large number of breeds that are characterized by marked phenotypic differences and thus constitute a valuable model to study genome evolution in response to processes such as selection and domestication. Detection of "signatures of selection" is a useful approach to study the evolutionary pressures experienced throughout history. In the present study, signatures of selection were investigated in five cattle breeds farmed in Italy using a multivariate approach.

METHODS

A total of 4094 bulls from five breeds with different production aptitudes (two dairy breeds: Italian Holstein and Italian Brown Swiss; two beef breeds: Piemontese and Marchigiana; and one dual purpose breed: Italian Simmental) were genotyped using the Illumina BovineSNP50 v.1 beadchip. Canonical discriminant analysis was carried out on the matrix of single nucleotide polymorphisms (SNP) genotyping data, separately for each chromosome. Scores for each canonical variable were calculated and then plotted in the canonical space to quantify the distance between breeds. SNPs for which the correlation with the canonical variable was in the 99th percentile for a specific chromosome were considered to be significantly associated with that variable. Results were compared with those obtained using an FST-based approach.

RESULTS

Based on the results of the canonical discriminant analysis, a large number of signatures of selection were detected, among which several had strong signals in genomic regions that harbour genes known to have an impact on production and morphological bovine traits, including MSTN, LCT, GHR, SCD, NCAPG, KIT, and ASIP. Moreover, new putative candidate genes were identified, such as GCK, B3GALNT1, MGAT1, GALNTL1, PRNP, and PRND. Similar results were obtained with the FST-based approach.

CONCLUSIONS

The use of canonical discriminant analysis on 50 K SNP genotypes allowed the extraction of new variables that maximize the separation between breeds. This approach is quite straightforward, it can compare more than two groups simultaneously, and relative distances between breeds can be visualized. The genes that were highlighted in the canonical discriminant analysis were in concordance with those obtained using the FST index.

Synergies between assisted reproduction technologies and functional genomics

This review, is a synopsis of advanced reproductive technologies in farm animals, including the discussion of their limiting factors as revealed by the study of offspring derived from embryos produced in vitro and through cloning. These studies show that the problems of epigenetic mis-programming, which were reported in the initial stages of assisted reproduction, still persist. The importance of whole-genome analyses, including the methylome and transcriptome, in improving embryo biotechnologies in farm animals, are discussed. Genome editing approaches for the improvement of economically-relevant traits in farm animals are also described. Efficient farm animal embryo biotechnologies, including cloning and the most recent technologies such as genome editing, will effectively complement the latest strategies to accelerate genetic improvement of farm animals.

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability

Background

The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA).

Results

We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12–10 ka ago.

Conclusions

Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds.

Electronic supplementary material

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

Adaptations to climate-mediated selective pressures in sheep

Following domestication, sheep (Ovis aries) have become essential farmed animals across the world through adaptation to a diverse range of environments and varied production systems. Climate-mediated selective pressure has shaped phenotypic variation and has left genetic "footprints" in the genome of breeds raised in different agroecological zones. Unlike numerous studies that have searched for evidence of selection using only population genetics data, here, we conducted an integrated coanalysis of environmental data with single nucleotide polymorphism (SNP) variation. By examining 49,034 SNPs from 32 old, autochthonous sheep breeds that are adapted to a spectrum of different regional climates, we identified 230 SNPs with evidence for selection that is likely due to climate-mediated pressure. Among them, 189 (82%) showed significant correlation (P ≤ 0.05) between allele frequency and climatic variables in a larger set of native populations from a worldwide range of geographic areas and climates. Gene ontology analysis of genes colocated with significant SNPs identified 17 candidates related to GTPase regulator and peptide receptor activities in the biological processes of energy metabolism and endocrine and autoimmune regulation. We also observed high linkage disequilibrium and significant extended haplotype homozygosity for the core haplotype TBC1D12-CH1 of TBC1D12. The global frequency distribution of the core haplotype and allele OAR22_18929579-A showed an apparent geographic pattern and significant (P ≤ 0.05) correlations with climatic variation. Our results imply that adaptations to local climates have shaped the spatial distribution of some variants that are candidates to underpin adaptive variation in sheep.

Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel

In this study we used a medium density panel of SNP markers to perform population genetic analysis in five Italian cattle breeds. The BovineSNP50 BeadChip was used to genotype a total of 2,935 bulls of Piedmontese, Marchigiana, Italian Holstein, Italian Brown and Italian Pezzata Rossa breeds. To determine a genome-wide pattern of positive selection we mapped the F st values against genome location. The highest F st peaks were obtained on BTA6 and BTA13 where some candidate genes are located. We identified selection signatures peculiar of each breed which suggest selection for genes involved in milk or meat traits. The genetic structure was investigated by using a multidimensional scaling of the genetic distance matrix and a Bayesian approach implemented in the STRUCTURE software. The genotyping data showed a clear partitioning of the cattle genetic diversity into distinct breeds if a number of clusters equal to the number of populations were given. Assuming a lower number of clusters beef breeds group together. Both methods showed all five breeds separated in well defined clusters and the Bayesian approach assigned individuals to the breed of origin. The work is of interest not only because it enriches the knowledge on the process of evolution but also because the results generated could have implications for selective breeding programs.

Food for Healthy Living and Active Ageing

The link between diet and health has been recognized since the Grecian period; as Hippocrates said, "Let food be your medicine and medicine be your food". Although the primary goals of diet are meeting nutritional requirements and providing energy, there is increasing awareness that a correct and balanced diet may prevent the insurgence of diet-related pathologies and/or improve well-being and life expectancy, also reflecting on the ageing process. Research on the interaction among nutrients, gut microbiota and host metabolism is presently unravelling the molecular mechanisms underlying the positive and negative effects of traditional diets on health and ageing, providing useful information for the design of innovative foods targeting specific needs and segments of the population. The food supply chain plays a key role in ensuring quality and safety through both comprehensive quality management and inspection systems and a focused innovation process mainly devoted to the creation of functional foods. However, innovation and scientific development pose a problem of information asymmetry towards final consumers; thus, regulatory aspects and private and public communication strategies must be efficiently developed.

Nutrition and Ageing

The world elderly population is rapidly increasing. This demographic change represents a new challenge for the society and demands for a multisectorial intervention to promote a long, healthy, and active life span. Between the factors that contribute in fostering a long healthy life, the nutritional regime plays a central role and is recognized as a major factor in the onset of chronic diseases. A better understanding of the interaction between nutrition and ageing is essential to unravel the mechanisms responsible for these positive/negative effects and to identify diet components promoting the quality of life in the old age and to contribute to the prevention of late-life disabilities. At Università Cattolica del Sacro Cuore, the research activity in food science is focusing on four main objectives: food quality, food safety, functional foods and diet balancing. These objectives are the target of multidisciplinary ongoing and future research activities for a better understanding of the link between diet and ageing. Briefly, the different activities are addressed to the study of the following subjects: the most relevant factors affecting food choices and habits of old aged persons; the effects of long term low dose supplementation of conjugated linoleic acid in mouse; the use of low glycemic index and high resistant starch foods to prevent diabetes and obesity; the adjuvant effect of food bacteria for vaccination; the role of food ingredients in disease; the immunosuppression effect of mycotoxins, and its relevance in ageing people; the production of sustainable and natural antioxidant ingredients to encourage a healthy diet. Our research projects emphasize an holistic and integrated approach that, by bringing together complementary research groups, can combine the collective expertise and thus provide a comprehensive assessment of the role of nutrition in healthy ageing people.

Genetic diversity of sheep breeds from Albania, Greece, and Italy assessed by mitochondrial DNA and nuclear polymorphisms (SNPs)

We employed mtDNA and nuclear SNPs to investigate the genetic diversity of sheep breeds of three countries of the Mediterranean basin: Albania, Greece, and Italy. In total, 154 unique mtDNA haplotypes were detected by means of D-loop sequence analysis. The major nucleotide diversity was observed in Albania. We identified haplogroups, A, B, and C in Albanian and Greek samples, while Italian individuals clustered in groups A and B. In general, the data show a pattern reflecting old migrations that occurred in postneolithic and historical times. PCA analysis on SNP data differentiated breeds with good correspondence to geographical locations. This could reflect geographical isolation, selection operated by local sheep farmers, and different flock management and breed admixture that occurred in the last centuries.

Landscape genomics and biased FST approaches reveal single nucleotide polymorphisms under selection in goat breeds of North-East Mediterranean

Background

In this study we compare outlier loci detected using a FST based method with those identified by a recently described method based on spatial analysis (SAM). We tested a panel of single nucleotide polymorphisms (SNPs) previously genotyped in individuals of goat breeds of southern areas of the Mediterranean basin (Italy, Greece and Albania). We evaluate how the SAM method performs with SNPs, which are increasingly employed due to their high number, low cost and easy of scoring.

Results

The combined use of the two outlier detection approaches, never tested before using SNP polymorphisms, resulted in the identification of the same three loci involved in milk and meat quality data by using the two methods, while the FST based method identified 3 more loci as under selection sweep in the breeds examined.

Conclusion

Data appear congruent by using the two methods for FST values exceeding the 99% confidence limits. The methods of FST and SAM can independently detect signatures of selection and therefore can reduce the probability of finding false positives if employed together. The outlier loci identified in this study could indicate adaptive variation in the analysed species, characterized by a large range of climatic conditions in the rearing areas and by a history of intense trade, that implies plasticity in adapting to new environments.

Double muscling in Marchigiana beef breed is caused by a stop codon in the third exon of myostatin gene

Double muscling is a partially recessive trait present in some beef breeds. It shows a high frequency in some breeds, while in others the frequency is low, and double-muscled individuals are rare. The double muscling is caused by an allelic series of mutations that cause a loss of function of the myostatin gene ( GDF8). We describe here a new mutation in the myostatin gene in Marchigiana breed, a typical beef breed of Central Italy, in which rare double-muscling individuals have been described. A PCR product of the third exon was sequenced in subjects phenotypically showing double muscling, and a G > T transversion was discovered that introduces a premature stop codon. The variant found adds to the large series of mutations present in cattle, and particularly to the only two causative of double muscling in the third exon. A PCR-RFLP test is described for the rapid and effective identification of both heterozygous and homozygous subjects. It was applied to a larger survey carried on the same and also in two other beef breeds, Chianina and Romagnola. Further individuals carrying the new variant were found in Marchigiana, but none in the other breeds. The results may be important for a better comprehension of the role of myostatin in muscular development, for commercial use and for the inference of phylogeny of this gene.