ddRAD sequencing based genotyping of six indigenous dairy cattle breeds of India to infer existing genetic diversity and population structure

Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle

The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from -0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075-35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.

Analysing the genetic diversity and population structure of five native Turkish cattle breeds using SNP data

The conservation and sustainable utilization of cattle genetic resources necessitate a comprehensive understanding of their genetic diversity and population structure. This study provides an analysis of five native Turkish cattle breeds: Anatolian Black (ANB), Turkish Grey (TUR), Anatolian Southern Yellow (ASY), East Anatolian Red (EAR), and South Anatolian Red (SAN) using 50 K SNP data. These breeds were compared with three European breeds, Simmental (SIM), Holstein (HOL), and Jersey (JER), and three Asian Zebu breeds: Arabic Zebu (ZAR), Nelore (NEL), and Red Sindhi (RSI). Genetic diversity indices demonstrated moderate heterogeneity among the breeds, with TUR exhibiting the highest observed heterozygosity (Ho = 0.35). Wright's Fst values indicated significant genetic differentiation, particularly between Turkish breeds and both European (Fst = 0.035-0.071) and Asian breeds (Fst = 0.025-0.150). Principal component analysis distinguished the unique genetic profiles of each breed cluster. Admixture analysis revealed degrees of shared genetic ancestry, suggesting historical gene flow between Turkish, European, and Asian breeds. Analysis of molecular variance (AMOVA) attributed approximately 58% of the variation to population differences. Nei's genetic distances highlighted the closer genetic relatedness within Turkish breeds (distance ranges between 0.032 and 0.046) and suggested a more relative affinity of TUR with European breeds. The study's phylogenetic assessments elucidate the nuanced genetic relationships among these breeds, with runs of homozygosity (ROH) analysis indicating patterns of ancestral relatedness and moderate levels of inbreeding, particularly evident in Turkish breeds. Our findings provide valuable insights into the genetic landscape of Turkish cattle, offering a crucial foundation for informed conservation and breeding strategies aimed at preserving these breeds' genetic integrity and heritage.

Genetic diversity, population structure analysis and codon substitutions of Indicine Badri cattle using ddRAD sequencing

The present study was carried out on 96 animals representing three distinct colour variants of Badri cattle to investigate the genetic diversity, population structure and substitution mutations in the genetic codons due to single nucleotide variations. The DNA samples of 96 Badri cows were genotyped using a double digestion restriction associated DNA (ddRAD) sequencing approach. A standardized bioinformatics pipeline was employed to identify single nucleotide polymorphisms (SNPs), initially detecting 7,168,552 SNPs through alignment with the Bos indicus reference genome assembly. Subsequent stringent quality filtration yielded 65,483 high-confidence SNPs for downstream analysis. Genetic diversity analysis of the Badri cattle population resulted in average values of 0.145, 0.088, and 0.091 for Shannon's diversity Index (I), Simpson's Diversity (h), and Simpson's Unbiased Diversity (uh), respectively. Genetic similarities between the black and brown, black and grey, and brown and grey Badri variants were found to be 0.9972, 0.9980 and 0.9970, respectively. Tajima's D diversity value was observed to be significant and positive for 99.29% of high-confidence SNPs (65,483). STRUCTURE analysis showed admixture among the three Badri colour variants, suggesting a lack of genetic differentiation. Annotation of high-confidence SNPs regarding genetic codon changes indicated maximum substitutions in the GGC with GGT (22 occurrences), followed by AAC to AGC (20 occurrences), GAA to TAA (19 occurrences) and CAA to CAG (19 occurrences). The study concludes there are genetic similarities among colour variants, lack of rare alleles, balancing selection, sudden population contraction and genetic codon substitutions within the Badri cattle population. Insights derived from SNP data analysis hold potential significance for conservation initiatives and breed improvement programs for indicine cattle.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

The Identification of Functional Genes Affecting Fat-Related Meat Traits in Meat-Type Pigeons Using Double-Digest Restriction-Associated DNA Sequencing and Molecular Docking Analysis

The Chinese indigenous Shiqi (SQ) pigeon and the imported White King (WK) pigeon are two meat-type pigeon breeds of economical and nutritional importance in China. They displayed significant differences in such meat quality traits as intramuscular fat (IMF) content and fatty acid (FA) compositions in the breast muscles. In this study, we aimed to screen candidate genes that could affect fat-related meat quality traits in meat-type pigeons. We investigated the polymorphic variations at the genomic level using double-digest restriction-associated DNA (ddRAD) sequencing in 12 squabs of SQ and WK pigeons that exhibited significant inter-breed differences in IMF content as well as FA and amino acid compositions in the breast muscles, and screened candidate genes influencing fat-related traits in squabs through gene ontology analysis and pathway analysis. By focusing on 6019 SNPs, which were located in genes with correct annotations and had the potential to induce changes in the encoded proteins, we identified 19 genes (ACAA1, ACAA2, ACACB, ACADS, ACAT1, ACOX3, ACSBG1, ACSBG2, ACSL1, ACSL4, ELOVL6, FADS1, FADS2, HACD4, HADH, HADHA, HADHB, MECR, OXSM) as candidate genes that could affect fat-related traits in squabs. They were significantly enriched in the pathways of FA metabolism, degradation, and biosynthesis (p < 0.05). Results from molecular docking analysis further revealed that three non-synonymous amino acid alterations, ACAA1(S357N), ACAA2(T234I), and ACACB(H1418N), could alter the non-bonding interactions between the enzymatic proteins and their substrates. Since ACAA1, ACAA2, and ACACB encode rate-limiting enzymes in FA synthesis and degradation, alterations in the enzyme-substrate binding affinity may subsequently affect the catalytic efficiency of enzymes. We suggested that SNPs in these three genes were worthy of further investigation into their roles in explaining the disparities in fat-related traits in squabs.

A genome-wide association study of coat color in Chinese Rex rabbits

Coat color is an important phenotypic characteristic of the domestic rabbit (Oryctolagus cuniculus) and has specific economic importance in the Rex rabbit industry. Coat color varies considerably among different populations of rabbits, and several causal genes for this variation have been thoroughly studied. Nevertheless, the candidate genes affecting coat color variation in Chinese Rex rabbits remained to be investigated. In this study, we collected blood samples from 250 Chinese Rex rabbits with six different coat colors. We performed genome sequencing using a restriction site-associated DNA sequencing approach. A total of 91,546 single nucleotide polymorphisms (SNPs), evenly distributed among 21 autosomes, were identified. Genome-wide association studies (GWAS) were performed using a mixed linear model, in which the individual polygenic effect was fitted as a random effect. We detected a total of 24 significant SNPs that were located within a genomic region on chromosome 4 (OCU4). After re-fitting the most significant SNP (OCU4:13,434,448, p = 1.31e-12) as a covariate, another near-significant SNP (OCU4:11,344,946, p = 7.03e-07) was still present. Hence, we conclude that the 2.1-Mb genomic region located between these two significant SNPs is significantly associated with coat color in Chinese Rex rabbits. The well-studied coat-color-associated agouti signaling protein (ASIP) gene is located within this region. Furthermore, low genetic differentiation was also observed among the six coat color varieties. In conclusion, our results confirmed that ASIP is a putative causal gene affecting coat color variation in Chinese Rex rabbits.