Consumers have shown more and more interest in high-quality and healthy dairy products and buffalo milk is commercially more viable than other milks in producing superior dairy products due to its higher contents of fat, crude protein, and total solids. Metabolomics is one of the most powerful strategies in molecular mechanism research however, little study has been focused on the milk metabolites in different buffalo species. Therefore, the aim of this study was to explore the underlying molecular mechanism of the fatty synthesis and candidate biomarkers by analyzing the metabolomic profiles. Milk of three groups of buffaloes, including 10 Mediterranean, 12 Murrah, and 10 crossbred buffaloes (Murrah × local swamp buffalo), were collected and UPLC-Q-Orbitrap HRMS was used to obtain the metabolomic profiles. Results showed that milk fatty acid in Mediterranean buffalo was significantly higher than Murrah buffalo and crossbred buffalo. A total of 1837/726 metabolites was identified in both positive and negative electrospray ionization (ESI±) mode, including 19 significantly different metabolites between Mediterranean and Murrah buffalo, and 18 different metabolites between Mediterranean and crossbred buffalo. We found 11 of the different metabolites were both significantly different between Mediterranean vs. Murrah group and Mediterranean vs crossbred group, indicating that they can be used as candidate biomarkers of Mediterranean buffalo milk. Further analysis found that the different metabolites were mainly enriched in fat synthesis related pathways such as fatty acid biosynthesis, unsaturated fatty acid biosynthesis, and linoleic acid metabolism, indicating that the priority of different pathways affected the milk fat content in different buffalo species. These specific metabolites may be used as biomarkers in the identification of milk quality and molecular breeding of high milk fat buffalo.

The giant freshwater prawn (Macrobrachium rosenbergii) exhibits sex dimorphism between the male and female individuals. To date, the molecular mechanism governing gonadal development was unclear, and limited data were available on the gonad transcriptome of M. rosenbergii. Here, we conducted comprehensive gonadal transcriptomic analysis of female (ZW), super female (WW), and male (ZZ) M. rosenbergii for gene discovery. A total of 70.33 gigabases (Gb) of sequences were generated. There were 115,338 unigenes assembled with a mean size of 1196 base pair (bp) and N50 of 2195 bp. Alignment against the National Center for Biotechnology Information (NCBI) non-redundant nucleotide/protein sequence database (NR and NT), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, SwissProt database, Protein family (Pfam), Gene ontology (GO), and the eukaryotic orthologous group (KOG) database, 36,282 unigenes were annotated at least in one database. Comparative transcriptome analysis observed that 10,641, 16,903, and 3393 genes were significantly differentially expressed in ZW vs. ZZ, WW vs. ZZ, and WW vs. ZW samples, respectively. Enrichment analysis of differentially expressed genes (DEGs) resulted in 268, 153, and 42 significantly enriched GO terms, respectively, and a total of 56 significantly enriched KEGG pathways. Additionally, 23 putative sex-related genes, including Gtsf1, IR, HSP21, MRPINK, Mrr, and other potentially promising candidate genes were identified. Moreover, 56,241 simple sequence repeats (SSRs) were identified. Our findings provide a valuable archive for further functional analyses of sex-related genes and future discoveries of underlying molecular mechanisms of gonadal development and sex determination.