A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes.

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5' untranslated region (5' UTR) of HOXB13 is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

circUBE3C modulates myoblast development by binding to miR-191 and upregulating the expression of p27.

Noncoding RNAs, including miRNAs (microRNAs) and circRNAs (circular RNA), are crucial regulators of myoblast proliferation and differentiation during muscle development. However, the specific roles and molecular mechanisms of circRNAs in muscle development remain poorly understood. Based on the existing circRNA-miRNA-mRNA network, our study focuses on circUBE3C, exploring its differential expression in fetal and adult muscle tissue of the cattle and investigating its impact on myoblast proliferation, apoptosis, and differentiation. The functional analysis of overexpression plasmids and siRNAs (small interfering RNAs) targeting circUBE3C was comprehensively evaluated by employing an array of advanced assays, encompassing CCK-8 (cell counting kit-8), EdU (5-ethynyl-20-deoxyuridine), flow cytometry, western blot analysis, and RT-qPCR. In vivo investigations indicated that overexpression of circUBE3C impedes the process of skeletal muscle regeneration. Mechanistically, we demonstrated that circUBE3C interacts with miR-191 and alleviates the suppression of p27 through cytoplasmic separation, bioinformatics prediction, dual-luciferase reporter assay, and RIP (RNA immunoprecipitation). Our findings indicate that the novel circRNA circUBE3C competitively binds to miR-191, thereby inhibiting proliferation and promoting apoptosis in bovine primary myoblasts and unveiling a regulatory pathway in bovine skeletal muscle development. These findings expand our understanding of circRNA functions in mammals and provide a basis for further exploration of their role in myogenesis and muscle diseases.

Genetic variant of the sheep E2F8 gene and its associations with litter size.

The economic efficiency of sheep breeding, aiming to enhance productivity, is a focal point for improvement of sheep breeding. Recent studies highlight the involvement of the Early Region 2 Binding Factor transcription factor 8 (E2F8) gene in female reproduction. Our group's recent genome-wide association study (GWAS) emphasizes the potential impact of the E2F8 gene on prolificacy traits in Australian White sheep (AUW). Herein, the purpose of this study was to assess the correlation of the E2F8 gene with litter size in AUW sheep breed. This work encompassed 659 AUW sheep, subject to genotyping through PCR-based genotyping technology. Furthermore, the results of PCR-based genotyping showed significant associations between the P1-del-32bp bp InDel and the fourth and fifth parities litter size in AUW sheep; the litter size of those with genotype ID were superior compared to those with DD and II genotypes. Thus, these results indicate that the P1-del-32bp InDel within the E2F8 gene can be useful in marker-assisted selection (MAS) in sheep.

Comprehensive comparison of two types of algorithm for circRNA detection from short-read RNA-Seq.

MOTIVATION: Circular RNA is generally formed by the 'back-splicing' process between the upstream splice acceptor and the downstream donor in/not in the regulation of the corresponding RNA-binding proteins or cis-elements. Therefore, more and more software packages have been developed and they are mostly based on the identification of the back-spliced junction reads. However, recent studies developed two software tools that can detect circRNA candidates by constructing k-mer table or/and de Bruijn graph rather than reads mapping. RESULTS: Here, we compared the precision, sensitivity and detection efficiency between software tools based on different algorithms. Eleven representative detection tools with two types of algorithm were selected for the overall pipeline analysis of RNA-seq datasets with/without RNase R treatment in two cell lines. Precision, sensitivity, AUC, F1 score and detection efficiency metrics were assessed to compare the prediction tools. Meanwhile, the sensitivity and distribution of highly expressed circRNAs before and after RNase R treatment were also revealed by their enrichment, unaffected and depleted candidate frequencies. Eventually, we found that compared to the k-mer based tools, CIRI2 and KNIFE based on reads mapping had relatively superior and more balanced detection performance regardless of the cell line or RNase R (-/+) datasets. AVAILABILITY AND IMPLEMENTATION: All predicted results and source codes can be retrieved from https://github.com/luffy563/circRNA_tools_comparison. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An upstream deletion polymorphism within the goat Period circadian regulator 1 (PER1) gene was associated with growth traits.

The Period circadian regulator (PER) gene family, including PER1, PER2 and PER3, codes transcriptional repressors which could accurately control biological rhythms. PER1/2 gene was proved to be associated with bone mass and PER1 gene was associated with insulin-like growth factor binding proteins 3 (IGFBP3) levels in serum. However, it was few studies reported genetic effects of PER gene on growth traits at the individual level. In this study, we identified the potential insertion/deletion (indel) loci in PER1/2/3 gene, and then explored the relationship between goat growth traits and the frequency of genotype in Shaanbei white cashmere goats (n = 827). As a result, a 9 bp indel within PER1 gene (g.27528003-27528011 del.TGCTGCTGC; rs642467689) was identified using molecular biology techniques. In addition, there existed significant correlation between the 9 bp indel and body height, height at hip cross, chest depth, body length index and cannon circumference index of goats. These results suggested that the 9 bp indel variation in PER1 gene was associated with goat growth traits, providing the theoretical basic for the role of PER1 gene in goat breeding.

Three novel simple sequence repeats (SSRs) identified by MALDI-TOF-MS method were associated with backfat in pig.

Backfat trait is an important economic trait and highly heritable, but difficult to evaluate. Thus, it is of great significance to explore optimal backfat thickness of pigs by using marker-assisted selection (MAS) to speed up its breeding process and improve economic efficiency. This study aimed to investigate the relationship between genetic variations (e.g., SSRs) and backfat of Qinghai Bamei pigs using MALDI-TOF Mass Spectrometry (MALDI-TOF-MS). Herein, five alternative SSR loci (namely V1, V2, V3, V4 and V5) were selected for subsequent detection. The results suggested that 3 (141-, 143- and 145-), 3 (128-, 130- and 132-), 2 (160- and 162-), 2 (136- and 139-) and 3 (170-, 184- and 192-) alleles of V1, V2, V3, V4 and V5 were found, respectively. Subsequent analysis showed that there was linkage equilibrium among five SSRs and Hap19 (13.1%) (141-/132-/160-/139-/192-) had the highest haplotype frequency. Among these five SSR loci, V1, V2 and V3 loci were significantly associated to the backfat of Qinghai Bamei sows. These findings enriched the study of SSRs in Qinghai Bamei pigs, and (AC)n (Chr15:85485851-85485995), (AC)n (Chr10:52724583-52724713) and (TG)n (Chr4:90732644-90732802) could be utilized as the candidate locus for MAS in pig industry.HIGHLIGHTSFive novel SSR loci was identified in pigs through MALDI-TOF MS.V1, V2 and V3 loci was were significantly associated to the backfat of pigs.

The InDel variants of sheep IGF2BP1 gene are associated with growth traits.

Insulin-like growth factor 2 mRNA binding protein 1 (IGF2BP1) plays positive roles in the growth, proliferation of cells and early embryos development by binding mRNA targets. Recently, it had been shown that some polymorphic loci within IGF2BP1 gene were associated with growth traits in animals, especially in goats. Therefore, it has been hypothesized that some variants within IGF2BP1 gene may be also involved in growth traits of sheep. Nine insertion/deletion (InDel) mutations within IGF2BP1 were identified and three loci were polymorphic. Meanwhile, the association analyses between three InDels and growth traits were carried out in 745 sheep. The results showed that all InDels included 5 bp InDel in downstream region, 9 bp InDel in intron 4 and 15 bp InDel in intron 2 within IGF2BP1 were significantly associated with growth traits (p<.05). Furthermore, at 5 and 9 bp InDel loci, the individuals of heterozygous genotype (ID) had superior growing performance especially at body weight (BW). In all, three InDels were crucial variants correlated with growth traits and could be applied in marker-assisted selection (MAS) in sheep.

Survey of the relationship between polymorphisms within the BMPR1B gene and sheep reproductive traits.

The BMPRIB gene is one of the main genes that can be used as a molecular genetic marker for the early selection of highly productive ewes. It is well-documented that the p.Q249R (g.746A > G) is the first mutation in the kinase domain of the BMPR1B gene that is highly related to increased ovulation rate and litter size. It is likely that the presence of the p.Q249R mutation in the sheep population is one of the factors contributing to the outstanding productivity of the sheep. Moreover, in recent years, researchers have been explored other polymorphisms in the BMPR1B gene with respect to reproductive traits in sheep. Therefore, we carried out the current study to evaluate the association between polymorphisms in this gene and sheep litter size from all appropriate studies. As a result, among 41 polymorphisms in the ovine BMPRIB gene, eight variants, including p.Q249R (g.746A > G), g.29362047T > C, g.29427689G > A, BMPR1B-2 (ss:1960972599), g.29382337G > A, g.29382340G > A, rs1092293287 (10 bp insertion/deletion) and g.29380965A > G were found to be associated with litter size in sheep. This systematic analysis presents the most current data evidence for BMPRIB polymorphisms, highlighting the need for further large-scale studies to determine more important variants.

Screen of small fragment mutations within the sheep thyroid stimulating hormone receptor gene associated with litter size.

The thyroid stimulating hormone receptor (TSHR), a glycoprotein hormone receptor, plays an important role in metabolic regulation and photoperiod control in the time of reproduction in birds and mammals. Previous genome-wide association studies revealed that the TSHR gene was related to reproduction and its function was identified in female reproduction, but rare studies reported the polymorphism of TSHR gene. However, the molecular mutations of the TSHR gene in sheep have not been reported so far. Herein, we explored potential polymorphisms of the sheep TSHR gene, and a 29 bp nucleotide sequence variant (rs1089565492) was identified in the AUW sheep. There were three genotypes of the 29 bp variant locus detected which named 'II' 'DD' and 'ID' been identified. Association analysis results showed the 29 bp variant was significantly associated with the litter size of the AUW sheep (p < 0.05). This finding suggests that the 29 bp nucleotide sequence variant within TSHR gene could be a candidate marker of reproduction traits for sheep breeding improving through the marker-assisted selection (MAS).

A novel indel within the bovine SEPT7 gene is associated with ovary length.

The ovary can generate oocytes and secrete female hormones and thus is of great significance to animal fertility. In turn, the functioning of this organ has an effect on the profit margins of the livestock breeding industry. As the development-regulating gene and target gene of miR-202, SEPT7 might play an important role in ovarian growth. Therefore, we hypothesized that SEPT7 is related to ovarian traits owing to the regulation of gonad-specific miR-202. To further investigate the connection between bovine SEPT7 and ovarian development, we analyzed data from 408 samples. After genotyping and analyzing three selected loci, we found that two out of the three loci (L1 and L5) were polymorphic, of which the minimum allelic frequencies were 0.417 (L1) and 0.094 (L5). Moreover, one novel indel L1 of SEPT7 was associated with ovarian length (p < 0.05). More specifically, individuals with II and ID genotypes have longer ovaries than those with the DD genotype. Our work shows that SEPT7 can be selected as a testing marker gene for animal fertility. Our findings contribute to improving the prospects of the cattle industry and the wider use of genetic techniques in breeding.

Investigation of Genetic Effects of Nucleotide Variants Within the Goat PRNT Gene on Growth Performance.

Our previous study has firstly pointed that three nucleotide variants (g.-11C > T, g.117A > G, and g.149C > T) of the goat PRNT gene can significantly influence litter size. Given litter size is positively correlated with growth performance, we consider whether the PRNT gene also acts on the growth performance in goats. In this work, a correlation analysis among different litter size types and growth traits of Shaanbei white cashmere (SBWC) goats was performed, and results showed that a positive correlation did exist in our detected population (P < 0.01). Then, the association among different genotypes of three variations and goat growth performance was measured. Our results pointed to g.117A > G being significantly associated with the cannon circumference (P = 4.60E-05) while no significant effect was found between another two SNPs and growth traits after the Bonferroni's correction (P*n < 0.05). Together, this is the first report about the influence of the PRNT gene on the growth of goat and g.117A > G can be regarded as a possible DNA marker applying for MAS breeding.

A deletion mutation within the goat AKAP13 gene is significantly associated with litter size.

A-kinase anchoring protein 13 (AKAP13) is one of the AKAP protein family members, which is correlated with estrogen receptors (ERs) and progesterone receptor (PR) activity. Consequently, the AKAP13 gene is considered to be one of the candidate genes for regulating female fertility. Hence, the objectives of this study were to discover the potential insertion/deletion (indel) variants within the AKAP13 gene and evaluate their associations with litter size of Shaanbei white cashmere goats (SBWC) to screen candidate genes for the molecular marker-assisted selection (MAS). Ultimately, we found the 16-bp deletion of AKAP13 gene which displayed three genotypes (II, ID and DD). However, it was not confirmed to Hardy-Weinberg equilibrium (HWE) in the tested population. Statistical analysis demonstrated that this 16-bp indel locus was significantly associated with litter size in goats (p < 0.05), in which the ID genotype was a key genotype for increasing litter size in goats. Besides, independent χ2 tests between different genotypes and litter size showed that high-prolific groups had higher frequency of the 'D' allele (p < 0.05). Briefly, AKAP13 gene is a candidate gene for improving fertility, and its 16-bp indel locus can be used as a valid DNA molecular marker for the MAS in goat breeding.

Insertion/deletion (Indel) variant of the goat RORA gene is associated with growth traits.

RAR related orphan receptor A (RORA), which encodes the retinoid-acid-related orphan receptor alpha (RORα), is a clock gene found in skeletal muscle. Several studies have shown that RORα plays an important role in bone formation, suggesting that RORA gene may take part in the regulation of growth and development. The purpose of this research is to study the insertion/deletion (indel) variations of the RORA gene and investigate the relationship with the growth traits of Shaanbei white cashmere (SBWC) goats. Herein, the current study identified that the P4-11-bp and P11-28-bp deletion sites are polymorphic among 12 pairs of primers within the RORA gene in the SBWC goats (n = 641). Moreover, the P11-28-bp deletion locus was significantly related to the body height (p = 0.046), height at hip cross (p = 0.012), and body length (p = 0.003). Both of P4-11-bp and P11-28-bp indels showed the moderate genetic diversity (0.25

Investigation on mRNA expression and genetic variation within goat SMAD2 gene and its association with litter size.

The SMAD family member 2 (SMAD2), a member of the TGF-beta superfamily, executes a significant part in the oogenesis and ovulation process. A genome-wide selective sweep analysis also found SMAD2 was different in the fertility groups of Laoshan dairy goats; whether this gene was linked to litter size was unknown. Therefore, SMAD2 was chosen to study its effects on Shaanbei white cashmere goat reproduction and mRNA expression profile. Herein, the mRNA expression level of SMAD2 was firstly determined in female goat tissues, revealing significant differences in mRNA levels of different tissues (p < 0.05), including ovary tissue, indicating a potential role for SMAD2 in goat prolificacy. Then, using six pairs of primers, only one indel locus (P3-Del-12-bp) was found to be polymorphic in goat SMAD2 (n = 501). ANOVA also revealed that a P3-Del-12-bp deletion was significantly related to first-born litter size (p = 0.037). The Chi-square (χ2) test revealed that the ID genotype was significantly more prevalent in mothers with multiple lambs (p = 0.01), indicating that heterozygous individuals (ID) are more likely to produce multiple lambs. Our findings suggest that the SMAD2 gene's P3-Del-12-bp deletion could be used to improve goat breeds by assisting with litter size selection.

Investigation of the association between goat DNMT3B gene polymorphism and growth traits.

The DNA methyltransferase 3 beta (DNMT3B) gene is key for DNA methylation and has been well recognized in regulating growth and development. A previous observation indicated that an 11-bp indel of DNMT3B affected the reproductive traits in goats, yet the effect of this polymorphism on body measurement traits in goats has not been reported. This study aims to investigate the associations between DNMT3B gene polymorphism and goat growth traits. We investigated this 11-bp indel in 2184 goats and three genotypes have been found in Shaanbei white cashmere goat (SBWC): insertion/insertion (II), deletion/deletion (DD) and insertion/deletion (ID). Only ID and DD genotypes were detected in Nubian goats and Guizhou heima goat (GZHM). The allele frequencies analyzed revealed that the 'D' allele frequencies were higher in all three goat breeds. Further association analysis demonstrated that this indel is markedly associated with the cannon circumference (CC) and cannon circumference index (CCI) of SBWC and cannon circumference (CC) of Nubian goats (p < .05). The CC and CCI are essential indicators to measure the growth status of goats. In summary, our study sheds some light on the potential impact of the 11-bp indel polymorphism of the DNMT3B gene on improving the growth traits in goats.

Relationships between novel nucleotide variants within the colony-stimulating factor 1 receptor (CSF1R) gene and mastitis indicators in sheep.

Colony-stimulating factor 1 receptor (CSF1R) plays an important role in the process of innate immunity and inflammation, thus it was hypothesized that the CSF1R gene might affect the occurrence of mammalian mastitis. The purpose of this study was to investigate the association between nucleotide variations of CSF1R gene and mastitis in Australian white sheep (AUWs). Two indel variants (Intron5-27 bp and Intron5-22 bp) within the CSF1R gene have been found in AUWs. The Chi-square test for different mastitis symptoms demonstrated that individuals without symptoms of mastitis had higher 'I' allele frequencies and 'II' genotype frequencies (p < 0.01). We found strong correlation between mastitis and lactation score through Pearson correlation analysis. Therefore, we also analyzed the relationship between the two indel loci and lactation, we found that the lactation ability of individuals with type II was stronger than that of DD genotype at the Intron5-22 bp (p < 0.05). Additionally, we found that the combined genotype of the two loci was significantly associated with mastitis (p < 0.01). These findings indicated that CSF1R mutations were significantly associated with mastitis, and could affect lactation performance, suggesting that two deletion sites could be used as the effective molecular markers against mastitis in sheep breeding.

Polymorphism within the GATA binding protein 4 gene is significantly associated with goat litter size.

GATA binding protein 4 (GATA4) is a typical transcription binding factor, and its main functions include regulating the proliferation, differentiation and apoptosis of ovarian granulosa cells, promoting spermatogenesis and sex differentiation, implying that this gene have possibly roles in animal reproduction. This study aims to detect five potential insertion/deletions (indels) of the GATA4 gene in 606 healthy unrelated Shaanbei white cashmere (SBWC) goats and analyze its association with the litter size. The electrophoresis and DNA sequencing identified two polymorphic indels (e.g., P4-Del-8bp and P5-Ins-9bp indel). Then T-test analysis showed that P4-Del-8bp was significantly correlated with litter size (p = 0.022) because of two different genotypes detected, e.g., insertion-deletion (ID) and deletion-deletion (DD), and the average litter size of individuals with DD genotype goats was higher than that of others. However, there was no correlation between P5-Ins-9bp and lambing of goats. Chi-square (X2) test found that the distribution of and P4-Del-8bp genotypes (X2 = 6.475, p = 0.011) was significantly different between single and multiple-lamb groups, while P5-Ins-9bp (X2 = 0.030, p = 0.862) was not. Therefore, these findings revealed that P4-Del-8bp polymorphism of goat GATA4 gene was a potential molecular marker significantly associated with litter size, which can be used for the marker-assisted selection (MAS) breeding to improve goat industry.

A 24-bp indel within the sheep AHR gene is associated with litter size.

Aryl Hydrocarbon Receptor (AHR) is a member of the PER-ARNT-SIM (PAS) family, which could mediate various biological processes, for instance, the balance of the immune system, cell proliferation, differentiation, vascular tissue remodeling and reproduction ability regulation. A previous research showed that the AHR gene exerted important functions on the pig reproduction, implying that it could serve as a candidate gene related to animal reproductive traits. Here, the aim of this work was to identify potential insertion/deletion (indel) mutations of the AHR gene in three sheep breeds and analyze the associations between these mutations and reproductive traits. Results showed that a 24-bp indel was uncovered three genotypes (II, ID and DD) in the Australian White sheep (AuW) and Lanzhou fat-tail sheep (LZFT) population, while there were only two genotypes (ID and DD) in Luxi black-headed sheep (LXBH). Moreover, the Fisher's exact test showed that the 24-bp indel mutation was significantly associated with litter size and live litter size in AuW sheep (Fisher's p < 0.05). Therefore, the 24-bp indel of sheep AHR gene can contribute to sheep marker-assisted selection breeding and further improve the sheep reproductive performance.

The detection of mutation within goat cell division cycle 25 A and its effect on kidding number.

Cell division cycle 25 A (CDC25A) accounts for an essential function on early folliculogenesis of female mammals, especially regulating the function of intra-ovarian, thus this gene is pinpointed as a candidate gene that influences the kidding number of goat. On this ground, the purpose of this study was to investigate whether the reported 20-nt nucleotide variants locus (rs639467625) of the CDC25A gene influences kidding number in Shaanbei white cashmere goat (SBWC). The χ2-test showed that there were more ID genotypes in mothers of multiple lambs than in mothers of single lambs. Interestingly, this indel locus was related to the first-born kidding number in the group of SBWC goats (p < 0.05). Similarly, the result of the t-test was consistent with the result of the χ2-test, showed the kidding number of ID genotype individuals was large than that of II individuals (p < 0.05). These findings proved that the different genotypes of CDC25A have impacts on goat kidding numbers. Thus, the results led us to speculate that the ID genotype of CDC25A was one of the main indel influencing goat kidding numbers. Simultaneously, this study was expected to provide useful DNA markers for superior individuals selection by marker-assisted selection (MAS) and make a contribution to goats breeding.

A 7-nt nucleotide sequence variant within the sheep KDM3B gene affects female reproduction traits.

Lysine demethylase 3B (KDM3B) gene is a histone demethylase, demonstrating specific demethylation of the histone H3 lysine 9. It was detected as a sheep reproductive candidate gene by genome-wide scans, and related studies also showed its significance in female reproductive process. However, rare study researched its polymorphism. Herein, we hypothesized that the polymorphisms of KDM3B gene were associated with sheep reproduction traits. A 7-nt nucleotide sequence variant (rs1088697156) within KDM3B gene was identified in a total of 888 individuals, including the Australian White (AUW) sheep and Lanzhou Fat-tailed (LFT) sheep. II (insertion/insertion) and ID (insertion/deletion) genotypes of 7-nt variant were detected, which were at Hardy-Weinberg equilibrium (HWE) in detected breeds. Association analysis illustrated the 7-nt variant was significantly associated with the litter size, duration of pregnancy, live lamb number, live lamb rate, stillbirth number, stillbirth rate of average and different parity (P < 0.05) in AUW sheep. Moreover, 'ID' was the dominant genotype with excellent consistency in reproductive traits. It is instrumental to select individuals with ID genotype for improving the sheep reproduction traits. These findings suggest that the 7-nt variant within KDM3B gene can be used as a candidate marker of reproduction traits for sheep breeding improvement by marker-assisted selection.