Comparative gene expression profiling of milk somatic cells of Sahiwal cattle and Murrah buffaloes

Adipose gene expression profiles in Northern Finncattle, Mirandesa cattle, Yakutian cattle and commercial Holstein cattle

The drastic change in global climate has led to in-depth studies of the geneticresources of native cattle adapted to challenging environments. Native cattle breeds may harbor unique genetic mechanisms that have enabled them adapt to their given environmental conditions. Adipose tissues are key factors in the regulation of metabolism and energy balance and are crucial for the molecular switches needed to adapt to rapid environmental and nutritional changes. The transcriptome landscape of four adipose tissues was used in this study to investigate the differential gene expression profiles in three local breeds, Yakutian cattle (Sakha Republic), Northern Finncattle (Finland), Mirandesa cattle (Portugal) and commercial Holstein cattle. A total of 26 animals (12 cows, 14 bulls) yielded 81 samples of perirenal adipose tissue (n = 26), metacarpal adipose tissue (n = 26), tailhead adipose tissue (n = 26) and prescapular adipose tissue (n = 3). More than 17,000 genes were expressed in our dataset. Principal component analysis of the normalized expression profiles revealed a differential expression profile of the metacarpal adipose tissue. We found that the genes upregulated in the metacarpal adipose tissue of Yakutian cattle, such as NR4A3, TEKT3, and FGGY, were associated with energy metabolism and response to cold temperatures. In Mirandesa cattle, the upregulated genes in perirenal adipose tissue were related to immune response and inflammation (AVPR2, CCN1, and IL6), while in Northern Finncattle, the upregulated genes appeared to be involved in various physiological processes, including energy metabolism (IGFBP2). According to the sex-based comparisons, the most interesting result was the upregulation of the TPRG1 gene in three tissues of Yakutian cattle females, suggesting that adaptation is related to feed efficiency. The highest number of differentially expressed genes was found between Yakutian cattle and Holstein, several of which were associated with immunity in Yakutian cattle, indicating potential differences in disease resistance and immunity between the two breeds. This study highlights the vast difference in gene expression profiles in adipose tissues among breeds from different climatic environments, most likely highlighting selective pressure and the potential significance of the uniquely important regulatory functions of metacarpal adipose tissue.

Yak DEFB123 alleviates lung injury caused by Klebsiella pneumoniae through MAPKs signaling pathway

Beta-defensins, such as β-defensin 123 (DEFB123), are vital components of the immune system's defense against infections due to their strong antimicrobial properties and capacity for modulating the body's immunological responses. In this study, we successfully cloned and analyzed the yak DEFB123 gene sequence. Subsequently, we obtained recombinant protein DEFB123 (rDEFB123) through prokaryotic expression. Our results demonstrate that rDEFB123 effectively inhibits the growth of Escherichia coli, Klebsiella pneumoniae, and Staphylococcus aureus. Furthermore, rDEFB123 enhances the phagocytic activity of macrophages by regulating specific factors. In a mouse model infected with Klebsiella pneumoniae, the administration of rDEFB123 showed significantly lower levels of serum ALT and AST compared to the control group. Moreover, IFN-γ and IgG were significantly increased in the rDEFB123-treated groups, indicating an enhanced immune response. In the MAPKs signaling pathway of the infected mouse lungs, the expressions of JNK, TRAF2, TRAF6, MIF, and IL-1β genes were downregulated in the rDEFB123-treated groups. Moreover, the levels of p-JNK protein were significantly decreased in these groups as well. Klebsiella pneumoniae caused systemic infection with organ damage in mice. However, the administration of rDEFB123 suppressed the expressions of inflammatory factors, thereby mitigating organ injury and regulating the activity of apoptosis-related factors to enhance immunity. These findings provide valuable theoretical data for future exploration of the functionality and potential applications of DEFB123 in yak.

Identification of candidate genes associated with primary feathers of tianfu nonghua ducks based on Genome-wide association studies

The primary feathers of ducks have important economic value in the poultry industry. This study quantified the primary feather phenotype of Nonghua ducks, including the primary feathers' length, area, distribution of black spots, and feather symmetry. And genome-wide association analysis was used to screen candidate genes that affect the primary feather traits. The genome-wide association study (GWAS) results identified the genetic region related to feather length (FL) on chromosome 2. Through Linkage disequilibrium (LD) analysis, candidate regions (chr2: 115,246,393-116,501,448 bp) were identified and were further annotated to 5 genes: MRS2, GPLD1, ALDH5A1, KIAA0319, and ATP9B. Secondly, candidate regions related to feather black spots were identified on chromosome 21. Through LD analysis, the candidate regions (chr21: 163,552-2,183,853 bp) were screened and further annotated to 47 genes. Among them, STK4, CCN5, and YWHAB genes were related to melanin-related pathways or pigment deposition, which may be key genes affecting the distribution of black spots on feathers. In addition, we also screened 125 genes on multiple chromosomes that may be related to feather symmetry. Among them, significant SNPs on chromosome 1 were further identified as candidate regions (chr1: 142,118,209-142,223,605 bp) through LD analysis and annotated into 2 genes, TGFBRAP1 and LOC113839965. These results reported the genetic basis of the primary feather from multiple phenotypes, and offered valuable insights into the genetic basis for the growth and development of duck feathers and feather color pattern.

Mammary fat globules as a source of mRNA to model alterations in the expression of some milk component genes during lactation in bovines

BACKGROUND

The milk's nutritional value is determined by its constituents, including fat, protein, carbohydrates, and minerals. The mammary gland's ability to produce milk is controlled by a complex network of genes. Thereby, the fat, protein, and lactose synthesis must be boost in milk to increase milk production efficiency. This can be accomplished by fusing genetic advancements with proper management practices. Therefore, this study aimed to investigate the association between the Lipoprotein lipase (LPL), kappa casein CSN3, and Glucose transporter 1 (GLUT1) genes expression levels and such milk components as fat, protein, and lactose in different dairy breeds during different stages of lactation.

METHODS

To achieve such a purpose, 94 milk samples were collected (72 samples from 36 multiparous black-white and red-white Holstein-Friesian (HF) cows and 22 milk samples from 11 Egyptian buffaloes) during the early and peak lactation stages. The milk samples were utilized for milk analysis and genes expressions analyses using non- invasive approach in obtaining milk fat globules (MFGs) as a source of Ribonucleic acid (RNA).

RESULTS

LPL and CSN3 genes expressions levels were found to be significantly higher in Egyptian buffalo than Holstein-Friesian (HF) cows as well as fat and protein percentages. On the other hand, GLUT1 gene expression level was shown to be significantly higher during peak lactation than early lactation. Moreover, lactose % showed a significant difference in peak lactation phase compared to early lactation phase. Also, fat and protein percentages were significantly higher in early lactation period than peak lactation period but lactose% showed the opposite pattern of Egyptian buffalo.

CONCLUSION

Total RNA can be successfully obtained from MFGs. The results suggest that these genes play a role in glucose absorption and lactose synthesis in bovine mammary epithelial cells during lactation. Also, these results provide light on the differential expression of these genes among distinct Holstein-Friesian cow breeds and Egyptian buffalo subspecies throughout various lactation phases.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

Revelation of genes associated with energy generating metabolic pathways in the fighter type Aseel chicken of India through skeletal muscle transcriptome sequencing

In this study, comparative analysis of skeletal muscle transcriptome was carried out for four biological replicates of Aseel, a fighter type breed and Punjab Brown, a meat type breed of India. The profusely expressed genes in both breeds were related to muscle contraction and motor activity. Differential expression analysis identified 961 up-regulated and 979 down-regulated genes in Aseel at a threshold of log2 fold change ≥ ±2.0 (padj<0.05). Significantly enriched KEGG pathways in Aseel included metabolic pathways and oxidative phosphorylation, with higher expression of genes associated with fatty acid beta-oxidation, formation of ATP by chemiosmotic coupling, response to oxidative stress, and muscle contraction. The highly connected hub genes identified through gene network analysis in the Aseel gamecocks were HNF4A, APOA2, APOB, APOC3, AMBP, and ACOT13, which are primarily associated with energy generating metabolic pathways. The up-regulated genes in Punjab Brown chicken were found to be related to muscle growth and differentiation. There was enrichment of pathways such as focal adhesion, insulin signaling pathway and ECM receptor interaction in these birds. The results presented in this study help to improve our understanding of the molecular mechanisms associated with fighting ability and muscle growth in Aseel and Punjab Brown chicken, respectively.

Invited review: From heat stress to disease-Immune response and candidate genes involved in cattle thermotolerance

Heat stress implies unfavorable effects on primary and functional traits in dairy cattle and, in consequence, on the profitability of the whole production system. The increasing number of days with extreme hot temperatures suggests that it is imperative to detect the heat stress status of animals based on adequate measures. However, confirming the heat stress status of an individual is still challenging, and, in consequence, the identification of novel heat stress biomarkers, including molecular biomarkers, remains a very relevant issue. Currently, it is known that heat stress seems to have unfavorable effects on immune system mechanisms, but this information is of limited use in the context of heat stress phenotyping. In addition, there is a lack of knowledge addressing the molecular mechanisms linking the relevant genes to the observed phenotype. In this review, we explored the potential molecular mechanisms explaining how heat stress affects the immune system and, therefore, increases the occurrence of immune-related diseases in cattle. In this regard, 2 relatively opposite hypotheses are under focus: the immunosuppressive action of cortisol, and the proinflammatory effect of heat stress. In both hypotheses, the modulation of the immune response during heat stress is highlighted. Moreover, it is possible to link candidate genes to these potential mechanisms. In this context, immune markers are very valuable indicators for the detection of heat stress in dairy cattle, broadening the portfolio of potential biomarkers for heat stress.

Expression profile of different classes of proteases in milk derived somatic cells across different lactation stages of indigenous cows (Bos indicus) and riverine buffaloes (Bubalus bubalis)

Proteases play a significant role in milk and its products by affecting flavor, texture and longevity. The expression of endogenous proteases varies across different stages of lactation. The study was conducted to understand the transcriptional pattern of different classes of protease-pathways associated genes (CTSB, CTSD, CTSH, CTSL, CTSK, CTSS, CTSZ, PLAU, PLAT) and potential protease inhibitors (SERPIN E2 and SERPIN F2) in 40 milk somatic cells (MSC) samples isolated during early, peak, mid and late lactation stages of Sahiwal cows and Murrah buffaloes - the two most important dairy breeds of India. In Sahiwal cows, except CTSK and PLAU, the expression of other proteases class was not affected significantly (p > 0.05) across lactation stages. However, in Murrah buffaloes, the expression of different proteases increased as the lactation progressed. Most of the proteases showed lower expression during early and peak lactation stages while their expression tends to increase during mid to late lactation stages. The overall trend was somewhat similar in both the dairy species albeit the level of expression was higher in buffalo MSC as compared to cow MSC. The study has provided valuable information on expression kinetics of different proteases in milk somatic cells of two major dairy breeds of India.

Identification of Key Genes and Potential Pathways Associated with Mastitis Induced by E. coli

Mastitis is one of the most important diseases of dairy cattle. It is an infectious disease leading to an inflammatory reaction in the cow's mammary gland. Escherichia coli is one of the common bacteria which induce mastitis in cows. The aim of this study was to identify key genes and potential pathways associated with mastitis induced by E. coli in dairy cattle using bioinformatics analysis. The gene expression profile of ten samples including five adjacent tissues from a quarter infected with Escherichia coli and five tissues from a healthy quarter of dairy cattle was assessed using GEO2R. Gene ontology and pathway analysis were performed using bioinformatics tools. A total of 156 differentially expressed genes were detected which 95 genes were upregulated and 61 genes were downregulated in adjacent tissue of quarter infected compared with healthy tissue. Cellular oxidant detoxification and oxidation-reduction process were the most significant biological process terms in gene ontology analysis. The most important pathways of DEGs were the biosynthesis of amino acids, p53 signaling pathway, and Metabolic pathways. Three important modules were identified and their path enrichment analysis was performed. There are 10 core genes, among which SOD2, COL1A2, COL3A1, POSTN, ALDH18A1, and CBS may be the main genes associated with mastitis, which can be considered as candidate genes in the prevention and carly diagnosis program of mastitis.

First report on delineation of differentially expressed genes and pathways in milk somatic cells of mastitic and healthy Murrah buffaloes

Despite immense contribution of buffaloes as dairy species, limited studies have addressed the bubaline mastitis as compared to cattle. This was the first differential transcriptomic study investigating the alterations induced by clinical mastitis in buffalo milk relative to healthy controls. Comparative gene expression profiling of three biological replicates of each group identified 1014 up-regulated and 999 down-regulated genes in the diseased buffaloes (Fold change > 2, FDR < 0.05). Activation of immune and inflammatory responses were the most enriched GO terms in the mastitic animals, with higher transcript abundance of many genes coding for anti-microbial proteins such as β-defensins, perforin, granzymes, granulysin, cathelicidins etc. Analysis of the gene regulatory interactions of the up-regulated DEGs identified many hub genes that govern the cellular and macromolecular metabolic processes (E2F4, E2F1, RBL2, FOXM1, IRF1 and MYB). This study contributes to an insightful understanding of molecular mechanisms governing immune response of buffaloes to mastitis.

Changes in the lipidome of water buffalo milk during intramammary infection by non-aureus Staphylococci

This study aimed to determine the lipidome of water buffalo milk with intramammary infection (IMI) by non-aureus staphylococci (NAS), also defined as coagulase-negative staphylococci, using an untargeted lipidomic approach. Non-aureus Staphylococci are the most frequently isolated pathogens from dairy water buffalo milk during mastitis. A total of 17 milk samples from quarters affected by NAS-IMI were collected, and the lipidome was determined by liquid chromatography-quadrupole time-of-flight mass spectrometry. The results were compared with the lipidome determined on samples collected from 16 healthy quarters. The study identified 1934 different lipids, which were classified into 15 classes. The abundance of 72 lipids changed in NAS-IMI milk compared to healthy quarters. Significant changes occurred primarily in the class of free fatty acids. The results of this study provided first-time insight into the lipidome of dairy water buffalo milk. Moreover, the present findings provide evidence that NAS-IMI induces changes in water buffalo milk's lipidome.

Transcriptional profiling of buffalo mammary gland with different milk fat contents

Milk fat is the most important energy substance in milk and contributes to its quality and health benefits. Water buffalo milk is well known for its high milk quality with higher fat contents compared with cattle milk. Dehong buffalo is a unique local swamp breed in Yunnan Province with higher milk fat and excellent milk quality which provides a good model for the investigation of the molecular mechanisms of milk fat deposition. In this study, we used RNA-seq to obtain mammary tissue transcriptomics of buffalo with different milk fat phenotypes including high(H), medium (M)and low (L) fat content groups. Comparative analyses of buffalo among three groups yielded differentially expressed genes (DEGs). Analyzing the number of different genes among H_VS_L, H_VS_M, and M_VS_L showed the same expression pattern between H_VS_M. The increasing expression levels of genes including CSN1S1, BTN1A1, LALBA, ALDH1L2, SCD and MUC15, and down-regulated expression levels of genes containing CCL2, CRABP2, ADTRP, CLU and C4A in H_VS_L and M_VS_L were found. GO and KEGG enriched pathways revealed these DEGs involved in milk protein and fat as well as immune response. The findings would enhance the understanding of the interplay between the milk composition and immune response, which suggests that the immune capacity should be considered when we tried to improve the milk quality.