Hepatic transcript profiling in beef cattle: Effects of feeding endophyte-infected tall fescue seeds.

The objective of our study was to evaluate the effect of endophyte-infected tall fescue (E+) seeds intake on liver tissue transcriptome in growing Angus × Simmental steers and heifers through RNA-seq analysis. Normal weaned calves (~8 months old) received either endophyte-free tall fescue (E-; n = 3) or infected tall fescue (E+; n = 6) seeds for a 30-d period. The diet offered was ad libitum bermudagrass (Cynodon dactylon) hay combined with a nutritional supplement of 1.61 kg (DM basis) of E+ or E- tall fescue seeds, and 1.61 kg (DM basis) of energy/protein supplement pellets for a 30-d period. Dietary E+ tall fescue seeds were included in a rate of 20 µg of ergovaline/kg BW/day. Liver tissue was individually obtained through biopsy at d 30. After preparation and processing of the liver samples for RNA sequencing, we detected that several metabolic pathways were activated (i.e., upregulated) by the consumption of E+ tall fescue. Among them, oxidative phosphorylation, ribosome biogenesis, protein processing in endoplasmic reticulum and apoptosis, suggesting an active mechanism to cope against impairment in normal liver function. Interestingly, hepatic protein synthesis might increase due to E+ consumption. In addition, there was upregulation of "thermogenesis" KEGG pathway, showing a possible increase in energy expenditure in liver tissue due to consumption of E+ diet. Therefore, results from our study expand the current knowledge related to liver metabolism of growing beef cattle under tall fescue toxicosis.

Natural clines and human management impact the genetic structure of Algerian honey bee populations.

BACKGROUND: The Algerian honey bee population is composed of two described subspecies A. m. intermissa and A. m. sahariensis, of which little is known regarding population genomics, both in terms of genetic differentiation and of possible contamination by exogenous stock. Moreover, the phenotypic differences between the two subspecies are expected to translate into genetic differences and possible adaptation to heat and drought in A. m. sahariensis. To shed light on the structure of this population and to integrate these two subspecies in the growing dataset of available haploid drone sequences, we performed whole-genome sequencing of 151 haploid drones. RESULTS: Integrated analysis of our drone sequences with a similar dataset of European reference populations did not detect any significant admixture in the Algerian honey bees. Interestingly, most of the genetic variation was not found between the A. m. intermissa and A. m. sahariensis subspecies; instead, two main genetic clusters were found along an East-West axis. We found that the correlation between genetic and geographic distances was higher in the Western cluster and that close-family relationships were mostly detected in the Eastern cluster, sometimes at long distances. In addition, we selected a panel of 96 ancestry-informative markers to decide whether a sampled bee is Algerian or not, and tested this panel in simulated cases of admixture. CONCLUSIONS: The differences between the two main genetic clusters suggest differential breeding management between eastern and western Algeria, with greater exchange of genetic material over long distances in the east. The lack of detected admixture events suggests that, unlike what is seen in many places worldwide, imports of queens from foreign countries do not seem to have occurred on a large scale in Algeria, a finding that is relevant for conservation purposes. In addition, the proposed panel of 96 markers was found effective to distinguish Algerian from European honey bees. Therefore, we conclude that applying this approach to other taxa is promising, in particular when genetic differentiation is difficult to capture.

Alterations in ileal transcriptomics during an intestinal barrier challenge in lactating Holstein cows fed a Saccharomyces cerevisiae fermentation product identify potential regulatory processes.

Stressors such as lack of access to feed, hot temperatures, transportation, and pen changes can cause impairment of ruminal and intestinal barrier function, also known as "leaky gut". Despite the known benefits of some nutritional approaches during periods of stress, little is understood regarding the underlying mechanisms, especially in dairy cows. We evaluated the effect of feeding a Saccharomyces cerevisiae fermentation product (SCFP; NutriTek, Diamond V, Cedar Rapids, IA) on the ileal transcriptome in response to feed restriction (FR), an established model to induce intestinal barrier dysfunction. Multiparous cows [97.1 ±â€…7.6 days in milk (DIM); n = 5/group] fed a control diet or control plus 19 g/d SCFP for 9 wk were subjected to an FR challenge for 5 d during which they were fed 40% of their ad libitum intake from the 7 d before FR. All cows were slaughtered at the end of FR, and ileal scrapping RNA was used for RNAseq (NovaSeq 6000, 100 bp read length). Statistical analysis was performed in R and bioinformatics using the KEGG (Kyoto Encyclopedia of Genes and Genomes) and GO databases. One thousand six hundred and ninety-six differentially expressed genes (DEG; FDR-adjusted P ≤ 0.10) were detected in SCFP vs. control, with 451 upregulated and 1,245 downregulated. "Mucin type O-glycan biosynthesis" was the top downregulated KEGG pathway due to downregulation of genes catalyzing glycosylation of mucins (GCNT3, GALNT5, B3GNT3, GALNT18, and GALNT14). An overall downregulation of cell and tissue structure genes (e.g., extracellular matrix proteins) associated with collagen (COL6A1, COL1A1, COL4A1, COL1A2, and COL6A2), laminin (LAMB2), and integrins (ITGA8, ITGA2, and ITGA5) also were detected with SCFP. A subset of DEG enriched in the GO term "extracellular exosome" and "extracellular space". Chemokines within "Cytokine-cytokine receptor interaction pathways" such as CCL16, CCL21, CCL14, CXCL12, and CXCL14 were downregulated by SCFP. The "Glutathione metabolism" pathway was upregulated by SCFP, including GSTA1 and RRM2B among the top upregulated genes, and GSTM1 and GPX8 as top downregulated genes. There were 9 homeobox transcription factors among the top 50 predicted transcription factors using the RNAseq DEG dataset, underscoring the importance of cell differentiation as a potential target of dietary SCFP. Taken together, SCFP downregulated immune-, ECM-, and mucin synthesis-related genes during FR. Homeobox transcription factors appear important for the transcriptional response of SCFP.

Stressors such as lack of access to feed, hot temperatures, transportation, and disease contribute to diminished gut epithelial barrier integrity in livestock. RNA-sequencing technology and bioinformatics were used to evaluate genome-wide mRNA abundance profiles in ileal tissue from dairy cows fed Saccharomyces cerevisiae fermentation product (SCFP) or an unsupplemented control diet during an intestinal challenge induced by feed restriction. Molecular responses were characterized according to metabolic pathways and other biological categories. Genes associated with "Mucin type O-glycan biosynthesis" and "Extracellular matrix-receptor interaction" were downregulated due to SCFP relative to controls. Alterations in cytokine and chemokine mRNA profiles induced by SCFP underscored differences in tissue immune response. Overall, SCFP altered the transcriptome of ileal tissue damaged by feed restriction.

Hepatic transcript profiling in beef cattle: Effects of rumen-protected niacin supplementation.

The objective of our study was to assess the effect of rumen-protected niacin supplementation on the transcriptome of liver tissue in growing Angus × Simmental steers and heifers through RNA-seq analysis. Consequently, we wanted to assess the known role of niacin in the physiological processes of vasodilation, detoxification, and immune function in beef hepatic tissue. Normal weaned calves (~8 months old) were provided either a control diet or a diet supplemented with rumen-protected niacin (6 g/hd/d) for a 30-day period, followed by a liver biopsy. We observed a significant list of changes at the transcriptome level due to rumen-protected niacin supplementation. Several metabolic pathways revealed potential positive effects to the animal's liver metabolism due to administration of rumen-protected niacin; for example, a decrease in lipolysis, apoptosis, inflammatory responses, atherosclerosis, oxidative stress, fibrosis, and vasodilation-related pathways. Therefore, results from our study showed that the liver transcriptional machinery switched several metabolic pathways to a condition that could potentially benefit the health status of animals supplemented with rumen-protected niacin. In conclusion, based on the results of our study, we can suggest the utilization of rumen-protected niacin supplementation as a nutritional strategy could improve the health status of growing beef cattle in different beef production stages, such as backgrounding operations or new arrivals to a feedlot.

Identification of Ancestry Informative Markers in Mediterranean Trout Populations of Molise (Italy): A Multi-Methodological Approach with Machine Learning.

Brown trout (Salmo trutta), like many other freshwater species, is threated by the release in its natural environment of alien species and the restocking with allochthonous conspecific stocks. Many conservation projects are ongoing and several morphological and genetic tools have been proposed to support activities aimed to restore genetic integrity status of native populations. Nevertheless, due to the complexity of degree of introgression reached up after many generations of crossing, the use of dichotomous key and molecular markers, such as mtDNA, LDH-C1* and microsatellites, are often not sufficient to discriminate native and admixed specimens at individual level. Here we propose a reduced panel of ancestry-informative SNP markers (AIMs) to support on field activities for Mediterranean trout management and conservation purpose. Starting from the genotypes data obtained on specimens sampled in the main two Molise's rivers (Central-Southern Italy), a 47 AIMs panel was identified and validated on simulated and real hybrid population datasets, mainly through a Machine Learning approach based on Random Forest classifier. The AIMs panel proposed may represent an interesting and cost-effective tool for monitoring the level of introgression between native and allochthonous trout population for conservation purpose and this methodology could be also applied in other species.

Study of the Fatty Acid Profile of Milk in Different Sheep Breeds: Evaluation by Multivariate Factorial Analysis.

A multivariate analysis was used to investigate the fatty acid (FA) profile in three different Italian sheep breeds: Comisana, Massese, and Sarda. A sample of 852 animals was considered: 118 Massese, 303 Comisana, 431 Sarda. Sarda sheep were divided into two groups, based on their breeding origin (298 and 133 reared in Sardinia and Tuscany, respectively). Sarda sheep, bred both in Sardinia and in Tuscany, were considered in different groups, both because in these two regions most of the sheep of this breed are reared, and because they differ in geographical characteristics and in the farming system. The individual milk FA composition of dairy ewes was analyzed with multivariate factor analysis. The extracted factors were representative of the following eight groups of fatty acids or functions: factor 1 (odd branched fatty acids and long-chain fatty acids), factor 2 (sn3_position), factor 3 (alternative biohydrogenation), factor 4 (SCD_1), factor 5 (SCD_2), factor 6 (SCD_3), factor 7 (fat secretion) and factor 8 (omega-3). A factor analysis suggested the presence of different metabolic pathways for de novo short- and medium-chain fatty acids and Δ9-desaturase products. The ANOVA of factor scores highlighted the significant effects of the breed. The results of the present study showed that breed is an important factor in defining the fatty acid profile of milk, combined with the effect of the diet. Breeds reared in the same farming system (Comisana, Massese and Sarda reared in Tuscany) showed significant differences for all the factors extracted. At the same time, we found differences between the Sarda sheep reared in Sardinia and Tuscany, two different regions of Italy.

Presence of Lactic Acid Bacteria in the Intestinal Tract of the Mediterranean Trout (Salmo macrostigma) in Its Natural Environment.

Knowledge of the composition of the gut microbiota in freshwater fish living in their natural habitat has taxonomic and ecological importance. Few reports have been produced on the composition of the gut microbiota and on the presence of LAB in the intestines of freshwater fish that inhabit river environments. In this study, we investigated the LAB community that was present in the gastrointestinal tract (GIT) of Mediterranean trout (Salmo macrostigma) that colonized the Biferno and Volturno rivers of the Molise region (Italy). The partial 16S rRNA gene sequences of these strains were determined for the species-level taxonomic placement. The phylogenetic analysis revealed that the isolated LABs belonged to seven genera (Carnobacterium, Enterococcus, Lactobacillus, Lactiplantibacillus, Vagococcus, Lactococcus, and Weissella). The study of the enzymatic activities showed that these LABs could contribute to the breakdown of polysaccharides, proteins, and lipids. In future studies, a greater understanding of how the LABs act against pathogens and trigger the fish immune response may provide practical means to engineer the indigenous fish microbiome and enhance disease control and fish health.

Genotyping of Two Mediterranean Trout Populations in Central-Southern Italy for Conservation Purposes Using a Rainbow-Trout-Derived SNP Array.

Mediterranean trout is a freshwater fish of particular interest with economic significance for fishery management, aquaculture and conservation biology. Unfortunately, native trout populations' abundance is significantly threatened by anthropogenic disturbance. The introduction of commercial hatchery strains for recreation activities has compromised the genetic integrity status of native populations. This work assessed the fine-scale genetic structure of Mediterranean trout in the two main rivers of Molise region (Italy) to support conservation actions. In total, 288 specimens were caught in 28 different sites (14 per basins) and genotyped using the Affymetrix 57 K rainbow-trout-derived SNP array. Population differentiation was analyzed using pairwise weighted FST and overall F-statistic estimated by locus-by-locus analysis of molecular variance. Furthermore, an SNP data set was processed through principal coordinates analysis, discriminant analysis of principal components and admixture Bayesian clustering analysis. Firstly, our results demonstrated that rainbow trout SNP array can be successfully used for Mediterranean trout genotyping. In fact, despite an overwhelming number of loci that resulted as monomorphic in our populations, it must be emphasized that the resulted number of polymorphic loci (i.e., ~900 SNPs) has been sufficient to reveal a fine-scale genetic structure in the investigated populations, which is useful in supporting conservation and management actions. In particular, our findings allowed us to select candidate sites for the collection of adults, needed for the production of genetically pure juvenile trout, and sites to carry out the eradication of alien trout and successive re-introduction of native trout.

Single-Step Genome Wide Association Study Identifies QTL Signals for Untrimmed and Trimmed Thigh Weight in Italian Crossbred Pigs for Dry-Cured Ham Production.

Protected Designation of Origin (PDO) dry-cured ham is the most important product in the Italian pig breeding industry, mainly oriented to produce heavy pig carcasses to obtain hams of the right weight and maturity. Recently, along with the traditional traits swine breeding programs have aimed to include novel carcass traits. The identification at the genome level of quantitative trait loci (QTLs) affecting such new traits helps to reveal their genetic determinism and may provide information to be integrated in prediction models in order to improve prediction accuracy as well as to identify candidate genes underlying such traits. This study aimed to estimate genetic parameters and perform a single step genome wide association studies (ssGWAS) on novel carcass traits such as untrimmed (UTW) and trimmed thigh weight (TTW) in two pig crossbred lines approved for the ham production of the Italian PDO. With this purpose, phenotypes were collected from ~1800 animals and 240 pigs were genotyped with Illumina PorcineSNP60 Beadchip. The single-step genomic BLUP procedure was used for the heritability estimation and to implement the ssGWAS. QTL were characterized based on the variance of 10-SNP sliding window genomic estimated breeding values. Moderate heritabilities were detected and QTL signals were identified on chromosome 1, 4, 6, 7, 11 and 15 for both traits. As expected, the genetic correlation among the two traits was very high (~0.99). The QTL regions encompassed a total of 249 unique candidate genes, some of which were already reported in association with growth, carcass or ham weight traits in pigs. Although independent studies are required to further verify our findings and disentangle the possible effects of specific linkage disequilibrium in our population, our results support the potential use of such new QTL information in future breeding programs to improve the reliability of genomic prediction.

Genetic Regulation of Biomarkers as Stress Proxies in Dairy Cows.

Stress in livestock reduces productivity and is a welfare concern. At a physiological level, stress is associated with the activation of inflammatory responses and increased levels of harmful reactive oxygen species. Biomarkers that are indicative of stress could facilitate the identification of more stress-resilient animals. We examined twenty-one metabolic, immune response, and liver function biomarkers that have been associated with stress in 416 Italian Simmental and 436 Italian Holstein cows which were genotyped for 150K SNPs. Single-SNP and haplotype-based genome-wide association studies were carried out to assess whether the variation in the levels in these biomarkers is under genetic control and to identify the genomic loci involved. Significant associations were found for the plasma levels of ceruloplasmin (Bos taurus chromosome 1-BTA1), paraoxonase (BTA4) and γ-glutamyl transferase (BTA17) in the individual breed analysis that coincided with the position of the genes coding for these proteins, suggesting that their expression is under cis-regulation. A meta-analysis of both breeds identified additional significant associations with paraoxonase on BTA 16 and 26. Finding genetic associations with variations in the levels of these biomarkers suggests that the selection for high or low levels of expression could be achieved rapidly. Whether the level of expression of the biomarkers correlates with the response to stressful situations has yet to be determined.

Unique adaptations in neonatal hepatic transcriptome, nutrient signaling, and one-carbon metabolism in response to feeding ethyl cellulose rumen-protected methionine during late-gestation in Holstein cows.

BACKGROUND: Methionine (Met) supply during late-pregnancy enhances fetal development in utero and leads to greater rates of growth during the neonatal period. Due to its central role in coordinating nutrient and one-carbon metabolism along with immune responses of the newborn, the liver could be a key target of the programming effects induced by dietary methyl donors such as Met. To address this hypothesis, liver biopsies from 4-day old calves (n = 6/group) born to Holstein cows fed a control or the control plus ethyl-cellulose rumen-protected Met for the last 28 days prepartum were used for DNA methylation, transcriptome, metabolome, proteome, and one-carbon metabolism enzyme activities. RESULTS: Although greater withers and hip height at birth in Met calves indicated better development in utero, there were no differences in plasma systemic physiological indicators. RNA-seq along with bioinformatics and transcription factor regulator analyses revealed broad alterations in 'Glucose metabolism', 'Lipid metabolism, 'Glutathione', and 'Immune System' metabolism due to enhanced maternal Met supply. Greater insulin sensitivity assessed via proteomics, and efficiency of transsulfuration pathway activity suggested beneficial effects on nutrient metabolism and metabolic-related stress. Maternal Met supply contributed to greater phosphatidylcholine synthesis in calf liver, with a role in very low density lipoprotein secretion as a mechanism to balance metabolic fates of fatty acids arising from the diet or adipose-depot lipolysis. Despite a lack of effect on hepatic amino acid (AA) transport, a reduction in metabolism of essential AA within the liver indicated an AA 'sparing effect' induced by maternal Met. CONCLUSIONS: Despite greater global DNA methylation, maternal Met supply resulted in distinct alterations of hepatic transcriptome, proteome, and metabolome profiles after birth. Data underscored an effect on maintenance of calf hepatic Met homeostasis, glutathione, phosphatidylcholine and taurine synthesis along with greater efficiency of nutrient metabolism and immune responses. Transcription regulators such as FOXO1, PPARG, E2F1, and CREB1 appeared central in the coordination of effects induced by maternal Met. Overall, maternal Met supply induced better immunometabolic status of the newborn liver, conferring the calf a physiologic advantage during a period of metabolic stress and suboptimal immunocompetence.

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 (h2 ) values showed relevant variations across different factors in each breed (0.03 ≤ h2  ≤ 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.

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.

PANEV: an R package for a pathway-based network visualization.

BACKGROUND: During the last decade, with the aim to solve the challenge of post-genomic and transcriptomic data mining, a plethora of tools have been developed to create, edit and analyze metabolic pathways. In particular, when a complex phenomenon is considered, the creation of a network of multiple interconnected pathways of interest could be useful to investigate the underlying biology and ultimately identify functional candidate genes affecting the trait under investigation. RESULTS: PANEV (PAthway NEtwork Visualizer) is an R package set for gene/pathway-based network visualization. Based on information available on KEGG, it visualizes genes within a network of multiple levels (from 1 to n) of interconnected upstream and downstream pathways. The network graph visualization helps to interpret functional profiles of a cluster of genes. CONCLUSIONS: The suite has no species constraints and it is ready to analyze genomic or transcriptomic outcomes. Users need to supply the list of candidate genes, specify the target pathway(s) and the number of interconnected downstream and upstream pathways (levels) required for the investigation. The package is available at https://github.com/vpalombo/PANEV.

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.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds.

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.

Breed and adaptive response modulate bovine peripheral blood cells' transcriptome.

BACKGROUND: Adaptive response includes a variety of physiological modifications to face changes in external or internal conditions and adapt to a new situation. The acute phase proteins (APPs) are reactants synthesized against environmental stimuli like stress, infection, inflammation. METHODS: To delineate the differences in molecular constituents of adaptive response to the environment we performed the whole-blood transcriptome analysis in Italian Holstein (IH) and Italian Simmental (IS) breeds. For this, 663 IH and IS cows from six commercial farms were clustered according to the blood level of APPs. Ten extreme individuals (five APP+ and APP- variants) from each farm were selected for the RNA-seq using the Illumina sequencing technology. Differentially expressed (DE) genes were analyzed using dynamic impact approach (DIA) and DAVID annotation clustering. Milk production data were statistically elaborated to assess the association of APP+ and APP- gene expression patterns with variations in milk parameters. RESULTS: The overall de novo assembly of cDNA sequence data generated 13,665 genes expressed in bovine blood cells. Comparative genomic analysis revealed 1,152 DE genes in the comparison of all APP+ vs. all APP- variants; 531 and 217 DE genes specific for IH and IS comparison respectively. In all comparisons overexpressed genes were more represented than underexpressed ones. DAVID analysis revealed 369 DE genes across breeds, 173 and 73 DE genes in IH and IS comparison respectively. Among the most impacted pathways for both breeds were vitamin B6 metabolism, folate biosynthesis, nitrogen metabolism and linoleic acid metabolism. CONCLUSIONS: Both DIA and DAVID approaches produced a high number of significantly impacted genes and pathways with a narrow connection to adaptive response in cows with high level of blood APPs. A similar variation in gene expression and impacted pathways between APP+ and APP- variants was found between two studied breeds. Such similarity was also confirmed by annotation clustering of the DE genes. However, IH breed showed higher and more differentiated impacts compared to IS breed and such particular features in the IH adaptive response could be explained by its higher metabolic activity. Variations of milk production data were significantly associated with APP+ and APP- gene expression patterns.