Novel InDels of GHR, GHRH, GHRHR and Their Association with Growth Traits in Seven Chinese Sheep Breeds

Identifying candidate genes and biological pathways in muscle development through multi-tissue transcriptome comparisons between male and female geese

Males and females have long shown disparities in body weight and height; yet, the underlying mechanisms influencing growth and development remain unclear. Male and female Zhedong White Geese (ZDW) geese have long been selected for large body size and egg production, respectively. This led to a large difference in body weight between males and females, making them a unique model for studying the effects of sex on growth and development. This study aimed to elucidate these mechanisms by comparing the transcriptomes of muscle and pituitary tissues in male and female ZDW geese to identify the critical genes responsible for the effects of sex on growth performance. Our analysis revealed 1101 differentially expressed genes (DEGs) in leg musculature (507 upregulated, 594 downregulated), 773 DEGs in breast musculature (311 upregulated, 462 downregulated), and 517 DEGs in the pituitary gland (281 upregulated, 236 downregulated) between male and female geese. These DEGs were significantly enriched in gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with endocrine metabolism (e.g., hormonal activities), muscle formation (e.g., sarcomere and myofibril), and bone formation (e.g., bone morphogenesis and cartilage formation). The upregulated genes in males were enriched in KEGG pathways involving nutrient digestion and absorption (vitamin and protein), as well as the secretion of digestive juices (gastric acid and bile). Through protein-protein interaction analyses, we also observed high-density gene networks related to muscle fiber development, calcium ion metabolism, mitochondrial respiratory chain, and bone development. Therefore, our multi-tissue transcriptome analysis provides a deeper understanding of the complex and systematic gender-driven effects on growth and development in geese. IGF1, GHRHR, and NCAPG-LCORL and pathways related to myogenesis might play vital roles in gender differences before hormones exert their effect.

Genome-wide unraveling SNP pairwise epistatic effects associated with sheep body weight

Epistatic effects are an important part of the genetic effect of complex traits in livestock. In this study, we used 218 synthetic ewes from the Xinjiang Academy of Agricultural Reclamation in China to identify interacting paired with genome-wide single nucleotide polymorphisms (SNPs) associated with birth weight, weaning weight, and one-yearling weight. We detected 2 and 66 SNP-SNP interactions of sheep birth weight and weaning weight, respectively. No significant epistatic interaction of one-year-old body weight was detected. The genetic interaction of sheep body weight is dynamic and time-dependent. Most significant interactions of weaning body weight contributed 1% or higher. In the weaning weight trait, 66 significant SNP pairs consisted of 98 single SNPs covering 23 chromosomes, 5 of which were nonsynonymous SNPs (nsSNPs), resulting in single amino acid substitution. We found that genes that interact with transcription factors (TFs) are target genes for the corresponding TFs. Four epitron networks affecting weaning weight, including subnetworks of HIVEP3 and BACH2 transcription factors, constructed using significant SNP pairs, were also analyzed and annotated. These results suggest that transcription factors may play an important role in explaining epistatic effects. It provides a new idea to study the genetic mechanism of weight developing.

Exploring the Impact of Insertion/Deletion in FTO and PLIN1 Genes on Morphometric Traits in Sheep

This study aimed to identify InDels from the FTO and PLIN1 genes and to analyze their association with morphometric traits in Hu sheep (HS), Dupor sheep (DS), and Small Tail Han sheep (STHS). The FTO and PLIN1 genes were genotyped using the insertion/deletion (InDel) method. A one-way ANOVA with SPSS 26.0 software (IBM Corp, Armonk, NY, USA) was used to assess the effect of the InDel FTO and PLIN1 genes on morphometric traits. The results revealed significant associations between certain InDels and the morphometric traits in different breeds of sheep. Specifically, FTO-2 was significantly associated with cannon circumference (CaC) in HS rams and body height (BoH) in HS ewes (p < 0.05). FTO-2 was also significantly associated with chest width (ChW), CaC, head length (HeL), and coccyx length (CoL) in the STHS breed (p < 0.05). FTO-3 showed significant associations with BoH in HS rams and BoH, back height (BaH), ChW, and chest depth (ChD) in HS ewes (p < 0.05). FTO-3 was also significantly associated with ChW in the DS and STHS breeds (p < 0.05). FTO-5 was significantly associated with body weight (BoW) in the DS breed and BoH in the STHS breed (p < 0.05). Furthermore, PLIN1 was significantly related to BoW in the DS breed and was significantly associated with CoL and forehead width (FoW) in the STHS breed (p < 0.05). In conclusion, the study suggested that InDels in the FTO and PLIN1 genes could provide practical information to improve morphometric traits in sheep breeding.

Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers

This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.

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

Rapid detection of InDel within the KDM3B gene in five sheep breeds using the mathematical expectation (ME) method

Lysine demethylase 3B (KDM3B), a candidate gene associated with bone formation and growth, and differentiation of osteoblast, might affect the animal growth traits. Herein, the insertion/deletion (InDel) of the KDM3B gene was quickly detected in 882 sheep from five breeds using the mathematical expectation (ME) method. The results showed that there were two genotypes of 7-bp variation in KDM3B, including II (insertion/insertion) and ID (insertion/deletion), and the frequency of two genotypes varied among the five sheep breeds. Association analysis results demonstrated that the 7-bp indel was significantly associated with chest depth of LFT sheep (P = 0.012), and body weight (P = 0.006), body height (P = 0.030), chest depth (P = 0.043), chest circumference (P = 0.016), abdominal width (P = 0.035) and height at hip cross (P = 0.022) in LXBH sheep. Moreover, II genotype was the predominant genotype with excellent consistency in sheep growth traits (P < 0.05). Collectively, the above results suggest that this locus can be used as an effective molecular marker to improve the sheep growth traits and provide a scientific basis for the development of sheep breeding.

A comparative investigation on H3K27ac enhancer activities in the brain and liver tissues between wild boars and domesticated pigs

Dramatic phenotypic differences between domestic pigs and wild boars (Sus scrofa) provide opportunities to investigate molecular mechanisms underlying the formation of complex traits, including morphology, physiology and behaviour. Most studies comparing domestic pigs and wild boars have focused on variations in DNA sequences and mRNA expression, but not on epigenetic changes. Here, we present a genome‐wide comparative study on H3K27ac enhancer activities and the corresponding mRNA profiling in the brain and liver tissues of adult Bama Xiang pigs (BMXs) and Chinese wild boars (CWBs). We identified a total of 1,29,487 potential regulatory elements, among which 11,241 H3K27ac peaks showed differential activity between CWBs and BMXs in at least one tissue. These peaks were overrepresented by binding motifs of FOXA1, JunB, ATF3 and BATF, and overlapped with differentially expressed genes that are involved in female mating behaviour, response to growth factors and hormones, and lipid metabolism. We also identified 4118 nonredundant super‐enhancers from ChIP‐Seq data on H3K27ac. Notably, we identified differentially active peaks located close to or within candidate genes, including TBX19, MSTN, AHR and P2RY1, which were identified in DNA sequence‐based population differentiation studies. This study generates a valuable dataset on H3K27ac profiles of the brain and liver from domestic pigs and wild boars, which helps gain insights into the changes in enhancer activities from wild boars to domestic pigs.

Polymorphisms of the Growth Hormone Releasing Hormone Receptor Gene Affect Body Conformation Traits in Chinese Dabieshan Cattle

This study was performed to expose the polymorphisms of the growth hormone-releasing hormone receptor gene in Chinese Dabieshan cattle, evaluate its effect on body conformation traits, and find potential molecular markers in Chinese cattle. The GHRHR structure and the phylogenetic tree were analyzed using bioinformatics software. The polymorphism of the GHRHR gene in 486 female cattle was genotyped by PCR-RFLP and DNA sequencing, and the association between SNPs and body conformation traits of Chinese Dabieshan cattle was analyzed by one-way ANOVA in SPSS software. GHRHR was often conserved in nine species, and its sequence of cattle was closest to sheep and goats. Six polymorphic SNPs were identified, g.10667A > C and g.10670A > C were missense mutation. The association analysis indicated that the six SNPs significantly influenced the body conformation traits of Chinese Dabieshan cattle (p < 0.05). Six haplotypes were identified and Hap1 (-CAACGA-) had the highest frequency (36.10%). The Hap3/5 (-GCCCCCGGAAGG-) exhibited a significantly greater wither height (WH), hip height (HH), heart girth (HG), and hip width (HW) (p < 0.05). Overall, the polymorphisms of GHRHR affected the body conformation traits of Chinese Dabieshan cattle, and the GHRHR gene could be used as a molecular marker in Dabieshan cattle breeding programs.

Detection of selection signatures in South African Mutton Merino sheep using whole-genome sequencing data

The South African Mutton Merino (SAMM), a dual-purpose (meat and wool) sheep breed, is characterized by its excellent performance on growth, carcass traits and meat quality compared to other fine-wool Merino breeds. Nowadays, the SAMM breed has been widely used to cross with commercial and indigenous fine-wool or coarse-wool breeds to improve the growth and meat performance in many countries. To date, however, little is known about the genetic basis for its prominent characteristics. In this study, whole-genome sequences of 10 SAMM were sequenced and the selection signatures were analyzed together with those of 39 Australian Merino and Chinese Merino (wool-type Merino) by F_ST , iHS, and XP-EHH methods. In total, 313 genes in 277 regions were identified by at least 2 methods with the signal of selection and 21 of them were identified by all three methods. We highlighted a list of interesting genes, including GHR, LCORL, SMO, NCAPG, DCC, IBSP, PPARGC1A, PACRGL, PRDM5, XYLB, AHCYL2, TEFM, AFG1L, and FAM184B, which have been shown to be involved in growth, carcass traits, and meat quality by previous studies. Herein, GHR, encoding a transmembrane receptor for growth hormone, is the most notable one. We report the first study on selection signatures analysis of SAMM at whole-genome sequence level. These results provide new insights into the genetic mechanisms underlying the growth and carcass traits in SAMM.

Polymorphisms in BMPRIB gene affect litter size in Chinese indigenous sheep breed

The BMPRIB gene belongs to the TGF-β superfamily and is considered to be a regulator of sheep reproductive performance. Single nucleotide polymorphisms (SNPs) of BMPRIB gene in the Small Tail Han, Hu, Mongolian, Oula, Gansu Alpine Fine-wool, Dorper and Australian White sheep were detected by Sanger sequencing. Five SNPs (rs427897187 G > A, rs418841713 A > G, rs159952533 T > C, rs429416173 C > A and rs403555643 A > G) of BMPRIB gene were identified. For rs427897187 G > A, further analysis revealed that genotype GG and GA had 0.26 (p < 0.05) and 0.33 (p < 0.05) litter size less than those with genotype AA in Oula sheep. For rs403555643 A > G, further analysis revealed that genotype GG and AG had 0.65 (p < 0.05) and 0.38 (p < 0.05) litter size more than those with genotype AA in Oula sheep, and genotype GG had 0.56 (p < 0.05) litter size more than those with genotype AA in Mongolian sheep. The results showed that rs427897187 G > A and rs403555643 A > G are potential molecular markers wich could improve litter size of Chinese indigenous sheep and be used in Chinese indigenous sheep breeding.