Two Insertion/Deletion Variants within SPAG17 Gene Are Associated with Goat Body Measurement Traits

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Detection of small sequence variations within the goat GHR gene and its effects on growth traits

Growth hormone receptor (GHR) gene is considered to be an important candidate gene in growth traits. Therefore, the purpose of this study was to detect whether there were potential indel variations in the GHR gene that were related to the growth traits of the Shaanbei white cashmere goats (SBWC). In this study, genomic DNA from 931 healthy SBWC individuals were used to verify the relationship between the indel of the GHR gene and growth traits. Two indel variants, P49-bp indel in intron 1 and P1410-bp indel in 3'-UTR, were confirmed. Association analyses demonstrated that these two indel polymorphism loci were associated with the chest circumference and chest width of SBWC. Additionally, for the P49-bp and P1410-bp indel loci, the ID and II genotypes were dominant genotypes, respectively. Moreover, the genotypic distributions of these two indel loci in SBWC were significantly different from those in three other Chinese indigenous goat breeds (HNBG, GZDG and IMWC) (p < 0.05). Taken together, two indel loci (P49-bp indel and P1410-bp indel) both significantly affected the growth traits of goats. This illustrated that these two indel loci might be the potential DNA marker for use in improving the selection and breeding of goats.

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 Collapse

Key Words

• RAR related orphan receptor A (RORA)
• clock gene
• goat
• growth traits
• insertion/deletion

Collapse

MESH Headings

• Pregnancy
• Female
• Animals
• Litter Size/genetics
• Goats/physiology
• INDEL Mutation/genetics
• Phenotype

Collapse

Grants Collapse

Affiliation(s)

Qian Zhou College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Huina Hu College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuta Yang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Yuxin Kang College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianyong Lan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
Xianfeng Wu Institute of Animal Husbandry and Veterinary, Fujian Academy of Agricultural Sciences, Fuzhou, Fujian, China
Zhengang Guo College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China Animal Husbandry and Veterinary Science Institute of Bijie city, Bijie, Guizhou, China
Chuanying Pan College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China

Collapse

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.

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.

Association analysis of indels in the VISFATIN gene with five cattle breeds and their growth traits

VISFATIN is an adipose cytokine that has been proved to correlate with growth and development traits. In a previous study from our lab, two insertion/deletions (indels; including a 35-bp insertion at its intron 4 and a 6-bp deletion in intron 5) were identified within the VISFATIN gene. To validate these indels and evaluate their association with growth traits in Chinese cattle, a total of 413 samples from four Chinese indigenous breeds and 217 samples from Chinese breeds were detected. Three genotypes (WW, WI and II) at intron 4 were detected based on the 35-bp insertion (allele I) or deletion (allele W) and showed moderate polymorphism in all samples. Two genotypes (WW and WD) at intron 5 were detected based on the 6-bp deletion (allele D) or insertion (allele W) in Xianan (XN) cattle and Jinnan (JN) cattle population but showed poor polymorphisms. Association analysis illustrated that the indel at intron 4 is significantly associated with chest girth, rump length and body weight in Ji'an (JA) cattle and the indel at intron 5 can cause a significant difference in rump length in JN cattle. To our knowledge, it is the first time it has been shown that indels within the VISFATIN gene are associated with growth traits in the two Chinese indigenous cattle breeds. These findings suggest that the VISFATIN gene can be used as a molecular marker for JN and JA cattle 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.

Polymorphism of GTF2A1 gene is associated with litter size in sheep

General transcription factor IIA subunit 1 (GTF2A1) is required for transcriptional activation from most promoters in eukaryotic transcription. Previous whole-genome association analyses articles have predicted the effect of this gene on lambing in sheep. In the study, a total of nine insertion/deletion (indel) variants named L1 to L9 in the gene were selected to detect in 550 Australian White sheep (AuW) adult ewes. The polymorphisms were detected in four loci (L1, L2, L3 and L8locus) and the polymorphism information content (PIC) values were 0.270, 0.375, 0.372 and 0.314. Further, our study revealed that L1, L2 and L3 loci of GTF2A1 gene were significantly related to the first parity litter size and the polymorphism of L8 was significantly correlated with litter size in the second parity. In details, for the first parity, individuals with the II genotype of the L1 locus had higher little size than that with the ID genotype, individuals with the ID and DD genotype of the L2 locus had higher little size than that with the II genotype, and individuals with the DD genotype of the L3 locus had higher little size than that with the II genotype. All four loci do not follow to Hardy-Weinberg equilibrium, and have no linkage between them. In conclusion, the polymorphisms of GTF2A1 were confirmed and analysis results demonstrated that there are some relevance between difference genotypes and litter size, and these findings may provide new insights for accelerating sheep molecular breeding through molecular marker-assisted selection (MAS).

InDel and CNV within the AKAP13 Gene Revealing Strong Associations with Growth Traits in Goat

A-kinase-anchoring protein 13 (AKAP13) is a member of the AKAP protein family that has been found to be associated with bone formation. Thus, we investigated the AKAP13 gene as a potential candidate gene for molecular-marker-assisted selection (MAS) in breeding. Our aim was to explore genetic variations (InDel and CNV) within the AKAP13 gene of Shaanbei white cashmere (SBWC) goats and analyze their relationship with growth traits. Ultimately, we identified three InDel loci (16-bp deletion, 15-bp insertion, and 25-bp deletion) and three CNVs, and the 16-bp and 15-bp loci were significantly associated with goat body length (p < 0.05). Both the 16-bp deletion variant and the 15-bp insertion variant facilitated an increase in body length in goats. In addition to this, there was a certain superposition effect between 16-bp and 15-bp loci, although there was no linkage. Additionally, the CNV1 locus was significantly correlated with body height and body length of goats (p < 0.05), and CNV2 was significantly correlated with chest depth, chest circumference, and cannon circumference of goats (p < 0.05). Individuals with gain type showed excellent growth performance. In conclusion, the InDel and CNV loci that we have identified could possibly serve as effective molecular markers in goat breeding, which is very essential for improving efficiency and success of breeding. Moreover, our findings provide a new avenue for further research into the function of the AKAP13 gene.

Screening Indels from the whole genome to identify the candidates and their association with economic traits in several goat breeds

In the present study, the re-sequencing of our previous whole-genome sequencing (WGS) for selected individuals of Dazu-black goat (DBG) and Inner-Mongolia Cashmere goat (IMCG) breeds was used to detect and compare the differentiation in Indels depending on the reference genome of goat. Then, three selected candidate Indels rs668795676, rs657996810, and rs669452874 of the three genes SUFU, SYCP2L and GLIPR1L1, respectively, have been chosen, based on the results of prior GWAS across the genome, and examined for their association with body weight and dimensions (body height, body length, heart girth, chest width, cannon circumference, and chest depth) by kompetitive allele specific PCR assay for 342 goats from the three studied goat breeds (DBG, n = 203; ♂99, ♀104), IMCG (n = 65; 15♂, 50♀), and Hechuan white goat (HWG, n = 74; 34♂, 40♀) breeds. The analysis of 192.747 Gb WGS revealed an average 334,151 Indels in the whole genome of DBG and IMCG breeds. Chromosome 1 had a maximum number of mutations (Indels) of 58,497 and 55,527 for IMCG and DBG, respectively, while chromosome 25 had the least number of mutations of 15,680 and 16,103 for IMCG and DBG, respectively. The majority of Indels were either Ins or Del of short fragments of 1-5 bp, which covered 79.06 and 71.78% of the total number of Indels mutations in IMCG and DBG, respectively. Comparing the differences of Indels between the studied goat breeds revealed 100 and 110 unique Indels for IMCG and DBG, respectively. The Indels loci in the intron region were unique for both studied goat breeds which were related to 30 and 38 candidate genes in IMCG and DBG, respectively, including SUFU, SYCP2L, and GLIPR1L1 genes. Concerning rs669452874 locus, body height and body length of Del/Del genotype in DBG were significantly higher (P < 0.05) than that of Ins/Del genotype, while body height and body length of Del/Del genotype in IMCG were significantly higher (P < 0.01) than those of Ins/Ins genotype, whereas body height and body length and heart girth of Del/Del genotype in HWG were significantly higher (P < 0.01) than those of the Ins/Del and Ins/Ins genotypes. Thus, Del/Del genotype of rs669452874 locus can be used as a candidate molecular marker related to the body dimensions in the studied goat breeds.

Effects of Genetic Variation of the Sorting Nexin 29 (SNX29) Gene on Growth Traits of Xiangdong Black Goat

Previous studies have found that the copy number variation (CNV) and insertion/deletion (indels) located in the sorting nexin 29 (SNX29) gene, which is an important candidate gene related to meat production and quality, are associated with growth traits of African goats and Shaanbei white cashmere goats. However, the genetic effects of SNX29 genetic variation on growth traits of Xiangdong black (XDB) goat (a representative meat goat breed in China) are still unclear. The purpose of this study was to detect the mRNA expression level of SNX29 and to explore the genetic effects of CNV and indel within SNX29 on growth traits and gene expression in XDB goat. The SNX29 mRNA expression profile showed that the SNX29 was highly expressed in adipose tissues, indicating that the SNX29 gene could play a key role in subcutaneous adipose deposition of XDB goat. 17 bp indel (g.10559298-10559314), 21 bp indel (g.10918982-10919002) and CNV were detected in 516 individuals of XDB goat by PCR or qPCR. The association analysis of SNX29 CNV with growth traits in XDB goats showed that SNX29 CNV was significantly correlated with chest circumference and abdominal circumference (p < 0.01), and the normal type of SNX29 CNV goat individuals were more advantageous. For the mRNA expression of SNX29 gene, individuals with SNX29 copy number normal type had a higher trend than that of SNX29 gene with copy number gain type in longissimus dorsi muscle (p = 0.07), whereas individuals with SNX29 copy number gain type had a higher trend in abdominal adipose (p = 0.09). Overall, these results suggested that the SNX29 gene could play an important role in growth and development of XDB goats and could be used for marker-assisted selection (MAS) in XDB goats.

Relationships between the Mini-InDel Variants within the Goat CFAP43 Gene and Body Traits

The cilia- and flagella-associated protein 43 (CFAP43) gene encodes a member of the cilia- and flagellum-associated protein family. Cilia on the cell surface influence intercellular signaling and are involved in biological processes such as osteogenesis and energy metabolism in animals. Previous studies have shown that insertion/deletion (InDel) variants in the CFAP43 gene affect litter size in Shaanbei white cashmere (SBWC) goats, and that litter size and body traits are correlated in this breed. Therefore, we hypothesized that there is a significant relationship between InDel variants within the CFAP43 gene and body traits in SBWC goats. Herein, we first investigated the association between three InDel variant loci (L-13, L-16, and L-19 loci) within CFAP43 and body traits in SBWC goats (n = 1827). Analyses revealed that the L-13, L-16, and L-19 loci were significantly associated with chest depth, four body traits, and three body traits, respectively. The results of this study are in good agreement with those previously reported and could provide useful molecular markers for the selection and breeding of goats for body traits.

Goat SNX29: mRNA expression, InDel and CNV detection, and their associations with litter size

The sorting nexin 29 (SNX29) gene, a member of the SNX family, is associated with material transport and lipid metabolism. Previous studies have shown that lipid metabolism affects reproductive function in animals. Thus, we hypothesized there is a correlation between the SNX29 gene and reproductive trait. To date, studies on the relationship between the SNX29 gene and reproductive traits are limited. Therefore, the purpose of this study was to examine the polymorphism in the SNX29 gene and its correlation with litter size. Herein, the mRNA expression levels of SNX29 were assayed in various goat tissue. Surprisingly, we found that SNX29 was highly expressed in the corpus luteum, large and small follicles. This result led us to suggest that the SNX29 gene has a critical role in reproduction. We further detected potential polymorphisms in Shaanbei white cashmere (SBWC) goats, including insertion/deletion (InDel, n = 2,057) and copy number variation (CNV, n = 1,402), which were related to fertility. The 17 bp deletion (n = 1004) and the 20 bp deletion (n = 1,053) within the SNX29 gene were discovered to be significantly associated with litter size (P < 0.05), and individuals the ID genotype of P1-Del-17 bp and the DD genotype of P2-Del-20bp had larger litter size. Additionally, the four CNV loci had significant correlations with litter size (P < 0.01) in our detected population. In CNV5, individuals with the median genotype were superior compared to those with loss or gain genotype in term of litter size, and in other three CNVs showed better reproductive trait in the gain genotype. Briefly, these findings suggest that SNX29 could be used as a candidate gene for litter size in goat breeding through marker-assisted selection (MAS).

KAT2B Gene Polymorphisms Are Associated with Body Measure Traits in Four Chinese Cattle Breeds

Simple Summary

Genetic improvement is one of the most important keys to overcoming the shortcomings of beef production. Identifying molecular markers related to growth characteristics and meat quality is significant in improving beef cattle breeds. Studies have shown that KAT2B, a transcriptional co-activator regulating the acetylation modification of histones, may be involved in the development and metabolism of muscle and adipose. However, there are no reports on investigating KAT2B genetic variation in Chinese native cattle. Firstly, this manuscript reports the initial bioinformatics analysis of KAT2B, finding that KAT2B protein is highly conserved among ruminants. The KAT2B gene expression profile in Qinchuan cattle was characterized, showing the spatiotemporal specificity of KAT2B gene expression in tissues such as adipose and liver during fetal and adult periods. Then, the investigation of KAT2B gene polymorphisms was carried out. Three SNPs of the KAT2B gene were identified and were found to be correlated with multiple body measurements in Fu, Qinchuan, Yak, and Chaidam cattle. These findings suggest that these three SNPs of KAT2B can serve as the molecular markers to select individuals for beef cattle breed improvement.

Abstract

Identifying molecular markers related to growth characteristics or meat quality is significant for improving beef cattle breeds. K(lysine) acetyltransferase 2B (KAT2B) is a transcriptional co-activator regulating the acetylation modification of histones, which may be involved in the development and metabolism of muscle and adipose. However, investigations of KAT2B genetic variations in Chinese native cattle are still limited. This study aimed to identify crucial single nucleotide polymorphisms (SNPs) influencing the body measurements of Chinese native cattle. Biological evolution and conservation analysis showed that KAT2B was highly conserved among the ruminants. By qPCR assay, KAT2B gene expression was found to be spatiotemporally specific in bovine tissues such as adipose and liver. By the RFLP-PCR method, three SNPs of KAT2B (g.T61908C, g.T62131C, and g.C73406T) were identified in 827 individuals of four Chinese cattle breeds, including Qinchuan (n = 658), Fu (n = 52), Yak (n = 48), and Chaidam (n = 69) cattle. Association analysis between these KAT2B polymorphisms and the body measurements of Chinese native cattle revealed significant observations. The genetic effects of g.T61908C, g.T62131C, and g.C73406T on the associated phenotypes were illustrated in each breed. In Qinchuan cattle, g.T62131C was significantly associated with better body height, chest width, hip width, and withers height, for which TC and/or TT were the advantageous genotype. In Fu cattle, TT genotype of g.T61908C was associated to better body length, while individuals with TT or CC of g.T62131C showed higher circumference of cannon bone than those with TC genotype. In Yak, individuals with TT genotype of g.C73406T had heavier body weight. In Chaidam cattle, TC genotype of g.C73406T was associated to superior body weight, while CC genotype of g.C73406T was associated to superior chest girth and circumference of cannon bone. These findings suggest that KAT2B gene polymorphisms can be used as the molecular markers for the early molecular marker-assisted selection in beef cattle breeding programs.

Detection of InDel and CNV of SPAG17 gene and their associations with bovine growth traits

Sperm-associated antigen 17 (SPAG17) gene encodes a central pair protein, which is involved in flagellar motility, male fertility and skeletal growth in ruminants. The insertions/deletions (indels) and copy number variations (CNVs) influence phenotypic traits by altering the sequences and copy numbers of functional genes, respectively. This study identified a novel 8-bp indel of SPAG17 gene in 1520 individuals from eight different cattle breeds, as well as a novel CNV region in 355 animals. The correlation analysis of indel showed that the individuals of ID genotype had superior performance traits such as body height (p = 0.038) and body slanting length (p = 0.041) as compared to other genotypes in Xianan cattle. For the CNV, different copy numbers were closely related to the body height in Qinchuan (p = 0.045) and body weight in Xianan (p = 0.036) breeds. Importantly, significant difference was observed between the 8-bp indel and the copy number loss in Xianan breed (p < 0.01). These findings indicated that the variations within the bovine SPAG17 gene can be considered as an effective DNA molecular marker for beef cattle 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 χ² 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.

Genetic variations of bovine PCOS-related DENND1A gene identified in GWAS significantly affect female reproductive traits

Genome-wide association studies (GWAS) have identified DENND1A as a potential candidate gene linked to the fertility-related phenotypes in dairy cows. However, to date, no studies have examined the association of the DENND1A insertion/deletions (indels) to bovine fertility on a large scale. Herein, two indel sites, including P4-del-26-bp and P8-ins-15-bp were identified in 1064 Holstein cows. The values of the minor allelic frequency (MAF) ranged between 0.471 (deletion) and 0.230 (deletion), respectively, and combined four different haplotypes by analyzing the haplotype combination. It is noteworthy that P4-del-26-bp is associated with the ovarian width (P = 0.0004) and corpus luteum diameter (P = 0.004). Meanwhile, P8-ins-15-bp was found to have a significant association with the ovarian width (P = 0.020), ovarian weight (P = 0.004), the number of mature follicles (P = 0.020), and diameter of the mature follicles (P = 0.016). Furthermore, the combinatorial analysis showed that the two indel combined-genotypes were significantly related to several reproductive traits (ovarian width, ovarian weight, etc.). Collectively, our findings indicated that these two novel indels and their combinations are correlated with the reproductive traits, and hence, they can serve in the marker-assisted selection (MAS) in cattle breeding. Nevertheless, further functional experiments are needed for understanding the mechanisms of these indels in cattle reproduction in a better way.

Exploration of Genetic Variants within the Goat A-Kinase Anchoring Protein 12 (AKAP12) Gene and Their Effects on Growth Traits

Simple Summary

AKAP12, the family of A-kinase anchoring proteins (AKAPs), plays an important role in the regulation of growth and development. There have been no corresponding studies of the effect of the AKAP12 gene on growth traits in goats. In our previous study, 7 bp (intron 3) and 13 bp (3′UTR) indels within the AKAP12 gene significantly influenced AKAP12 gene expression. This study expected to identify the association between these two genetic variations and growth-related traits in 1405 Shaanbei white cashmere (SBWC) goats. The P1–7 bp indel locus was significantly correlated with height at hip cross (HHC; p < 0.05) and the P2–13 bp indel locus was associated with body weight, body length, chest depth, chest width, hip width, chest circumference and cannon (bone) circumference in SBWC goats (p < 0.05). These results prove that the AKAP12 gene plays an important role in the growth and development of goats.

Abstract

The A-kinase anchoring protein 12 gene (AKAP12) is a scaffold protein, which can target multiple signal transduction effectors, can promote mitosis and cytokinesis and plays an important role in the regulation of growth and development. In our previous study, P1–7 bp (intron 3) and P2–13 bp (3′UTR) indels within the AKAP12 gene significantly influenced AKAP12 gene expression. Therefore, this study aimed to identify the association between these two genetic variations and growth-related traits in Shaanbei white cashmere goats (SBWC) (n = 1405). Herein, we identified two non-linkage insertions/deletions (indels). Notably, we found that the P1–7 bp indel mutation was related to the height at hip cross (HHC; p < 0.05) and the P2–13 bp indel was associated with body weight, body length, chest depth, chest width, hip width, chest circumference and cannon (bone) circumference in SBWC goats (p < 0.05). Overall, the two indels’ mutations of AKAP12 affected growth traits in goats. Compared to the P1–7 bp indel, the P2–13 bp indel is more suitable for the breeding of goat growth traits.

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.

Integrating genetics, physiology and morphology to study desert adaptation in a lizard species

Integration of multiple approaches is key to understand the evolutionary processes of local adaptation and speciation. Reptiles have successfully colonized desert environments, that is, extreme and arid conditions that constitute a strong selective pressure on organisms. Here, we studied genomic, physiological and morphological variations of the lizard Liolaemus fuscus to detect adaptations to the Atacama Desert. By comparing populations of L. fuscus inhabiting the Atacama Desert with populations from the Mediterranean forests from central Chile, we aimed at characterizing features related to desert adaptation. We combined ddRAD sequencing with physiological (evaporative water loss, metabolic rate and selected temperature) and morphological (linear and geometric morphometrics) measurements. We integrated the genomic and phenotypic data using redundancy analyses. Results showed strong genetic divergence, along with a high number of fixed loci between desert and forest populations. Analyses detected 110 fixed and 30 outlier loci located within genes, from which 43 were in coding regions, and 12 presented non-synonymous mutations. The candidate genes were associated with cellular membrane and development. Desert lizards presented lower evaporative water loss than those from the forest. Morphological data showed that desert lizards had smaller body size, different allometry, larger eyeballs and more dorsoventrally compressed heads. Our results suggest incipient speciation between desert and forest populations. The adaptive signal must be cautiously interpreted since genetic drift could also contribute to the divergence pattern. Nonetheless, we propose water and resource availability, and changes in habitat structure, as the most relevant challenges for desert reptiles. This study provides insights of the mechanisms that allow speciation as well as desert adaptation in reptiles at multiple levels, and highlights the benefit of integrating independent evidence.

An insertion/deletion within the CREB1 gene identified using the RNA-sequencing is associated with sheep body morphometric traits

In a previous study, the cyclic AMP response element-binding protein 1 (CREB1) gene, which is likely involved in the regulation of fat metabolism in sheep adipose tissue, was identified using RNA sequencing. CREB1 is a transcription factor that participates in the regulation of cell proliferation, differentiation, and survival as well as energy metabolism. Therefore, based on preliminary studies, this study aimed to reveal the correlation between the insertion/deletion (indel) polymorphism of the CREB1 gene and sheep growth traits. One insertion variation of the ovine CREB1 gene, C3-ins-26 bp, was investigated in 1847 Chinese and Mongolian sheep breeds. The minor allele frequencies in the CREB1 gene varied from 0.021 to 0.938. Further, statistical analyses indicated that the C3-ins-26 bp indel in the CREB1 gene was significantly related to various body measurements (body length, height, and index; chest width, depth, and width index; cannon circumference index; and height at the hip cross) in a Tan sheep population (p < 0.05). Collectively, these findings may provide important insights into marker-assisted selection of sheep.

Identification of three new microsatellites and their effects on body measurement traits in pigs using time of flight-mass spectrometry (TOF-MS)

The body status of livestock affects their physiological function and productive performances. Microsatellites, one of the most used DNA markers, have been found to be associated with pig productive traits. However, their identifications and effects on body measurement traits of the Chinese Qinghai Bamei pig still uncovered. According to our previous sequencing data, in this study, three novel microsatellites were found in this breed. Using time of flight-mass spectrometry (TOF-MS) method, these microsatellites were further identified in a large Bamei pig population. TOF-MS spectra showed that there are three microsatellites loci, named P1, P2 and P3. These microsatellites were linkage equilibrium based on the values of D' and r² tests. Association results demonstrated that P1 locus was associated with the body length, body height and chest width and the beneficial genotype was 150-/150-bp (p < 0.05); and P2 locus was associated with the body height (p < 0.05), and the 145-/145-bp, 145-/147-bp and 145-/149-bp were claimed as favorable genotypes and 145-bp allele was considered as the favorable allele. These findings suggested that P1 and P2 microsatellites might be considered as the candidate genetic markers to select pigs with superior body sizes, especially in local breed.

A novel hypomorphic allele of Spag17 causes primary ciliary dyskinesia phenotypes in mice

Primary ciliary dyskinesia (PCD) is a human condition of dysfunctional motile cilia characterized by recurrent lung infection, infertility, organ laterality defects and partially penetrant hydrocephalus. We recovered a mouse mutant from a forward genetic screen that developed many of the hallmark phenotypes of PCD. Whole-exome sequencing identified this primary ciliary dyskinesia only (Pcdo) allele to be a nonsense mutation (c.5236A>T) in the Spag17 coding sequence creating a premature stop codon (K1746*). The Pcdo variant abolished several isoforms of SPAG17 in the Pcdo mutant testis but not in the brain. Our data indicate differential requirements for SPAG17 in different types of motile cilia. SPAG17 is essential for proper development of the sperm flagellum and is required for either development or stability of the C1 microtubule structure within the central pair apparatus of the respiratory motile cilia, but not the brain ependymal cilia. We identified changes in ependymal ciliary beating frequency, but these did not appear to alter lateral ventricle cerebrospinal fluid flow. Aqueductal stenosis resulted in significantly slower and abnormally directed cerebrospinal fluid flow, and we suggest that this is the root cause of the hydrocephalus. The Spag17Pcdo homozygous mutant mice are generally viable to adulthood but have a significantly shortened lifespan, with chronic morbidity. Our data indicate that the c.5236A>T Pcdo variant is a hypomorphic allele of Spag17 that causes phenotypes related to motile, but not primary, cilia. Spag17Pcdo is a useful new model for elucidating the molecular mechanisms underlying central pair PCD pathogenesis in the mouse.This article has an associated First Person interview with the first author of the paper.

A 17-bp InDel (rs668420586) within goat CHCHD7 gene located in growth-related QTL affecting body measurement traits

The Coiled-Coil-Helix-Coiled-Coil-Helix Domain Containing 7 (CHCHD7) gene was located in a growth-related major QTL that participated in the process of bone cells metabolism in animals by regulating mitochondrial copper homeostasis and cytochrome C oxidase assembly. Therefore, we speculated that CHCHD7 gene might be involved in animal growth and body size. Herein, we discovered a 17-bp insertion/deletion (indel) within the goat CHCHD7 gene. Then, we detected this variation in Shaanbei White Cashmere (SBWC, n = 1055) goats and Inner Mongolia White Cashmere (IMWC, n = 743) goats (Alathai type) using the mathematical expectation (ME) method. We then analyzed the correlation between these genotypes and goat body measurement traits. The results showed that the minor allelic frequency (MAF) was 0.011 in SBWC, and 0.048 in IMWC. In SBWC and IMWC, the reaction time by ME method was reduced by 36.78% and 27.59%, respectively, compared to the traditional method of screening samples one by one. Moreover, in SBWC goats, the 17-bp indel was significantly associated with body measurement traits (e.g. body height, and body length) in adults. In IMWC goats, the 17-bp indel was correlated with body measurement traits (e.g. body height) in weaners. In SBWC and IMWC goat populations, the body measurement traits of the individuals homozygous for 17-bp indel were higher than those in heterozygous individuals, except for the case of cannon circumference in IMWC weaners. These findings showed that the 17-bp insertion mutation within the goat CHCHD7 gene significantly affected body morphometric traits, and could provide a basis for marker-assisted selection (MAS) breeding of cashmere goats.

Detecting novel Indel variants within the GHR gene and their associations with growth traits in Luxi Blackhead sheep

The growth hormone is important in the regulation of metabolism and energy homeostasis and acts through a growth hormone receptor (GHR). In this work, genetic variations within the ovine GHR gene were identified and tested for associations with body morphometric traits in Chinese Luxi Blackhead (LXBH) sheep. Novel deletion loci in the LXBH GHR gene included P2-del-23 bp and P8-del-23 bp indel variants. The polymorphic information content (PIC) was 0.329 in P2-del-23 bp and 0.257 in P8-del-23 bp. Moreover, both indel polymorphisms were not at Hardy-Weinberg equilibrium (p < 0.05) in the LXBH population. Statistical analyses revealed that the P2-del-23 bp and P8-del-23 bp indels were significantly associated (p < 0.05) with several growth traits in rams and ewes, including body weight, body height, chest depth, chest width, chest circumference, cannon circumference, paunch girth and hip width. Among the tested sheep, the body traits of those with genotype DD were superior to those with II and ID genotypes, suggesting that the 'D' allele was responsible for the positive effects on growth traits. Thus, these results indicate that the P2-del-23 bp and P8-del-23 bp indel sites and the DD genotype can be useful in marker-assisted selection in sheep.

Detection of a 4 bp Mutation in the 3'UTR Region of Goat Sox9 Gene and Its Effect on the Growth Traits

Simple Summary

The sex determining region Y (SRY)-type high mobility group (HMG) box 9 (Sox9) gene is critically important in the formation and development of cartilage and is considered the “main regulator” of chondrogenesis. Additionally, a large number of studies have shown that mutations in a single allele of human Sox9 can lead to campomelic dysplasia syndrome. Therefore, the mutations of Sox9 have been the subject of increasing interest among researchers. However, no studies to date have examined the association between Sox9 gene variants and growth traits in goats. Here, we detected a 4 bp indel in the 3′Untranslated Regions (3′UTR) region of Sox9 in Shaanbei white cashmere (SBWC) goats (n = 1109) and studied the association between this indel and growth traits. The 4 bp indel of Sox9 was significantly associated with body length, heart girth, hip width, and all body measurement indexes (p < 0.05) in SBWC goats. Thus, this deletion could be used as an effective molecular marker for maximizing the growth traits of goats in breeding programs.

Abstract

The SRY-type HMG box 9 (Sox9) gene plays an important role in chondrocyte development as well as changes in hypertrophic chondrocytes, indicating that Sox9 can regulate growth in animals. However, no studies to date have examined the correlation between variations in Sox9 and growth traits in goats. Here, we found a 4 bp indel in the 3′UTR of Sox9 and verified its association with growth traits in Shaanbei white cashmere goats (n = 1109). The frequencies of two genotypes (ID and II) were 0.397 and 0.603, respectively, and polymorphic information content (PIC) values showed that the indel had a medium PIC (PIC > 0.25). The 4 bp indel was significantly correlated with body length (p = 0.006), heart girth (p = 0.001), and hip width (p = 4.37 × 10 ⁻⁴). Notably, individuals with the ID genotype had significantly superior phenotypic traits compared with individuals bearing the II genotype. Hence, we speculated that the 4 bp indel is an important mutation affecting growth traits in goat, and may serve as an effective DNA molecular marker for marker-assisted selection in goat breeding programs.

Goat CMTM2: mRNA expression profiles of different alternative spliced variants and associations analyses with growth traits

CKLF like MARVEL transmembrane domain containing 2 (CMTM2) plays crucial roles in spermiogenesis, skeletogenous, growth, and development through PI3K/Akt and other pathways. The purpose of this study was to explore the expression profile and variation of different spliced CMTM2 gene in Shaanbei white cashmere goats, as well as to find the relationships between a CMTM2 promoter region 14 bp genetic variant and growth traits in 1366 Shaanbei white cashmere goats. In this study, we identified alternative CMTM2 splicing and detected the effects of the spliced variants on mRNA expression levels in tissues. Meanwhile, an unreported spliced variant of CMTM2 in goat was identified using in CDS cloning and RT-PCR, namely, CMTM2-AS2. Compared with the normal transcript (CMTM2-AS1), the novel variant had the higher expression level in muscle and liver tissues, indicating that it plays an effective role in growth traits. Furthermore, a 14 bp deletion was detected within CMTM2 promoter region, and the different genotypes were significantly associated with growth traits (e.g., body length, circumference of cannon bone) in the large group of 1366 individuals in Shaanbei white cashmere goats. We found that the body length of the individuals with II (n = 571) genotype had better phenotypes than those with DD (n = 118) and ID (n = 650) genotypes. These results have direct guiding significance for goat breeding in the future and provide a new idea for studying the characteristics and functions of CMTM2 gene in goats.

Screening of Deletion Variants within the Goat PRDM6 Gene and Its Effects on Growth Traits

Simple Summary

Genome-wide association studies found that the PR/SET Domain 6 (PRDM6) gene mutation was associated with bone development, bone density, and body mass index. This study found a 12 bp deletion variation within the PRDM6 gene in Shaanbei white cashmere goats in a large sample size (n = 1044). This variation polymorphism was associated with multiple goat growth traits in the yearling period, including heart girth (p = 0.027), cannon circumference (p = 0.008), chest depth (p = 2.10 × 10⁻⁵), chest width (p = 0.004), body height (p = 0.032), body length (p = 0.044), and hip-width (p = 0.014). The effects of the 12 bp variation were found to make no difference on adult goat growth traits. Taken together, these results demonstrate that the 12 bp deletion variant plays an important role in the early growth and development of goats and could be considered as a useful and effective molecular marker for goat breeding selection in the growth stage.

Abstract

By genome-wide association studies, the PRDM6 gene has been shown to affect multiple, apparently unrelated inherited traits, including bone density and body mass index. Therefore, it is considered a potentially pleiotropic gene. In this study, we identified a 12 bp deletion variant (NC_030814.1:rs651603667, g: 79985625-79985636delTTGACTGATCCA) within the PRDM6 gene in a large sample (SBWC goats; n = 1044). All goat samples were collected in Shaanxi province in July 2018. The frequency of the wt allele was higher than the frequency of the del allele, and this mutation polymorphism confirmed to be consistent with the Hardy–Weinberg equilibrium (p > 0.05). Further results showed that in a group of goats in the yearling period (18 months old, n = 567), this deletion variant of the PRDM6 gene was associated with heart girth (p = 0.027), cannon circumference (p = 0.008), chest depth (p = 2.10 × 10⁻⁵), chest width (p = 0.004), body height (p = 0.032), body length (p = 0.044) and hip-width (p = 0.014). For adult SBWC goats (36 months old, n = 477), the effects of the 12 bp variation on growth-related traits were found to make no difference. These findings show that the 12 bp deletion within the goat PRDM6 gene plays an important role in the early growth and development of goats. Using the 12 bp mutation, breeders can quickly and effectively select excellent individual goats at an early stage.

An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats

The casein alpha s1 (CSN1S1) gene encodes α-s1 casein, one of the proteins constituting milk, which affects milk performance, as well as improving the absorption of calcium and bone development in mammals. A previous study found that an 11-bp insertion/deletion (indel) of this gene strongly affected litter size in goats. However, to our knowledge, the relationships between this polymorphism and the milk performance and body measurement traits of goats have not been reported. In this paper, the previously identified indel has been recognized in three Chinese goat breeds, namely the Guanzhong dairy goat (GZDG; n = 235), Shaanbei white cashmere goat (SBWC; n = 1092), and Hainan black goat (HNBG; n = 278), and the following three genotypes have been studied for all of the breeds: insertion/insertion (II), deletion/deletion (DD), and insertion/deletion (ID). The allele frequencies analyzed signified that the frequencies of the "D" allele were higher (47.8%-65.5%), similar to the previous report, which indicates that this polymorphism is genetically stable in different goat breeds. Further analysis showed that this indel was markedly associated with milk fat content, total solids content, solids-not-fat content, freezing point depression, and acidity in GZDG (p < 0.05), and also affected different body measurement traits in all three breeds (p < 0.05). The goats with II genotypes had superior milk performance, compared with the others; however, goats with DD genotypes had better body measurement sizes. Hence, it may be necessary to select goats with an II or DD genotype, based on the desired traits, while breeding. Our study provides information on the potential impact of the 11-bp indel polymorphism of the CSN1S1 gene for improving the milk performance and body measurement traits in goats.