Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing.

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

Global dispersal and adaptive evolution of domestic cattle: a genomic perspective.

Domestic cattle have spread across the globe and inhabit variable and unpredictable environments. They have been exposed to a plethora of selective pressures and have adapted to a variety of local ecological and management conditions, including UV exposure, diseases, and stall-feeding systems. These selective pressures have resulted in unique and important phenotypic and genetic differences among modern cattle breeds/populations. Ongoing efforts to sequence the genomes of local and commercial cattle breeds/populations, along with the growing availability of ancient bovid DNA data, have significantly advanced our understanding of the genomic architecture, recent evolution of complex traits, common diseases, and local adaptation in cattle. Here, we review the origin and spread of domestic cattle and illustrate the environmental adaptations of local cattle breeds/populations.

Whole genome sequencing of simmental cattle for SNP and CNV discovery.

BACKGROUD: The single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity. RESULTS: In this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 - 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility. CONCLUSIONS: In conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

Genome-wide association study reveals that the IBSP locus affects ear size in cattle.

Ear size is a classical model for hot climate adaptation following the evolution, but the genetic basis of the traits associated with ear size remains to be elucidated. Here, we performed a genome-wide association study on 158 cattle to explain the genetic mechanism of ear size. One region on BTA6 between 36.79 and 38.80 Mb included 50 suggestive SNPs and 4 significant SNPs that were significantly associated with ear size. The most significant locus (P = 1.30 × 10-8) was a missense mutation (T250I) on the seventh exon of integrin-binding sialoprotein (IBSP), which had an allele substitution effect of 23.46 cm2 for ear size. Furthermore, this mutation will cause changes in the three-dimensional structure of the protein. To further identify genes underlying this typical feature, we performed a genome scan among nine cattle breeds with different ear sizes by using SweeD. Results suggested that IBSP was under positive selection among four breeds with relatively large ear sizes. The expression levels of IBSP in ear tissues of large- and small-ear cattle were significantly different. A haplotype diversity survey of this missense mutation in worldwide cattle breeds strongly implied that the origin of this missense mutation event was Bos taurus. These findings have important theoretical importance for the exploration of major genes associated with ear size and provide important molecular markers for the identification of cattle germplasm resources.

Goat CLSTN2 gene: tissue expression profile, genetic variation, and its associations with litter size.

Calsyntenin-2 (CLSTN2) is involved in cell proliferation, differentiation, cell death, tumorigenesis, and follicular expression. Although CLSTN2 has been identified as a potential candidate gene for sheep prolificacy, no studies have been done on its effect on goat prolificacy. The purpose of this study was to identify mRNA expression and genetic variation within goat CLSTN2, and its association with prolificacy. Herein, we uncovered significant differences in mRNA levels of the CLSTN2 gene in different tissues in female goats (p < 0.01), including ovary tissue. Nine putative indels were designed to investigate their correlation to litter size, but only one 16-bp deletion was discovered in female Shaanbei white cashmere goats (n = 902). We discovered that a 16-bp deletion within the CLSTN2 gene was significantly correlated with first-born litter size (p = 0.0001). As shown by the chi-squared test, the genotypic II of single-lambs and multi-lambs was dramatically higher than with genotype ID (p = 0.005). Our findings suggest that indel within the CLSTN2 gene is a candidate gene affecting prolificacy in goats and may be used for Marker Assisted Selection (MAS) in goats.

Maternal genetic diversity, differentiation and phylogeny of three white yak breeds/populations in China.

White yak is a unique economic livestock animal on the Qilian Mountains of Qinghai-Tibet Plateau (QTP). Here, 369 mitochondrial D-loop sequences from three Chinese white yak breeds/populations (Tianzhu, Menyuan and Huzhu) were comprehensively analyzed to indicate their maternal genetic diversity, differentiation and phylogenetic relationship. Our results showed that the haplotype diversity (Hd) was found to be highest in Tianzhu white yak (Hd ± SD = 0.9501 ± 0.0058), while the lowest was recorded in Huzhu white yak (Hd ± SD = 0.7178 ± 0.0474). Totally, the haplotype and nucleotide diversities of white yak were 0.9407 ± 0.0069 and 0.0187 ± 0.0094, respectively, indicating an abundant maternal genetic diversity in white yak. Estimates of FST (Fixation Index) values showed a high genetic differentiation between Tianzhu and Menyuan populations (FST = 0.2928, p < 0.05) as well as that between Tianzhu and Huzhu populations (FST = 0.2721, p < 0.05), but a moderate genetic differentiation between Menyuan and Huzhu populations (FST = 0.1352, p > 0.05) was observed. Cluster analysis based on FST values among populations indicated that the genetic relationship between Menyuan and Huzhu white yak was closer, but they had a far genetic relationship with Tianzhu white yak. Maternal phylogenetic analysis revealed that white yak represented two maternal lineages (I and II), suggesting two maternal origins. In addition, it is notable that taurine mtDNA haplotypes were detected in Tianzhu and Huzhu white yak populations, indicating taurine genetic introgression to some extent. Our study would provide useful information for the conservation and utilization of white yak.

Muscle transcriptome analysis reveal candidate genes and pathways related to fat and lipid metabolism in Yunling cattle.

Yunling cattle (YL) is a recently developed beef breed harboring a quarter of Yunnan ancestral cattle genome, spanning over past 30 years. Compared with Diqing cattle (DQ), a Yunnan native cattle breed, YL presents various advantages, including rapid growth and exquisite meat quality. However, the molecular mechanisms underlying these phenotypic differences are not clearly understood. To further identify the candidate genes responsible for the quality of the meat in the muscle, longissimus dorsi (LD) muscle was used for RNA-Seq analysis. A total of 508 differentially expressed genes (DEGs) were identified in YL (adjusted p-value <0.01 and log2FoldChange >1), of which 243 were up-regulated and 265 were down-regulated. Functional association analysis showed that the identified DEGs mainly enriched the lipid and fat metabolism pathways. Moreover, it was also observed that several fat-related genes were differentially expressed in both cattle breeds, including three up-regulated genes (MOGAT1, ACSM3, PLPP2) and two down-regulated genes (ADIG, GPAT3). In addition, alternative splice analysis was also performed revealing an important 9-11 exon skipping variation of GPAM gene (crucial for beef marbling) in YL, which is three times higher than that in DQ, suggesting that this variation might have played the central role in the 'snow beef' effect in YL. We believe that our results will help in understanding the mechanism of muscle development and promote the further breeding programs in YL cattle.

A novel SNP of TECPR2 gene associated with heat tolerance in Chinese cattle.

Heat stress affects the animal production and causes serious economic losses to the husbandry. Tectonin beta-propeller repeat containing 2 (TECPR2) gene plays an important role in autophagy which may affect the temperature sensation in animals. A missense mutation (XM_024981840.1:c.3989 G > A p.Arg1330His) of the transcripts X4 in the bovine TECPR2 gene was identified. In this study, the c.3989 G > A variant in TECPR2 gene was genotyped in a total of 25 cattle breeds (520 individuals). Our results indicated that the frequency of A allele showed a decreasing pattern from southern cattle to northern cattle, while the frequency of G allele showed the opposite pattern, which was consistent with the climate distribution of China. Compared with the GG genotype, southern cattle carried more the AA and AG genotypes. Furthermore, the association results carried out that the frequencies of genotypes (GG, AG, AA) and the value of climate parameters (mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) were significantly correlated (p < 0.01). Hence, we speculated that the c.3989 G > A variant of TECPR2 gene was associated with the heat tolerance trait in Chinese cattle and the locus may be considered as a molecular marker for Chinese cattle breeding.

Whole-Genome Sequencing of Endangered Dengchuan Cattle Reveals Its Genomic Diversity and Selection Signatures.

Dengchuan cattle are the only dairy yellow cattle and endangered cattle among Yunnan native cattle breeds. However, its genetic background remains unclear. Here, we performed whole-genome sequencing of ten Dengchuan cattle. Integrating our data with the publicly available data, Dengchuan cattle were observed to be highly interbred than other cattle in the dataset. Furthermore, the positive selective signals were mainly manifested in candidate genes and pathways related to milk production, disease resistance, growth and development, and heat tolerance. Notably, five genes (KRT39, PGR, KRT40, ESR2, and PRKACB) were significantly enriched in the estrogen signaling pathway. Moreover, the missense mutation in the PGR gene (c.190T > C, p.Ser64Pro) showed a homozygous mutation pattern with higher frequency (83.3%) in Dengchuan cattle. In addition, a large number of strong candidate regions matched genes and QTLs related to milk yield and composition. Our research provides a theoretical basis for analyzing the genetic mechanism underlying Dengchuan cattle with excellent lactation and adaptability, crude feed tolerance, good immune performance, and small body size and also laid a foundation for genetic breeding research of Dengchuan cattle in the future.

Whole-Genome Analyses Reveal Genomic Characteristics and Selection Signatures of Lincang Humped Cattle at the China-Myanmar Border.

The location on the Yunnan border with Myanmar and its unique cultural landscape has shaped Lincang humped cattle over time. In the current study, we investigated the genetic characteristics of 22 Lincang humped cattle using whole-genome resequencing data. We found that Lincang humped cattle derived from both Indian indicine and Chinese indicine cattle depicted higher levels of genomic diversity. Based on genome-wide scans, candidate genomic regions were identified that were potentially involved in local thermal and humid environmental adaptions, including genes associated with the body size (TCF12, SENP2, KIF1C, and PFN1), immunity (LIPH, IRAK3, GZMM, and ELANE), and heat tolerance (MED16, DNAJC8, HSPA4, FILIP1L, HELB, BCL2L1, and TPX2). Missense mutations were detected in candidate genes IRAK3, HSPA4, and HELB. Interestingly, eight missense mutations observed in the HELB gene were specific to the indicine cattle pedigree. These mutations may reveal differences between indicine and taurine cattle adapted to variable climatic conditions. Our research provides new insights into the genetic characteristics of Lincang humped cattle representing Lincang and Pu'er areas as an important channel for the migration of Indian indicine from domestication centers toward southwestern China.

Rich maternal and paternal genetic diversity and divergent lineage composition in wild yak (Bos mutus).

Wild yak (Bos mutus) is a vulnerable bovine species on the Qinghai-Tibetan Plateau (QTP). So far, most studies on molecular genetic diversity of wild yak have focused on autosomal and mtDNA variations based on small number of samples. In this study, we analyzed 84 D-loop and 24 whole mitogenome sequences of wild yak to further comprehensively explore its maternal genetic diversity and lineage composition. Meanwhile, using six yak Y-specific polymorphic markers (i.e., SRY4, USP9Y, UTY19, AMELY3, OFD1Y10 and INRA189), we assessed the paternal genetic diversity and lineage composition based on eight wild yak. Our results showed that wild yak exhibited abundant maternal genetic diversity with haplotype diversities of 0.9621 ± 0.0078 and 0.9928 ± 0.0144 in the D-loop and whole mitogenome sequences, respectively. Maternal phylogenetic analysis of wild yak uncovered three defined lineages (mt-I, mt-II and mt-III). Similarly, profuse paternal genetic diversity was observed in wild yak with Y-haplotype diversity (Hd) at 0.8214 ± 0.1007. Two Y-haplogroups (Y1 and Y2) with four Y-haplotypes (yH1-yH4) were identified in paternal phylogenetic analysis, indicating wild yak to be of two paternal lineages. This study of genetic diversity and lineage composition of wild yak would provide useful information for the genetic resource conservation and utilization of this vulnerable wild species.

Genome-wide association study identifies quantitative trait loci affecting cattle temperament.

Cattle temperament is an interesting trait due to its correlation with production efficiency, labor safety, and animal welfare. To date, however, its genetic basis is not clearly understood. Here, we performed a genome-wide association study for a series of temperament traits in cattle, assessed with via open field and novel object tests, using autosomal single nucleotide polymorphisms (SNPs) derived from the whole-genome sequence. We identified 37 and 29 genome-wide significant loci in the open field and novel object tests, respectively. Gene set analysis revealed the most significant pathway was the neuroactive ligand-receptor interaction pathway, which may be essential for emotional control in cattle. Analysis of the expression levels of 18 tissue-specific genes based on transcriptomic data showed enrichment in the brain, with some candidate genes involved in psychiatric and neurodegenerative diseases in humans. Based on principal component analysis, the first principal component explained the largest variance in the open field and novel object test data, and the most significant loci were assigned to SORCS3 and SESTD1, respectively. Our findings should help facilitate cattle breeding for sound temperament by pyramiding favorable alleles to further improve cattle production.

Rich maternal and paternal genetic diversity and divergent lineage composition in wild yak (Bos mutus).

Wild yak (Bos mutus) is a vulnerable bovine species on the Qinghai-Tibetan Plateau. So far, most studies on the molecular genetic diversity of wild yak have focused on autosomal and mtDNA variations based on the small number of samples. In this study, we analyzed 84 D-loop and 24 whole mitogenome sequences of wild yak to further comprehensively explore its maternal genetic diversity and lineage composition. Meanwhile, using six yak Y-specific polymorphic markers (i.e., SRY4, USP9Y, UTY19, AMELY3, OFD1Y10 and INRA189), we assessed the paternal genetic diversity and lineage composition based on eight wild yak. Our results showed that wild yak exhibited abundant maternal genetic diversity with haplotype diversities of 0.9621 ± 0.0078 and 0.9928 ± 0.0144 in the D-loop and whole mitogenome sequences, respectively. Maternal phylogenetic analysis of wild yak uncovered three defined lineages (mt-I, mt-II and mt-III). Similarly, profuse paternal genetic diversity was observed in wild yak with Y-haplotype diversity at 0.8214 ± 0.1007. Two Y-haplogroups (Y1 and Y2) and four Y-haplotypes (yH1-yH4) were identified in paternal phylogenetic analysis, indicating wild yak to be of two paternal lineages. The present study of genetic diversity and lineage composition of wild yak would provide useful information for the genetic resource conservation and utilization of this vulnerable wild species.

Distribution of the variant at AKIRIN2: c.*188G > A in Chinese cattle.

Marbling score (MS), is an economically important trait in cattle. Previous results showed that a SNP (c.*188G > A) of akirin 2 (AKIRIN2) gene was associated with MS in Japanese Black cattle and Korean cattle. However, the distribution of the genotypic frequency of the single nucleotide polymorphism (SNP) has not been explored in Chinese cattle. In this study, we used polymerase chain reaction (PCR) and DNA sequencing to detect the variation in 1296 individuals from 39 Chinese cattle breeds, one semi-wild bovine species (Dulong) and three introduced breeds (Angus, Holstein and Brahman). Our study found the frequency of the A allele at this locus roughly diminished from north to south in Chinese cattle, and we detected statistically significant differences between Angus and Brahman (p < 0.05), Dulong and another two breeds (Angus and Holstein; p < 0.01) using Chi-Square Independence Test. Our results reflected the variation of AKIRIN2: c.*188G > A in Chinese cattle, which would help us better understand Chinese cattle genetic resources and provide reference for further research.

Characterization of whole mitogenome sequence of the Tongde yak (Bos grunniens).

Tongde County is located in the southeast of Qinghai Province, China, harboring rich yak genetic resources. In the present study, the complete mitochondrial genome (mitogenome) of the Tongde yak (Bos grunniens) was firstly sequenced using Illumina sequencing technique and the corresponding sequence characterization was identified. Our results showed that the mitogenome of Tongde yak is a circular molecule with 16,323 bp length consisting of 37 genes (13 protein-coding genes, 2 rRNA genes, 22 tRNA genes) and a non-coding control region (D-loop), which is consistent with most bovine species. The overall nucleotide composition was found as: A (33.72%), T (27.27%), C (25.80%), and G (13.21%), respectively, yielding a higher AT content (60.99%). The complete mitogenome sequence of Tongde yak would provide useful information for further studies on its genetic resource conservation and molecular breeding programmes in the future.

Antiproliferative activity of 3,5-disubstituted tetrahydro-2H-1,3,5-thiadiazine thione (THTT) derivatives and evaluation as potential prodrug.

Four series of tetrahydro-2H-1,3,5-thiadiazine thione derivatives were screened for their in vitro antiproliferative activities against two human cancerous PC3 and HeLa cell lines. The cytotoxicity of all the compounds (series A-D) was also determined on mammalian mouse fibroblast 3T3 cells. Most of the compounds showed significant anticancer potential against both cancer cell lines within the range of IC50 = 6.4-29.9 and 2.4-23.8 M respectively when compared with standard doxorubicin (IC50 = 0.3 M). All compounds demonstrated a notable selectivity for Hela cells and found either non-toxic or relatively less toxic for 3T3 cell lines model. The structure-activity relationship indicated that antiproliferative activity mainly influenced by the nature and position of substituents at thidiazine nucleus. In general, the presence of aryl groups for example 3,4-(OMe) 2.Bzl and CH(Ph)Me at N-3 position resulted in a significant activity. Under enzymatic hydrolysis, complete conversion (100%) of ester derivative of thiadiazine thione (10a) into its acidic counterpart (7c) was achieved during 20 min which indicated that these types of THTT ester derivatives can be a possible lead for future investigations as prodrug anticancer probes.

An atlas of CNV maps in cattle, goat and sheep.

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

DGAT1 K232A polymorphism is associated with milk production traits in Chinese cattle.

The production traits of cattle, especially milk trait, are of great significance to human life. A quantitative trait loci (QTL) associated with milk fat content was detected in the centromeric region of cattle chromosome 14. This QTL harbors a strong candidate gene called DGAT1 responsible for the milk quality. A non-conservative substitution of lysine by alanine (K232A) was found in DGAT1 gene producing a strong effect on milk composition and yield. The lysine (K allele) is associated with increased milk fat content, while the decreased milk fat content is linked to the alanine (A allele) amino acid. To estimate the frequencies of the DGAT1 K232A polymorphism in Chinese cattle breeds, PCR and DNA sequencing methods were used to investigate the polymorphism of DGAT1 K232A in a total of 682 individuals, including 655 Chinese cattle and 27 Holstein cattle. The results demonstrated that the frequency of K allele gradually elevated from the northern group to the southern group of native Chinese cattle, whereas the frequency of A allele showed a contrary pattern, displaying a significant geographical difference across native Chinese cattle breeds. Our results confirm that the southern cattle group has higher milk fat content than that of the northern group.

Whole genome analyses revealed genomic difference between European taurine and East Asian taurine.

European taurine and East Asian taurine are two main clades in Bos taurus, but their genomic differences are not clearly elucidated. Here, we sequenced 16 Mongolian cattle genomes and compared them to the 92 genomes of 10 representative breeds worldwide. We found the highest LD level in Mishima cattle and the fastest LD decay in European taurine. Phylogenetic analysis revealed that Mongolian, Hanwoo and Mishima cattle were clustered into East Asian taurine. From selective sweep, gene annotation, functional enrichment and differential expression analysis, we identified selective signals including genes and/or pathways related to rapid growth and large body size in European taurine, and superior meat quality in East Asian taurine. Our findings will help us understand the evolutionary history and formation process of the breeds and provide theoretical materials regarding the genetic mechanism underlying breed characteristics and molecular breeding programmes of the taurine clades in the future.

Whole-genome resequencing reveals diversity and selective signals in Longlin goat.

The Longlin goat is one of the most valuable livestock species in Guangxi Autonomous Region of China, but its genomic diversity and selective signals are not clearly elucidated. Here we compared 20 genomes of Longlin goat to 66 genomes of other seven goat breeds worldwide to analyze patterns of Longlin goat genetic variation. We found the lowest linkage disequilibrium at the large distances between SNPs associated with the highest effective population size in the recent generations ago in Longlin goat. The eight goat breeds could be divided into Euro-African and East Asian goat population. Interestingly, like East Asian taurine, the same two migration phases might have occurred in the history of East Asian goat. More importantly, we identified selective signals implicated in immune resistance to disease, especially for skin disease, in Longlin goat. Our findings will not only help understand the evolutionary history and breed characteristic but can provide valuable resources for conservation of germplasm resources and implementation of crossbreeding programs.