Combined approaches to reveal genes associated with litter size in Yunshang black goats

Whole-genome resequencing reveals new mutations in candidate genes for Beichuan-white goat prolificacya

In this study, we evaluated the copy number variation in the genomes of two groups of Beichuan-white goat populations with large differences in litter size by FST method, and identified 1739 genes and 485 missense mutations in the genes subject to positive selection. Through functional enrichment, ITGAV, LRP4, CDH23, TPRN, RYR2 and CELSR1 genes, involved in embryonic morphogenesis, were essential for litter size trait, which received intensive attention. In addition, some mutation sites of these genes have been proposed (ITGAV: c.38C > T; TPRN: c.133A > T, c.1192A > G, c.1250A > C; CELSR1: c.7640T > C), whose allele frequencies were significantly changed in the high fecundity goat group. Besides, we found that new mutations at these sites altered the hydrophilicity and 3D structure of the protein. Candidate genes related to litter size in this study and their missense mutation sites were identified. These candidate genes are helpful to understand the genetic mechanism of fecundity in Beichuan white goat, and have important significance for future goat breeding.

Multiomic meta-analysis suggests a correlation between steroid hormone-related genes and litter size in goats

Litter size is a key indicator of production performance in livestock. However, its genetic basis in goats remains poorly understood. In this work, a genome-wide selection sweep analysis (GWSA) on 100 published goat genomes with different litter rates was performed for the first time to identify candidate genes related to kidding rate. This analysis was combined with the public RNA-sequencing data of ovary tissues (follicular phase) from high- and low-yielding goats. A total of 2278 genes were identified by GWSA. Most of these genes were enriched in signaling pathways related to ovarian follicle development and hormone secretion. Moreover, 208 differentially expressed genes between groups were obtained from the ovaries of goats with different litter sizes. These genes were substantially enriched in the cholesterol and steroid synthesis signaling pathways. Meanwhile, the weighted gene co-expression network was used to perform modular analysis of differentially expressed genes. The results showed that seven modules were reconstructed, of which one module showed a very strong correlation with litter size (r = -0.51 and p-value <0.001). There were 51 genes in this module, and 39 hub genes were screened by Pearson's correlation coefficient between core genes > 0.4, correlation coefficient between module members > 0.80 and intra-module connectivity ≥5. Finally, based on the results of GWSA and hub gene Venn analysis, seven key genes (ACSS2, HECW2, KDR, LHCGR, NAMPT, PTGFR and TFPI) were found to be associated with steroid synthesis and follicle growth development. This work contributes to understanding of the genetic basis of goat litter size and provides theoretical support for goat molecular breeding.

Association of litter size with the ruminal microbiome structure and metabolomic profile in goats

The Yunshang black goat is a renowned mutton specialist breed mainly originating from China that has excellent breeding ability with varying litter sizes. Litter size is an important factor in the economics of goat farming. However, ruminal microbiome structure might be directly or indirectly regulated by pregnancy-associated factors, including litter sizes. Therefore, the current experiment aimed to evaluate the association of different litter sizes (low versus high) with ruminal microbiome structure by 16S rRNA gene sequencing and metabolomic profiling of Yunshang black does. A total of twenty does of the Yunshang Black breed, approximately aged between 3 and 4 years, were grouped (n = 10 goats/group) into low (D-l) and high (D-h) litter groups according to their litter size (the lower group has ≤ 2 kids/litter and the high group has ≧ 3 kids/litter, respectively). All goats were sacrificed, and collected ruminal fluid samples were subjected to 16S rRNA sequencing and LC-MS/MC Analysis for ruminal microbiome and metabolomic profiling respectively. According to PCoA analysis, the ruminal microbiota was not significantly changed by the litter sizes among the groups. The Firmicutes and Bacteroidetes were the most dominant phyla, with an abundance of 55.34% and 39.62%, respectively. However, Ruminococcaceae_UCG-009, Sediminispirochaeta, and Paraprevotella were significantly increased in the D-h group, whereas Ruminococcaceae_UCG-010 and Howardella were found to be significantly decreased in the D-l group. The metabolic profiling analysis revealed that litter size impacts metabolites as 29 and 50 metabolites in positive and negative ionic modes respectively had significant differences in their regulation. From them, 16 and 24 metabolites of the D-h group were significantly down-regulated in the positive ionic mode, while 26 metabolites were up-regulated in the negative ionic mode for the same group. The most vibrant identified metabolites, including methyl linoleate, acetylursolic acid, O-desmethyl venlafaxine glucuronide, melanostatin, and arginyl-hydroxyproline, are involved in multiple biochemical processes relevant to rumen roles. The identified differential metabolites were significantly enriched in 12 different pathways including protein digestion and absorption, glycerophospholipid metabolism, regulation of lipolysis in adipocytes, and the mTOR signaling pathway. Spearman's correlation coefficient analysis indicated that metabolites and microbial communities were tightly correlated and had significant differences between the D-l and D-h groups. Based on the results, the present study provides novel insights into the regulation mechanisms of the rumen microbiota and metabolomic profiles leading to different fertility in goats, which can give breeders some enlightenments to further improve the fertility of Yunshang Black goats.

Temporal Transcriptome Dynamics of Longissimus dorsi Reveals the Mechanism of the Differences in Muscle Development and IMF Deposition between Fuqing Goats and Nubian Goats

In this study, we measured the growth performance and intramuscular fat (IMF) content of the Longissimus dorsi (LD) of Fuqing goats (FQs) and Nubian goats (NBYs), which exhibit extreme phenotypic differences in terms of their production and meat quality traits. RNA-Seq analysis was performed, and transcriptome data were obtained from the LD tissue of 3-month fetuses (E3), 0-month lambs (0M), 3-month lambs (3M), and 12-month lambs (12M) to reveal the differences in the molecular mechanisms regulating the muscle development and IMF deposition between FQs and NBYs. The results showed that a higher body weight and average daily gain were observed in the NBYs at three developmental stages after birth, whereas a higher IMF content was registered in the FQs at 12M. Additionally, transcriptome profiles during the embryonic period and after birth were completely different for both FQs and NBYs. Moreover, DEGs (KIF23, CCDC69, CCNA2, MKI67, KIF11, RACGAP1, NUSAP1, SKP2, ZBTB18, NES, LOC102180034, CAPN6, TUBA1A, LOC102178700, and PEG10) significantly enriched in the cell cycle (ko04110) at E3 (FQs vs. NBYs), and DEGs (MRPS7, RPS8, RPL6, RPL4, RPS11, RPS10, RPL5, RPS6, RPL8, RPS13, RPS24, RPS15, RPL23) significantly enriched in ribosomes (ko03010) at 0M (FQs vs. NBYs) related to myogenic differentiation and fusion were identified. Meanwhile, the differences in glucose and lipid metabolism began at the E3 timepoint and continued to strengthen as growth proceeded in FQs vs. NBYs. DEGs (CD36, ADIROQR2, ACACA, ACACB, CPT1A, IGF1R, IRS2, LDH-A, PKM, HK2, PFKP, PCK1, GPI, FASN, FADS1, ELOVL6, HADHB, ACOK1, ACAA2, and ACSL4) at 3M (FQs vs. NBYs) and 12M (FQs vs. NBYs) significantly enriched in the AMPK signaling pathway (ko04152), insulin resistance (ko04931), the insulin signaling pathway (ko04910), fatty acid metabolism (ko01212), and glycolysis/gluconeogenesis (ko00010) related to IMF deposition were identified. Further, the results from this study provide the basis for future studies on the mechanisms regulating muscle development and IMF deposition in different breeds of goats, and the candidate genes identified could be used in the selection process.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Genetic characterization of the Latvian local goat breed and genetic traits associated with somatic cell count

The Latvian local goat (LVK) breed represents the only native domestic goat breed in Latvia, but its limited population places it within the endangered category. However, the LVK breed has not yet undergone a comprehensive genetic characterization. Therefore, we completed whole genome sequencing to reveal the genetic foundation of the LVK breed while identifying genetic traits linked to the somatic cell count (SCC) levels. The study included 40 genomes of LVK goats sequenced to acquire at least 35x or 10x coverage. A Principal component analysis, a genetic distance tree, and an admixture analysis showed LVK's similarity to some European breeds, such as Finnish Landrace, Alpine, and Saanen, which aligns with the breed's history. An analysis of genome-wide heterozygosity, nucleotide diversity, and LD analysis indicated that the LVK population exhibits substantial levels of genetic diversity. LVK genome was dominated by short runs of homozygosity (ROHs, ≤ 500 kb) with a median length of 25 kb. With FROH 2.49%, average inbreeding levels were low; however, FROH ranged broadly from 0.13 to 12.2%. With the exception of one pure-blood breeding buck exhibiting FROH of 9.3% and FSNP of 8.5%, animals with at least 66% LVK ancestry showed moderate or no inbreeding. Overall, this study demonstrated that the LVK goats can be differentiated from imported breeds, although the population has a complex genetic structure. We were able to identify potential genetic traits associated with SCC levels, although the kinship of the animals and the heterogenic substructure of the population might have largely influenced the association analysis. We identified 26 genetic variants associated with SCC levels, which included the potentially relevant SNP rs662053371 in the OSBPL8 gene, indicating a potential signal linked to lipid metabolism in goats. To conclude, these findings present valuable insight into the genetic structure of the LVK breed for the conservation of local genetic resources.

Association analysis of polymorphisms in SLK, ARHGEF9, WWC2, GAB3, and FSHR genes with reproductive traits in different sheep breeds

The aim was to investigate the relationship between polymorphisms of gene mutation loci and reproductive traits in local sheep breeds (Duolang Sheep) and introduced sheep breeds (Suffolk, Hu Sheep) in Xinjiang to provide new molecular markers for the selection and breeding of high fecundity sheep. The expression pattern of typing successful genes in sheep tissues was investigated by RT-qPCR technology, providing primary data for subsequent verification of gene function. The 26 mutation loci of WWC2, ARHGEF9, SLK, GAB3, and FSHR genes were typed using KASP. Association analyses were performed using SPSS 25.0, and the typing results showed that five genes with six loci, WWC2 (g.14962207 C>T), ARHGEF9 (g.48271079 C>A), SLK (g.27107842 T>C, g.27108855 G>A), GAB3 (g.86134602 G>A), and FSHR (g.80789180 T>G) were successfully typed. The results of the association analyses showed that WWC2 (g.14962207 C>T), SLK (g.27108855 G>A), ARHGEF9 (g.48271079 C>A), and FSHR (g.80789180 T>G) caused significant or extremely significant effects on the litter size in Duolang, Suffolk and Hu Sheep populations. The expression distribution pattern of the five genes in 12 sheep reproduction-related tissues was examined by RT-qPCR. The results showed that the expression of the SLK gene in the uterus, the FSHR gene in the ovary, and the ARHGEF9 gene in hypothalamic-pituitary-gonadal axis-related tissues were significantly higher than in the tissues of other parts of the sheep. WWC2 and GAB3 genes were highly expressed both in reproductive organs and visceral tissues. In summary, the WWC2 (g.14962207 C>T), SLK (g.27108855 G>A), ARHGEF9 (g.48271079 C>A), and FSHR (g.80789180 T>G) loci can be used as potential molecular markers for detecting differences in reproductive performance in sheep. Due to variations in typing results, the SLK (g.27107842 T>C) and GAB3 (g.86134602 G>A) loci need further validation.

Identification of genomic characteristics and selective signals in Guizhou black goat

BACKGROUND

Guizhou black goat is one of the indigenous black goat breeds in the southwest region of Guizhou, China, which is an ordinary goat for mutton production. They are characterized by moderate body size, black coat, favorite meat quality with tender meat and lower odor, and tolerance for cold and crude feed. However, little is known about the genetic characteristics or variations underlying their important economic traits.

RESULTS

Here, we resequenced the whole genome of Guizhou black goat from 30 unrelated individuals breeding in the five core farms. A total of 9,835,610 SNPs were detected, and 2,178,818 SNPs were identified specifically in this breed. The population structure analysis revealed that Guizhou black goat shared a common ancestry with Shaanbei white cashmere goat (0.146), Yunshang black goat (0.103), Iran indigenous goat (0.054), and Moroccan goat (0.002). However, Guizhou black goat showed relatively higher genetic diversity and a lower level of linkage disequilibrium than the other seven goat breeds by the analysis of the nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity. Based on FST and θπ values, we identified 645, 813, and 804 selected regions between Guizhou black goat and Yunshang black goat, Iran indigenous goat, and cashmere goats. Combined with the results of XP-EHH, there were 286, 322, and 359 candidate genes, respectively. Functional annotation analysis revealed that these genes are potentially responsible for the immune response (e.g., CD28, CD274, IL1A, TLR2, and SLC25A31), humility-cold resistance (e.g., HBEGF, SOSTDC1, ARNT, COL4A1/2, and EP300), meat quality traits (e.g., CHUK, GAB2, PLAAT3, and EP300), growth (e.g., GAB2, DPYD, and CSF1), fertility (e.g., METTL15 and MEI1), and visual function (e.g., PANK2 and NMNAT2) in Guizhou black goat.

CONCLUSION

Our results indicated that Guizhou black goat had a high level of genomic diversity and a low level of linkage disequilibrium in the whole genome. Selection signatures were detected in the genomic regions that were mainly related to growth and development, meat quality, reproduction, disease resistance, and humidity-cold resistance in Guizhou black goat. These results would provide a basis for further resource protection and breeding improvement of this very local breed.

Genome-Wide Association Study of Milk Composition in Karachai Goats

This study is first to perform a genome-wide association study (GWAS) to investigate the milk quality traits in Karachai goats. The objective of the study was to identify candidate genes associated with milk composition traits based on the identification and subsequent analysis of all possible SNPs, both genome-wide (high-confidence) and suggestive (subthreshold significance). To estimate the milk components, 22 traits were determined, including several types of fatty acids. DNA was extracted from ear tissue or blood samples. A total of 167 Karachai goats were genotyped using an Illumina GoatSNP53K BeadChip panel (Illumina Inc., San Diego, CA, USA). Overall, we identified 167 highly significant and subthreshold SNPs associated with the milk components of Karachai goats. A total of 10 SNPs were located within protein-coding genes and 33 SNPs in close proximity to them (±0.2 Mb). The largest number of genome-wide significant SNPs was found on chromosomes 2 and 8 and some of them were associated with several traits. The greatest number of genome-wide significant SNPs was identified for crude protein and lactose (6), and the smallest number-only 1 SNP-for freezing point depression. No SNPs were identified for monounsaturated and polyunsaturated fatty acids. Functional annotation of all 43 SNPs allowed us to identify 66 significant candidate genes on chromosomes 1, 2, 3, 4, 5, 8, 10, 13, 16, 18, 21, 23, 25, 26, and 27. We considered these genes potential DNA markers of the fatty acid composition of Karachai goat milk. Also, we found 12 genes that had a polygenic effect: most of them were simultaneously associated with the dry matter content and fatty acids (METTL, SLC1A 8, PHACTR1, FMO2, ECI1, PGP, ABCA3, AMDHD2). Our results suggest that the genes identified in our study affecting the milk components in Karachai goats differed from those identified in other breeds of dairy goats.

Polymorphism and expression of GLUD1 in relation to reproductive performance in Jining Grey goats

Understanding the molecular mechanism of mammalian reproduction (puberty and prolificacy) will play a part in improving animal reproductive performance. GLUD1 (glutamate dehydrogenase 1) is important for mammalian reproduction, as shown in previous studies; however, its roles in puberty and prolificacy have rarely been reported. In this study, we designed seven pairs of primers (P1 to P7) for cloning and sequencing genomic DNA of Jining Grey goats and Liaoning Cashmere goats. Primer 8 (P8) was designed to detect single nucleotide polymorphism (SNP) of the GLUD1 in both sexually precocious and high-fecundity breeds (Jining Grey, Nanjiang Brown and Matou goats) and sexually late-maturing and low-fecundity breeds (Liaoning Cashmere, Inner Mongolia Cashmere and Taihang goats) by PCR-RFLP (restriction fragment length polymorphism). The real-time quantitative polymerase chain reaction (RT-qPCR) technique was used to detect the expression of GLUD1 in a variety of tissues. The results showed that the A197C mutation was only found in the amplification product of P6. For this SNP locus, only two genotypes (AA and AC) were detected in Nanjiang Brown goats, while three genotypes (AA, AC and CC) were detected in the other five breeds. In Jining Grey goats, the frequency of genotypes AA, AC and CC was 0.69, 0.26 and 0.05, respectively. In Jining Grey goats, AA genotype had 0.54 (P < 0.05 ) and 0.3 (P < 0.05 ) more kids than the CC and AC genotype, respectively, and no significant difference (P > 0.05 ) was found in kidding number between the AC and CC genotype. GLUD1 was expressed in five tissues of different developmental stages. The expression level of GLUD1 in the hypothalamus was higher than that in the other four tissues except during puberty of Liaoning Cashmere goats. In puberty in goats, GLUD1 expression was significantly higher in ovaries than that in the juvenile period (P < 0.01 ). RT-qPCR results showed that the expression of GLUD1 in ovaries may relate to the puberty of goats. The present study preliminarily indicated that there might be an association between the 197 locus of GLUD1 and sexual precocity in goats, and allele A of GLUD1 was a potential DNA marker for improving kidding number in Jining Grey goats.

Uterus proliferative period ceRNA network of Yunshang black goat reveals candidate genes on different kidding number trait

Pregnancy loss that occurs in the uterus is an important and widespread problem in humans and farm animals and is also a key factor affecting the fecundity of livestock. Understanding the differences in the fecundity of goats may be helpful in guiding the breeding of goats with high fecundity. In this study, we performed RNA sequencing and bioinformatics analysis to study the uterus of Yunshang black goats with high and low fecundity in the proliferative period. We identified mRNAs, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) by analyzing the uterine transcriptomes. The target genes of the identified miRNAs and lncRNAs were predicted, and miRNA-mRNA interaction and competitive endogenous RNA (ceRNA) networks were constructed. By comparisons between low- and high-fecundity groups, we identified 1,674 differentially expressed mRNAs (914 were upregulated, and 760 were downregulated), 288 differentially expressed lncRNAs (149 were upregulated, and 139 were downregulated), and 17 differentially expressed miRNAs (4 were upregulated, and 13 were downregulated). In addition, 49 miRNA-mRNA pairs and 45 miRNA-lncRNA pairs were predicted in the interaction networks. We successfully constructed a ceRNA interaction network with 108 edges that contained 19 miRNAs, 11 mRNAs, and 73 lncRNAs. Five candidate genes (PLEKHA7, FAT2, FN1, SYK, and ITPR2) that were annotated as cell adhesion or calcium membrane channel protein were identified. Our results provide the overall expression profiles of mRNAs, lncRNAs, and miRNAs in the goat uterus during the proliferative period and are a valuable reference for studies into the mechanisms associated with the high fecundity, which may be helpful to guide goat to reduce pregnancy loss.

Integrated Transcriptome Analysis Reveals the Crucial mRNAs and miRNAs Related to Fecundity in the Hypothalamus of Yunshang Black Goats during the Luteal Phase

A normal estrus cycle is essential for the breeding of goats, and the luteal phase accounts for most of the estrus cycle. The corpus luteum (CL) formed during the luteal phase is a transient endocrine gland that is crucial for the reproductive cycle and pregnancy maintenance, and is controlled by many regulatory factors. However, the molecular mechanism of the hypothalamus effect on the reproductive performance of different litter sizes during the luteal phase of goats has not been elucidated. In this study, RNA-sequencing was used to analyze the mRNA and miRNA expression profiles of the hypothalamic tissues with the high-fecundity goats during the luteal phase (LP-HF) and low-fecundity goats during the luteal phase (LP-LF). The RNA-seq results found that there were 1963 differentially expressed genes (DEGs) (890 up-regulated and 1073 down-regulated). The miRNA-seq identified 57 differentially expressed miRNAs (DEMs), including 11 up-regulated and 46 down-regulated, of which 199 DEGs were predicted to be potential target genes of DEMs. Meanwhile, the functional enrichment analysis identified several mRNA-miRNA pairs involved in the regulation of the hypothalamic activity, such as the common target gene MEA1 of novel-miR-972, novel-miR-125 and novel-miR-403, which can play a certain role as a related gene of the reproductive development in the hypothalamic-pituitary-gonadal (HPG) axis and its regulated network, by regulating the androgen secretion. While another target gene ADIPOR2 of the novel-miR-403, is distributed in the hypothalamus and affects the reproductive system through a central role on the HPG axis and a peripheral role in the gonadal tissue. An annotation analysis of the DE miRNA-mRNA pairs identified targets related to biological processes, such as anion binding (GO:0043168) and small molecule binding (GO: 0036094). Subsequently, the KEGG(Kyoto Encyclopedia of Genes and Genomes) pathways were performed to analyze the miRNA-mRNA pairs with negatively correlated miRNAs. We found that the GnRH signaling pathway (ko04912), the estrogen signaling pathway (ko04915), the Fc gamma R-mediated phagocytosis (ko04666), and the IL-17 signaling pathway (ko04657), etc., were directly and indirectly associated with the reproductive process. These targeting interactions may be closely related to the reproductive performance of goats. The results of this study provide a reference for further research on the molecular regulation mechanism for the high fertility in goats.

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.

Oviduct Transcriptomic Reveals the Regulation of mRNAs and lncRNAs Related to Goat Prolificacy in the Luteal Phase

Simple Summary

The kidding number is an important reproductive trait in domestic goats. The oviduct, as one of the most major organs, is directly involved in the reproductive process, providing nutrition and a location for early embryonic development. The current study provides genome-wide expression profiles of mRNA and long noncoding RNAs (lncRNAs) expression in Yunshang black goat, a new breed of meat goat bred in China with a high kidding number. During the luteal phases, oviduct mRNAs and lncRNAs associated with high- and low-fecundity Yunshang black goats were identified, and their potential biological functions were predicted using GO, KEGG, and GSEA enrichment analysis. These findings shed light on the oviduct-based prolificacy mechanism in goats.

Abstract

The oviduct is associated with embryo development and transportation and regulates the pregnancy success of mammals. Previous studies have indicated a molecular mechanism of lncRNAs in gene regulation and reproduction. However, little is known about the function of lncRNAs in the oviduct in modulating goat kidding numbers. Therefore, we combined RNA sequencing (RNA-seq) to map the expression profiles of the oviduct at the luteal phase from high- and low-fecundity goats. The results showed that 2023 differentially expressed mRNAs (DEGs) and 377 differentially expressed lncRNAs (DELs) transcripts were screened, and 2109 regulated lncRNA-mRNA pairs were identified. Subsequently, the genes related to reproduction (IGF1, FGFRL1, and CREB1) and those associated with embryonic development and maturation (DHX34, LHX6) were identified. KEGG analysis of the DEGs revealed that the GnRH- and prolactin-signaling pathways, progesterone-mediated oocyte maturation, and oocyte meiosis were related to reproduction. GSEA and KEGG analyses of the target genes of DELs demonstrated that several biological processes and pathways might interact with oviduct functions and the prolificacy of goats. Furthermore, the co-expression network analysis showed that XLOC_029185, XLOC_040647, and XLOC_090025 were the cis-regulatory elements of the DEGs MUC1, PPP1R9A, and ALDOB, respectively; these factors might be associated with the success of pregnancy and glucolipid metabolism. In addition, the GATA4, LAMA2, SLC39A5, and S100G were trans-regulated by lncRNAs, predominantly mediating oviductal transport to the embryo and energy metabolism. Our findings could pave the way for a better understanding of the roles of mRNAs and lncRNAs in fecundity-related oviduct function in goats.

Analysis of the Association of Two SNPs in the Promoter Regions of the PPP2R5C and SLC39A5 Genes with Litter Size in Yunshang Black Goats

Screening for candidate genes and genetic variants associated with litter size is important for goat breeding. The aim of this study was to analyze the relationship between single nucleotide polymorphisms (SNPs) in PPP2R5C and SLC39A5 and litter size in Yunshang black goats. KASP genotyping was used to detect the SNP genetic markers in the PPP2R5C and SLC39A5 in a population of 569 Yunshang black goats. The results show that there were two SNPs in the PPP2R5C and SLC39A5 promoter regions. Association analysis revealed that the polymorphisms PPP2R5C g.65977743C>T and SLC39A5 g.50676693T>C were significantly associated with the litter size of the third parity of Yunshang black goats (p < 0.05). To further explore the regulatory mechanism of the two genes, the expression of different genotypes of PPP2R5C and SLC39A5 was validated by RT-qPCR and Western blotting. The expression of PPP2R5C was significantly higher in individuals with the TT genotype than in those with the TC and CC genotypes (p < 0.05). The expression of SLC39A5 was also significantly higher in individuals with the TT genotype than in TC and CC genotypes (p < 0.05). Dual luciferase reporter analysis showed that the luciferase activity of PPP2R5C-C variant was significantly higher than that of PPP2R5C-T variant (p < 0.05). The luciferase activity of SLC39A5-T variant was significantly higher than that of SLC39A5-C variant (p < 0.05). Software was used to predict the binding of transcription factors to the polymorphic sites, and the results show that SOX18, ZNF418, and ZNF667 and NKX2-4 and TBX6 might bind to PPP2R5C g.65977743C>T and SLC39A5 g.50676693T>C, respectively. These results provide new insights into the identification of candidate genes for marker-assisted selection (MAS) in goats.

Identification of Genomic Regions and Candidate Genes Associated with Body Weight and Body Conformation Traits in Karachai Goats

The objective of this study was to identify the SNPs and candidate genes related to body weight and seven body conformation traits at the age of 8 months in the Russian aboriginal Karachai goats (n = 269) by conducting genome-wide association studies (GWAS), using genotypes generated by Goat SNP BeadChip (Illumina Inc., USA). We identified 241 SNPs, which were significantly associated with the studied traits, including 47 genome-wide SNPs (p < 10−5) and 194 suggestive SNPs (p < 10−4), distributed among all goat autosomes except for autosome 23. Fifty-six SNPs were common for two and more traits (1 SNP for six traits, 2 SNPs for five traits, 12 SNPs for four traits, 20 SNPs for three traits, and 21 SNPs for two traits), while 185 SNPs were associated with single traits. Structural annotation within a window of 0.4 Mb (±0.2 Mb from causal SNPs) revealed 238 candidate genes. The largest number of candidate genes was identified at Chr13 (33 candidate genes for the five traits). The genes identified in our study were previously reported to be associated with growth-related traits in different livestock species. The most significant genes for body weight were CRADD, HMGA2, MSRB3, MAX, HACL1 and RAB15, which regulate growth processes, body sizes, fat deposition, and average daily gains. Among them, the HMGA2 gene is a well-known candidate for prenatal and early postnatal development, and the MSRB3 gene is proposed as a candidate gene affecting the growth performance. APOB, PTPRK, BCAR1, AOAH and ASAH1 genes associated with withers height, rump height and body length, are involved in various metabolic processes, including fatty acid metabolism and lipopolysaccharide catabolism. In addition, WDR70, ZBTB24, ADIPOQ, and SORCS3 genes were linked to chest width. KCNG4 was associated with rump height, body length and chest perimeter. The identified candidate genes can be proposed as molecular markers for growth trait selection for genetic improvement in Karachai goats.

Genetic Diversity and Selection Signatures in Jianchang Black Goats Revealed by Whole-Genome Sequencing Data

Understanding the genetic composition of indigenous goats is essential to promote the scientific conservation and sustainable utilization of these breeds. The Jianchang Black (JC) goat, a Chinese native breed, is solid black and exhibits crude feed tolerance, but is characterized by a low growth rate and small body size. Based on the whole-genome sequencing data for 30 JC, 41 Jintang Black (JT), and 40 Yunshang Black (YS) goats, and 21 Bezoar ibexes, here, we investigated the genetic composition of JC goats by conducting analyses of the population structure, runs of homozygosity (ROH), genomic inbreeding, and selection signature. Our results revealed that JT and YS showed a close genetic relationship with a non-negligible amount of gene flows but were genetically distant from JC, apart from Bezoars. An average of 2039 ROHs were present in the autosomal genome per individual. The ROH-based inbreeding estimates in JC goats generally showed moderate values ranging from 0.134 to 0.264, mainly due to rapid declines in the effective population size during recent generations. The annotated genes (e.g., IL2, IL7, and KIT) overlapping with ROH islands were significantly enriched in immune-related biological processes. Further, we found 61 genes (e.g., STIM1, MYO9A, and KHDRBS2) under positive selection in JC goats via three complementary approaches, which may underly genetic adaptations to local environmental conditions. Our findings provided references for the conservation and sustainable utilization of JC goats.

Effect of Upregulation of Transcription Factor TFDP1 Binding Promoter Activity Due to RBP4 g.36491960G>C Mutation on the Proliferation of Goat Granulosa Cells

Retinol-binding protein 4 (RBP4), a member of the lipocalin family, is a specific carrier of retinol (vitamin A) in the blood. Numerous studies have shown that RBP4 plays an important role in mammalian embryonic development and that mutations in RBP4 can be used for the marker-assisted selection of animal reproductive traits. However, there are few studies on the regulation of reproduction and high-prolificacy traits by RBP4 in goats. In this study, the 5′ flanking sequence of RBP4 was amplified, and a G>C polymorphism in the promoter region -211 bp (g.36491960) was detected. An association analysis revealed that the respective first, second and third kidding number and mean kidding number of nanny goats with CC and GC genotypes (2.167 ± 0.085, 2.341 ± 0.104, 2.529 ± 0.107 and 2.189 ± 0.070 for CC and 2.052 ± 0.047, 2.206 ± 0.057, 2.341 ± 0.056 and 2.160 ± 0.039 for GC) were significantly higher (p < 0.05) than those with the GG genotype (1.893 ± 0.051, 2.027 ± 0.064, 2.107 ± 0.061 and 1.74 ± 0.05). The luciferase assay showed that luciferase activity was increased in C allele individuals compared with that in G allele individuals. A competitive electrophoretic mobility shift assay (EMSA) showed that individuals with the CC genotype had a stronger promoter region binding capacity than those with the GG genotype. In addition, transcription factor prediction software showed that the RBP4 g.36491960G>C mutation added a novel binding site for transcription factor DP-1 (TFDP1). RT−qPCR results showed that the expression of TFDP1 was significantly higher in the high-prolificacy group than in the low-prolificacy group, and the expression of RBP4 was higher in both the CC and GC genotypes than that in the GG genotype. TFDP1 overexpression significantly increased the expression of RBP4 mRNA (p < 0.05) and the expression of the cell proliferation factors cyclin-D1, cyclin-D2 and CDK4 (p < 0.05). The opposite trend was observed after interference with TFDP1. Both the EdU and CCK-8 results showed that TFDP1 expression could regulate the proliferation of goat ovarian granulosa cells. In summary, our results showed that RBP4 g.36491960G>C was significantly associated with fecundity traits in goats. The g.36491960G>C mutation enhanced the transcriptional activity of RBP4 and increased the expression of RBP4, thus improving the fertility of Yunshang black goats.

Single nucleotide polymorphisms in the 3′ UTR of follistatin-like 4 and scavenger receptor class B member 1 are associated with Dazu black goat litter size

The untranslated regions (UTR) of genes play crucial roles in regulating gene expression at the post-transcriptional level such as affecting mRNA stabilization. In this study, 26 single nucleotide polymorphisms (SNPs) and one deletion located in UTR were genotyped from 186 Dazu black goats via SNaPshot, and the correlation between genotype and litter size was analyzed. Results indicated that two SNP loci, SNP_chr17-20182525 and SNP_chr7-65652612, which are located at the 3′UTR of scavenger receptor class B member 1 and follistatin-like 4, are significantly (p<0.05) correlated with the litter size of first parity goats. SNP_chr7-65652612 is also significantly associated with the total litter size of first and second parity offspring (p<0.05). In conclusion, SNP_chr7-65652612 and SNP_chr17-20182525 have correlation with the litter size of Dazu black goat and they are potential genetic markers for litter size breeding.

Genome-Wide Association Study of Body Conformation Traits by Whole Genome Sequencing in Dazu Black Goats

Simple Summary

Body conformation traits are economically important in the goat meat industry. Good growth performance in goats, including an accelerated growth rate, can improve carcass weight and meat yield. The identification of genetic variants associated with these traits provides a basis for the genetic improvement of growth performance. In this study, we measured six body conformation traits, including body height, body length, cannon circumference, chest depth, chest width, and heart girth. By a genome-wide association study of a Chinese meat goat breed, 53 significant single nucleotide polymorphisms and 42 candidate genes associated with these traits were detected. These findings improve our understanding of the genetic basis of body conformation traits in goats.

Abstract

Identifying associations between genetic markers and economic traits has practical benefits for the meat goat industry. To better understand the genomic regions and biological pathways contributing to body conformation traits of meat goats, a genome-wide association study was performed using Dazu black goats (DBGs), a Chinese indigenous goat breed. In particular, 150 DBGs were genotyped by whole-genome sequencing, and six body conformation traits, including body height (BH), body length (BL), cannon circumference (CC), chest depth (CD), chest width (CW), and heart girth (HG), were examined. In total, 53 potential SNPs were associated with these body conformation traits. A bioinformatics analysis was performed to evaluate the genes located close to the significant SNPs. Finally, 42 candidate genes (e.g., PSTPIP2, C7orf57, CCL19, FGF9, SGCG, FIGN, and SIPA1L) were identified as components of the genetic architecture underlying body conformation traits. Our results provide useful biological information for the improvement of growth performance and have practical applications for genomic selection in goats.

Genomic characterization of the world's longest selection experiment in mouse reveals the complexity of polygenic traits

BACKGROUND

Long-term selection experiments are a powerful tool to understand the genetic background of complex traits. The longest of such experiments has been conducted in the Research Institute for Farm Animal Biology (FBN), generating extreme mouse lines with increased fertility, body mass, protein mass and endurance. For >140 generations, these lines have been maintained alongside an unselected control line, representing a valuable resource for understanding the genetic basis of polygenic traits. However, their history and genomes have not been reported in a comprehensive manner yet. Therefore, the aim of this study is to provide a summary of the breeding history and phenotypic traits of these lines along with their genomic characteristics. We further attempt to decipher the effects of the observed line-specific patterns of genetic variation on each of the selected traits.

RESULTS

Over the course of >140 generations, selection on the control line has given rise to two extremely fertile lines (>20 pups per litter each), two giant growth lines (one lean, one obese) and one long-distance running line. Whole genome sequencing analysis on 25 animals per line revealed line-specific patterns of genetic variation among lines, as well as high levels of homozygosity within lines. This high degree of distinctiveness results from the combined effects of long-term continuous selection, genetic drift, population bottleneck and isolation. Detection of line-specific patterns of genetic differentiation and structural variation revealed multiple candidate genes behind the improvement of the selected traits.

CONCLUSIONS

The genomes of the Dummerstorf trait-selected mouse lines display distinct patterns of genomic variation harbouring multiple trait-relevant genes. Low levels of within-line genetic diversity indicate that many of the beneficial alleles have arrived to fixation alongside with neutral alleles. This study represents the first step in deciphering the influence of selection and neutral evolutionary forces on the genomes of these extreme mouse lines and depicts the genetic complexity underlying polygenic traits.

Hypothalamic Transcriptome Analysis Reveals the Crucial MicroRNAs and mRNAs Affecting Litter Size in Goats

The hypothalamus was the coordination center of the endocrine system, which played an important role in goat reproduction. However, the molecular mechanism of hypothalamus regulating litter size in goats was still poorly understood. This study aims to investigate the key functional genes associated with prolificacy by hypothalamus transcriptome analysis of goats. In this research, an integrated analysis of microRNAs (miRNAs)-mRNA was conducted using the hypothalamic tissue of Yunshang black goats in the follicular stage. A total of 72,220 transcripts were detected in RNA-seq. Besides, 1,836 differentially expressed genes (DEGs) were identified between high fecundity goats at the follicular phase (FP-HY) and low fecundity goats at the follicular phase (FP-LY). DEGs were significantly enriched in 71 Gene Ontology (GO) terms and 8 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The transcriptome data suggested that DEGs such as BMPR1B, FGFR1, IGF1 and CREB1 are directly or indirectly involved in many processes like hypothalamic gonadal hormone secretion. The miRNA-seq identified 1,837 miRNAs, of which 28 differentially expressed miRNAs (DEMs). These DEMs may affect the nerve cells survival of goat hypothalamic regulating the function of target genes and further affect the hormone secretion activities related to reproduction. They were enriched in prolactin signaling pathway, Jak-STAT signaling pathway and GnRH signaling pathway, as well as various metabolic pathways. Integrated analysis of DEMs and DEGs showed that 87 DEGs were potential target genes of 28 DEMs. After constructing a miRNA-mRNA pathway network, we identified several mRNA-miRNAs pairs by functional enrichment analysis, which was involved in hypothalamic nerve apoptosis. For example, NTRK3 was co-regulated by Novel-1187 and Novel-566, as well as another target PPP1R13L regulated by Novel-566. These results indicated that these key genes and miRNAs may play an important role in the development of goat hypothalamus and represent candidate targets for further research. This study provides a basis for further explanation of the basic molecular mechanism of hypothalamus, but also provides a new idea for a comprehensive understanding of prolificacy characteristics in Yunshang black goats.

Identification of genes associated with litter size combining genomic approaches in Luzhong mutton sheep

Litter size is one of the most important reproductive traits of sheep, which has pronounced effects on the profit of husbandry enterprises and enthusiasm of breeders. Despite the importance of litter size, the underlying genetic mechanisms have not been entirely elucidated. Therefore, based on a high-density SNP chip, genome-wide comparative analysis was performed between two groups with different fecundity to reveal candidate genes linked to litter size via detection of homozygosity and selection signatures in Luzhong mutton sheep. Consequently, nine promising genes were identified from six runs of homozygosity islands, and functionally linked to reproduction (ACTL7A, ACTL7B, and ELP1), embryonic development (KLF5 and PIBF1), and cell cycle (DACH1, BORA, DIS3, and MZT1). A total of 128 genes were observed under selection, of which HECW1 and HTR1E were related to total lambs born, GABRG3, LRP1B, and MACROD2 to teat number, and AGBL1 to reproductive seasonality. Additionally, the presence of inbreeding depression implies the urgency of reasonable mating system to increase litter size in the present herd. These findings provide a comprehensive insight to the genetic makeup of litter size, and also contribute to implementation of marker-assisted selection in sheep.

Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep

The litter size of domestic goats and sheep is an economically important trait that shows variation within breeds. Strenuous efforts have been made to understand the genetic mechanisms underlying prolificacy in goats and sheep. However, there has been a paucity of research on the genetic convergence of prolificacy between goats and sheep, which likely arose because of similar natural and artificial selection forces. Here, we performed comparative genomic and transcriptomic analyses to identify the genetic convergence of prolificacy between goats and sheep. By combining genomic and transcriptomic data for the first time, we identified this genetic convergence in (1) positively selected genes (CHST11 and SDCCAG8), (2) differentially expressed genes (SERPINA14, RSAD2, and PPIG at follicular phase, and IGF1, GPRIN3, LIPG, SLC7A11, and CHST15 at luteal phase), and (3) biological pathways (genomic level: osteoclast differentiation, ErbB signaling pathway, and relaxin signaling pathway; transcriptomic level: the regulation of viral genome replication at follicular phase, and protein kinase B signaling and antigen processing and presentation at luteal phase). These results indicated the potential physiological convergence and enhanced our understanding of the overlapping genetic makeup underlying litter size in goats and sheep.

Whole-genome sequencing to identify candidate genes for litter size and to uncover the variant function in goats (Capra hircus)

To select candidate genes for goat prolificacy, we managed six multi- and six single-kid female goats at the same feeding level and in the same management mode over a 4-year period. These goats showed stable differences in litter size over five continuous parturition records. Whole-genome re-sequencing was used in all 12 to select candidate genes, namely, AURKA, ENDOG, SOX2, RORA, GJA10, RXFP2, CDC25C, and NANOS3, by the strength of their differentiation signals. Most of the selected genes were enriched in the coiled coil process and ovarian development, which suggests that the coiled coil process has a potential regulatory effect on fecundity. Detection of the distribution of variants and association analyses with litter size in 400 goats showed that NANOS3 exon mutations may lead to a transformation of the protein structure. The variation in CDC25C, ENDOG, and NANOS3 showed a significant association with litter size. These results can contribute to the improvement of reproduction traits in the artificial breeding of goats.