Jifeng Zhong

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-21014-9

Current knowledge about the evolutionary history of donkeys is still incomplete due to the lack of archeological and whole-genome diversity data. To fill this gap, we have de novo assembled a chromosome-level reference genome of one male Dezhou donkey and analyzed the genomes of 126 domestic donkeys and seven wild asses. Population genomics analyses indicate that donkeys were domesticated in Africa and conclusively show reduced levels of Y chromosome variability and discordant paternal and maternal histories, possibly reflecting the consequences of reproductive management. We also investigate the genetic basis of coat color. While wild asses show diluted gray pigmentation (Dun phenotype), domestic donkeys display non-diluted black or chestnut coat colors (non-Dun) that were probably established during domestication. Here, we show that the non-Dun phenotype is caused by a 1 bp deletion downstream of the TBX3 gene, which decreases the expression of this gene and its inhibitory effect on pigment deposition.

Peroxisome proliferator-activated receptor alpha (PPARα) regulates responses to chemical or physical stress in part by altering expression of genes involved in proteome maintenance. In this research, polymerase chain reaction (PCR) technique was used to amplify 766 and 589 bp fragments of intron 3 and 7 of PPARα gene in Chinese Holstein (n = 771). Sequencing results showed that three novel single nucleotide polymorphisms (SNPs) were identified at position 44087 (G/A), 65550 (G/A), and 65676(G/A) in the PPARα gene. PCR–restriction fragment length polymorphism technology was used to genotype the three SNPs. Association analysis showed that cows with H1H8 (P < 0.05), H2H8 (P < 0.01), H5H7 (P < 0.05), H5H8 (P < 0.05), and H8H8 (P < 0.05) haplotype combinations had lower potassium content in erythrocytes than those with H2H6 haplotype combination. Cows with H1H8, and H8H8 haplotype combinations had lower decrease rate of milk yield than those with H2H6 and H6H8 haplotype combinations (P < 0.05). Cows with H2H8 and H8H8 haplotype combinations had lower rectal temperature than those with H5H8 and H7H7 haplotype combinations (P < 0.05). In conclusion, H8H8 haplotype combination may be advantageous for heat resistance traits in Chinese Holstein cattle.

Molecular chaperones have been understood to be preferentially transcribed to prevent perturbations in response to various stresses. In this study, three single nucleotide polymorphisms (SNPs), g.324G>C, g.589C>T and g.651C>G in Heat shock factor binding protein 1 (HSBP1) gene were found and genotyped in 930 Chinese Holstein cattle. The results indicated that only g.589C>T polymorphism locus met Hardy–Weinberg equilibrium (P > 0.05). Pair linkage disequilibrium analysis and haplotype construction of HSBP1 gene were performed using SHEsis software. Seven haplotypes were constructed and fourteen haplotype combinations were found. Association analysis showed that H2H2 haplotype combination was advantageous for thermo tolerance breeding in Chinese Holstein. The cows with H2H2 haplotype combination have lower decrease rate of milk yield than those with H2H3 haplotype combination (P < 0.05) and lower potassium content in erythrocytes (PCE) than those with H2H5 (P < 0.05), H4H4 (P < 0.05) and H4H5 (P < 0.01) haplotype combination. The association between SNP and thermo tolerance traits showed that PCE of cows with GG genotype was lower than those with CG genotype at g.651C>G locus (P < 0.01). Pair linkage disequilibrium analysis revealed that the three loci were at a strong disequilibrium state. So we presumed that the effect of H2H2 haplotype combination on thermo tolerance traits major due to the SNP of g.651C>G.

Understanding bovine metabolism and its relationship with milk products is important in cow breeding. In the present work, the metabolic network in the mammary gland tissue of cattle was reconstructed with the available bovine genome information using several public datasets from NCBI, Uniprot, and KEGG. The network consisted of 1,743 metabolites named by KEGG compound numbers as nodes and 657 enzymes that catalyzed the corresponding reactions as edges. The characteristics of the network were analyzed. The top 20 hub metabolites were determined, and the mean path length was identified to be 6.52. Moreover, 11 key enzymes with significant changes in expression under the condition of mastitis were identified and analyzed by integrating the microarray expression data of normal and clinical mastitis. Aside from the GATM gene, 10 downregulated enzymes were detected in bovine with mastitis. In addition, many of the identified enzymes were involved in amino acid metabolisms or had a direct connection to amino acid metabolisms. These results indicate that mastitis could affect the expression of enzymes, which is vital in some amino acid metabolisms, resulting in the reduction of milk proteins. The present work provides information that may improve the understanding on bovine milk production and mastitis.

The signaling lymphocyte-activating molecule family 7 (SLAMF7) proteins serve as adhesion molecules on the surface of a variety of mature hematopoietic cells, and also partially control certain innate and adaptive immune responses. We characterized three novel bovine SLAMF7 splice variants, designated as SLAMF7-AS1, AS2, and AS3. All three novel SLAMF7 isoforms are derived from the complete transcripts (SLAMF7-complete) via alternative splicing (AS). The patterns of the three splice variants are exon skipping and alternative 5′ splice sites. Bovine SLAMF7 transcripts are expressed in mammary tissue, as demonstrated by real-time PCR. The levels of the complete transcript expression in the normal mammary tissues were higher than that in Staphylococcus aureus (Staph. aureus)-induced mastitis mammary tissues. However, it was not significant for the mRNA expression level comparison between these two kinds of mammary. The SLAMF-AS2 isoforms are expressed the lowest levels among the three transcripts in both normal and infected mammary tissues. This study provides clues for a better understanding of bovine SLAMF7 gene function.

Leucine aminopeptidase 3 (LAP3) is an aminopeptidase which catalyses the removal of N-terminal amino acids and is involved in protein maturation and degradation. In this study, we detected the polymorphisms of LAP3 gene by polymerase chain reaction-single-strand conformation polymorphism (PCR-SSCP) and DNA sequencing methods in 916 individuals from three Chinese cattle breeds including Chinese Holstein, Luxi Yellow and Bohai Black. One novel single nucleotide polymorphism (SNP) (g.24564G>A ss196003366) and four previously deposited SNPs in the GenBank database (g.24794T>G, g.24803T>C, g.24846T>C, g.25415T>C) were detected. Three of the SNPs (g.24794T>G, g.24803T>C, g.24846T>C) were firstly found to be linked completely and regarded as a SNP g.24794M>N by PCR-SSCP and DNA sequencing in the tested breeds. The allelic frequencies and genetic indices of the SNPs were different in three Chinese cattle populations. The SNPs and their genetic effects on milk production traits in Chinese Holsteins were evaluated. Least squares analysis showed that cows with genotype MM had higher fat percentage and protein percentage than genotype NN (P < 0.05); and the cows with g.25415T>C-CC genotype had higher protein rate than ones with TT genotype (P < 0.05). In addition, eight haplotypes and 23 combined genotypes were identified based on the nine genotypes and the association between combined genotypes and milk production traits were analyzed. Statistic results showed that the cows with genotype combination MAT/MGC have higher protein and fat rate and lower SCS. Our finding demonstrated that the LAP3 gene possibly contributed to conducting association analysis and can be used as molecular marker in milk production traits and other performance for animal breeding.

Heat stress induces heat shock proteins (HSPs) expression and HSP70 family is one of them that have been reported to involve in cellular protection against heat stress. But whether there is any association of genetic variation in the HSP70A1A gene with thermo tolerance is unknown. PCR-SSCP and DNA sequencing were used to detect possible SNPs in HSP70A1A gene in 890 Chinese Holstein cattle. Three fragments were amplified and five novel mutations were found in HSP70A1A gene in Chinese Holstein cattle. G/A mutation was found at nucleotide 1524 in coding region that resulting two genotypes of AA and AB. T/C mutation was found at nucleotide 3494 in 3′-UTR resulting three genotypes of CC, CD and DD. The other three point mutations, G/C at nucleotide 6400, C/T at nucleotide 6600 and G/A at nucleotide 6601 were also found in 3′-UTR resulting six genotypes of EE, EF, FF, EG, FG and GG. Linkage disequilibrium (LD) analysis showed that the nucleotide 6400, 6600 and 6601 were in strong LD (D > 0.75). Association analysis indicated that AB, DD and FF genotype were the thermo tolerance genotype. SBYE Green I was used to quantify HSP70A1A mRNA expression in different tissues through quantitative real-time PCR assay. The results of the real-time PCR showed that the expression of HSP70A1A mRNA in the heart was significantly higher than that in the other tissues (P < 0.05). We presume that these mutations could be used in marker assisted selection for anti-heat stress cows in our breeding program.

As an essential signaling modulator, Src gene appears to be necessary for increased expression of the prolactin receptor, normal downstream signaling, and alveolar cell organization. In this study, we detected the polymorphism of Src gene by polymerase chain reaction-restriction fragment length polymorphism assay (PCR-RFLP) and DNA sequencing methods in 985 individuals from three Chinese cattle breeds. Three novel single nucleotide polymorphisms (SNPs) (g.14062C>T ss161151834, g.17302G>A ss161151835, g.18107T>C ss161151836) were detected. Least squares analysis showed that cows with g.14062C>T-CC genotypes and g.18107T>C-TT genotypes had the highest protein rate, while the cows with g.17302G>A-GG genotype had higher 305 d milk yield (p < 0.05), fat yield (p < 0.01) and protein yield (p < 0.01) than the ones with genotypes g.17302G>A-GA. These results revealed the statistical significant effects of three SNPs of the Src gene on the milk production traits in Chinese Holstein. In addition, based on the nine genotypes constructed from 27 combined haplotypes, the association analysis between combined haplotypes and milk production traits was carried out. Statistic results showed that the cows with combined haplotypes H2H4(CCGATT) had the highest fat rate and the highest protein rate and the cows with combined haplotypes H1H8(TCGATC) and H3H7(TCGGCC) had greater 305 d mild yield than H1H2(CCAATC)(P < 0.05). Our finding demonstrated that the Src gene possibly contributed to conducting association analysis and can be recognized as genetic marker in milk production traits and other performance for animal breeding and genetics.