De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Neuron-derived neurotrophic factor promotes the differentiation of intramuscular and subcutaneous adipocytes in goat.

Adipocyte play an important role in human health and meat quality by influencing the tenderness, flavor, and juiciness of mutton It has been shown that neuron-derived neurotrophic factor (NENF) is closely related to energy metabolism and adipocyte differentiation in bovine. However, the role of NENF in the goats remains unclear. The aim of this study was to detect the expression of NENF in goat subcutaneous and intramuscular adipocytes, temporal expression profiles of the NENF, and overexpressed NENF on the differentiation of different adipocytes. In this study, PCR amplification successfully cloned the goat NENF gene with a fragment length of 521 bp. In addition, the time point of highest expression of NENF differed between these two adipocytes differentiation processes. Overexpression of NENF in intramuscular and subcutaneous adipocytes promoted the expression levels of differentiation markers CEBPß and SREBP, which in turn promoted the differentiation of intramuscular and subcutaneous adipocytes. This study will provide basic data for further study of the role of goats in goat adipocyte differentiation and for the final elucidation of its molecular mechanisms in regulating goat adipocyte deposition.

Comparison of genomic prediction methods for early growth traits of Inner Mongolia cashmere goats based on multi trait models.

Inner Mongolia cashmere goat is an excellent livestock breed formed through long-term natural selection and artificial breeding, and is currently a world-class dual-purpose breed producing cashmere and meat. Multi trait animal model is considered to significantly improve the accuracy of genetic evaluation in livestock and poultry, enabling indirect selection between traits. In this study, the pedigree, genotype, environment, and phenotypic records of early growth traits of Inner Mongolia cashmere goats were used to build multi trait animal model., Then three methods including ABLUP, GBLUP, and ssGBLUP wereused to estimate the genetic parameters and genomic breeding values of early growth traits (birth weight, weaning weight, average daily weight gain before weaning, and yearling weight). The accuracy and reliability of genomic estimated breeding value are further evaluated using the five fold cross validation method. The results showed that the heritability of birth weight estimated by three methods was 0.13-0.15, the heritability of weaning weight was 0.13-0.20, heritability of daily weight gain before weaning was 0.11-0.14, and the heritability of yearling weight was 0.09-0.14, all of which belonged to moderate to low heritability. There is a strong positive genetic correlation between weaning weight and daily weight gain before weaning, daily weight gain before weaning and yearling weight, with correlation coefficients of 0.77-0.79 and 0.56-0.67, respectively. The same pattern was found in phenotype correlation among traits. The accuracy of the estimated breeding values by ABLUP, GBLUP, and ssGBLUP methods for birth weight is 0.5047, 0.6694, and 0.7156, respectively; the weaning weight is 0.6207, 0.6456, and 0.7254, respectively; the daily weight gain before weaning was 0.6110, 0.6855, and 0.7357 respectively; and the yearling weight was 0.6209, 0.7155, and 0.7756, respectively. In summary, the early growth traits of Inner Mongolia cashmere goats belong to moderate to low heritability, and the speed of genetic improvement is relatively slow. The genetic improvement of other growth traits can be achieved through the selection of weaning weight. The ssGBLUP method has the highest accuracy and reliability in estimating genomic breeding value of early growth traits in Inner Mongolia cashmere goats, and is significantly higher than that from ABLUP method, indicating that it is the best method for genomic breeding of early growth weight in Inner Mongolia cashmere goats.

Genetic parameters for reproductive traits and somatic cell score in U.S. dairy goats.

There are no studies regarding the estimation of genetic parameters and genetic trends for reproductive traits and somatic cells in goats. Their knowledge allows optimization of selection schemes. The objective of this study was to estimate genetic parameters and genetic and phenotypic trends for age at first kidding (AFK), kidding interval (KIN) and somatic cell score (SCS). Analyses were conducted within and across seven US goat breeds, namely, Nubian (NU), Alpine (AL), LaMancha (LM), Toggenburg (TO), Saanen (SA), Nigerian Dwarf (ND) and Oberhasli (OB), and a set of all of these breeds (AB). The restricted maximum likelihood methodology and trivariate animal models were used. Genetic and phenotypic trends were estimated using regression models. The average and standard deviation of AFK, KIN and SCS for AB were 573.6 ± 178.5 days, 418.8 ± 125.5 days and 4.67 ± 2.23 Log2, respectively. The heritabilities (h2) and standard errors of AFK, KIN and SCS for AB were 0.28 ± 0.02, 0.04 ± 0.02 and 0.22 ± 0.01, respectively. The h2 ranged from 0.15 (SA) to 0.37 (NU) for AFK, from 0.04 (AB) to 0.10 (AL) for KIN, and from 0.11 (TO) to 0.26 (LM and ND) for SCS. Genetic correlations between AFK and KIN and between AFK and SCS for AB were positive and weak (0.07 and 0.12, respectively) but significant (P < 0.01). Genetic correlations between SCS and KIN were significant (P < 0.01) for all the breeds and ranged from -0.15 (NU) to 0.44 (AL). Genetic correlations between AFK and SCS in the NU and AL breeds were similar (approximately 0.21). A positive genetic trend was found for KIN in the SA breed, which caused an increase in the number of days between consecutive kiddings. The genetic trend of SCS for the NU, AL and ND breeds was negative and decreased annually, which is beneficial for producers. These first results show the intensity and direction of some favorable/unfavorable relationships between AFK or KIN and SCS Log2 in some U.S. goat genetic groups.

Double-Tube Multiplex TaqMan Real-Time PCR for the Detection of Eight Animal-Derived Dairy Ingredients.

Milk and dairy products represent important sources of nutrition in our daily lives. The identification of species within dairy products holds importance for monitoring food adulteration and ensuring traceability. This study presented a method that integrated double-tube and duplex real-time polymerase chain reaction (PCR) with multiplex TaqMan probes to enable the high-throughput detection of animal-derived ingredients in milk and dairy products. The detection system utilized one pair of universal primers, two pairs of specific primers, and eight animal-derived specific probes for cow, buffalo, goat, sheep, camel, yak, horse, and donkey. These components were optimized within a double-tube and four-probe PCR multiplex system. The developed double-tube detection system could simultaneously identify the above eight targets with a detection limit of 10-0.1 pg/µL. Validation using simulated adulterated milk samples demonstrated a detection limit of 0.1%. The primary advantage of this method lies in the simplification of the multiplex quantitative real-time PCR (qPCR) system through the use of universal primers. This method provides an efficient approach for detecting ingredients in dairy products, providing powerful technical support for market supervision.

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.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats.

BACKGROUND: Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate. RESULTS: To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds. CONCLUSIONS: The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Development and verification of a 10K liquid chip for Hainan black goat based on genotyping by pinpoint sequencing of liquid captured targets.

BACKGROUND: China has thousands years of goat breeding and abundant goat genetic resources. Additionally, the Hainan black goat is one of the high-quality local goat breeds in China. In order to conserve the germplasm resources of the Hainan black goat, facilitate its genetic improvement and further protect the genetic diversity of goats, it is urgent to develop a single nucleotide polymorphism (SNP) chip for Hainan black goat. RESULTS: In this study, we aimed to design a 10K liquid chip for Hainan black goat based on genotyping by pinpoint sequencing of liquid captured targets (cGPS). A total of 45,588 candidate SNP sites were obtained, 10,677 of which representative SNP sites were selected to design probes, which finally covered 9,993 intervals and formed a 10K cGPS liquid chip for Hainan black goat. To verify the 10K cGPS liquid chip, some southern Chinese goat breeds and a sheep breed with similar phenotype to the Hainan black goat were selected. A total of 104 samples were used to verify the clustering ability of the 10K cGPS liquid chip for Hainan black goat. The results showed that the detection rate of sites was 97.34% -99.93%. 84.5% of SNP sites were polymorphic. The heterozygosity rate was 3.08%-36.80%. The depth of more than 99.4% sites was above 10X. The repetition rate was 99.66%-99.82%. The average consistency between cGPS liquid chip results and resequencing results was 85.58%. In addition, the phylogenetic tree clustering analysis verified that the SNP sites on the chip had better clustering ability. CONCLUSION: These results indicate that we have successfully realized the development and verification of the 10K cGPS liquid chip for Hainan black goat, which provides a useful tool for the genome analysis of Hainan black goat. Moreover, the 10K cGPS liquid chip is conducive to the research and protection of Hainan black goat germplasm resources and lays a solid foundation for its subsequent breeding work.

Transcriptomic and metabolomic dissection of skeletal muscle of crossbred Chongming white goats with different meat production performance.

BACKGROUND: The transcriptome and metabolome dissection of the skeletal muscle of high- and low- growing individuals from a crossbred population of the indigenous Chongming white goat and the Boer goat were performed to discover the potential functional differentially expressed genes (DEGs) and differential expression metabolites (DEMs). RESULTS: A total of 2812 DEGs were detected in 6 groups at three time stages (3,6,12 Month) in skeletal muscle using the RNA-seq method. A DEGs set containing seven muscle function related genes (TNNT1, TNNC1, TNNI1, MYBPC2, MYL2, MHY7, and CSRP3) was discovered, and their expression tended to increase as goat muscle development progressed. Seven DEGs (TNNT1, FABP3, TPM3, DES, PPP1R27, RCAN1, LMOD2) in the skeletal muscle of goats in the fast-growing and slow-growing groups was verified their expression difference by reverse transcription-quantitative polymerase chain reaction. Further, through the Liquid chromatography-mass spectrometry (LC-MS) approach, a total of 183 DEMs in various groups of the muscle samples and these DEMs such as Queuine and Keto-PGF1α, which demonstrated different abundance between the goat fast-growing group and slow-growing group. Through weighted correlation network analysis (WGCNA), the study correlated the DEGs with the DEMs and identified 4 DEGs modules associated with 18 metabolites. CONCLUSION: This study benefits to dissection candidate genes and regulatory networks related to goat meat production performance, and the joint analysis of transcriptomic and metabolomic data provided insights into the study of goat muscle development.

MiR-103-5p deficiency suppresses lipid accumulation via upregulating PLSCR4 and its host gene PANK3 in goat mammary epithelial cells.

Lipids are intimately related to the unique flavor and nutritional values of goat milk. MicroRNAs (miRNA) participate in the regulation of various biological functions, including the synthesis and degradation of lipids. Several studies have shown that miR-103 is involved in the regulation of lipid metabolism, however, the molecular mechanism by which miR-103 regulates lipid metabolism in goat mammary gland is poorly understood. In this study, miR-103 was knocked out in goat mammary epithelial cells (GMECs) by CRISPR/Cas9, and the accumulation of lipid droplets, triglycerides, and cholesterol in the cells was suppressed subsequently. Overexpression or knockdown of miR-103-5p and miR-103-3p in GMECs revealed that it was miR-103-5p that promoted lipid accumulation but not miR-103-3p. In addition, Pantothenate Kinase 3 (PANK3), the host gene of miR-103, and Phospholipid Scramblase 4 (PLSCR4) were identified as the target genes of miR-103-5p by dual fluorescein and miRNA pulldown. Furthermore, we identified that cellular lipid levels were negatively regulated by PANK3 and PLSCR4. Lastly, in miR-103 knockout GMECs, the knockdown of PANK and PLSCR4 rescued the lipid accumulation. These findings suggest that miR-103-5p promotes lipid accumulation by targeting PLSCR4 and the host gene PANK3 in GMECs, providing new insights for the regulation of goat milk lipids via miRNAs.

Accuracy of Genomic prediction for fleece traits in Inner Mongolia Cashmere goats.

The fleece traits are important economic traits of goats. With the reduction of sequencing and genotyping cost and the improvement of related technologies, genomic selection for goats has become possible. The research collect pedigree, phenotype and genotype information of 2299 Inner Mongolia Cashmere goats (IMCGs) individuals. We estimate fixed effects, and compare the estimates of variance components, heritability and genomic predictive ability of fleece traits in IMCGs when using the pedigree based Best Linear Unbiased Prediction (ABLUP), Genomic BLUP (GBLUP) or single-step GBLUP (ssGBLUP). The fleece traits considered are cashmere production (CP), cashmere diameter (CD), cashmere length (CL) and fiber length (FL). It was found that year of production, sex, herd and individual ages had highly significant effects on the four fleece traits (P < 0.01). All of these factors should be considered when the genetic parameters of fleece traits in IMCGs are evaluated. The heritabilities of FL, CL, CP and CD with ABLUP, GBLUP and ssGBLUP methods were 0.26 ~ 0.31, 0.05 ~ 0.08, 0.15 ~ 0.20 and 0.22 ~ 0.28, respectively. Therefore, it can be inferred that the genetic progress of CL is relatively slow. The predictive ability of fleece traits in IMCGs with GBLUP (56.18% to 69.06%) and ssGBLUP methods (66.82% to 73.70%) was significantly higher than that of ABLUP (36.73% to 41.25%). For the ssGBLUP method is significantly (29% ~ 33%) higher than that with ABLUP, and which is slightly (4% ~ 14%) higher than that of GBLUP. The ssGBLUP will be as an superiors method for using genomic selection of fleece traits in Inner Mongolia Cashmere goats.

Association of polymorphism in the promotor area of the caprine BMPR1B gene with litter size and body measurement traits in Damani goats.

This study aimed to explore polymorphisms in the promoter region of the caprine BMPR1B (Bone morphogenetic protein receptor 1 beta) gene and its association with body measurement and litter size traits in Damani does. A total of 53 blood samples were collected to analyze the association between the BMPR1B gene polymorphism and 11 phenotypic traits in Damani female goats. The results revealed that three novel SNPs were identified in the promoter region of the caprine BMPR1B gene, including g.67 A > C (SNP1), g.170 G > A(SNP2), and g.501A > T (SNP3), among which the SNP1 and SNP2 were significantly (p < 0.05) associated with litter size and body measurement traits in Damani goats. In SNP1 the AC genotype could be used as a marker for litter size, and the CC genotype for body weight in Damani goats. In SNP2, the genotype GG was significantly (p < 0.05) associated with ear and head length. Therefore, we can conclude from the present study, that genetic variants AC and CC of the caprine BMPR1B gene could be used as genetic markers for economic traits through marker-assisted selection for the breed improvement program of the Damani goat.

Transcriptome analysis of mRNA and miRNA in the development of LeiZhou goat muscles.

The progression of muscle development is a pivotal aspect of animal ontogenesis, where miRNA and mRNA exert substantial influence as prominent players. It is important to understand the molecular mechanisms involved in skeletal muscle development to enhance the quality and yield of meat produced by Leizhou goats. We employed RNA sequencing (RNA-SEQ) technology to generate miRNA-mRNA profiles in Leizhou goats, capturing their developmental progression at 0, 3, and 6 months of age. A total of 977 mRNAs and 174 miRNAs were found to be differentially expressed based on our analysis. Metabolic pathways, calcium signaling pathways, and amino acid synthesis and metabolism were found to be significantly enriched among the differentially expressed mRNA in the enrichment analysis. Meanwhile, we found that among these differentially expressed mRNA, some may be related to muscle development, such as MYL10, RYR3, and CSRP3. Additionally,, we identified five muscle-specific miRNAs (miR-127-3p, miR-133a-3p, miR-193b-3p, miR-365-3p, and miR-381) that consistently exhibited high expression levels across all three stages. These miRNAs work with their target genes (FHL3, SESN1, PACSIN3, LMCD1) to regulate muscle development. Taken together, our findings suggest that several miRNAs and mRNAs are involved in regulating muscle development and cell growth in goats. By uncovering the molecular mechanisms involved in muscle growth and development, these findings contribute valuable knowledge that can inform breeding strategies aimed at enhancing meat yield and quality in Leizhou goats.

Selection signatures and landscape genomics analysis to reveal climate adaptation of goat breeds.

Goats have achieved global prominence as essential livestock since their initial domestication, primarily owing to their remarkable adaptability to diverse environmental and production systems. Differential selection pressures influenced by climate have led to variations in their physical attributes, leaving genetic imprints within the genomes of goat breeds raised in diverse agroecological settings. In light of this, our study pursued a comprehensive analysis, merging environmental data with single nucleotide polymorphism (SNP) variations, to unearth indications of selection shaped by climate-mediated forces in goats. Through the examination of 43,300 SNPs from 51 indigenous goat breeds adapting to different climatic conditions using four analytical methods: latent factor mixed models (LFMM), F-statistics (Fst), Extended haplotype homozygosity across populations (XPEHH), and spatial analysis method (SAM), A total of 74 genes were revealed to display clear signs of selection, which are believed to be influenced by climatic conditions. Among these genes, 32 were consistently identified by at least two of the applied methods, and three genes (DENND1A, PLCB1, and ITPR2) were confirmed by all four approaches. Moreover, our investigation yielded 148 Gene Ontology (GO) terms based on these 74 genes, underlining pivotal biological pathways crucial for environmental adaptation. These pathways encompass functions like vascular smooth muscle contraction, cellular response to heat, GTPase regulator activity, rhythmic processes, and responses to temperature stimuli. Of significance, GO terms about endocrine regulation and energy metabolic responses, key for local adaptation were also uncovered, including biological processes, such as cell differentiation, regulation of peptide hormone secretion, and lipid metabolism. These findings contribute to our knowledge of the genetic structure of climate-triggered adaptation across the goat genome and have practical implications for marker-assisted breeding in goats.

Melatonin's effect on hair follicles in a goat (Capra hircus) animal model.

Introduction: Melatonin can treat androgenetic alopecia in males. Goats can be used as animal models to study melatonin treatment for human alopecia. In this study, a meta-analysis of melatonin's effects on goat hair follicles was pursued. Methods: Literature from the last 20 years was searched in Scopus, Science Direct, Web of Science and PubMed. Melatonin's effect on goat hair follicles and litter size were performed through a traditional meta-analysis and trial sequential analysis. A network meta-analysis used data from oocyte development to blastocyst. The hair follicle genes regulated by melatonin performed KEGG and PPI. We hypothesized that there are differences in melatonin receptors between different goats, and therefore completed melatonin receptor 1A homology modelling and molecular docking. Results: The results showed that melatonin did not affect goat primary follicle or litter size. However, there was a positive correlation with secondary follicle growth. The goat melatonin receptor 1A SNPs influence melatonin's functioning. The wild type gene defect MR1 is a very valuable animal model. Discussion: Future studies should focus on the relationship between goat SNPs and the effect of embedded melatonin. This study will provide theoretical guidance for the cashmere industry and will be informative for human alopecia research.

Comparative metabolomics reveals serum metabolites changes in goats during different developmental stages.

Goats can provide meat, milk and skins for humans and are livestock with high economic benefits. However, despite their economic significance, the comprehensive analysis of goats' serum metabolic profile and its intricate alterations throughout their developmental journey remains conspicuously absent. To investigate the stage-specificity and dynamic change characteristics of metabolites during the growth and development of goats, this study compared the alterations in serum hormone levels and serum biochemical markers across different developmental stages of female goats (1, 60, 120 and 180 days old; n = 5). Additionally, a serum untargeted LC-MS metabolomics analysis was conducted. A total of 504 DAMs were identified with age. The results indicated that PE, PC, Lyso-PE, Lyso-PC and FAFHA may play important roles in lipid metabolism in goats after birth. Weighted gene co-expression network analysis (WGCNA) identified two metabolite modules (Turquoise and Yellow) and key metabolites within these modules that were significantly associated with phenotypic features. L-carnitine may be a metabolite related to muscle development in goats. The findings of this study demonstrate notable variations in serum metabolites across distinct developmental phases in goats. Lipids and organic acids play important roles in different developmental stages of goats.

Genetic evaluations of dairy goats with few pedigree data: different approaches to use molecular information.

One of the limitations of implementing animal breeding programs in small-scale or extensive production systems is the lack of production records and genealogical records. In this context, molecular markers could help to gain information for the breeding program. This study addresses the inclusion of molecular data into traditional genetic evaluation models as a random effect by molecular pedigree reconstruction and as a fixed effect by Bayesian clustering. The methods were tested for lactation curve traits in 14 dairy goat herds with incomplete phenotypic data and pedigree information. The results showed an increment of 37.3% of the relationships regarding the originals with MOLCOAN and clustering into five genetic groups. Data leads to estimating additive variance, error variance, and heritability with four different models, including pedigree and molecular information. Deviance Information Criterion (DIC) values demonstrate a greater fitting of the models that include molecular information either as fixed (genetic clusters) or as random (molecular matrix) effects. The molecular information of simple markers can complement genetic improvement strategies in populations with little information.

A comparison of genetic and genomic breeding values in Saanen and Alpine goats.

Nowadays, several countries are developing or adopting genomic selection in the dairy goat sector. The most used method to estimate breeding values is Single-Step Genomic Best Linear Unbiased Prediction (ssGBLUP) which offers several advantages in terms of computational process and accuracy of the estimated breeding values (EBVs). Saanen and Alpine are the predominant dairy goat breeds in Italy, and both have similar breeding programs where EBVs for productive traits are currently calculated using BLUP. This work describes the implementation of genomic selection for these two breeds in Italy, aligning with the selection practices already carried out in the international landscape. The available dataset included 3 611 genotyped animals, 11 470 lactation records, five traits (milk, protein and fat yields, and fat and protein percentages), and three-generation pedigrees. EBVs were estimated using BLUP, GBLUP, and ssGBLUP both with single and multiple trait approaches. The methods were compared in terms of correlation between EBVs and genetic trends. Results were also validated with the linear regression method excluding part of the phenotypic data. In both breeds, EBVs and GEBVs were strongly correlated and the trend of each trait was similar comparing the three methods. The average increase in accuracy across traits and methods amounted to +13 and +10% from BLUP to ssGBLUP for Alpine and Saanen breeds, respectively. Results indicated higher prediction accuracy and correlation for GBLUP and ssGBLUP compared to BLUP, implying that the use of genotypes increases the accuracy of EBVs, particularly in the absence of phenotypic data. Therefore, ssGBLUP is likely to be the most effective method to enhance genetic gain in Italian Saanen and Alpine goats.

The SLC19A1-AS/miR-1343/WNT11 axis is a novel positive regulatory ceRNA network governing goat granulosa cell proliferation.

Long noncoding RNAs (lncRNAs), as competitive endogenous RNAs (ceRNAs), can directly or indirectly affect the proliferation and apoptosis of granulosa cells by regulating microRNA (miRNA) pathways. A ceRNA network of the SLC19A1-AS-miR-1343-WNT11 axis was constructed via comprehensive transcriptome sequencing of ovaries from goats with various fertility levels to further elucidate the function and regulatory mechanism of SLC19A1-AS in modulating miR-1343 and WNT11 during granulosa cell proliferation and apoptosis. Subsequent validation experiments were conducted in vitro using granulosa cells. In these experiments, we performed RNA immunoprecipitation (RIP) and identified SLC19A1-AS as a ceRNA in goat granulosa cells that promoted proliferation. Through bioinformatics prediction, luciferase reporter gene assays, and RNA pulldown assays, we confirmed that SLC19A1-AS acts as a sponge for miR-1343, preventing its binding to WNT11 mRNA and thereby increasing the expression of WNT11. This interaction also influenced the proliferation and apoptosis of granulosa cells. Our study systematically validated the biological function of the lncRNA-miRNA-mRNA ceRNA network in goat ovaries and revealed the potential regulatory mechanism by which SLC19A1-AS functions as a ceRNA in granulosa cells. These findings are expected to provide an important experimental foundation for further elucidating the physiological regulatory network of the ovary and contributing to reproductive health in goats.

Runs of Homozygosity Detection and Selection Signature Analysis for Local Goat Breeds in Yunnan, China.

Runs of Homozygosity (ROH) are continuous homozygous DNA segments in diploid genomes, which have been used to estimate the genetic diversity, inbreeding levels, and genes associated with specific traits in livestock. In this study, we analyzed the resequencing data from 10 local goat breeds in Yunnan province of China and five additional goat populations obtained from a public database. The ROH analysis revealed 21,029 ROH segments across the 15 populations, with an average length of 1.27 Mb, a pattern of ROH, and the assessment of the inbreeding coefficient indicating genetic diversity and varying levels of inbreeding. iHS (integrated haplotype score) was used to analyze high-frequency Single-Nucleotide Polymorphisms (SNPs) in ROH regions, specific genes related to economic traits such as coat color and weight variation. These candidate genes include OCA2 (OCA2 melanosomal transmembrane protein) and MLPH (melanophilin) associated with coat color, EPHA6 (EPH receptor A6) involved in litter size, CDKAL1 (CDK5 regulatory subunit associated protein 1 like 1) and POMC (proopiomelanocortin) linked to weight variation and some putative genes associated with high-altitude adaptability and immune. This study uncovers genetic diversity and inbreeding levels within local goat breeds in Yunnan province, China. The identification of specific genes associated with economic traits and adaptability provides actionable insights for utilization and conservation efforts.