SNV discovery and functional candidate gene identification for milk composition based on whole genome resequencing of Holstein bulls with extremely high and low breeding values

Exploration of the Polymorphism Distribution of Bovine HMGA2 Gene in Worldwide Breeds and Its Associations with Ovarian Traits

The high-mobility group AT-hook 2(HMGA2) gene has been widely studied in the context of cancer and animal growth. However, recently, several studies have uncovered its critical role in cell proliferation. A genome-wide association study (GWAS) further suggests that the HMGA2 gene is a candidate gene in fertility, indicating its connection not only to growth traits but also to reproduction, specifically ovarian traits. Thus, this study aimed to analyze the distribution of the HMGA2 gene in 54 bovine breeds worldwide, identify important short fragment variants (indels), and investigate the relationship between HMGA2 and ovarian development. The dataset included genotypic information from a bovine population of 634 individuals (n = 634). After genotyping and analyzing four selected loci, we found that one out of four loci, rs133750033 (P4-D22-bp), was polymorphic. Our results also reveal that this indel of HMGA2 is significantly associated with certain ovarian traits (p < 0.05). Specifically, it has connection with ovarian length (p = 0.004) and ovarian height (p = 0.026) during diestrus. Additionally, we discovered a higher expression of the HMGA2 gene in Asian cattle breeds. In summary, this study suggests that HMGA2 has the potential to serve as an animal fertility testing marker gene. Moreover, these findings contribute to a more promising outlook for the bovine industry.

Detection distribution of CNVs of SNX29 in three goat breeds and their associations with growth traits

As a member of the SNX family, the goat sorting nexin 29 (SNX29) is initially identified as a myogenesis gene. Therefore, this study aimed to examine the polymorphism in the SNX29 gene and its association with growth traits. In this study, we used an online platform to predict the structures of the SNX29 protein and used quantitative real-time PCR to detect potential copy number variation (CNV) in Shaanbei white cashmere (SBWC) goats (n = 541), Guizhou black (GB) goats (n = 48), and Nubian (NB) goats (n = 39). The results showed that goat SNX29 protein belonged to non-secretory protein. Then, five CNVs were detected, and their association with growth traits was analyzed. In SBWC goats, CNV1, CNV3, CNV4, and CNV5 were associated with chest width and body length (P < 0.05). Among them, the CNV1 individuals with gain and loss genotypes were superior to those individuals with a median genotype, but CNV4 and CNV5 of individuals with the median genotype were superior to those with the loss and gain genotypes. In addition, individuals with the gain genotype had superior growth traits in CNV3. In brief, this study suggests that the CNV of SNX29 can be used as a molecular marker in goat breeding.

Assessing Bos taurus introgression in the UOA Bos indicus assembly

Background

Reference genomes are essential in the analysis of genomic data. As the cost of sequencing decreases, multiple reference genomes are being produced within species to alleviate problems such as low mapping accuracy and reference allele bias in variant calling that can be associated with the alignment of divergent samples to a single reference individual. The latest reference sequence adopted by the scientific community for the analysis of cattle data is ARS_UCD1.2, built from the DNA of a Hereford cow (Bos taurus taurus—B. taurus). A complementary genome assembly, UOA_Brahman_1, was recently built to represent the other cattle subspecies (Bos taurus indicus—B. indicus) from a Brahman cow haplotype to further support analysis of B. indicus data. In this study, we aligned the sequence data of 15 B. taurus and B. indicus breeds to each of these references.

Results

The alignment of B. taurus individuals against UOA_Brahman_1 detected up to five million more single-nucleotide variants (SNVs) compared to that against ARS_UCD1.2. Similarly, the alignment of B. indicus individuals against ARS_UCD1.2 resulted in one and a half million more SNVs than that against UOA_Brahman_1. The number of SNVs with nearly fixed alternative alleles also increased in the alignments with cross-subspecies. Interestingly, the alignment of B. taurus cattle against UOA_Brahman_1 revealed regions with a smaller than expected number of counts of SNVs with nearly fixed alternative alleles. Since B. taurus introgression represents on average 10% of the genome of Brahman cattle, we suggest that these regions comprise taurine DNA as opposed to indicine DNA in the UOA_Brahman_1 reference genome. Principal component and admixture analyses using genotypes inferred from this region support these taurine-introgressed loci. Overall, the flagged taurine segments represent 13.7% of the UOA_Brahman_1 assembly. The genes located within these segments were previously reported to be under positive selection in Brahman cattle, and include functional candidate genes implicated in feed efficiency, development and immunity.

Conclusions

We report a list of taurine segments that are in the UOA_Brahman_1 assembly, which will be useful for the interpretation of interesting genomic features (e.g., signatures of selection, runs of homozygosity, increased mutation rate, etc.) that could appear in future re-sequencing analysis of indicine cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00688-1.

Genome-wide assessment of diversity and differentiation between original and modern Brown cattle populations

Identifying genomic regions involved in the differences between breeds can provide information on genes that are under the influence of both artificial and natural selection. The aim of this study was to assess the genetic diversity and differentiation among four different Brown cattle populations (two original vs. two modern populations) and to characterize the distribution of runs of homozygosity (ROH) islands using the Illumina Bovine SNP50 BeadChip genotyping data. After quality control, 34 735 SNPs and 106 animals were retained for the analyses. Larger heterogeneity was highlighted for the original populations. Patterns of genetic differentiation, multidimensional scaling, and the neighboring joining tree distinguished the modern from the original populations. The F_ST -outlier identified several genes putatively involved in the genetic differentiation between the two groups, such as stature and growth, behavior, and adaptability to local environments. The ROH islands within both the original and the modern populations overlapped with QTL associated with relevant traits. In modern Brown (Brown Swiss and Italian Brown), ROH islands harbored candidate genes associated with milk production traits, in evident agreement with the artificial selection conducted to improve this trait in these populations. In original Brown (Original Braunvieh and Braunvieh), we identified candidate genes related with fat deposition, confirming that breeding strategies for the original Brown populations aimed to produce dual-purpose animals. Our study highlighted the presence of several genomic regions that vary between Brown populations, in line with their different breeding histories.

Polymorphic variants of bovine ADCY5 gene identified in GWAS analysis were significantly associated with ovarian morphological related traits

The reproductive performance (e.g. fertility) of dairy cows, which declined over past few decades due to the intense and intensive selection, needs to be improved. Previous genome-wide association study (GWAS) of female Holstein screened the Adenylate cyclase 5 (ADCY5) as the candidate gene for cow fertility. As a member of the adenylyl cyclases family, adenylate cyclase 5 (ADCY5) is famous for regulating extrapyramidal motor system related various neuropsychiatric diseases, and its genetic variant is reported to associate with lower birth and placenta weight which leads to asymmetric fetal growth restriction. It was hypothesized that ADCY5 may affect the fertility of cows by regulating the processes of ovarian development. Herein, genomic DNA from 768 ovaries samples of healthy unrelated Holstein cow were used to screen potential insertion/deletion (indel) mutations using eight pairs of primers, and we found three novel polymorphic indel variants, namely, rs385624978 (P3-D_11-bp), rs433028962 (P5-I_19-bp) and rs382393457 (P8-D_19-bp). The minor allelic frequencies (MAF) of P3-D_11-bp, P5-I_19-bp and P8-D_19-bp loci were 0.188, 0.365 and 0.06, respectively, and there were 7 different haplotypes. Additionally, linkage disequilibrium analysis demonstrated no linkage among them. Importantly, P3-D_11-bp locus was significantly related to both ovarian width (P = 1.0E-6) and corpus luteum diameter (P = 0.015); P5-I_19-bp locus had a significant relation with corpus albicans diameter (P = 0.030) and ovaries with mutational homozygous genotype produced a superior corpus albicans diameter than those with other genotypes. Briefly, three novel indel mutations of bovine ADCY5 gene were identified and two of them were uncovered to be significantly correlated with ovarian phenotypic traits or corpus luteum or albicans traits. These findings contributed to the application of molecular marker-assisted selection (MAS) in improving female fertility in cattle, which could accelerate the development of the cattle industry.

The domesticated buffalo - An emerging model for experimental and therapeutic use of extraembryonic tissues

Large animals play important roles as model animals for biomedical sciences and translational research. The water buffalo (Bubalus bubalis) is an economically important, multipurpose livestock species. Important assisted reproduction techniques, such as in vitro fertilization, cryo-conservation of sperm and embryos, embryo transfer, somatic cell nuclear transfer, genetic engineering, and genome editing have been successfully applied to buffaloes. Recently, detailed whole genome data and transcriptome maps have been generated. In addition, rapid progress has been made in stem cell biology of the buffalo. Apart from embryonic stem cells, bubaline extra-embryonic stem cells have gained particular interest. The multipotency of non-embryonic stem cells has been revealed, and their utility in basic and applied research is currently investigated. In particular, success achieved in bubaline extra-embryonic stem cells may have important roles in experimental biology and therapeutic regenerative medicine. Progress in other farm animals in assisted reproduction techniques, stem cell biology and genetic engineering, which could be of importance for buffalo, will also be briefly summarized.