Genetic Variations and mRNA Expression of Goat DNAH1 and Their Associations with Litter Size

Unique feather color characteristics and transcriptome analysis of hair follicles in Liancheng White ducks

Avian feather color is a fascinating trait, and the genetic mechanism of duck plumage formation is still in the preliminary stage. In this study, feather color of Liancheng White ducks was analyzed by determination of melanin content and RNA-seq analysis. In this research, 9 ducks from Mallards (n = 3), Liancheng White (n = 3) and Pekin ducks (n = 3) were used by high performance liquid chromatography (HPLC) and Masson-Fontana staining to reveal the difference of feather melanin content. RNA-seq from 11 hair follicle tissues (1- and 8-wk-old) of Liancheng White ducks (n = 5) and Pekin ducks (n = 7) was used to analyze the candidate genes for the feather melanin synthesis, and Immunofluorescence experiment was used to show the protein expression in 6 black- and white-feathered ducks. Pectorale, skin, liver, fat, brain, heart, kidney, lung, spleen of an 8-wk-old black-feathered Mallard were collected for candidate gene expression. The results showed that the contents of feathers, beak, web melanin in Liancheng White ducks were higher than in Pekin ducks (p < 0.05). Melanin within hair follicles was located in the barb ridge and hair matrix of black feather duck, also we found that TYRP1, TYR, SOX10 genes were differentially expressed between Liancheng White and Pekin ducks (p < 0.05), and these genes were mainly expressed showed in duck skin tissues. This study revealed the unique feather color phenotype of Liancheng White duck shedding light on the transcriptome that underlies it.

Analysis of goat PPP6C mRNA profile, detection of genetic variations, and their associations with litter size

The Protein Phosphatase 6 Catalytic Subunit (PPP6C) is evolutionarily a conserved gene in eukaryotes known to play a significant role in mammalian reproduction. This study aimed to investigate expression patterns of PPP6C and explore its association with litter size in Shaanbei white cashmere (SBWC) goats. Initially, we determined the mRNA expression levels of PPP6C in both male and female goats across multiple tissues. The results showed that PPP6C mRNA was expressed in multiple tissues, with higher levels in the testis and fallopian tubes, suggesting its involvement in goat reproduction. Additionally, we identified a novel 19 bp InDel within the PPP6C gene in a population of 1030 SBWC goats, which exhibited polymorphism. Statistical analysis revealed a significant association between the19 bp InDel mutation and litter size (P < 0.05). Subsequent, bioinformatics analysis, including linkage disequilibrium (LD) block and selective scanning, highlighted the linkage tendency among most InDel loci did not stand out within B-8 block, there were still some InDel loci linked to the 19 bp within a relatively narrow region. Furthermore, comparative analysis with Bezoars, these selective signals all indicated that this gene was under higher selection pressure, implying that the 19 bp InDel locus within the PPP6C is potentially associated with domesticated traits, particularly in relation to litter size. The results of the present study suggest that the PPP6C is a vital candidate gene affecting prolificacy in goats, with implications for selective breeding programs for goat breeds.

Genome-Wide DNA Methylation Analysis and Functional Validation of Litter Size Traits in Jining Grey Goats

DNA methylation (DNAm) is associated with the reproductive system. However, the genetic mechanism through which DNAm regulates gene expression and thus affects litter size in goats is unclear. Therefore, in the present work, genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues were comprehensively analyzed via WGBS, and RNA-Seq data were combined to identify candidate genes associated with litter size traits in the Jining Grey goat. Finally, BSP and RT-qPCR were used to verify the sequencing results of the key genes. Notably, the DNMT genes were downregulated at the expression level in the HP group. Both groups exhibited comparable levels of methylation. A total of 976 differentially methylated regions (DMRs) (973 DMRs for CG and 3 DMRs for CHG) and 310 differentially methylated genes (DMGs) were identified in this study. Through integration of WGBS and RNA-Seq data, we identified 59 differentially methylated and differentially expressed genes (DEGs) and ultimately screened 8 key DMGs (9 DMRS) associated with litter size traits in Jining Grey goats (SERPINB2: chr24_62258801_62259000, NDRG4: chr18_27599201_27599400, CFAP43: chr26_27046601_27046800, LRP1B. chr2_79720201_79720400, EPHA6: chr1_40088601_40088800, TTC29: chr17_59385801_59386000, PDE11A: chr2_117418601_117418800 and PGF: chr10_ 16913801_16914000 and chr10_16916401_16916600). In summary, our research comprehensively analyzed the genome-wide DNAm profiles of HP and LP Jining Grey goat ovary tissues. The data findings suggest that DNAm in goat ovaries may play an important role in determining litter size.

ALAS1 associated with goat kidding number trait was regulated by the transcription factor ASCL2 to affect granulosa cell proliferation

ALAS1 is a member of the α-oxoamine synthase family, which is the first rate-limiting enzyme for heme synthesis and is important for maintaining intracellular heme levels. In the ovary, ALAS1 is associated with the regulation of ovulation-related mitochondrial P450 cytochromes, steroid metabolism, and steroid hormone production. However, there are few studies on the relationship between ALAS1 and reproductive traits in goats. In this study, a mutation located in the promoter region of ALAS1 (g.48791372C>A) was found to be significantly (p < 0.05) associated with the kidding number of Yunshang black goats. Specifically, the mean kidding number in the first three litters and the kidding numbers of all three litters were significantly (p < 0.05) higher in individuals with the CA genotype or AA genotype than in those with the CC genotype. To further investigate the regulatory mechanism of ALAS1, the expression of ALAS1 in goat ovarian tissues with different genotypes was verified by real-time quantitative PCR. The results showed that the expression of ALAS1 was significantly higher in the ovaries of individuals with AA genotype than those with AC and CC genotypes (p < 0.01), and the expression trend of transcription factor ASCL2 was consistent with ALAS1. Additionally, the ALAS1 g.48791372C>A mutation created a new binding site for the transcription factor ASCL2. The luciferase activity assay indicated that the mutation increased the promoter activity of ALAS1. Overexpression of the transcription factor ASCL2 induced increased expression of ALAS1 in goat granulosa cells (p < 0.05). The opposite trend was shown for the inhibition of ASCL2 expression. The results of real-time quantitative PCR, EdU and Cell Counting Kit-8 assays indicated that the transcription factor ASCL2 increased the proliferation of goat granulosa cells by mediating the expression of ALAS1. In conclusion, the transcription factor ASCL2 positively regulated the transcriptional activity and expression levels of ALAS1, altering granulosa cell proliferation and the kidding number in goats.

Goat MyoD1: mRNA expression, InDel and CNV detection and their associations with growth traits

The Myogenic differentiation 1 (MyoD1) gene is a crucial regulator of muscle formation and differentiation. However, there are few studies on the mRNA expression pattern of the goat MyoD1 gene and its effect on goat growth and development. To address this, we investigated the mRNA expression of the MyoD1 gene in several tissues of fetal and adult goats, containing heart, liver, spleen, lung, kidney and skeletal muscle. The results focused on the expression of the MyoD1 gene in skeletal muscle of fetal goats was much higher than adult goats, suggesting its important role in skeletal muscle formation and development. Following, a total of 619 Shaanbei White Cashmere goats (SBWCs) were used to monitor the InDel (Insertion/Deletion) and CNV (Copy Number Variation) variations of the MyoD1 gene. Three InDel loci were identified, and there was no significant correlation with goat growth traits. Furthermore, a CNV locus containing the MyoD1 gene exon with three types (Loss type, Normal type, Gain type) were identified. The association analysis results showed that the CNV locus was significantly associated with body weight, height at hip cross, heart girth and hip width in SBWCs (P < 0.05). Meanwhile, the Gain type of CNV exhibited the best growth traits and good consistency among three types in goats, suggesting its potential as a DNA marker for marker-assisted breeding of goats. Overall, our study provided a scientific basis for breeding goats with better growth and development traits.

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.

Association analysis of PLIN2 gene polymorphisms and lambing performance in Qianbei Ma goats

Qianbei Ma goats are one of the three most valued local goat breeds in Guizhou, China; furthermore, it has lower litter size performance. The purpose of this study was to explore the correlation between SNP (single-nucleotide polymorphisms) of PLIN2 gene and lambing performance. The bioinformatics analysis, DNA sequencing, RT-qPCR and correlation analysis methods were used to analyse the evolutionary relationship of PLIN2 protein in 13 species, to detect the expression pattern of PLIN2 gene in the gonad axis of Qianbei sheep, to explore the dominant genotype of PLIN2 related to lambing traits and to screen molecular markers related to lambing performance to guide the breeding of Qianbei Ma goats. Results showed that the Qianbei Ma goat PLIN2 protein had the closest genetic relationship with sheep and the furthest from mice, there were significant or extremely significant differences in the expression levels of the PLIN2 gene in the gonadal axis of the mothers of single- and multi-lamb groups. Compared with the reference sequence, four SNPs were found, which were g.1006 C → A and g.1171 A → G in the first and second intron regions of the PLIN2 gene, g.8514 C → T in the exon 8 region and g.9122 A → T in the 3'UTR. The correlation analysis showed that g.1006 C → A, g.8514 C → T and g.9122 A → T had significant indigenous effects on the lambing performance of Qianbei Ma goats (p < .05). The number of third births for diploid H2H5 was significantly higher than that of diploid H1H2, and the number of first to third births for diploid H2H5 was large and stable. The results showed that PLIN2 gene could be used as a candidate gene related to lambing traits of Qianbei Ma goat.