Flow Cytometry-Detected Immunological Markers and on Farm Recorded Parameters in Composite Cow Milk as Related to Udder Health Status.

Flow cytometry is a powerful technology used in many fields of cell biology. It is also used as a routine method to count somatic cells in milk and to characterize bovine milk leukocytes. In this study, we used flow cytometry to simultaneously assess the viability, the percentage of the single subsets of leukocytes and to quantify the expression of CD11b, an immunological marker of cell activation status. Immunological markers were then related with on farm recorded parameters as milk electrical conductivity (MEC) and average milk flow rate (MFR). Composite milk samples were collected from 43 cows, nine of which had naturally infected udders and 34 of which had no infected udders. First, the milk samples were classified according to bacteriological test in positive and negative. The results showed that the negative samples to bacteriological test had: (i) significantly higher percentages of live lymphocytes; (ii) significantly lower percentages of CD11b+ leukocytes; (iii) significantly lower MEC and higher MFR values. Then, samples were classified in three groups according to somatic cell count (SCC): Group A (n = 15), samples with SCC ≤ 100,000 cells/mL, all negative to bacteriological analysis; Group B (n = 11), samples with 100,000 < SCC < 300,000 cells/m, with four samples positive to bacteriological analysis; Group C (n = 17), samples with SCC ≥ 300,000 cell/mL with five samples positive to bacteriological analysis. Multivariate discriminant analysis was used to verify which flow cytometry immunological markers and on farm recorded parameters could better discriminate among the different groups of SCCs. Linear discriminant analysis (LDA) indicated that 5 of the 10 parameters could best be used to reveal the differences between positive and negative samples among the considered groups of SCCs. Furthermore, the Canonical discriminant analysis (CDA) indicated that composite milk samples with different SCC and infection status were clearly separated from each other in a two-dimensional space. In conclusion, the study highlighted that: (1) the conventional flow cytometry analysis applied on milk samples is a useful tool to investigate immunological parameters as potential indicators of udder health; (2) the combined evaluation of live milk leukocytes and recorded farm parameters could be useful to assess udder health status in dairy cows. The results obtained from this pilot study on few data require new and larger trials to be confirmed.

Identification of the complete coding cDNAs and expression analysis of B4GALT1, LALBA, ST3GAL5, ST6GAL1 in the colostrum and milk of the Garganica and Maltese goat breeds to reveal possible implications for oligosaccharide biosynthesis.

BACKGROUND: Milk sialylated oligosaccharides (SOS) play crucial roles in many biological processes. The most abundant free SOS in goat's milk are 3'sialyllactose (3'-SL), 6'sialyllactose (6'-SL) and disialyllactose (DSL). The production of these molecules is determined genetically by the expression of glycosyltransferases and by the availability of nucleotide sugar substrates, but the precise mechanisms regulating the differential patterns of milk oligosaccharides are not known. We aimed to identify the complete cDNAs of candidate genes implicated in SOS biosynthesis (B4GALT1, LALBA, ST3GAL5, ST6GAL1) and to analyse their expression during lactation in the Garganica and Maltese goat breeds. Moreover, we analysed the colostrum and milk contents of 3'-SL, 6'-SL and disialyllactose (DSL) and the possible correlations between expressed genes and SOS. RESULTS: We identified the complete coding cDNAs of B4GALT1 (HQ700335.1), ST3GAL5 (KF055858.2), and ST6GAL1 (HQ709167.1), the single nucleotide polymorphism (SNPs) of these genes and 2 splicing variants of the ST6GAL1 cDNA. RT-qPCR analysis showed that LALBA and ST6GAL1 were the genes with the highest and lowest expression in both breeds, respectively. The interaction effects of the breeds and sampling times were associated with higher levels of B4GALT1 and ST3GAL5 gene expression in Garganica than in Maltese goats at kidding. B4GALT1, LALBA, and ST3GAL5 gene expression changed from kidding to 60 and 120 days in Maltese goats, while in Garganica goats, a difference was observed only for the LALBA gene. Breed and lactation effects were also found for SOS contents. Positive correlations of B4GALT1, LALBA, ST3GAL5, and ST6GAL1 with 3'-SL/6'SL and DSL were found. CONCLUSIONS: The genetic effect on the oligosaccharide content of milk was previously highlighted in bovines, and this study is the first to investigate this effect in two goat breeds (Garganica and Maltese) during lactation. The genetic variability of candidate genes involved in SOS biosynthesis highlights their potential role in affecting gene expression and ultimately biological function. The investigation of gene regulatory regions as well as the examination of other sialyltransferase genes will be needed to identify the genetic pattern leading to a higher SOS content in the autochtonous Garganica breed and to protect it using a focused breeding strategy.

Novel and known signals of selection for fat deposition in domestic sheep breeds from Africa and Eurasia.

Genomic regions subjected to selection frequently show signatures such as within-population reduced nucleotide diversity and outlier values of differentiation among differentially selected populations. In this study, we analyzed 50K SNP genotype data of 373 animals belonging to 23 sheep breeds of different geographic origins using the Rsb (extended haplotype homozygosity) and FST statistical approaches, to identify loci associated with the fat-tail phenotype. We also checked if these putative selection signatures overlapped with regions of high-homozygosity (ROH). The analyses identified novel signals and confirmed the presence of selection signature in genomic regions that harbor candidate genes known to affect fat deposition. Several genomic regions that frequently appeared in ROH were also identified within each breed, but only two ROH islands overlapped with the putative selection signatures. The results reported herein provide the most complete genome-wide study of selection signatures for fat-tail in African and Eurasian sheep breeds; they also contribute insights into the genetic basis for the fat tail phenotype in sheep, and confirm the great complexity of the mechanisms that underlie quantitative traits, such as the fat-tail.

Assessment of Poly(ADP-ribose) Polymerase1 (PARP1) expression and activity in cells purified from blood and milk of dairy cattle.

Poly(ADP-ribosyl)ation (PAR) is a post-translational protein modification catalysed by enzyme member of the poly(ADP-ribose) polymerases (PARPs) family. The activation of several PARPs is triggered by DNA strand breakage and the main PARP enzyme involved in this process is PARP1. Besides its involvement in DNA repair, PARP1 is involved in several cellular processes including transcription, epigenetics, chromatin re-modelling as well as in the maintenance of genomic stability. Moreover, several studies in human and animal models showed PARP1 activation in various inflammatory disorders. The aims of the study were (1) to characterize PARP1 expression in bovine peripheral blood mononuclear cells (PBMC) and (2) to evaluate PAR levels as a potential inflammatory marker in cells isolated from blood and milk samples following different types of infection, including mastitis. Our results show that (i) bovine PBMC express PARP1; (ii) lymphocytes exhibit higher expression of PARP1 than monocytes; (iii) PARP1 and PAR levels were higher in circulating PBMCs of infected cows; (iv) PAR levels were higher in cells isolated from milk with higher Somatic Cell Counts (SCC > 100,000 cells/mL) than in cells from milk with low SCCs. In conclusion, these findings suggest that PARP1 is activated during mastitis, which may prove to be a useful biomarker of mastitis.

New Insights on Water Buffalo Genomic Diversity and Post-Domestication Migration Routes From Medium Density SNP Chip Data.

The domestic water buffalo is native to the Asian continent but through historical migrations and recent importations, nowadays has a worldwide distribution. The two types of water buffalo, i.e., river and swamp, display distinct morphological and behavioral traits, different karyotypes and also have different purposes and geographical distributions. River buffaloes from Pakistan, Iran, Turkey, Egypt, Romania, Bulgaria, Italy, Mozambique, Brazil and Colombia, and swamp buffaloes from China, Thailand, Philippines, Indonesia and Brazil were genotyped with a species-specific medium-density 90K SNP panel. We estimated the levels of molecular diversity and described population structure, which revealed historical relationships between populations and migration events. Three distinct gene pools were identified in pure river as well as in pure swamp buffalo populations. Genomic admixture was seen in the Philippines and in Brazil, resulting from importations of animals for breed improvement. Our results were largely consistent with previous archeological, historical and molecular-based evidence for two independent domestication events for river- and swamp-type buffaloes, which occurred in the Indo-Pakistani region and close to the China/Indochina border, respectively. Based on a geographical analysis of the distribution of diversity, our evidence also indicated that the water buffalo spread out of the domestication centers followed two major divergent migration directions: river buffaloes migrated west from the Indian sub-continent while swamp buffaloes migrated from northern Indochina via an east-south-eastern route. These data suggest that the current distribution of water buffalo diversity has been shaped by the combined effects of multiple migration events occurred at different stages of the post-domestication history of the species.

Genomic Regions in Local Endangered Sheep Encode Potentially Favorable Genes.

The economic evaluation of farm animal genetic resources plays a key role in developing conservation programs. However, to date, the link between diversity as assessed by neutral genetic markers and the functional diversity is not yet understood. Two genome-wide comparisons, using over 44,000 Single Nucleotide Polymorphisms, identified the markers with the highest difference in allele frequency between the Alpago endangered breed and two clusters, composed of four specialized dairy sheep, and four meat breeds respectively. The genes in proximity of these markers were mapped to known pathways of the Gene Ontology to determine which ones were most represented. Our results indicated that the differences of the Alpago breed from the more productive sheep rely upon genes involved in cellular defense and repair mechanisms. A higher number of different markers and genes were detected in the comparison with the specialized dairy sheep. These genes play a role in complex biological processes: metabolic, homeostatic, neurological system, and macromolecular organization; such processes may possibly explain the evolution of gene function as a result of selection to improve milk yield.

Genome assembly and transcriptome resource for river buffalo, Bubalus bubalis (2n = 50).

Water buffalo is a globally important species for agriculture and local economies. A de novo assembled, well-annotated reference sequence for the water buffalo is an important prerequisite for studying the biology of this species, and is necessary to manage genetic diversity and to use modern breeding and genomic selection techniques. However, no such genome assembly has been previously reported. There are 2 species of domestic water buffalo, the river (2 n = 50) and the swamp (2 n = 48) buffalo. Here we describe a draft quality reference sequence for the river buffalo created from Illumina GA and Roche 454 short read sequences using the MaSuRCA assembler. The assembled sequence is 2.83 Gb, consisting of 366 983 scaffolds with a scaffold N50 of 1.41 Mb and contig N50 of 21 398 bp. Annotation of the genome was supported by transcriptome data from 30 tissues and identified 21 711 predicted protein coding genes. Searches for complete mammalian BUSCO gene groups found 98.6% of curated single copy orthologs present among predicted genes, which suggests a high level of completeness of the genome. The annotated sequence is available from NCBI at accession GCA_000471725.1.

RNA-Sequencing for profiling goat milk transcriptome in colostrum and mature milk.

BACKGROUND: In this work we aimed at sequencing and assembling the goat milk transcriptome corresponding at colostrum and 120 days of lactation. To reconstruct transcripts we used both the genome as reference, and a de novo assembly approach. Additionally, we aimed at identifying the differentially expressed genes (DEGs) between the two lactation stages and at analyzing the expression of genes involved in oligosaccharides metabolism. RESULTS: A total of 44,635 different transcripts, organized in 33,757 tentative genes, were obtained using the goat genome as reference. A significant sequence similarity match was found for 40,353 transcripts (90%) against the NCBI NT and for 35,701 (80%) against the NR databases. 68% and 69% of the de novo assembled transcripts, in colostrum and 120 days of lactation samples respectively, have a significant match with the merged transcriptome obtained using Cufflinks/Cuffmerge. CSN2, PAEP, CSN1S2, CSN3, LALBA, TPT1, FTH1, M-SAA3, SPP1, GLYCAM1, EEF1A1, CTSD, FASN, RPS29, CSN1S1, KRT19 and CHEK1 were found between the top fifteen highly expressed genes. 418 loci were differentially expressed between lactation stages, among which 207 and 122 were significantly up- and down-regulated in colostrum, respectively. Functional annotation and pathway enrichment analysis showed that in goat colostrum somatic cells predominate biological processes involved in glycolysis, carbohydrate metabolism, defense response, cytokine activity, regulation of cell proliferation and cell death, vasculature development, while in mature milk, biological process associated with positive regulation of lymphocyte activation and anatomical structure morphogenesis are enriched. The analysis of 144 different oligosaccharide metabolism-related genes showed that most of these (64%) were more expressed in colostrum than in mature milk, with eight expressed at very high levels (SLCA3, GMSD, NME2, SLC2A1, B4GALT1, B3GNT2, NANS, HEXB). CONCLUSIONS: To our knowledge, this is the first study comparing goat transcriptome of two lactation stages: colostrum and 120 days. Our findings suggest putative differences of expression between stages and can be envisioned as a base for further research in the topic. Moreover because a higher expression of genes involved in immune defense response, carbohydrate metabolism and related to oligosaccharide metabolism was identified in colostrum we here corroborate the potential of goat milk as a natural source of lactose-derived oligosaccharides and for the development of functional foods.

Evaluation of leptin receptor expression on buffalo leukocytes.

Experimental evidences support a direct role for leptin in immunity. Besides controlling food intake and energy expenditure, leptin was reported to be involved in the regulation of the immune system in ruminants. The aim of this work was to highlight the expression of leptin receptor (LEPR) on Bubalus bubalis immune cells using a multi-approach assessment: flow cytometry, confocal microscopy and gene expression analysis. Flow cytometric analysis of LEPR expression showed that peripheral blood monocytes were the predominant cells expressing LEPR. This result was corroborated by confocal microscopy and RT-PCR analysis. Moreover, among lymphocytes, LEPR was mainly expressed by B lymphocytes and Natural Killer cells. Evidence of LEPR expression on buffalo blood leukocytes showed to be a good indicator of the responsivity of these cells to leptin, so confirming the involvement of leptin in buffalo immune response.

A selective genotyping approach identifies QTL in a simulated population.

BACKGROUND: Identification of QTLs for important phenotypic traits, through the use of medium-density genome-wide SNP panels, is one of the most challenging areas in animal genetics, for preventing the time-consuming direct sequencing of putative candidate genes, when searching for the mutations that affect the trait. Appropriate statistical analyses allow the identification of genomic regions associated with the investigated trait in the genotyped population. METHODS: The selective genotyping technique was applied to 1000 genotyped animals with known phenotype. Sliding windows composed of five consecutive SNPs were created for each chromosome; we assumed that the QTLs were encoded by the windows showing the highest difference in the frequency of the same alleles between the most divergent productive groups (the two tails of the distribution). RESULTS: Ten windows affected at least one trait. For five of these windows, the highest and significant effect was given by one only SNP, which could therefore be taken as the QTL itself. CONCLUSIONS: In this study we proposed a simple method to identify genomic regions associated to the phenotype under study. The identification of the DNA region is the first step to search for the mutation which is really responsible for the trait variability, through the direct sequencing of the genome regions that encode the QTL.

Signatures of selection identify loci associated with milk yield in sheep.

BACKGROUND: Identification of genomic regions that have been targets of selection for phenotypic traits is one of the most challenging areas of research in animal genetics, particularly in livestock where few annotated genes are available. In this study a genome-wide scan using the Illumina SNP50K Beadchip was performed in the attempt to identify genomic regions associated with milk productivity in sheep. The ovine genomic regions encoding putative candidate genes were compared with the corresponding areas in Bos taurus, as the taurine genome is better annotated. RESULTS: A total of 100 dairy sheep were genotyped on the Illumina OvineSNP50K Beadchip. The Fisher's exact test of significance of differences of allele frequency between each pair of the two tails of the distribution of top/worse milk yielders was performed for each marker. The genomic regions where highly divergent milk yielders showed different allele frequencies at consecutive markers was extracted from the OAR v3.1 Ovine (Texel) Genome Assembly, and was compared to the corresponding areas in Bos taurus, allowing the detection of two genes, the Palmdelphin and the Ring finger protein 145. These genes encoded non-synonymous mutations correlated with the marker alleles. CONCLUSION: The innovation of this study was to show that the DNA genotyping with the Illumina SNP50K Beadchip allowed to detect genes, and mutations in the genes, which have not yet been annotated in the livestock under investigation.

Splicing variants of SERPINA1 gene in ovine milk: characterization of cDNA and identification of polymorphisms.

The serine protease inhibitor, clade A, member 1 (SERPINA1) is the gene for a protein called alpha-1-antitrypsin (AAT), which is a member of the serine protease inhibitor (serpin) superfamily of proteins. By conformational change, serpins control several chemical reactions inhibiting the activity of proteases. AAT is the most abundant endogenous serpin in blood circulation and it is present in relatively high concentration in human milk as well as in bovine and porcine colostrum. Here we report for the first time the molecular characterization and sequence variability of the ovine SERPINA1 cDNA and gene. cDNAs from mammary gland and from milk were PCR amplified, and three different transcripts (1437, 1166 and 521bp) of the SERPINA1 gene were identified. We amplified and sequenced different regions of the gene (5' UTR, from exon 2 to exon 5 and 3' UTR), and we found that the exon-intron structure of the gene is similar to that of human and bovine. We detected a total of 97 SNPs in cDNAs and gene sequences from 10 sheep of three different breeds. In adult sheep tissues a SERPINA1 gene expression analysis indicated a differential expression of the three different transcripts. The finding reported in this paper will aid further studies on possible involvement of the SERPINA1 gene in different physiological states and its possible association with production traits.

Milk fatty acid variability: effect of some candidate genes involved in lipid synthesis.

In this work, the genetic variation of milk FA was investigated in three different bovine breeds, the Jersey, the Piedmontese and the Valdostana, and at different lactation stages. All animals were genotyped for 21 Single Nucleotide Polymorphisms located within nine candidate genes involved in lipid synthesis: diacylglycerol acyltransferase 1 and 2 (DGAT1, 2); stearoyl-CoA desaturase (SCD); growth hormone receptor (GHR); fatty acid synthase (FASN); acyl-CoA dehydrogenase (ACAD); fatty acid binding protein (FABP4); lipoprotein lipase (LPL); and leptin gene (LEP). The highest milk-fat Jersey breed also showed the highest content of saturated FA. Throughout lactation, the breeds showed a similar variation in the FA, with a decrease in the short-chain, this was accompanied by a general increase in the long chain FA at the end of lactation. The increase in long chain saturated FA was particularly evident in the case of the Jersey. The effect of SCD gene on the C14 desaturation index was confirmed; the DGAT1 gene was polymorphic only in the Jersey breed, but its effect was confirmed only on milk fat content; three further potential candidate genes were identified: first, the FABP4 gene, which was found to influence medium and long chain FA in all the breeds, but not the desaturation indices; second, the FASN gene, which was found to influence the amount of PUFA in the Piedmontese and the Valdostana, and third, the LPL gene, which was found to affect fat content in the Piedmontese.