Goats have become one of the most adaptive and important livestock species distributed in developing countries in recent years. The Hainan Black goat is a native goat breed of the Hainan region that is generally well-liked by the local population and is thus raised in large numbers. However, the genomic diversity and selective signals of the Hainan Black goat have not been clearly elucidated yet. Therefore, in this study, we performed whole-genome resequencing of 16 Hainan Black goats and compared the results with those of 71 goats of 6 other breeds from different geographic regions. Principal component analysis (PCA) and phylogenetic analysis identified seven lineages for all goats. Hainan Black goats showed the most similarity with Leizhou goats and the least similarity with Boer goats. Selective sweep analysis identified candidate genes associated with various functions, including immune resistance to disease (TNFAIP2 (TNF alpha induced protein 2) and EXOC3L4 (exocyst complex component 3 like 4)), melanin biosynthetic process (CDH15 (cadherin 15), ASIP (agouti signaling protein), and PARD3 (par-3 family cell polarity regulator)), and light sensitivity (CNGB3 (cyclic nucleotide gated channel subunit beta 3) and CNBD1 (cyclic nucleotide binding domain containing 1)), underlying strong selection signatures in Hainan Black goats. The melanin biosynthetic process, circadian entrainment, regulation of cyclic adenosine 3,5-monophosphate (cAMP)-mediated signaling, and the Rap-1 signaling pathway were significantly enriched in Hainan Black and Alashan Cashmere goats. This result may be important for understanding each trait. Selection signature analysis revealed candidate single nucleotide polymorphisms (SNPs) and genes correlated with the traits of Hainan Black goats. Collectively, our results provide valuable insights into the genetic basis of specific traits correlated with the Hainan island climate, artificial selection in certain local goat breeds, and the importance of protecting breed resources.

Background: China has thousands years of goat breeding and abundant goat genetic resources. Additionally, the Hainan black goat is one of the high-quality local goat breeds in China. In order to conserve the germplasm resources of the Hainan black goat, facilitate its genetic improvement and further protect the genetic diversity of goats, it is urgent to develop a single nucleotide polymorphism (SNP) chip for Hainan black goat. Results: In this study, we aimed to design a 10K liquid chip for Hainan black goat based on genotyping by pinpoint sequencing of liquid captured targets (cGPS). A total of 45,588 candidate SNP sites were obtained, 10,677 of which representative SNP sites were selected to design probes, which finally covered 9,993 intervals and formed a 10K cGPS liquid chip for Hainan black goat. To verify the 10K cGPS liquid chip, some southern Chinese goat breeds and a sheep breed with similar phenotype to the Hainan black goat were selected. A total of 104 samples were used to verify the clustering ability of the 10K cGPS liquid chip for Hainan black goat. The results showed that the detection rate of sites was 97.34% -99.93%. 84.5% of SNP sites were polymorphic. The heterozygosity rate was 3.08%-36.80%. The depth of more than 99.4% sites was above 10X. The repetition rate was 99.66%-99.82%. The average consistency between cGPS liquid chip results and resequencing results was 85.58%. In addition, the phylogenetic tree clustering analysis verified that the SNP sites on the chip had better clustering ability. Conclusion: These results indicate that we have successfully realized the development and verification of the 10K cGPS liquid chip for Hainan black goat, which provides a useful tool for the genome analysis of Hainan black goat. Moreover, the 10K cGPS liquid chip is conducive to the research and protection of Hainan black goat germplasm resources and lays a solid foundation for its subsequent breeding work.

Pasteurella multocida can cause goat hemorrhagic sepsis and endemic pneumonia. Respiratory epithelial cells are the first line of defense in the lungs during P. multocida infection. These cells act as a mechanical barrier and activate immune response to protect against invading pathogenic microorganisms. Upon infection, P. multocida adheres to the cells and causes changes in cell morphology and transcriptome. ATAC-seq was conducted to determine the changes in the chromatin open region of P. multocida-infected goat bronchial epithelial cells based on transcriptional regulation. A total of 13,079 and 28,722 peaks were identified in the control (CK) and treatment (T) groups (P. multocida infection group), respectively. The peaks significantly increased after P. multocida infection. The specific peaks for the CK and T groups were annotated to 545 and 6632 genes, respectively. KEGG pathway enrichment analysis revealed that the specific peak-related genes in the T group were enriched in immune reaction-related pathways, such as Fc gamma R-mediated phagocytosis, MAPK signaling pathway, bacterial invasion of epithelial cells, endocytosis, and autophagy pathways. Other cellular component pathways were also enriched, including the regulation of actin cytoskeleton, adherent junction, tight junction, and focal adhesion. The differential peaks between the two groups were subsequently analyzed. Compared to those in the CK group, 863 and 11 peaks were upregulated and downregulated, respectively, after the P. multocida infection. Fifty-six known transcription factor motifs were revealed in upregulated peaks in the P. multocida-infected group. By integrating ATAC-seq and RNA-seq, some candidate genes (SETBP1, RASGEF1B, CREB5, IRF5, TNF, CD70) that might be involved in the goat bronchial epithelial cell immune reaction to P. multocida infection were identified. Overall, P. multocida infection changed the structure of the cell and caused chromatin open regions to be upregulated. In addition, P. multocida infection actively mobilized the host immune response with the inflammatory phenotype. The findings provide valuable information for understanding the regulatory mechanisms of P. multocida-infected goat bronchial epithelial cells.

Buffalo meat is gaining popularity for its nutritional properties, such as its low fat and cholesterol content. However, it is often unsatisfactory to consumers due to its dark color and low tenderness. There is currently limited research on the regulatory mechanisms of buffalo meat quality. Xinglong buffalo are raised in the tropical Hainan region and are undergoing genetic improvement from draught to meat production. For the first time, we evaluated the meat quality traits of Xinglong buffalo using the longissimus dorsi muscle and compared them to Hainan cattle. Furthermore, we utilized a multi-omics approach combining transcriptomics and metabolomics to explore the underlying molecular mechanism regulating meat quality traits. We found that the Xinglong buffalo had significantly higher meat color redness but lower amino acid content and higher shear force compared to Hainan cattle. Differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were identified, with them being significantly enriched in nicotinic acid and nicotinamide metabolic and glycine, serine, and threonine metabolic pathways. The correlation analysis revealed that those genes and metabolites (such as: GAMT, GCSH, PNP, L-aspartic acid, NADP+, and glutathione) are significantly associated with meat color, tenderness, and amino acid content, indicating their potential as candidate genes and biological indicators associated with meat quality. This study contributes to the breed genetic improvement and enhancement of buffalo meat quality.

Pasteurella multocida (P. multocida) is a zoonotic bacterium that predominantly colonizes the respiratory tract and lungs of a variety of farmed and wild animals, and causes severe respiratory disease. To investigate the characteristics of the host immune response induced by P. multocida strains of serotype A and D, high-throughput mRNA-Seq and miRNA-Seq were performed to analyze the changes in goat bronchial epithelial cells stimulated by these two serotypes of P. multocida for 4 h. Quantitative RT-PCR was used to validate the randomly selected genes and miRNAs. The results revealed 204 and 117 differentially expressed mRNAs (|log₂(Fold-change)| ≥ 1, p-value < 0.05) in the P. multocida serotype A and D stimulated groups, respectively. Meanwhile, the number of differentially expressed miRNAs (|log₂(Fold-change)| > 0.1, p-value < 0.05) were 269 and 290, respectively. GO and KEGG enrichment analyses revealed 13 GO terms (p-value < 0.05) and four KEGG pathways (p-value < 0.05) associated with immunity. In the serotype A-stimulated group, the immune-related pathways were the GABAergic synapse and Toll-like receptor signaling pathways, while in the serotype D-stimulated group, the immune-related pathways were the phagosome and B cell receptor signaling pathways. Based on the predicted results of TargetScan and miRanda, the differentially expressed mRNA–miRNA network of immune-related GO terms and KEGG pathways was constructed. According to the cell morphological changes and the significant immune-related KEGG pathways, it was speculated that the P. multocida serotype D strain-stimulated goat bronchial epithelial cells may induce a cellular immune response earlier than serotype A-stimulated cells. Our study provides valuable insight into the host immune response mechanism induced by P. multocida strains of serotype A and D.

Melioidosis is a zoonotic disease that is caused by Burkholderia pseudomallei, which is a serious public health and safety risk. In order to explore the global animal melioidosis risk distribution and its dynamic response to future climate scenarios, we collected global data about reported animal incidence sites. Data regarding the density of Burkholderia pseudomallei in the environment were created by collecting and sorting information about the Burkholderia pseudomallei occurrence sites in contaminated air, soil, and water. Combined with bioclimatic variables, the maximum entropy (MaxEnt) niche was modeled for global animal melioidosis. The findings indicate that under current bioclimatic conditions, global animal melioidosis risk regions are concentrated between 30° S and 30° N, with high-risk areas being mainly in Central America, the northern part of South America, and eastern and southern India, among others. Most countries will expand their risk regions under future climatic scenarios. Melioidosis risk expanding towards higher northern latitudes has led to new epidemic areas. In addition, the melioidosis risk area will contract in some areas. Therefore, we have provided a basis for global melioidosis surveillance and propose feasible strategies for prevention and control in high-risk regions, which will help countries to carry out targeted surveillance and prevention to reduce risks and losses.