Whole genome detection of recent selection signatures in Sarabi cattle: a unique Iranian taurine breed

Copy number variants selected during pig domestication inferred from whole genome resequencing

Over extended periods of natural and artificial selection, China has developed numerous exceptional pig breeds. Deciphering the germplasm characteristics of these breeds is crucial for their preservation and utilization. While many studies have employed single nucleotide polymorphism (SNP) analysis to investigate the local pig germplasm characteristics, copy number variation (CNV), another significant type of genetic variation, has been less explored in understanding pig resources. In this study, we examined the CNVs of 18 Wanbei pigs (WBP) using whole genome resequencing data with an average depth of 12.61. We identified a total of 8,783 CNVs (~30.07 Mb, 1.20% of the pig genome) in WBP, including 8,427 deletions and 356 duplications. Utilizing fixation index (Fst), we determined that 164 CNVs were within the top 1% of the Fst value and defined as under selection. Functional enrichment analyses of the genes associated with these selected CNVs revealed genes linked to reproduction (SPATA6, CFAP43, CFTR, BPTF), growth and development (NR6A1, SMYD3, VIPR2), and immunity (PARD3, FYB2). This study enhances our understanding of the genomic characteristics of the Wanbei pig and offers a theoretical foundation for the future breeding of this breed.

CeRNA regulates network and expression and SNP effect on NFKBIA of cashmere fineness

In this study, a total of 1140 Liaoning Cashmere Goats (LCG) were genotyped for single nucleotide polymorphism (SNP) of NFKBIA gene. There are 15 SNPs and 7 genotypes have been found, and G1547A (GG) genotype has been associated with cashmere fineness and cashmere yield. An integrated ceRNA regulatory network of NFKBIA gene was made. To prove NFKBIA and these non-coding RNAs (ncRNAs) may be related to cashmere fineness, we performed qPCR on these ncRNA in LCG coarse type skin (CT-LCG) and LCG fine type skin (FT-LCG). The result of qPCR showed lncRNA XLOC_011060 and ciRNA452 are at high expression level in CT-LCG, all miRNAs appear high expressed in FT-LCG, and mir-93 was the most significant difference between CT-LCG and FT-LCG. In addition, five miRNAs were selected for qPCR in different genotypes. The qPCR results showed that mir-93 might negatively regulate cashmere fineness and mir-17-5p may play a positive role in regulating cashmere fineness of individuals with G1355A (AG) genotype. These results demonstrated that NFKBIA gene is associated with cashmere fineness of LCG and G1547A (GG) genotype is the preferred marker genotype for cashmere fineness.

Polymorphism in ovine ADCY8 gene and its association with residual feed intake in Hu sheep

Feed efficiency makes up a large part of sheep production, which also has a crucial impact on the economic benefits of producers. This study explores the polymorphism of Adenylyl cyclase 8 gene associated with feed efficiency and detects the expression characteristics of ADCY8 in ten tissues of Hu sheep. The polymorphism of ADCY8 was recognized by using PCR amplification and Sanger sequencing, KASPar technology was used for genotyping subsequently, and the relationship between SNP and RFI is also studied. The results indicated that an intronic mutation g.24799148 C > T (rs 423395741) was identified in ADCY8, and association analysis showed that the SNP g.24799148 C > T (rs 423395741) was significantly associated with RFI at 100-120, 100-140, 100-160, and 100-180 days (p < 0.05). The quantitative real-time PCR (qRT-PCR) result showed that ADCY8 was expressed in ten tissues, and the expression of ADCY8 gene in rumen tissue was significantly higher than in the other tested tissues. Therefore, these results indicated that the ADCY8 mutation locus may be used as a candidate molecular marker for evaluating the feed efficiency of Hu sheep.

Study on lactation performance and development of KASP marker for milk traits in Xinjiang donkey (Equus asinus)

Donkey milk has high nutritional and medicinal value, but there are few researches in donkey milk traits, especially on genome. The whole lactation of 89 donkeys was recorded and it was found that Xinjiang donkey had good lactation performance while great differences among individuals. In our previous study, four genes including LGALS2, NUMB, ADCY8 and CA8 were identified as milk-associated with Chinese Kazakh house, based on Equine 670k Chip genomic analysis. And then 15 SNPs of the four key genes were conducted for genotyping in Xinjiang donkey in this study, one of Chinese indigenous breed, 14 SNPs were successful classified. And those SNPs were correlation analysis with milk yield of Xinjiang donkeys. The results showed that NUMB g.46709914T > G was significantly correlated with daily milk yield of Xinjiang donkey in the early, middle, and late periods, while ADCY8 g.48366302T > C, CA8 g.89567442T > G and CA8 g.89598328T > A were significantly correlated with lactation in the late periods. These results indicate that NUMB g.46709914T > G can be as markers of candidate genes for lactating traits in donkeys, SNPs of ADCY8 and CA8 as potential. Our findings will not only help confirm key genes for donkey milk traits, but also provide future for genomic selection in donkeys.

Assessing genomic diversity and selective pressures in Bashan cattle by whole-genome sequencing data

Specific ecological environments and domestication have continuously influenced the physiological characteristics of Chinese indigenous cattle. Among them, Bashan cattle belongs to one of the indigenous breeds. However, the genomic diversity of Bashan cattle is still unknown. Published whole-genome sequencing (WGS) data of 13 Bashan cattle and 48 worldwide cattle were used to investigate the genetic composition and selection characteristics of Bashan cattle. The population structure analysis revealed that Bashan cattle harbored ancestries with East Asian taurine and Chinese indicine. Genetic diversity analysis implied the relatively high genomic diversity in Bashan cattle. Through the identification of containing >5 nsSNPs or frameshift mutations genes in Bashan cattle, a large number of pathways related to sensory perception were discovered. CLR, θπ ratio, FST, and XP-EHH methods were used to detect the candidate signatures of positive selection in Bashan cattle. Among the identified genes, most of the enriched signal pathways were related to environmental information processing, biological systems, and metabolism. We mainly reported genes related to the nervous system (HCN1, KATNA1, FSTL1, GRIK2, and CPLX2), immune (CD244, SLAMF1, LY9, and CD48), and reproduction (AKR1C1, AKR1C3, AKR1C4, and TUSC3). Our findings will be significant in understanding the molecular basis underlying phenotypic variation of breed-related traits and improving productivity in Bashan cattle.

Polymorphism and expression of the HMGA1 gene and association with tail fat deposition in Hu sheep

Hu sheep is an excellent short fat-tailed breed in China. Fat deposition in Hu sheep tail affects carcass quality and consumes a lot of energy, leading to an increase in feed cost. The objective of this study was to analyze the effects of HMGA1 polymorphism on tail fat weight in Hu sheep. Partial coding and non-coding sequences of HMGA1 were amplified with PCR and single nucleotide polymorphisms (SNP) of HMGA1 in 1163 Hu sheep were detected using DNA sequencing and KASPar technology. RT-qPCR analysis was performed to test HMGA1 expression in different tissues. The results showed that the expression of HMGA1 was higher in the duodenum, liver, spleen, kidney, and lung than in the heart, muscle, rumen, tail fat, and lymph. A mutation, g.5312 C > T, was detected in HMGA1; g.5312 C > T was significantly associated with tail fat weight, relative weight of tail fat (body weight), and relative weight of tail fat (carcass) (p < 0.05). The tail fat weight of the TT genotype was remarkably higher than that of the CC and TC genotypes. Therefore, HMGA1 can be used as a genetic marker for marker-assisted selection of tail fat weight in Hu sheep.

Effects of starter feeding and early weaning on developmental expressions of IGF-I gene in liver and IGF-IR gene in rumen of lambs from birth to eighty-four days of age

The early weaning and starter feeding have significant effects on lamb growth and digestive tract development. However, it is not clear whether the expression of IGF-1/IGF-1R genes could be affected by feeding starter and weaning. Therefore, a total of 102 Hu male neonatal lambs were randomly divided into 3 groups: fed with starter at 42d + weaned at 56d (group A), fed with starter at 7d + weaned at 28d (group B) and fed with starter at 7d + weaned at 56d (group C), to explore the effects of starter feeding and weaning age on developmental expressions of IGF-I gene in liver and IGF-IR gene in rumen of Hu sheep. The results showed that IGF-I and IGF-IR genes were expressed extensively in various tissues of lambs, the expression of IGF-I was significantly higher in liver (p < 0.01), while the expression of IGF-IR was higher in rumen among gastrointestinal tissues. The mRNA level of IGF-I of group C was strikingly higher than that of group A at 28, 70, and 84d (p < 0.01) in liver, respectively. Rumen mRNA level of IGF-IR of group C was prominently higher than that of group A at 70d and 84d (p < 0.05), respectively. However, IGF-IR expression of group A was higher than that of group C at 14 and 42d (p < 0.05). The mRNA level of IGF-I of group B was significantly higher than that of group C at 42, 56, 70, and 84d (p < 0.01) in liver, while IGF-IR expression of group B was higher than that of group C at 70d (p < 0.01). In conclusion, early weaning and starter feeding affected liver IGF-I mRNA and rumen IGF-IR mRNA expression.

Assessing Population Structure and Signatures of Selection in Wanbei Pigs Using Whole Genome Resequencing Data

Wanbei pig (WBP) is one of the indigenous pig resources in China and has many germplasm characteristics. However, research on its genome is lacking. To assess the genomic variation, population structure, and selection signatures, we resequenced 18 WBP for the first time and performed a comprehensive analysis with resequenced data of 10 Asian wild boars. In total, 590.03 Gb of data and approximately 41 million variants were obtained. Polymorphism level (θπ) ratio and genetic differentiation (fixation index)-based cross approaches were applied, and 539 regions, which harbored 176 genes, were selected. Functional analysis of the selected genes revealed that they were associated with lipid metabolism (SCP2, APOA1, APOA4, APOC3, CD36, BCL6, ADCY8), backfat thickness (PLAG1, CACNA2D1), muscle (MYOG), and reproduction (CABS1). Overall, our results provide a valuable resource for characterizing the uniqueness of WBP and a basis for future breeding.

Expression and polymorphisms of CD8B gene and its associations with body weight and size traits in sheep

The growth traits are economically important traits in sheep. Improving growth rates will increase the profitability of producers. The aim of this study was to identify alleles of CD8B (encoding T-cell surface glycoprotein CD8 beta chain) that are aberrantly expressed in different tissues and to assess the effects and associations of its different genotypes on weight and size traits in sheep. Using quantitative real-time reverse transcription PCR arrays, expression profiling of CD8B was performed in various organs and tissues. CD8B was ubiquitously expressed, with very high expression in the lung, spleen, lymph, duodenum, and liver. One intronic mutation (chr3:62,718,030 (Oar_rambouillet_v1.0, same below) G > A) was identified using pooled DNA sequencing. Subsequently, the variants (AA, AG, and GG) were genotyped using the KASPar^® PCR single nucleotide polymorphism (SNP) genotyping system. The results of association analysis with body weight and body size traits in 1304 sheep showed that increases in multiple phenotypic traits correlated with the AA genotype (body weight, p < 0.05; body length, p < 0.05). Thus, SNP chr3:62,718,030 G > A is a promising molecular marker for marker-assisted selection in sheep breeding.

Polymorphisms in the ovine GP5 gene associated with blood physiological indices

This study aimed to analyze the effects of polymorphisms in GP5 on blood physiological indices of 1065 sheep. The coefficients of variation of the red blood cell count (RBC), hemoglobin concentration (HGB), mean platelet volume (MPV), mean erythrocyte hemoglobin content (MCH), and red blood cell distribution-coefficient of variation (RDW-CV) were greater than 10%, and there was a very significant correlation between the main indices such as RBC, white blood cell, and platelet count (PLT) and most other indices (p < 0.01). qRT-PCR showed that GP5 was expressed in the heart, liver, spleen, lung, kidney, rumen, duodenum, muscle, tail fat, and lymph tissue, with significantly higher expression in the lymph. Subsequently, we detected single nucleotide polymorphisms (SNPs) in GP5 from group, which identified synonymous mutation g.657 T > C in the first exon of GP5. Association analysis showed significant correlations between the SNP and the physiological traits (p < 0.05), in which the RBC, neutrophilic granulocyte (NEUT) and RDW-CV values in sheep with the TC genotype and TT genotype were markedly lower than those in the CC genotype (p < 0.05). Thus, GP5 polymorphisms could be candidate biomarkers to screen blood physiological indices.

Associations of CALM1 and DRD1 polymorphisms, and their expression levels, with Taihang chicken egg-production traits

Egg production by hens is an important reproductive performance index in the poultry industry. To investigate the effects of the CALM1 and DRD1 genes on egg production in chicken, their mRNA expression and single nucleotide polymorphisms (SNP) levels were investigated, and bioinformatics and egg-production association analyses were performed. Three SNPs (g.44069941G > A and g.44069889A > G in CALM1 and g.10742639C > T in DRD1) were detected in the exons and introns of CALM1 and DRD1 in 400 Taihang chickens. Among them, g.44069941G > A was significantly associated with Taihang chicken egg production on the 500th day (p < 0.05), whereas g.10742639C > T was significantly associated with the 300th day (p < 0.05). The expression levels of CALM1 and DRD1 in ovarian tissues of a high-yielding Taihang group were greater than in a low-yielding group (p < 0.05). The bioinformatics analysis revealed that the mutations influenced the mRNA secondary structures of CALM1 and DRD1. This study provides new insights into the potential effects of CALM1 and DRD1 polymorphisms on chicken egg production. The two SNPs g.44069941G > A and g.10742639C > T are potential molecular markers for improving the reproductive traits of Taihang chicken.

Transcriptome profiling of LncRNAs in sheep tail fat deposition

LncRNAs have recently received special attention due to their critical role in many important biological processes. There are few reports on its regulatory function in sheep fat deposition. In this study, two sheep populations with different tail types in Xinjiang, Bashibai sheep (fat-tailed) and the hybrid population of Bashibai sheep and wild argali (small-tailed) were selected for whole transcriptome sequencing from their tail tissues. First, 728 differentially expressed LncRNAs of tail fat between Bashibai and F2 sheep were identified by RNA-seq. Second, the tissue expression profile and relative expression difference between Bashibai and F2 sheep of 2 of 728 DE LncRNAs were analyzed by RT-PCR. LncRNA-MSTRG.24995 was highly expressed in tail fat, while lncRNA-MSTRG.36913 was highly expressed in subcutaneous fat. In addition, the expressions of LncRNA-MSTRG.24995 and LncRNA-MSTRG.36913 in tail fat of F2 sheep were significantly lower than that of Bashibai sheep, while those patterns in longissimus dorsi, quadriceps femoris and rumen were reversed. Third, the expression pattern of target genes FASN and THRSP in each tissue was similar with that of corresponding LncRNAs. The LncRNA-MSTRG.24995 directly affects tail fat deposition by FASN gene, while the LncRNA-MSTRG.36913 indirectly affects that by THRSP gene. This will help us to understand molecular mechanism of fat tail deposition from transcriptomic perspectives.

Identification of polymorphisms in the MSTN and ADRB3 genes associated with growth and ultrasound carcass traits in Kajli sheep

The aim of this investigation was to find single nucleotide polymorphisms (SNPs) in the myostatin (MSTN) and the beta-3 adrenergic receptor (ADRB3) genes associated with growth and ultrasound carcass traits in Kajli sheep. The five growth traits were birth weight (BWT-EBV), 120-day weight (120DWT-EBV), 180-day weight (180DWT-EBV), 270-day weight (270DWT-EBV), and 365-day weight (365DWT-EBV). The three ultrasound carcass traits were width (WLD) and depth of longissimus dorsi (DLD) and back fat thickness (BFT). The analysis of the MSTN sequence revealed one non-synonymous substitution (c.197T > A) in exon 1, one single nucleotide substitution (c.373 + 18G > T) in intron 1, and one synonymous substitution (c.861T > A) in exon 3. However, there were four single nucleotide synonymous substitutions (c.130C > T, c.294C > G, c.579G > T, and c.654C > G) in exon 1 of the ADRB3 gene. All the SNPs in the MSTN gene, except for c.373 + 18G > T, were in Hardy-Weinberg Equilibrium (HWE). Conversely, none of the SNPs found in ADRB3 were in HWE. Two of the MSTN SNPs (c.197T > A and c.373 + 18G > T) had significant associations with all evaluated growth and ultrasound carcass traits. The SNPs c.130C > T and c.294C > G in ADRB3 were significantly associated with 180DWT-EBV. Collectively, these findings indicate that several SNPs in the studied genes were significantly related to growth and carcass traits in Kajli sheep.

Analysis of the eighth intron polymorphism of NR6A1 gene in sheep and its correlation with lumbar spine number

For revealing molecular markers related to the traits of multiple lumbar vertebrae in sheep, we analyze the relationship between NR6A1 gene polymorphism and lumbar vertebrae number traits in Xinjiang Kazakh sheep. Lumbar muscle tissues were collected from 6-lumbar spine (L6) Kazak sheep and 7-lumbar spine (L7) Kazak sheep and the intron-8 of NR6A1 gene was amplified by PCR. The SNP locus was detected by the PCR-SSCP method. One-Way ANOVA and an Independent Chi-square Test is adopted to analyze the genotype association with lumbar spine number variation. There were two SNP loci in the intron-8 of the NR6A1 gene: IVS8-188 and IVS8-281. One-Way ANOVA and Independent Chi-square Test indicated a significant association between IVS8-281 and lumbar spine number. The SNP locus of NR6A1 gene intron 8 (IVS8-281G > A) could play a certain role in the variation of lumbar spine number in Xinjiang Kazakh sheep and demonstrates potential to accelerate the sheep breeding of selection process.

Shared Ancestry and Signatures of Recent Selection in Gotland Sheep

Gotland sheep, a breed native to Gotland, Sweden (an island in the Baltic Sea), split from the Gute sheep breed approximately 100 years ago, and since, has probably been crossed with other breeds. This breed has recently gained popularity, due to its pelt quality. This study estimates the shared ancestors and identifies recent selection signatures in Gotland sheep using 600 K single nucleotide polymorphism (SNP) genotype data. Admixture analysis shows that the Gotland sheep is a distinct breed, but also has shared ancestral genomic components with Gute (~50%), Karakul (~30%), Romanov (~20%), and Fjällnäs (~10%) sheep breeds. Two complementary methods were applied to detect selection signatures: A Bayesian population differentiation F_ST and an integrated haplotype homozygosity score (iHS). Our results find that seven significant SNPs (q-value < 0.05) using the F_ST analysis and 55 significant SNPs (p-value < 0.0001) using the iHS analysis. Of the candidate genes that contain significant markers, or are in proximity to them, we identify several belongings to the keratin genes, RXFP2, ADCY1, ENOX1, USF2, COX7A1, ARHGAP28, CRYBB2, CAPNS1, FMO3, and GREB1. These genes are involved in wool quality, polled and horned phenotypes, fertility, twining rate, meat quality, and growth traits. In summary, our results provide shared founders of Gotland sheep and insight into genomic regions maintained under selection after the breed was formed. These results contribute to the detection of candidate genes and QTLs underlying economic traits in sheep.

Genomic scans for selection signatures revealed candidate genes for adaptation and production traits in a variety of cattle breeds

Domestication and selection are the major driving forces responsible for the determinative genetic variability in livestock. These selection patterns create unique genetic signatures within the genome. BovineSNP50 chip data from 236 animals (seven indicine and five taurine cattle breeds) were analyzed in the present study. We implemented three complementary approaches viz. iHS (Integrated haplotype score), ROH (Runs of homozygosity), and F_ST, to detect selection signatures. A total of 179, 56, and 231 regions revealed 518, 277, and 267 candidate genes identified by iHS, ROH, and F_ST methods, respectively. We found several candidate genes (e.g., NCR3, ARID5A, HIST1H2BN, DEFB4, DEFB7, HSPA1L, HSPA1B, and DNAJB4) related to production traits and the adaptation of indigenous breeds to local environmental constraints such as heat stress and disease susceptibility. However, further studies are warranted to refine the findings using a larger sample size, whole-genome sequencing, and/or high density genotyping.