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.

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.

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 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.

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.

Detection of the 23-bp nucleotide sequence mutation in retinoid acid receptor related orphan receptor alpha (RORA) gene and its effect on sheep litter size.

Retinoid acid receptor related orphan receptor alpha (RORA) transcribes steroid-related genes to regulate estrogen synthesis. As an important reproductive trait, litter size relates to estrogen synthesis. Therefore, it is important to investigate the association between RORA gene and sheep litter size. In this study, one 23-bp nucleotide sequence mutation was identified in intron 1 of RORA gene in 532 female Australian White Sheep. Moreover, the polymorphic information content (PIC) values of this locus was 0.219. The litter size of ID genotype was significantly better than II genotype and DD genotype in the second born litter size (p < 0.05). This loci was related to third born litter size and the ID is the dominant genotype (p < 0.05). The association between combined genotypes and average litter size showed that sheep with heterozygote (ID) genotypes had larger lamb than homozygous (DD and II) genotypes. To sum, this study provided theoretical references for the comprehensively research of the function of RORA gene and the breeding of Australian White Sheep. The 23-bp indel variants could be considered as molecular markers for the second and third born litter size of sheep for MAS breeding.

Goat SMAD family member 1 (SMAD1): mRNA expression, genetic variants, and their associations with litter size.

SMAD family member 1 (SMAD1) is phosphorylated and activated by the BMP receptors, which help regulate ovulation rate, cell growth, apoptosis, and development. Previously, the genome-wide association study revealed that it has been associated with fecundity in sheep. However, its effect on litter size has not been investigated in goats. Therefore, this study aimed to determine the level of SMAD1 mRNA expression in various tissues and to identify its polymorphisms and their association with litter size in Shaanbei white cashmere goat (SBWC). As a result, RT-qPCR analysis showed that SMAD1 was expressed in various tissues in female SBWC goats, including the ovary (P < 0.05). Importantly, the mRNA expression level in the ovaries of mothers of multi-lambs had a higher level than the mothers of single lambs (P < 0.05). Moreover, two InDels (18-bp and 7-bp) in intron 1 of SMAD1 were polymorphic among ten potential loci. Both 18-bp and 7-bp InDels were significantly correlated with litter size (P = 0.014) and (P = 0.0001), respectively. As shown by the chi-squared test, genotypic distributions of 18-bp and 7-bp were significantly distinct between single-lamb (P = 0.02) and multi-lamb mothers (P = 0.002). Our findings confirm that two InDels in SMAD1 were significantly associated with litter size and suggest that they could be used to improve fertility traits in goat breeding strategies.

Detection of mRNA Expression and Copy Number Variations Within the Goat Fec B Gene Associated With Litter Size.

The Booroola fecundity (Fec B ) gene, as the first major fecundity gene identified in Booroola sheep, has attracted careful attention. So far, previous research have uncovered the FecB mutation (Q249R) as the main mutation by virtue of which sheep exhibits multiple lambing phenomena. This mutation is now being intensively studied and widely used. However, such effect of the FecB mutation has not been applied to goats, and similar types of the Fec B gene in goats still need to be studied. Thus, the current study attempted to verify potential mutations in the goat Fec B gene as well as investigate their functions related to fecundity. First, Fec B expression was investigated in six different goat tissues, and we found that Fec B expression was highest in the mammary gland, followed by the ovary. Next, the influence of the Fec B gene was analyzed from a new perspective, where five potential copy number variations (CNVs) (CNV1-5) within the Fec B gene were identified for the first time, and then their effects on litter size were measured. Our results point out that CNV3 (P = 3.44E-4) and CNV5 (P = 0.034) could significantly influence the litter size of goats. Identically, the combination genotype of CNV3 and CNV5 which consisted of their dominant genotypes was also significantly associated with goat litter size (P = 7.80E-5). Hence, CNV3 and CNV5 could serve as potential DNA molecular markers applied to DNA editing and DNA microarray. Additionally, the abovementioned study has laid a theoretical foundation for the detection of potential fertility-related quantitative trait loci within the goat Fec B gene.