Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo.

Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (FROH) ranged from 0.0201 to 0.0371. Pearson correlations between FROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.

Genome-wide association study for buffalo mammary gland morphology.

This research communication describes a genome-wide association study for Italian buffalo mammary gland morphology. Three single nucleotide polymorphisms (AX-85117983, AX-8509475 and AX-85117518) were identified to be significantly associated with buffalo anterior teat length, posterior teat length and distance between anterior and posterior teat, respectively. Two significant signals for buffalo mammary gland morphology were observed in two genomic regions on the chromosome 10, and chromosome 20. One of the regions located on the chromosome 10 has the most likely candidate genes ACTC1 and GJD2, both of which have putative roles in the regulation of mammary gland development. This study provides new insights into the genetic variants of buffalo mammary gland morphology and may be beneficial for understanding of the genetic regulation of mammary growth.

N-O Reduction and ROS-Mediated AKT/FOXO1 and AKT/P53 Pathways Are Involved in Growth Promotion and Cytotoxicity of Cyadox.

Cyadox is a novel derivative of quinoxaline-1,4-dioxides (QdNOs) with the potential to be developed as a feed additive. However, the pharmacological and toxicological bioactive molecules of cyadox and the molecular mechanism of its pharmacological and toxic actions remain unclear. In the present study, cyadox and its main metabolites of cy1, cy4, cy6, and cy12 were selected; the growth promotion characteristic was indicated by the mRNA level of EGF; and the cytotoxicity of cyadox was determined by methylthiazol tetrazolium bromide (MTT) assay, lactate dehydrogenase (LDH) release, and Annexin V-FITC/PI apoptosis detection kit with flow cytometry. The intracellular ROS, cyclin D1, and Akt/P53/FOXO1 signaling pathway were also investigated. Our data suggested that cyadox showed relatively higher activity than its metabolites, and the ROS was generated from N-O reduction of cyadox. Moreover, cyadox (2 µM) activated the Akt and increased the EGF, cyclin D1, and FOXO1 expression levels. Cyadox (100 µM) induced cytotoxicity in L02 cells in a concentration- and time-dependent manner. Additionally, the activated P53 pathway, hyperactivated Akt, and apoptosis were found in L02 cells after incubation with 100 µM cyadox. Our data demonstrated that Akt promoted cell survival when it was mildly activated by cyadox at 2 µM, and Akt leads to apoptosis when it was severely activated by cyadox at 100 µM. Thus, the present study revealed that N-O reduction of cyadox and ROS-mediated AKT/FOXO1 and AKT/P53 pathways were involved in growth promotion and cytotoxicity of cyadox.

DGAT1 polymorphism in Riverine buffalo, Swamp buffalo and crossbred buffalo.

This Research Communication describes the polymorphisms in the coding region of DGAT1 gene in Riverine buffalo, Swamp buffalo and crossbred buffalo, and associations between polymorphisms and milk production performance in Riverine buffalo. Two polymorphisms of DGAT1were identified, located in exon 13 and exon 17, respectively. The distribution of the genotypes of the two SNP loci in different buffalo population varied, especially the polymorphism located in exon 13 which was not found in the Swamp buffalo. Moreover, SNP located in exon 17 was a nonsynonymous switch resulting in the animo acid sequence changed from an arginine (Arg) to a histidine (His) at position 484. Both SNPs were in Hardy-Weinberg equilibrium, and the polymorphism of g.8330T>C in the exon 13 was significantly associated with peak milk yield, total milk yield and protein percentage. The C variant was associated with an increase in milk yield and peak yield but less in protein percentage compared with the T variant. The polymorphisms of g.9046T>C in exon 17 were significantly associated with fat percentage, in that the buffaloes with TT genotype had a significantly higher fat percentage than those with CC genotype. These findings reveal the difference in the genetic evolution of the DGAT1 between Riverine buffalo and Swamp buffalo, and provide evidence that the DGAT1 gene has potential effects for Riverine buffalo milk production traits, which can be used as a candidate gene for marker-assisted selection in buffalo breeding.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing.

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

Heat Stress Induces Shifts in the Rumen Bacteria and Metabolome of Buffalo.

Exposure to the stress (HS) negatively affects physiology, performance, reproduction and welfare of buffalo. However, the mechanisms by which HS negatively affects rumen bacteria and its associated metabolism in buffalo are not well known yet. This study aimed to gain insight into the adaption of bacteria and the complexity of the metabolome in the rumen of six buffalo during HS using 16S rDNA and gas chromatography metabolomics analyses. HS increased respiratory rate (p < 0.05) and skin temperature (p < 0.01), and it decreased the content of acetic acid (p < 0.05) and butyric acid (p < 0.05) in the rumen. Omics sequencing revealed that the relative abundances of Lachnospirales, Lachnospiraceae, Lachnospiraceae_NK3A20_group and Clostridia_UCG-014 were significantly (p < 0.01) higher under HS than non-heat stress conditions. Several bacteria at different levels, such as Lactobacillales, Streptococcus, Leuconostocaceae and Leissella, were significantly (p < 0.05) more abundant in the rumen of the non-heat stress than HS condition. Thirty-two significantly different metabolites closely related to HS were identified (p < 0.05). Metabolic pathway analysis revealed four key pathways: D-Alanine metabolism; Lysine degradation, Tropane; piperidine and pyridine alkaloid biosynthesis; and Galactose metabolism. In summary, HS may negatively affected rumen fermentation efficiency and changed the composition of rumen community and metabolic function.

Dosage Regimen Formulation and Therapeutic Effect Evaluation of Cyadox Nanosuspension Against Dairy Cow Mastitis Caused by Staphylococcus aureus.

AIMS: In this study, the dosage regimen establishment of cyadox nanosuspension against dairy cow mastitis caused by Staphylococcus aureus (S. aureus) was used as example to provide a general reference for the other novel nanocrystal preparations. METHODS: The effect of cyadox against S. aureus isolates from dairy cows were firstly estimated and then the dosing regimen of nanosuspension after intramammary administration was optimized according to the model of ex vivo pharmacokinetic (PK) and pharmacodynamic (PD). The therapeutic efficacy of the predicted dosage regimen was evaluated. RESULTS: The results demonstrated that cyadox has a concentration-dependent effect on S. aureus. The smallest and highest values of minimum inhibitory concentration (MIC) against 80 isolates was 8 and 64 µg/mL, respectively. The corresponding MIC 50 and MIC 90 was 16 and 32 µg/mL, respectively. The MIC against the pathogenic S. aureus SAHZ156001 in broth and milk were 16 and 32 µg/mL, respectively. The AUC 0-last and C max of cyadox in milk were 4442.877 µg*h/mL and 753.052 µg/mL, respectively. According to the inhibitory sigmoid E max modeling and dosage equation, the daily doses were predicted 1.6, 6.6, and 12.2 mL/gland to achieve bacteriostatic, bactericidal, and elimination effects. The dosage internal was daily administration for continuous three days. CONCLUSION: The clinical experiment showed that the efficient rates were 100, 100, and 90.9%, and the curative rates were 100, 81.8, and 63.6% in 12.2, 6.6 and 1.6 ml/gland groups, respectively. These results showed that cyadox nanosuspension had a good prospect as intramammary infusion to cure dairy cow mastitis infected by S. aureus. This study will be helpful for providing reference for nanocrystal preparation dosage regimen formulation.

Metabolomic and Lipidomic Approaches to Evaluate the Effects of Eucommia ulmoides Leaves on Milk Quality and Biochemical Properties.

Eucommia ulmoides leaves (EUL) contain a variety of natural bioactive compounds including chlorogenic acid, geniposide acid, and aucubin. These bioactive chemicals improve immune function and regulate lipid metabolism. The aim of this study was to investigate the effects of EUL on the biochemical properties of milk. Twenty Holstein dairy cows were randomly allocated to two groups fed a control (CTR, diet without EUL, n = 10) or EUL (diet containing 3% EUL, dry matter, n = 10) diet for 55 d. At the end of the experimental period (d 55), milk samples were collected and analyzed to determine their composition. Though levels of milk fat, protein, lactose, and total milk solids were similar between the groups, small molecules, metabolites, lipids, and cytokines differed. Compared with the CTR group, the EUL group had an improved cluster of differentiation (CD)4/CD8 ratio (P < 0.05) and lower interleukin (IL)-8 and IL-6 content (P < 0.05). Metabolomics analysis identified 14 metabolites including 7Z, 10Z, 13Z, 16Z, 19Z-docosapentaenoic acid (FC = 3.129), adrenic acid (FC = 2.830), and eicosapentaenoic acid (FC=1.685) as having significantly increased in the EUL group (P < 0.05) while 11 metabolites, including indole-2-carboxylic acid (FC = 0.636), cholic acid (FC = 0.430), and creatine (FC = 0.784) had significantly decreased (P < 0.05). Based on a constructed metabolome map, linoleic acid metabolism had the highest impact value for EUL. A total of 21 lipid classes and 1,094 lipid species were detected in the milk by lipidomic analysis, among which 40 differed significantly between the CTR and EUL groups. The present findings showed that the EUL altered milk composition. Correlation analysis showed that 7Z, 10Z, 13Z, 16Z, 19Z-docosapentaenoic acid, adrenic acid, and eicosapentaenoic acid levels were negatively correlated with those of the inflammatory factors IL-6 and IL-8 (P < 0.05), indicating that EUL improved milk quality by reducing inflammatory factors and increasing the CD4/CD8 ratio. Overall, our data demonstrate that EUL had positive effects on milk antioxidant parameters, immune indices, and micro-composition metabolism, thereby improving milk quality.

Transcriptome-Wide m6A Analysis Provides Novel Insights Into Testicular Development and Spermatogenesis in Xia-Nan Cattle.

Testis is the primary organ of the male reproductive tract in mammals that plays a substantial role in spermatogenesis. Improvement of our knowledge regarding the molecular mechanisms in testicular development and spermatogenesis will be reflected in producing spermatozoa of superior fertility. Evidence showed that N6-Methyladenosine (m6A) plays a dynamic role in post-transcription gene expression regulation and is strongly associated with production traits. However, the role of m6A in bovine testis has not been investigated yet. In this study, we conducted MeRIP-Seq analysis to explore the expression profiles of the m6A and its potential mechanism underlying spermatogenesis in nine bovine testes at three developmental stages (prepuberty, puberty and postpuberty). The experimental animals with triplicate in each stage were chosen based on their semen volume and sperm motility except for the prepuberty bulls and used for testes collection. By applying MeRIP-Seq analysis, a total of 8,774 m6A peaks and 6,206 m6A genes among the studied groups were identified. All the detected peaks were found to be mainly enriched in the coding region and 3'- untranslated regions. The cross-analysis of m6A and mRNA expression exhibited 502 genes with concomitant changes in the mRNA expression and m6A modification. Notably, 30 candidate genes were located in the largest network of protein-protein interactions. Interestingly, four key node genes (PLK4, PTEN, EGR1, and PSME4) were associated with the regulation of mammal testis development and spermatogenesis. This study is the first to present a map of RNA m6A modification in bovine testes at distinct ages, and provides new insights into m6A topology and related molecular mechanisms underlying bovine spermatogenesis, and establishes a basis for further studies on spermatogenesis in mammals.

Transcriptome Analysis Reveals Potential Regulatory Genes Related to Heat Tolerance in Holstein Dairy Cattle.

Heat stress affects the physiology and production performance of Chinese Holstein dairy cows. As such, the selection of heat tolerance in cows and elucidating its underlying mechanisms are vital to the dairy industry. This study aimed to investigate the heat tolerance associated genes and molecular mechanisms in Chinese Holstein dairy cows using a high-throughput sequencing approach and bioinformatics analysis. Heat-induced physiological indicators and milk yield changes were assessed to determine heat tolerance levels in Chinese Holstein dairy cows by Principal Component Analysis method following Membership Function Value Analysis. Results indicated that rectal temperature (RT), respiratory rate (RR), and decline in milk production were significantly lower (p < 0.05) in heat tolerant (HT) cows while plasma levels of heat shock protein (HSP: HSP70, HSP90), and cortisol were significantly higher (p < 0.05) when compared to non-heat tolerant (NHT) Chinese Holstein dairy cows. By applying RNA-Seq analysis, we identified 200 (81 down-regulated and 119 up-regulated) significantly (|log2fold change| ≥ 1.4 and p ≤ 0.05) differentially expressed genes (DEGs) in HT versus NHT Chinese Holstein dairy cows. In addition, 14 of which were involved in protein-protein interaction (PPI) network. Importantly, several hub genes (OAS2, MX2, IFIT5 and TGFB2) were significantly enriched in immune effector process. These findings might be helpful to expedite the understanding for the mechanism of heat tolerance in Chinese Holstein dairy cows.

Identifying Hub Genes for Heat Tolerance in Water Buffalo (Bubalus bubalis) Using Transcriptome Data.

Heat stress has a detrimental effect on the physiological and production performance of buffaloes. Elucidating the underlying mechanisms of heat stress is challenging, therefore identifying candidate genes is urgent and necessary. We evaluated the response of buffaloes (n = 30) to heat stress using the physiological parameters, ELISA indexes, and hematological parameters. We then performed mRNA and microRNA (miRNA) expression profiles analysis between heat tolerant (HT, n = 4) and non-heat tolerant (NHT, n = 4) buffaloes, as well as the specific modules, significant genes, and miRNAs related to the heat tolerance identified using the weighted gene co-expression network analysis (WGCNA). The results indicated that the buffaloes in HT had a significantly lower rectal temperature (RT) and respiratory rate (RR) and displayed a higher plasma heat shock protein (HSP70 and HSP90) and cortisol (COR) levels than those of NHT buffaloes. Differentially expressed analysis revealed a total of 753 differentially expressed genes (DEGs) and 16 differentially expressed miRNAs (DEmiRNAs) were identified between HT and NHT. Using the WGCNA analysis, these DEGs assigned into 5 modules, 4 of which were significantly correlation with the heat stress indexes. Interestingly, 158 DEGs associated with heat tolerance in the turquoise module were identified, 35 of which were found within the protein-protein interaction network. Several hub genes (IL18RAP, IL6R, CCR1, PPBP, IL1B, and IL1R1) were identified that significantly enriched in the Cytokine-cytokine receptor interaction. The findings may help further elucidate the underlying mechanisms of heat tolerance in buffaloes.

Integrative Analysis of Transcriptome and GWAS Data to Identify the Hub Genes Associated With Milk Yield Trait in Buffalo.

The mammary gland is the production organ in mammals that is of great importance for milk production and quality. However, characterization of the buffalo mammary gland transcriptome and identification of the valuable candidate genes that affect milk production is limited. Here, we performed the differential expressed genes (DEGs) analysis of mammary gland tissue on day 7, 50, 140, and 280 after calving and conducted gene-based genome-wide association studies (GWAS) of milk yield in 935 Mediterranean buffaloes. We then employed weighted gene co-expression network analysis (WGCNA) to identify specific modules and hub genes related to milk yield based on gene expression profiles and GWAS data. The results of the DEGs analysis showed that a total of 1,420 DEGs were detected across different lactation points. In the gene-based analysis, 976 genes were found to have genome-wide association (P ≤ 0.05) that could be defined as the nominally significant GWAS geneset (NSGG), 9 of which were suggestively associated with milk yield (P < 10-4). Using the WGCNA analysis, 544 and 225 genes associated with milk yield in the turquoise module were identified from DEGs and NSGG datasets, respectively. Several genes (including BNIPL, TUBA1C, C2CD4B, DCP1B, MAP3K5, PDCD11, SRGAP1, GDPD5, BARX2, SCARA3, CTU2, and RPL27A) were identified and considered as the hub genes because they were involved in multiple pathways related to milk production. Our findings provide an insight into the dynamic characterization of the buffalo mammary gland transcriptome, and these potential candidate genes may be valuable for future functional characterization of the buffalo mammary gland.

Transcriptome profiling of anti-müllerian hormone treated preantral/small antral mouse ovary follicles.

The predisposition for the initiation of folliculogenesis in mammals including humans is programmed to start at fetal life and continues until reproductive capacity. The follicles grow from a pool of primordial follicles which retain the major functions in the entire reproductive life of a female. Anti-müllerian hormone (AMH), a glycoprotein belonging to the transforming growth factor-beta family, has an inhibitory effect on ovarian follicle development. The key regulatory target genes in primordial follicle development are of paramount importance in reproductive biology of female. A systems biology method was used to find regulatory genes performing critical role in primordial follicle development. A complete in-depth bioinformatics analysis was performed to investigate the changes in transcriptome of preantral to small antral mouse follicles treated for 12 h and 24 h with two different concentrations; 50 and 200 ng/ml of AMH, and thereby identify candidate genes in time and concentration manner. Firstly, we found differentially expressed genes that were time and concentration dependent in response to AMH. The network analysis of these differentially expressed genes provided new candidate genes and pathways associated with inhibitory action of AMH on the primordial follicle development. To further emphasize the function of AMH, the key identified genes' protein-protein docking was analyzed and found the intracellular and extracellular protein-protein interaction. This study elucidates one of the novel mechanisms of AMH involvement in inhibition of ovarian follicle development which may lead to prolong productive life in female.

Genome-Wide SNP Data Revealed the Extent of Linkage Disequilibrium, Persistence of Phase and Effective Population Size in Purebred and Crossbred Buffalo Populations.

Linkage disequilibrium (LD) is a useful parameter for guiding the accuracy and power of both genome-wide association studies (GWAS) and genomic selection (GS) among different livestock species. The present study evaluated the extent of LD, persistence of phase and effective population size (Ne) for the purebred (Mediterranean buffalo; n = 411) and crossbred [Mediterranean × Jianghan × Nili-Ravi buffalo, n = 9; Murrah × Nili-Ravi × local (Xilin or Fuzhong) buffalo, n = 36] buffalo populations using the 90K Buffalo SNP genotyping array. The results showed that the average square of correlation coefficient (r 2) between adjacent SNP was 0.13 ± 0.19 across all autosomes for purebred and 0.09 ± 0.13 for crossbred, and the most rapid decline in LD was observed over the first 200 kb. Estimated r 2 ≥ 0.2 extended up to ~50 kb in crossbred and 170 kb in purebred populations, while average r 2 values ≥0.3 were respectively observed in the ~10 and 60 kb in the crossbred and purebred populations. The largest phase correlation (R P, C = 0.47) was observed at the distance of 100 kb, suggesting that this phase was not actively preserved between the two populations. Estimated Ne for the purebred and crossbred population at the current generation was 387 and 113 individuals, respectively. These findings may provide useful information to guide the GS and GWAS in buffaloes.