Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

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.

The aim of the present study was to assess the association of single nucleotide polymorphisms (SNPs) of Calpain (CAPN) gene with birth weight (BW), final weight (FW) and average daily gain (ADG) in three Egyptian sheep breeds: Barki, Rahmani and Ossimi. Blood samples were collected from 108 animals representing the three breeds. DNA was isolated using salting out procedure and then the quality and quantity of DNA extracted were measured. A 190 bp of CAPN was amplified by PCR using specific primers. The allele and genotype frequencies for all the identified SNPs were calculated. The PCR products corresponding to each genotype were sequenced to identify SNPs associated with the traits in question. Two SNPs (C→T) were detected in the nucleotides 44 and 154. For each SNP, the two mentioned alleles were named C and T, respectively. The sequenced CAPN segments were subjected to nucleotide blast at NCBI, which revealed 99% identity with that reported for sheep in Genbank. The TT was the least common genotype, whereas frequencies of CT and CC genotypes were fluctuated in the three sheep breeds under study. Animal carrier TT genotype had higher BW, FW and ADG than those with CT genotype, while the lowest values were associated with CC genotype. For the three traits under study, Rahmani had the highest estimates followed by Ossimi and Barki. Males exhibited heavier BW and FW as well as higher ADG compared with females. The results generated provide preliminary indication of the functional diversity present in Barki, Rahmani and Ossimi sheep and the possibility of using this polymorphism in Egyptian sheep genetic improvement.

The water buffalo is an important dual-purpose livestock that is widespread throughout central and southern China. However, there has been no characterization of the population genetics of Chinese buffalo. Using an Axiom buffalo genotyping array (Thermo Fisher Scientific, Wilmington, DE), we analyzed the genetic diversity, linkage disequilibrium pattern, and signature of selection in 176 Chinese buffaloes from 13 breeds. A total of 35,547 SNP passed quality control and were used for further analyses. Population genetic analysis revealed a clear separation between swamp and river types. Ten Chinese indigenous breeds were clustered into the swamp group, the Murrah and Nili-Ravi breeds were clustered into the river group, and the crossbred breed was closer to the river group. Genetic diversity analysis showed that the swamp group had a lower average expected heterozygosity. Linkage disequilibrium decay distance was much shorter in the swamp group compared with the river group, with an average square of correlation coefficient value of 0.2 of approximately 50 kb. Analysis of runs of homozygosity indicated extensive remote and recent inbreeding within swamp and river groups, respectively. Moreover, one genomic region under selection was detected between the river and swamp groups. Our findings contribute to our understanding of the characterization of population genetics in Chinese buffaloes, which in turn may be used in buffalo breeding programs.

Aim: The present study was performed to assess the association of single-nucleotide polymorphisms (SNPs) in the fatty acid-binding protein 4 (FABP4) gene with birth weight (BW), final weight (FW), and average daily gain (ADG) in three Egyptian sheep breeds. Materials and Methods: Genomic DNA was extracted from the blood samples of 50 male and female individuals representing Ossimi, Rahmani, and Barki sheep breeds. A 407 bp nucleotide (nt) segment from the first intron of FABP4 was amplified by polymerase chain reaction, sequenced, and analyzed in the different samples. Results: Sequence analysis of the determined segment (407 bp) revealed four SNPs (all transition types) at nt position 372 (CP011894.1:g.57605471) A>G, nt position 211 (CP011894.1:g.57605632) A>G, nt position 143 (CP011894.1:g.57605700) T>C, and nt position 111 (CP011894.1:g.57605732) T>C. The allelic and genotypic frequencies for the identified SNPs in the sheep breeds were calculated. At nt positions 372 and 211, two alleles were identified (A and G). Only two genotypes were present at nt position 372 (AA and AG), while three genotypes were present at nt position 211 (AA, AG, and GG). Two alleles (T and C) and three identified genotypes (TT, TC, and CC) were detected at nt positions 143 and 111. Analysis of the results revealed that AA genotype at nt position 372 is associated with higher estimates for BW, FW, and ADG when compared to all the other genotypes. Very high correlation coefficients were found between the genotypes 143-TT and 111-TT and also between 143-TC and 111-TC. The genotypes 372-AG, 211-GG, 211-AA, 143-TT, 143-CC, 111-TT, 111-TC, and 111-CC were associated with negative effects on BW, FW, and ADG. Conclusion: The detection of four SNPs in a partial sequence of the Egyptian ovine FABP4 gene intron 1 reflected that this gene harbors substantial diversity. In addition, a novel SNP at nt position 372 (CP011894.1:g.57605471) A>G was associated with higher estimates for BW, FW, and ADG.

The organic anion transporter (OAT) family is the subfamily of the solute carrier (SLC) superfamily, which plays a vital role in regulating essential nutrients in milk. However, little is known about the members’ identification, evolutionary basis, and function characteristics of OAT genes associated with milk performance in buffalo. Comparative genomic analyses were performed to identify the potential role of buffalo OAT genes in milk performance in this study. The results showed that a total of 10 and 7 OAT genes were identified in river buffalo and swamp buffalo, respectively. These sequences clustered into three groups based on their phylogenetic relationship and had similar motif patterns and gene structures in the same groups. Moreover, the river-specific expansions and homologous loss of OAT genes occurred in the two buffalo subspecies during the evolutionary process. Notably, the duplicated SLCO3A1 gene specific to river buffalo showed higher expression level in mammary gland tissue than that of swamp buffalo. These findings highlight some promising candidate genes that could be potentially utilized to accelerate the genetic progress in buffalo breeding programs. However, the identified candidate genes require further validation in a larger cohort for use in the genomic selection of buffalo for milk production.