Sequencing and automated whole-genome optical mapping of the genome of a domestic goat (Capra hircus)

Dataset of the rumen microbiota and epithelial transcriptomics and proteomics in goat affected by solid diets

Although early solid diet supplementation is a common practice to improve the growth and development in goat kids, its biological mechanism how solid diet induces rumen microbiota and epithelial development is still unknow. In this study, rumen fermentation parameters, 16S rRNA sequencing for rumen content and epithelial microbiota, transcriptomics and proteomics of epithelium were determined to classify the effects of solid diet supplementation. Here, we classified the changes of goat phenotypes (i.e., growth performance, rumen fermentation and development) and linked them to the changes of rumen microbiota, transcriptome and expressed proteins. The mechanism of solid diet improving rumen development was elucidated preliminarily. Moreover, different roles between the rumen content and epithelial microbiota were identified. Thess datasets expands our understanding of the association between the early diet intervention and rumen development, providing the useful information how nutrient strategy affects rumen function and subsequently improves the host growth. The generated data provides insights in the importance of rumen niche microbiota and microbe-host interactions, which benefits future studies.

Anthropogenic events and responses to environmental stress are shaping the genomes of Ethiopian indigenous goats

Anthropological and biophysical processes have shaped livestock genomes over Millenia and can explain their current geographic distribution and genetic divergence. We analyzed 57 Ethiopian indigenous domestic goat genomes alongside 67 equivalents of east, west, and north-west African, European, South Asian, Middle East, and wild Bezoar goats. Cluster, ADMIXTURE (K = 4) and phylogenetic analysis revealed four genetic groups comprising African, European, South Asian, and wild Bezoar goats. The Middle Eastern goats had an admixed genome of these four genetic groups. At K = 5, the West African Dwarf and Moroccan goats were separated from East African goats demonstrating a likely historical legacy of goat arrival and dispersal into Africa via the coastal Mediterranean Sea and the Horn of Africa. FST, XP-EHH, and Hp analysis revealed signatures of selection in Ethiopian goats overlaying genes for thermo-sensitivity, oxidative stress response, high-altitude hypoxic adaptation, reproductive fitness, pathogen defence, immunity, pigmentation, DNA repair, modulation of renal function and integrated fluid and electrolyte homeostasis. Notable examples include TRPV1 (a nociception gene); PTPMT1 (a critical hypoxia survival gene); RETREG (a regulator of reticulophagy during starvation), and WNK4 (a molecular switch for osmoregulation). These results suggest that human-mediated translocations and adaptation to contrasting environments are shaping indigenous African goat genomes.

Genome-wide comparative analysis reveals selection signatures for reproduction traits in prolific Suffolk sheep

The identification of genome-wide selection signatures can reveal the potential genetic mechanisms involved in the generation of new breeds through natural or artificial selection. In this study, we screened the genome-wide selection signatures of prolific Suffolk sheep, a new strain of multiparous mutton sheep, to identify candidate genes for reproduction traits and unravel the germplasm characteristics and population genetic evolution of this new strain of Suffolk sheep. Whole-genome resequencing was performed at an effective sequencing depth of 20× for genomic diversity and population structure analysis. Additionally, selection signatures were investigated in prolific Suffolk sheep, Suffolk sheep, and Hu sheep using fixation index (F ST) and heterozygosity H) analysis. A total of 5,236.338 Gb of high-quality genomic data and 28,767,952 SNPs were obtained for prolific Suffolk sheep. Moreover, 99 selection signals spanning candidate genes were identified. Twenty-three genes were significantly associated with KEGG pathway and Gene Ontology terms related to reproduction, growth, immunity, and metabolism. Through selective signal analysis, genes such as ARHGEF4, CATIP, and CCDC115 were found to be significantly correlated with reproductive traits in prolific Suffolk sheep and were highly associated with the mTOR signaling pathway, the melanogenic pathway, and the Hippo signaling pathways, among others. These results contribute to the understanding of the evolution of artificial selection in prolific Suffolk sheep and provide candidate reproduction-related genes that may be beneficial for the establishment of new sheep breeds.

The Role of BMP7 in the Proliferation of Hu Sheep Dermal Papilla Cells Is Influenced by DNA Methylation

Previous studies have shown that the BMP7 gene is differentially expressed in Hu sheep lamb skin of different pattern types, and its expression level is significantly correlated with hair follicle indices of different pattern types, but the molecular mechanism of the differential expression of the BMP7 gene remains unclear. This study investigated the effect of DNA methylation on the transcriptional expression of BMP7. Firstly, we found that the mRNA expression of the BMP7 gene and the activity of the core promoter of the BMP7 gene were upregulated after 5-Aza-Deoxycytidine-induced demethylation treatment using qRT-PCR and double luciferase reporter assay. Then, we found that the proliferation of Hu sheep DPCs in vitro was promoted after 5-Aza-Deoxycytidine-induced demethylation treatment through qRT-PCR, CCK-8, and EdU assay, and that the overexpression of DNMT1 in DPCs induced the opposite effect. In addition, the results of the cell cycle assay reveal that the percentage of cells in the S phase was increased after 5-Aza-Deoxycytidine-induced demethylation treatment, and that the percentage of cells in the S phase was decreased after overexpression of DNMT1 in DPCs. This study indicated that the differential expression of the BMP7 gene in different patterns of Hu sheep lamb skin may be regulated by DNA methylation modification. In addition, DNA methylation can regulate the proliferation and cell cycle of DPCs in Hu sheep.

Integrated analysis of differently expressed microRNAs and mRNAs at different postnatal stages reveals intramuscular fat deposition regulation in goats (Capra hircus)

Intramuscular fat refers to the adipose tissue distributed in the muscle. It is an important indicator that affects the quality of goat meat, and can directly affect the tenderness and flavor of goat meat. Our previous study revealed the mRNA that may be crucial for intramuscular fat deposition during goat growth; however, how the microRNAs (miRNAs) are involved in the process is largely unclear. In the present study, a total of 401 known miRNAs and 120 goat novel miRNAs, including 110 differentially expressed (DE) miRNAs, were identified among longissimus dorsi from three growth stages (2, 9, and 24 months) by miRNA sequencing. Combining analysis of the DE mRNAs and DE miRNAs was then performed by miRDB and miRwalk, and miR-145-5p and FOXO1, miR-487b-3p, and PPARG coactivator 1 α (PPARGC1A), miR-345-3p, and solute carrier family 2 member 4 (SLC2A4), etc. were shown to closely associate with lipid metabolism, which was then validated by a correlation analysis. The final DE mRNAs were significantly enriched in fatty acid transmembrane transport, fatty acid homeostasis, apelin signaling pathway, glucagon signaling pathway, insulin signaling pathway, and AMPK signaling pathway by gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis. Besides, miR-145-5p showed a certain effect on goat intramuscular fat metabolism by acting on the possible target gene Forkhead Box O1 (FOXO1). These data provide some theoretical support for improving the quality of goat meat.

Livestock species as emerging models for genomic imprinting

Genomic imprinting is an epigenetically-regulated process of central importance in mammalian development and evolution. It involves multiple levels of regulation, with spatio-temporal heterogeneity, leading to the context-dependent and parent-of-origin specific expression of a small fraction of the genome. Genomic imprinting studies have therefore been essential to increase basic knowledge in functional genomics, evolution biology and developmental biology, as well as with regard to potential clinical and agrigenomic perspectives. Here we offer an overview on the contribution of livestock research, which features attractive resources in several respects, for better understanding genomic imprinting and its functional impacts. Given the related broad implications and complexity, we promote the use of such resources for studying genomic imprinting in a holistic and integrative view. We hope this mini-review will draw attention to the relevance of livestock genomic imprinting studies and stimulate research in this area.

Genetic background of juniper (Juniperus spp.) consumption predicted by fecal near-infrared spectroscopy in divergently selected goats raised in harsh rangeland environments

BACKGROUND

Junipers (Juniperus spp.) are woody native, invasive plants that have caused encroachment problems in the U.S. western rangelands, decreasing forage productivity and biodiversity. A potential solution to this issue is using goats in targeted grazing programs. However, junipers, which grow in dry and harsh environmental conditions, use chemical defense mechanisms to deter herbivores. Therefore, genetically selecting goats for increased juniper consumption is of great interest for regenerative rangeland management. In this context, the primary objectives of this study were to: 1) estimate variance components and genetic parameters for predicted juniper consumption in divergently selected Angora (ANG) and composite Boer x Spanish (BS) goat populations grazing on Western U.S. rangelands; and 2) to identify genomic regions, candidate genes, and biological pathways associated with juniper consumption in these goat populations.

RESULTS

The average juniper consumption was 22.4% (± 18.7%) and 7.01% (± 12.1%) in the BS and ANG populations, respectively. The heritability estimates (realized heritability within parenthesis) for juniper consumption were 0.43 ± 0.02 (0.34 ± 0.06) and 0.19 ± 0.03 (0.13 ± 0.03) in BS and ANG, respectively, indicating that juniper consumption can be increased through genetic selection. The repeatability values of predicted juniper consumption were 0.45 for BS and 0.28 for ANG. A total of 571 significant SNP located within or close to 231 genes in BS, and 116 SNP related to 183 genes in ANG were identified based on the genome-wide association analyses. These genes are primarily associated with biological pathways and gene ontology terms related to olfactory receptors, intestinal absorption, and immunity response.

CONCLUSIONS

These findings suggest that juniper consumption is a heritable trait of polygenic inheritance influenced by multiple genes of small effects. The genetic parameters calculated indicate that juniper consumption can be genetically improved in both goat populations.

Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

High-throughput analysis of lncRNA in cows with naturally infected Staphylococcus aureus mammary gland

LncRNA (long non-coding RNA) is an RNA molecule with a length between 200 and 100,000 nt. It does not encode proteins and is involved in a variety of intracellular processes, becoming a research hotspot of genetics. To identify key lncRNAs associated with dairy mastitis, we collected mammary epithelial tissue samples of Normal disease-free Holstein cows (HCN) and unhealthy Holstein cows with Staphylococcus aureus (HCU) and performed RNA sequencing (RNA-seq) on the samples. A total of 270 differentially expressed lncRNAs and 500 differentially expressed mRNAs were identified by high-throughput sequencing and bioinformatics analysis. Furthermore, Hydrolase activity is the most enriched in GO, and ErbB signaling pathway is significantly enriched in KEGG. In addition, through qPCR validation of 5 candidate lncRNAs in HCN and HCU, four differentially expressed lncRNAs MSTRG.498, MSTRG57.1, MSTRG.41.1 and MSTRG 124.1 were confirmed to have significant differentially expressed in cow mastitis. Also, lncRNA MSTRG.498 and its target gene, SMC4, might directly or indirectly play a role in cow mastitis. The regulatory network of lncRNA-miRNA-mRNA has been inferred from a bioinformatics perspective, which may assist understand the underlying molecular mechanism of lncRNAs involved in regulating mastitis in cows. Our findings will provide meaningful resources for further research on the regulatory function of lncRNAs in cow mastitis.

High-throughput analysis of CircRNA in cows with naturally infected Staphylococcus aureus mammary gland

Circular RNAs (CircRNA) are a special type of non-coding RNA molecule with a closed ring structure and are not affected by RNA exonucases. It has stable expression, is not easy to degrade, and exists in most eukaryotes. However, circRNA regulation of cow mastitis has not been widely recognized. Mammary epithelial tissues were collected from healthy Holstein cows (HCN) and mastitis Holstein cows (HCU). RNA sequencing (RNA SEQ) was performed for the differentially expressed circRNAs, and analysis results showed that 19 differentially expressed circRNAs were identified in HCN and HCU, among which 6 circRNAs were up-regulated and 13 circRNAs were down-regulated. We randomly selected nine circRNAs for Q-PCR verification, and the results showed consistent expression. Three circRNAs: circRNA2860, circRNA5323 and circRNA4027 were confirmed to be significantly differentially expressed circRNAs in cow mastitis. Also, their host genes TRPS1, SLC12A2 and MYH11 might be directly or indirectly play a role in cow mastitis. Furthermore, RNA polymerase transcription factor binding and tight junction are most enriched in GO and KEGG pathways, respectively. In addition, the regulatory network of circRNA-miRNA has been inferred from a bioinformatics perspective, which may help to understand the underlying molecular mechanism of circRNAs involved in regulating mastitis in cows.

Identification of genetic associations and functional SNPs of bovine KLF6 gene on milk production traits in Chinese holstein

BACKGROUND

Our previous research identified the Kruppel like factor 6 (KLF6) gene as a prospective candidate for milk production traits in dairy cattle. The expression of KLF6 in the livers of Holstein cows during the peak of lactation was significantly higher than that during the dry and early lactation periods. Notably, it plays an essential role in activating peroxisome proliferator-activated receptor α (PPARα) signaling pathways. The primary aim of this study was to further substantiate whether the KLF6 gene has significant genetic effects on milk traits in dairy cattle.

RESULTS

Through direct sequencing of PCR products with pooled DNA, we totally identified 12 single nucleotide polymorphisms (SNPs) within the KLF6 gene. The set of SNPs encompasses 7 located in 5' flanking region, 2 located in exon 2 and 3 located in 3' untranslated region (UTR). Of these, the g.44601035G > A is a missense mutation that resulting in the replacement of arginine (CGG) with glutamine (CAG), consequently leading to alterations in the secondary structure of the KLF6 protein, as predicted by SOPMA. The remaining 7 regulatory SNPs significantly impacted the transcriptional activity of KLF6 following mutation (P < 0.005), manifesting as changes in transcription factor binding sites. Additionally, 4 SNPs located in both the UTR and exons were predicted to influence the secondary structure of KLF6 mRNA using the RNAfold web server. Furthermore, we performed the genotype-phenotype association analysis using SAS 9.2 which found all the 12 SNPs were significantly correlated to milk yield, fat yield, fat percentage, protein yield and protein percentage within both the first and second lactations (P < 0.0001 ~ 0.0441). Also, with Haploview 4.2 software, we found the 12 SNPs linked closely and formed a haplotype block, which was strongly associated with five milk traits (P < 0.0001 ~ 0.0203).

CONCLUSIONS

In summary, our study represented the KLF6 gene has significant impacts on milk yield and composition traits in dairy cattle. Among the identified SNPs, 7 were implicated in modulating milk traits by impacting transcriptional activity, 4 by altering mRNA secondary structure, and 1 by affecting the protein secondary structure of KLF6. These findings provided valuable molecular insights for genomic selection program of dairy cattle.

Progress and future perspectives of livestock genomics in India: a mini review

The field of genetics has evolved a lot after the emergence of molecular and advanced genomic technologies. The advent of Next Generation Sequencing, SNP genotyping platforms and simultaneous reduction in the cost of sequencing had opened the door to genomic research in farm animals. There are various applications of genomics in livestock, such as the use of genomic data: (i) to investigate genetic diversity and breed composition/population structure (ii) to identify genetic variants and QTLs related to economically important and ecological traits, genome-wide association studies (GWAS) and genomic signatures of selection; (iii) to enhance breeding programs by genomic selection. Compared to traditional methods, genomic selection is expected to improve selection response by increasing selection accuracy and reducing the generation interval due to early selection. Genomic selection (GS) in developed countries has led to rapid genetic gains, especially in dairy cattle, due to a well-established genetic evaluation system. Indian livestock system is still lagging behind developed nations in adopting these technologies. This review discusses the current status, challenges, and future perspectives of livestock genomics in India.

Genome-Wide Association and Pathway Analysis of Carcass and Meat Quality Traits in Karachai Young Goats

Goats with diverse economic phenotypic traits play an important role in animal husbandry. However, the genetic mechanisms underlying complex phenotypic traits are unclear in goats. Genomic studies of variations provided a lens to identify functional genes. The work aimed to search for candidate genes related to body measurements and body weight of Karachai goats and develop an experimental PCR-RV test system for genotyping significant SNPs. Comparison of GWAS results for ages 4 and 8 months revealed 58 common SNPs for significant genotypes. 11 common SNPs were identified for body weight, 4 SNPs-for group of traits withers height, rump height, body length, 2 SNPs-for withers height and rump height, 1 SNP-for body length and chest depth. Structural annotation of genomic regions covering a window of ±0.20 Mb showed the presence of 288 genes; 52 of them had the described functions in accordance with gene ontology. The main molecular functions of proteins encoded by these genes are the regulation of transcription, cell proliferation, angiogenesis, body growth, fatty acid and lipid metabolism, nervous system development, and spermatogenesis. SNPs common to body weight and localized within a window of ±200 kb from the structural genes CRADD, HMGA2, MSRB3, FUT8, MAX, and RAB15 were selected to create a test system. The study of meat productivity after slaughter and chemical analysis of muscle tissue in Karachai goats at the age of 8 months of different genotypes according to the identified SNPs revealed that rs268269710 is the most promising for further research and use in breeding. The GG genotype is associated with a larger live weight of animals, a larger carcass yield, the content of the boneless part in it, and the ratio of protein and adipose tissue in meat preferred for dietary nutrition. These results will contribute to the genetic improvement of Karachai goats.

The multi-kingdom microbiome of the goat gastrointestinal tract

BACKGROUND

Goat is an important livestock worldwide, which plays an indispensable role in human life by providing meat, milk, fiber, and pelts. Despite recent significant advances in microbiome studies, a comprehensive survey on the goat microbiomes covering gastrointestinal tract (GIT) sites, developmental stages, feeding styles, and geographical factors is still unavailable. Here, we surveyed its multi-kingdom microbial communities using 497 samples from ten sites along the goat GIT.

RESULTS

We reconstructed a goat multi-kingdom microbiome catalog (GMMC) including 4004 bacterial, 71 archaeal, and 7204 viral genomes and annotated over 4,817,256 non-redundant protein-coding genes. We revealed patterns of feeding-driven microbial community dynamics along the goat GIT sites which were likely associated with gastrointestinal food digestion and absorption capabilities and disease risks, and identified an abundance of large intestine-enriched genera involved in plant fiber digestion. We quantified the effects of various factors affecting the distribution and abundance of methane-producing microbes including the GIT site, age, feeding style, and geography, and identified 68 virulent viruses targeting the methane producers via a comprehensive virus-bacterium/archaea interaction network.

CONCLUSIONS

Together, our GMMC catalog provides functional insights of the goat GIT microbiota through microbiome-host interactions and paves the way to microbial interventions for better goat and eco-environmental qualities. Video Abstract.

The African Goat Improvement Network: a scientific group empowering smallholder farmers

The African Goat Improvement Network (AGIN) is a collaborative group of scientists focused on genetic improvement of goats in small holder communities across the African continent. The group emerged from a series of workshops focused on enhancing goat productivity and sustainability. Discussions began in 2011 at the inaugural workshop held in Nairobi, Kenya. The goals of this diverse group were to: improve indigenous goat production in Africa; characterize existing goat populations and to facilitate germplasm preservation where appropriate; and to genomic approaches to better understand adaptation. The long-term goal was to develop cost-effective strategies to apply genomics to improve productivity of small holder farmers without sacrificing adaptation. Genome-wide information on genetic variation enabled genetic diversity studies, facilitated improved germplasm preservation decisions, and provided information necessary to initiate large scale genetic improvement programs. These improvements were partially implemented through a series of community-based breeding programs that engaged and empowered local small farmers, especially women, to promote sustainability of the production system. As with many international collaborative efforts, the AGIN work serves as a platform for human capacity development. This paper chronicles the evolution of the collaborative approach leading to the current AGIN organization and describes how it builds capacity for sustained research and development long after the initial program funds are gone. It is unique in its effectiveness for simultaneous, multi-level capacity building for researchers, students, farmers and communities, and local and regional government officials. The positive impact of AGIN capacity building has been felt by participants from developing, as well as developed country partners.

Evolutionary activation of acidic chitinase in herbivores through the H128R mutation in ruminant livestock

Placental mammals' ancestors were insectivores, suggesting that modern mammals may have inherited the ability to digest insects. Acidic chitinase (Chia) is a crucial enzyme hydrolyzing significant component of insects' exoskeleton in many species. On the other hand, herbivorous animal groups, such as cattle, have extremely low chitinase activity compared to omnivorous species, e.g., mice. The low activity of cattle Chia has been attributed to R128H mutation. The presence of either of these amino acids correlates with the feeding behavior of different bovid species with R and H determining the high and low enzymatic activity, respectively. Evolutionary analysis indicated that selective constraints were relaxed in 67 herbivorous Chia in Cetartiodactyla. Despite searching for another Chia paralog that could compensate for the reduced chitinase activity, no active paralogs were found in this order. Herbivorous animals' Chia underwent genetic alterations and evolved into a molecule with low activity due to the chitin-free diet.

Taming transposable elements in livestock and poultry: a review of their roles and applications

Livestock and poultry play a significant role in human nutrition by converting agricultural by-products into high-quality proteins. To meet the growing demand for safe animal protein, genetic improvement of livestock must be done sustainably while minimizing negative environmental impacts. Transposable elements (TE) are important components of livestock and poultry genomes, contributing to their genetic diversity, chromatin states, gene regulatory networks, and complex traits of economic value. However, compared to other species, research on TE in livestock and poultry is still in its early stages. In this review, we analyze 72 studies published in the past 20 years, summarize the TE composition in livestock and poultry genomes, and focus on their potential roles in functional genomics. We also discuss bioinformatic tools and strategies for integrating multi-omics data with TE, and explore future directions, feasibility, and challenges of TE research in livestock and poultry. In addition, we suggest strategies to apply TE in basic biological research and animal breeding. Our goal is to provide a new perspective on the importance of TE in livestock and poultry genomes.

RUNX1T1 modulates myogenic differentiation by regulating the calcium signaling pathway and the alternative splicing of ROCK2

RUNX1T1 (Runt-related transcription factor 1, translocated to 1) plays a wide-ranging and diverse role in cellular development, including hematopoiesis and adipogenesis. However, little is known about the function of RUNX1T1 in the skeletal muscle development. Here, we assessed the impact of RUNX1T1 on the proliferation and myogenic differentiation of goat primary myoblasts (GPMs). It was observed that RUNX1T1 is highly expressed during the early stages of myogenic differentiation and the fetal stage. Moreover, the knockdown of RUNX1T1 promotes the proliferation and inhibits myogenic differentiation and mitochondrial biogenesis of GPMs. RNA sequencing analysis revealed that significantly differentially expressed genes in RUNX1T1 knockdown cells were enriched in the calcium signaling pathway. Additionally, we discovered that RUNX1T1 regulates alternative splicing (AS) events involved in myogenesis. We also show that silencing RUNX1T1 blocked the Ca²⁺ -CAMK signaling pathway and reduced the expression levels of muscle-specific isoforms of recombinant rho associated coiled coil containing crotein kinase 2 (ROCK2) during myogenic differentiation, partially explaining why RUNX1T1 deficiency leads to the impairment of myotube formation. These findings suggest that RUNX1T1 is a novel regulator of myogenic differentiation that regulates the calcium signaling pathway and AS of ROCK2. Overall, our results highlight the critical role of RUNX1T1 in myogenesis and broaden our understanding of myogenic differentiation.

The orf virus 129 protein can inhibit immune responses by interacting with host complement C1q binding protein in goat turbinate bone cells

OBJECTIVE

Contagious ecthyma is a severe and highly contagious disease caused by an orf virus (ORFV). The virus is responsible for substantial economic losses in the goat industry and threatens humans. We previously determined the role of ORFV129 protein, one of the five ankyrin-repeat proteins coded by the orf genome, in suppressing the transcription of pro-inflammatory cytokines IL-6, IL-1β and IFN-γ. In the present study, we identified 14 cellular proteins (complement C1q binding protein [C1QBP], MCM7, EIF5A, PKM, SLC6A, TSPAN6, ATP6AP2, GPS1, MMADHC, HSPB6, SLC35B1, MTF1, P3H4, and IL15RA) that interact with ORFV129 using a yeast two-hybrid system in goat turbinate bone cells (GFTCs). The interaction between ORFV129 and (C1QBP), an immune-related protein, was confirmed using immunofluorescence co-localization and co-immunoprecipitation assays. C1QBP overexpression inhibited ORFV replication, whereas the knockdown of C1QBP promoted ORFV replication in GFTCs. Furthermore, ORFV or ORFV129 increased C1QBP expression in GFTCs, indicated that ORFV129-C1QBP interaction might contribute to the ORFV-induced host immune process. In addition, our research showed that ORFV increased the expression of ORFV129, cytokine IL-6, IL-1β and IFN-γ. C1QBP overexpression induced IFN-γ production and reduced IL-6 and IL-1β production. Conversely, C1QBP knockdown induced IL-1β production and reduced IFN-γ and IL-1β production. Moreover, augmentation of ORFV129 expression enhanced the inhibition of the secretion of cytokines IL-6, IL-1β, and IFN-γ induced by the altered expression of C1QBP. These findings suggest different downstream pathways might be involved in regulating different cytokines induced by ORFV129 expression in GFTCs.

Whole-genome sequencing and genomic analysis of Norduz goat (Capra hircus)

Artificial and natural selective breeding of goats has resulted in many different goat breeds all around the world. Norduz goat is one of these breeds, and it is a local goat breed of Turkey. The goats are favorable due to pre-weaning viability and reproduction values compared to the regional breeds. Development in sequencing technologies has let to understand huge genomic structures and complex phenotypes. Until now, such a comprehensive study has not been carried out to understand the genomic structure of the Norduz goats, yet. In the study, the next-generation sequencing was carried out to understand the genomic structure of Norduz goat. Real-time PCR was used to evaluate prominent CNVs in the Norduz goat individuals. Whole genome of the goat was constructed with an average of 33.1X coverage level. In the stringent filtering condition, 9,757,980 SNPs, 1,536,715 InDels, and 290 CNVs were detected in the Norduz goat genome. Functional analysis of high-impact SNP variations showed that the classical complement activation biological process was affected significantly in the goat. CNVs in the goat genome were found in genes related to defense against viruses, immune response, and cell membrane transporters. It was shown that GBP2, GBP5, and mammalian ortholog GBP1, which are INF-stimulated GTPases, were found to be high copy numbers in the goats. To conclude, genetic variations mainly in immunological response processes suggest that Norduz goat is an immunologically improved goat breed and natural selection could take an important role in the genetical improvements of the goats.

Hair Follicle Transcriptome Analysis Reveals Differentially Expressed Genes That Regulate Wool Fiber Diameter in Angora Rabbits

Wool fiber diameter (WFD) is an important index of wool traits and the main determinant of wool quality and value. However, the genetic determinants of fiber diameter have not yet been fully elucidated. Here, coarse and fine wool of Wan strain Angora rabbits and their hair follicle traits were characterized. The results indicated significant differences in the diameters of wool fibers and their hair follicles. The RNA sequencing (RNA-Seq) technique was used to identify differences in gene expression in hair follicles between coarse and fine wool. In total, 2574 differentially expressed genes (DEGs) were found between the two hair follicle groups. Transcription factors, keratin-associated protein (KAP) and keratin (KRT) families, and ECM-related genes may control the structure of fine fibers in rabbits. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that skin development, epidermal cell and keratinocyte differentiation, epithelium development, and Notch and ribosome signaling pathways were significantly enriched, respectively. GSEA further filtered six important pathways and related core genes. PPI analysis also mined functional DEGs associated with hair structure, including LEF1, FZD3, SMAD3, ITGB6, and BMP4. Our findings provide valuable information for researching the molecular mechanisms regulating wool fiber and could facilitate enhanced selection of super-fine wool rabbits through gene-assisted selection in the future.

Molecular Cytogenetics in Domestic Bovids: A Review

The discovery of the Robertsonian translocation (rob) involving cattle chromosomes 1 and 29 and the demonstration of its deleterious effects on fertility focused the interest of many scientific groups on using chromosome banding techniques to reveal chromosome abnormalities and verify their effects on fertility in domestic animals. At the same time, comparative banding studies among various species of domestic or wild animals were found useful for delineating chromosome evolution among species. The advent of molecular cytogenetics, particularly the use of fluorescence in situ hybridization (FISH), has allowed a deeper investigation of the chromosomes of domestic animals through: (a) the physical mapping of specific DNA sequences on chromosome regions; (b) the use of specific chromosome markers for the identification of the chromosomes or chromosome regions involved in chromosome abnormalities, especially when poor banding patterns are produced; (c) better anchoring of radiation hybrid and genetic maps to specific chromosome regions; (d) better comparisons of related and unrelated species by comparative FISH mapping and/or Zoo-FISH techniques; (e) the study of meiotic segregation, especially by sperm-FISH, in some chromosome abnormalities; (f) better demonstration of conserved or lost DNA sequences in chromosome abnormalities; (g) the use of informatic and genomic reconstructions, in addition to CGH arrays, to predict conserved or lost chromosome regions in related species; and (h) the study of some chromosome abnormalities and genomic stability using PCR applications. This review summarizes the most important applications of molecular cytogenetics in domestic bovids, with an emphasis on FISH mapping applications.

Milk yield variation partially attributed to blood oxygen-mediated neutrophil activation in lactating dairy goats

Blood oxygen is an essential component for numerous biological processes of mammalian animals. Milk production of ruminants largely relies on the supply of nutrients, such as glucose, amino acids and fatty acids. To define the regulatory role of blood oxygen availability in regard to milk production, seventy-five healthy Guanzhong dairy goats with similar body weight, days in milk and parities were selected. For each animal, milk yield was recorded and milk sample was collected to determine compositions. Milk vein blood was collected to determine parameters including blood gas, physio-biochemistry and haematology. Another blood sample was prepared for transcriptome and RT-qPCR. Results showed that both pressure of oxygen (pO₂) in the milk vein (positively) and numbers of neutrophils in mammary vein (negatively) were associated with milk yield of the animals. To learn the role of pO₂ in blood cell functionality, twelve animals (six with higher yield (H-group) and six with lower yield (L-group)) from seventy-five goats were selected. Compared with animals in L-group, goats in H-group were higher in pO₂ but lower in pCO₂, lactate, lactate dehydrogenase activity and neutrophil abundance in milk vein, compared with L-group. The blood transcriptome analysis suggested that compared with L-group, animals in H-group were depressed in functionality including neutrophil activation and metabolic pathways including glycolysis, NF-κB and HIF-1. Our result revealed that lower milk production could be associated with neutrophil activation responding to low pO₂ in the mammary vein. In the meantime, we highlighted the potential importance of blood oxygen as a milk yield regulator.

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

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

Hair Follicle Development and Cashmere Traits in Albas Goat Kids

The objectives of this trial were to study the growth and development of hair follicles and cashmere traits in cashmere goats and to provide a theoretical basis for the regulation of secondary hair follicle development and the scientific breeding selection of cashmere goats. Twelve single-fetal female kids were selected as research objects. A long-term tracking plan was created to regularly determine their growth performance, cashmere performance, and hair follicle traits. The results showed no significant difference in live weight after the first and second combing. The cashmere yield and unit yield of the first combing were significantly higher than those of the second combing (p < 0.05). Sections of hair follicles showed that the primary hair follicles are almost fully developed by 1 month, and the secondary hair follicles are fully developed by 5-6 months after birth. The primary hair follicle density (PFD) and secondary hair follicle density (SFD) were highest at birth and decreased within 1 month; and SFD was stable at 5-6 months of age. The change of MSFD took a maximum time of 2 to 3 months. The S:P increase reached its peak at 6 months. BMP4 expression increased with time. FGF2, FGF21 and BMP7 were higher at 3 months old than at the other two-time points. In conclusion, this study determined the total development time of primary and secondary hair follicles from morphology and speculated that FGF2, FGF21, and BMP7 may play a regulatory role in developing secondary hair follicles. Therefore, the period from birth to 6 months of age was the best time to regulate secondary hair follicle development in cashmere goats kids. The traits of the hair follicle and cashmere at 6 months of age could be breeding selection indicators for cashmere goats.

Do non-bovine domestic animals produce A2 milk?: an in silico analysis

A2 milk is an easily digestible product since it has only A2 beta-casein. In cattle, the A1 and A2 alleles are found in the population and the A2 milk is produced from A2A2 animals. Little is known about these alleles in other domestic dairy species. The present study aims to analyze sequence of genetic material available on public databases and quantify the animals genotyped. Eight domestic species were analyzed. There is strong evidence that domestic non-bovine species only carry A2 beta-casein. The data reported here for goats already confirm it due to the large number of animals genotyped as well as buffaloes. It means that they naturally produce A2 milk and no selection must be done. Thus, the fact that A2 milk is easier to digest can be used to add value to dairy product of these species. It helps to conquer new markets. It also improves people's health and breeder profitability.

Comparative and Functional Analysis of miRNAs and mRNAs Involved in Muscle Fiber Hypertrophy of Juvenile and Adult Goats

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate several pathway intermediates and affect the skeletal muscle development in mice, pigs, sheep, and cattle. However, to date, only a small number of miRNAs have been reported in the muscle development of goats. In this report, the longissimus dorsi transcripts of one- and ten-month-old goats were analyzed by sequencing RNAs and miRNAs. The results showed that the ten-month-old Longlin goats had 327 up- and 419 down-regulated differentially expressed genes (DEGs) compared with the one-month-old. In addition, 20 co-up-regulated and 55 co-down-regulated miRNAs involved in the muscle fiber hypertrophy of goats were identified in ten-month-old Longlin and Nubian goats compared with one-month-old. Five miRNA-mRNA pairs (chi-let-7b-3p-MIRLET7A, chi-miR193b-3p-MMP14, chi-miR-355-5p-DGAT2, novel_128-LOC102178119, novel_140-SOD3) involved in the goat skeletal muscle development were identified by miRNA-mRNA negative correlation network analysis. Our results provided new insight into the functional roles of goat muscle-associated miRNAs, allowing a deeper understanding of the transformation of miRNA roles during mammalian muscle development.

MicroRNAs in Ruminants and Their Potential Role in Nutrition and Physiology

The knowledge of how diet choices, dietary supplements, and feed intake influence molecular mechanisms in ruminant nutrition and physiology to maintain ruminant health, is essential to attain. In the present review, we focus on the role of microRNAs in ruminant health and disease; additionally, we discuss the potential of circulating microRNAs as biomarkers of disease in ruminants and the state of technology for their detection, also considering the major difficulties in the transition of biomarker development from bench to clinical practice. MicroRNAs are an inexhaustible class of endogenous non-protein coding small RNAs of 18 to 25 nucleotides that target either the 3' untranslated (UTR) or coding region of genes, ensuring a tight post-transcriptionally controlled regulation of gene expression. The development of new "omics" technologies facilitated a fresh perspective on the nutrition-to-gene relationship, incorporating more extensive data from molecular genetics, animal nutrition, and veterinary sciences. MicroRNAs might serve as important regulators of metabolic processes and may present the inter-phase between nutrition and gene regulation, controlled by the diet. The development of biomarkers holds the potential to revolutionize veterinary practice through faster disease detection, more accurate ruminant health monitoring, enhanced welfare, and increased productivity. Finally, we summarize the latest findings on how microRNAs function as biomarkers, how technological paradigms are reshaping this field of research, and how platforms are being used to identify novel biomarkers. Numerous studies have demonstrated a connection between circulating microRNAs and ruminant diseases such as mastitis, tuberculosis, foot-and-mouth disease, fasciolosis, and metabolic disorders. Therefore, the identification and analysis of a small number of microRNAs can provide crucial information about the stage of a disease, etiology, and prognosis.

Transcriptome sequencing reveals the effects of circRNA on testicular development and spermatogenesis in Qianbei Ma goats

Circular RNAs (circRNAs) play an important role in regulating the mammalian reproductive system, especially testicular development and spermatogenesis. However, their functions in testicular development and spermatogenesis in the Qianbei Ma goat, the Guizhou endemic breed are still unclear. In this study, tissue sectioning and circRNAs transcriptome analysis were conducted to compare the changes of morphology and circular RNAs gene expression profile at four different developmental stages (0Y, 0-month-old; 6Y, 6-month-old; 12Y, 12-month-old; 18Y, 18-month-old). The results showed that the circumferences and area of the seminiferous tubule gradually increased with age, and the lumen of the seminiferous tubule in the testis differentiated significantly. 12,784 circRNAs were detected from testicular tissues at four different developmental stages by RNA sequencing, and 8,140 DEcircRNAs (differentially expressed circRNAs) were found in 0Y vs. 6Y, 6Y vs. 12Y, 12Y vs. 18Y and 0Y vs. 18Y, 0Y vs. 12Y, 6Y vs. 18Y Functional enrichment analysis of the source genes showed that they were mainly enriched in testicular development and spermatogenesis. In addition, the miRNAs and mRNAs associated with DECircRNAs in 6 control groups were predicted by bioinformatics, and 81 highly expressed DECircRNAs and their associated miRNAs and mRNAs were selected to construct the ceRNA network. Through functional enrichment analysis of the target genes of circRNAs in the network, some candidate circRNAs related to testicular development and spermatogenesis were obtained. Such as circRNA_07172, circRNA_04859, circRNA_07832, circRNA_00032 and circRNA_07510. These results will help to reveal the mechanism of circRNAs in testicular development and spermatogenesis, and also provide some guidance for goat reproduction.

Why Is It So Difficult To Have Competent Oocytes from In vitro Cultured Preantral Follicles?

The developmental competence of oocytes is acquired gradually during follicular development, mainly through oocyte accumulation of RNA molecules and proteins that will be used during fertilization and early embryonic development. Several attempts to develop in vitro culture systems to support preantral follicle development up to maturation are reported in the literature, but oocyte competence has not yet been achieved in human and domestic animals. The difficulties to have fertilizable oocytes are related to thousands of mRNAs and proteins that need to be synthesized, long-term duration of follicular development, size of preovulatory follicles, composition of in vitro culture medium, and the need of multi-step culture systems. The development of a culture system that maintains bidirectional communication between the oocyte and granulosa cells and that meets the metabolic demands of each stage of follicle growth is the key to sustain an extended culture period. This review discusses the physiological and molecular mechanisms that determine acquisition of oocyte competence in vitro, like oocyte transcriptional activity, follicle and oocyte sizes, and length and regulation of follicular development in murine, human, and domestic animal species. The state of art of in vitro follicular development and the challenges to have complete follicular development in vitro are also highlighted.

Estimates of genomic heritability and genome-wide association studies for blood parameters in Akkaraman sheep

The aim of this study was to estimate genomic heritability and the impact that genetic backgrounds have on blood parameters in Akkaraman sheep by conducting genome-wide association studies and regional heritability mapping analysis. Genomic heritability estimates for blood parameters ranged from 0.00 to 0.55, indicating that measured phenotypes have a low to moderate heritability. A total of 7 genome- and 13 chromosome-wide significant SNPs were associated with phenotypic changes in 15 blood parameters tested. Accordingly, SCN7A, SCN9A, MYADM-like, CCDC67, ITGA9, MGAT5, SLC19A1, AMPH, NTRK2, MSRA, SLC35F3, SIRT6, CREB3L3, and NAV3 genes as well as three undefined regions (LOC101117887, LOC106991526 and LOC105608461) were suggested as candidates. Most of the identified genes were involved in basic biological processes that are essential to immune system function and cellular growth; specific functions include cellular transport, histone deacetylation, cell differentiation, erythropoiesis, and endocytosis. The top significant SNP for HCT, MCH, and MCHC was found within a genomic region mainly populated by the MYADM-like gene family. This region was previously suggested to be under historical selection pressure in many sheep breeds from various parts of the world. These results have implications on animal breeding program studies due to the effect that the genetic background has on blood parameters, which underlying many productive and wellness related traits.

Effects of exogenous melatonin on expressional differences of immune-related genes in cashmere goats

The interplay between melatonin and immune system is well recognized in humans. The true integration of research on cashmere goat is still far from clear, especially for cashmere goat maintained in wool and cashmere growth. In this study, we applied various approaches to identify the complex regulated network between the immune-related genes and transcription factors (TFs) and to explore the relationship between melatonin and gene expression in cashmere goats. In total, 1,599 and 1756 immune-related genes were found in the blood and skin of cashmere goats, respectively, and 24 differentially expressed immune-related GO terms were highly expressed in blood after melatonin implantation. We studied the melatonin-dependent networks between the TFs and immune-related genes in cashmere goat. The 3 major regulatory networks were interconnected through TFs. The TFs, such as PHF5A, REXO4, STRAP, JUNB, GATAD2A, ZNF710, and VDR, were also expressed in the blood and skin tissue of cashmere goat. In addition, most genes in these networks, such as VDR, JUNB, and Trib3, were involved in WNT pathway, which is related to cashmere wool growth regulation. On the network basis, we developed a knockout mouse model to identify the network interaction. We observed that 8 high-sulfur protein genes, 12 keratin (KRT) genes, and 19 keratin associated protein (KRTAP) genes related to the growth of cashmere wool were almost not expressed in Trib3^−/− rat skin. Our results suggested that the expression of genes related to wool and cashmere growth may be regulated by the interaction network between genes affected by melatonin and immune-related genes. In summary, we outlined some particularly promising ways for future research on immune-related genes of cashmere goats and the role of melatonin in wool and cashmere growth.

De Novo Transcriptome Dataset Generation of the Swamp Buffalo Brain and Non-Brain Tissues

The sequenced data availability opened new horizons related to buffalo genetic control of economic traits and genomic diversity. The visceral organs (brain, liver, etc.) significantly involved in energy metabolism, docility, or social interactions. We performed swamp buffalo transcriptomic profiling of 24 different tissues (brain and non-brain) to identify novel transcripts and analyzed the differentially expressed genes (DEGs) of brain vs. non-brain tissues with their functional annotation. We obtained 178.57 Gb clean transcriptomic data with GC contents 52.77%, reference genome alignment 95.36%, exonic coverage 88.49%. Totally, 26363 mRNAs transcripts including 5574 novel genes were obtained. Further, 7194 transcripts were detected as DEGs by comparing brain vs. non-brain tissues group, of which 3,999 were upregulated and 3,195 downregulated. These DEGs were functionally associated with cellular metabolic activities, signal transduction, cytoprotection, and structural and binding activities. The related functional pathways included cancer pathway, PI3k-Akt signaling, axon guidance, JAK-STAT signaling, basic cellular metabolism, thermogenesis, and oxidative phosphorylation. Our study provides an in-depth understanding of swamp buffalo transcriptomic data including DEGs potentially involved in basic cellular activities and development that helped to maintain their working capacity and social interaction with humans, and also, helpful to disclose the genetic architecture of different phenotypic traits and their gene expression regulation.

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.

Whole-Genome Sequencing Reveals the Genomic Characteristics and Selection Signatures of Hainan Black Goat

Goats have become one of the most adaptive and important livestock species distributed in developing countries in recent years. The Hainan Black goat is a native goat breed of the Hainan region that is generally well-liked by the local population and is thus raised in large numbers. However, the genomic diversity and selective signals of the Hainan Black goat have not been clearly elucidated yet. Therefore, in this study, we performed whole-genome resequencing of 16 Hainan Black goats and compared the results with those of 71 goats of 6 other breeds from different geographic regions. Principal component analysis (PCA) and phylogenetic analysis identified seven lineages for all goats. Hainan Black goats showed the most similarity with Leizhou goats and the least similarity with Boer goats. Selective sweep analysis identified candidate genes associated with various functions, including immune resistance to disease (TNFAIP2 (TNF alpha induced protein 2) and EXOC3L4 (exocyst complex component 3 like 4)), melanin biosynthetic process (CDH15 (cadherin 15), ASIP (agouti signaling protein), and PARD3 (par-3 family cell polarity regulator)), and light sensitivity (CNGB3 (cyclic nucleotide gated channel subunit beta 3) and CNBD1 (cyclic nucleotide binding domain containing 1)), underlying strong selection signatures in Hainan Black goats. The melanin biosynthetic process, circadian entrainment, regulation of cyclic adenosine 3,5-monophosphate (cAMP)-mediated signaling, and the Rap-1 signaling pathway were significantly enriched in Hainan Black and Alashan Cashmere goats. This result may be important for understanding each trait. Selection signature analysis revealed candidate single nucleotide polymorphisms (SNPs) and genes correlated with the traits of Hainan Black goats. Collectively, our results provide valuable insights into the genetic basis of specific traits correlated with the Hainan island climate, artificial selection in certain local goat breeds, and the importance of protecting breed resources.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

Conservation status and historical relatedness of South African communal indigenous goat populations using a genome-wide single-nucleotide polymorphism marker

Indigenous goats form the majority of populations in smallholder, low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support of communal farming. Effective population size (N_e), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Across populations, the current N_e of Gauteng was the lowest at 371 animals, while the historical N_e across populations suggests that the ancestor N_e has decreased by 53.86%, 44.58%, 42.16%, and 41.16% in Free State (FS), North West (NW), Limpopo (LP), and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (F_ROH > 5 Mb) were highest in FS (0.08 ± 0.09) and lowest in the Eastern Cape (EC) (0.02 ± 0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR, and MAPK3 were identified (occurring in over 20% of the samples) in FS (n = 37), GP (n = 42), and NW (n = 2) populations only. The mean length of ROH across populations was 7.76 Mb and ranged from 1.61 Mb in KwaZulu-Natal (KZN) to 98.05 Mb (GP and NW). The distribution of ROH according to their size showed that the majority (n = 1949) of the detected ROH were > 5 Mb in length compared to the other categories. Assuming two hypothetical ancestral populations, the populations from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP share a common origin but have substantial admixture from the EC and NW populations. The findings revealed that the occurrence of high N_e and autozygosity varied largely across breeds in communal indigenous goat populations at recent and ancient events when a genome-wide single-nucleotide polymorphism (SNP) marker was used. The use of Illumina goat SNP50K BeadChip shows that there was a migration route of communal indigenous goat populations from the northern part (LP) of South Africa to the eastern areas of the KZN that confirmed their historical relatedness and coincides with the migration periods of the Bantu nation.

Genome wide association study identifies novel candidate genes for growth and body conformation traits in goats

Pakistan is third largest country in term of goat population with distinct characteristics of breeds and estimated population of 78.2 million. Punjab province has 37% of country's total population with seven important documented goat breeds namely Beetal, Daira Din Pannah, Nachi, Barbari, Teddi, Pahari and Pothwari. There is paucity of literature on GWAS for economically important traits i.e., body weight and morphometric measurements. Therefore, we performed GWAS using 50 K SNP Chip for growth in term of age adjusted body weight and morphometric measurements in order to identify genomic regions influencing these traits among Punjab goat breeds. Blood samples were collected from 879 unrelated animals of seven goat breeds along with data for body weight and morphometric measurements including body length, body height, pubic bone length, heart girth and chest length. Genomic DNA was extracted and genotyped using 50 K SNP bead chip. Association of genotypic data with the phenotypic data was performed using Plink 1.9 software. Linear mixed model was used for the association study. Genes were annotated from Capra hircus genome using assembly ARS1. We have identified a number of highly significant SNPs and respective candidate genes associated with growth and body conformation traits. The functional aspects of these candidate genes suggested their potential role in body growth. Moreover, pleiotropic effects were observed for some SNPs for body weight and conformation traits. The results of current study contributed to a better understanding of genes influencing growth and body conformation traits in goat.

Identification and Regulatory Network Analysis of Genes Related to Reproductive Performance in the Hypothalamus and Pituitary of Angus Cattle

In this study, we explored the gene expression patterns of the pituitary gland and hypothalamus of Angus cows at different growth and developmental stages by deep sequencing and we identified genes that affect bovine reproductive performance to provide new ideas for improving bovine fertility in production practice. We selected three 6-month-old (weaning period), three 18-month-old (first mating period), and three 30-month-old (early postpartum) Angus cattle. The physiological status of the cows in each group was the same, and their body conformations were similar. After quality control of the sequencing, the transcriptome analyses of 18 samples yielded 129.18 GB of clean data. We detected 13,280 and 13,318 expressed genes in the pituitary gland and hypothalamus, respectively, and screened 35 and 50 differentially expressed genes (DEGs) for each, respectively. The differentially expressed genes in both tissues were mainly engaged in metabolism, lipid synthesis, and immune-related pathways in the 18-month-old cows as compared with the 6-month-old cows. The 30-month-old cows presented more regulated reproductive behavior, and pituitary CAMK4 was the main factor regulating the reproductive behavior during this period via the pathways for calcium signaling, longevity, oxytocin, and aldosterone synthesis and secretion. A variant calling analysis also was performed. The SNP inversions and conversions in each sample were counted according to the different base substitution methods. In all samples, most base substitutions were represented by substitutions between bases A and G, and the probability of base conversion exceeded 70%, far exceeding the transversion. Heterozygous SNP sites exceeded 37.68%.

Analysis of World-Scale Mitochondrial DNA Reveals the Origin and Migration Route of East Asia Goats

Despite much attention on the history of goat evolution, information on origin, demographic history, and expansion route remains controversial. To address these questions, we collected 4,189 published goat DNA sequences including 1,228 sequences from 57 breeds in China and 2,961 sequences including 193 goat breeds from 71 other countries and carried out an integrated analysis. We found goat breeds from South China had the highest genetic diversity of lineage B, and subclades B2 only were found in Southwest China, suggesting that lineage B (particularly, subclade B2) probably originated from Southwest China and its surrounding areas. In addition, in this study, we found that lineage A from South China also presented higher genetic diversity and earlier expansion time (10, 606 years ago), even earlier than breeds from the Middle East. Hence, we speculated that South China and surrounding areas were the origin of lineage B and also the transportation hub for lineage A spreading to North China and Southwest Asia. Furthermore, according to the analysis of correlation between genetic differentiation value λ1 and λ2 and geographical distance, we further confirmed two phases of migration in goat breeds of North China. These results will contribute to a better understanding of the origin and migration history of domestic goat.

Ancient mitochondrial and modern whole genomes unravel massive genetic diversity loss during near extinction of Alpine ibex

Population bottlenecks can have dramatic consequences for the health and long-term survival of a species. Understanding of historic population size and standing genetic variation prior to a contraction allows estimating the impact of a bottleneck on the species genetic diversity. Although historic population sizes can be modelled based on extant genomics, uncertainty is high for the last 10-20 millenia. Hence, integrating ancient genomes provides a powerful complement to retrace the evolution of genetic diversity through population fluctuations. Here, we recover 15 high-quality mitogenomes of the once nearly extinct Alpine ibex spanning 8601 BP to 1919 CE and combine these with 60 published modern whole genomes. Coalescent demography simulations based on modern whole genomes indicate population fluctuations coinciding with the last major glaciation period. Using our ancient and historic mitogenomes, we investigate the more recent demographic history of the species and show that mitochondrial haplotype diversity was reduced to a fifth of the pre-bottleneck diversity with several highly differentiated mitochondrial lineages having co-existed historically. The main collapse of mitochondrial diversity coincides with elevated human population growth during the last 1-2 kya. After recovery, one lineage was spread and nearly fixed across the Alps due to recolonization efforts. Our study highlights that a combined approach integrating genomic data of ancient, historic and extant populations unravels major long-term population fluctuations from the emergence of a species through its near extinction up to the recent past.

A draft genome of Drung cattle reveals clues to its chromosomal fusion and environmental adaptation

Drung cattle (Bos frontalis) have 58 chromosomes, differing from the Bos taurus 2n = 60 karyotype. To date, its origin and evolution history have not been proven conclusively, and the mechanisms of chromosome fusion and environmental adaptation have not been clearly elucidated. Here, we assembled a high integrity and good contiguity genome of Drung cattle with 13.7-fold contig N50 and 4.1-fold scaffold N50 improvements over the recently published Indian mithun assembly, respectively. Speciation time estimation and phylogenetic analysis showed that Drung cattle diverged from Bos taurus into an independent evolutionary clade. Sequence evidence of centromere regions provides clues to the breakpoints in BTA2 and BTA28 centromere satellites. We furthermore integrated a circulation and contraction-related biological process involving 43 evolutionary genes that participated in pathways associated with the evolution of the cardiovascular system. These findings may have important implications for understanding the molecular mechanisms of chromosome fusion, alpine valleys adaptability and cardiovascular function.

Candidate Genes and Their Expressions Involved in the Regulation of Milk and Meat Production and Quality in Goats (Capra hircus)

Simple Summary

During the present decade, highly selected caprine farming has increased in popularity due to the hardiness and adaptability inherent to goats. Recent advances in genetics have enabled the improvement in goat selection efficiency. The present review explores how genetic technologies have been applied to the goat-farming sector in the last century. The main candidate genes related to economically relevant traits are reported. The major source of income in goat farming derives from the sale of milk and meat. Consequently, yield and quality must be specially considered. Meat-related traits were evaluated considering three functional groups (weight gain, carcass quality and fat profile). Milk traits were assessed in three additional functional groups (milk production, protein and fat content).

Abstract

Despite their pivotal position as relevant sources for high-quality proteins in particularly hard environmental contexts, the domestic goat has not benefited from the advances made in genomics compared to other livestock species. Genetic analysis based on the study of candidate genes is considered an appropriate approach to elucidate the physiological mechanisms involved in the regulation of the expression of functional traits. This is especially relevant when such functional traits are linked to economic interest. The knowledge of candidate genes, their location on the goat genetic map and the specific phenotypic outcomes that may arise due to the regulation of their expression act as a catalyzer for the efficiency and accuracy of goat-breeding policies, which in turn translates into a greater competitiveness and sustainable profit for goats worldwide. To this aim, this review presents a chronological comprehensive analysis of caprine genetics and genomics through the evaluation of the available literature regarding the main candidate genes involved in meat and milk production and quality in the domestic goat. Additionally, this review aims to serve as a guide for future research, given that the assessment, determination and characterization of the genes associated with desirable phenotypes may provide information that may, in turn, enhance the implementation of goat-breeding programs in future and ensure their sustainability.

Knockdown of CDR1as Decreases Differentiation of Goat Skeletal Muscle Satellite Cells via Upregulating miR-27a-3p to Inhibit ANGPT1

Myogenesis is a complex process controlled by several coding and non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) that are known to function as endogenous microRNAs (miRNAs) sponges. Cerebellar Degeneration-Related protein 1 antisense (CDR1as) is the most spotlighted circRNA that is known as an miR-7 sponge, which has bloomed circRNAs’ research in animal disease and physiology. Here, we screened for miRNAs and mRNA associated with CDR1as and further characterized their regulatory function during muscle differentiation. We found that a total of 43 miRNAs (including miR-107-3p, miR-125b-5p, miR-140-5p, miR-29a-3p, and miR-27a-3p upregulated) and 789 mRNAs (including ANGPT1, E2F2, CCN1, FGFR1, and MEF2C downregulated) were differentially expressed in goat skeletal muscle satellite cells (SMSCs). Further, knockdown of CDR1as and ANGPT1 inhibited SMSCs differentiation. miR-27a-3p was differentially upregulated after the knockdown of CDR1as in SMSCs. Overexpressed miR-27a-3p decreased SMSCs differentiation. Via RNAhybrid and luciferase, miR-27a-3p was identified to regulate ANGPT1. We discovered that miR-27a-3p has an inverse relationship with CDR1as and decreases the expression level of ANGPT1 during SMSCs differentiation. In summary, our study demonstrates that siCDR1as inhibits myoblast differentiation by downregulating ANGPT1 mRNA via miR-27a-3p in SMSCs.

Accuracy of Genomic Selection for Important Economic Traits of Cashmere and Meat Goats Assessed by Simulation Study

Genomic selection in plants and animals has become a standard tool for breeding because of the advantages of high accuracy and short generation intervals. Implementation of this technology is hindered by the high cost of genotyping and other factors. The aim of this study was to determine an optional marker density panel and reference population size for using genomic selection of goats, with speculation on the number of QTLs that affect the important economic traits of goats. In addition, the effect of buck population size in the reference population on the accuracy of genomic estimated breeding value (GEBV) was discussed. Based on the previous genetic evaluation results of Inner Mongolia White Cashmere Goats, live body weight (LBW, h2 = 0.11) and fiber diameter (FD, h2 = 0.34) were chosen to perform genomic selection in this study. Reasonable genome parameters and generation transmission processes were set, and phenotypic and genotype data of the two traits were simulated. Then, different sizes of the reference population and validation population were selected from progeny. The GEBVs were obtained by six methods, including GBLUP (Genomic Best Linear Unbiased Prediction), ssGBLUP (Single Step Genomic Best Linear Unbiased Prediction), BayesA, BayesB, Bayesian ridge regression, and Bayesian LASSO. The correlation coefficient between the predicted and realized phenotypes from simulation was calculated and used as a measure of the accuracy of GEBV in each trait. The results showed that the medium marker density Panel (45 K) could be used for genomic selection in goats, which can ensure the accuracy of the GEBV. The reference population size of 1,500 can achieve greater genetic progress in genomic selection for fiber diameter and live body weight in goats by comparing with the population size below this level. The accuracy of the GEBV for live body weight and fiber diameter was better when the number of QTLs was 100 and 50, respectively. Additionally, the accuracy of GEBV was discovered to be good when the buck population size was up to 200. Meanwhile, the accuracy of the GEBV for medium heritability traits (FDs) was found to be higher than the accuracy of the GEBV for low heritability traits (LBWs). These findings will provide theoretical guidance for genomic selection in goats by using real data.

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.

Defining the caprine γδ T cell WC1 multigenic array and evaluation of its expressed sequences and gene structure conservation among goat breeds and relative to cattle

Workshop cluster 1 (WC1) molecules are part of the scavenger receptor cysteine-rich (SRCR) superfamily and act as hybrid co-receptors for the γδ T cell receptor and as pattern recognition receptors for binding pathogens. These members of the CD163 gene family are expressed on γδ T cells in the blood of ruminants. While the presence of WC1⁺ γδ T cells in the blood of goats has been demonstrated using monoclonal antibodies, there was no information available about the goat WC1 gene family. The caprine WC1 multigenic array was characterized here for number, structure and expression of genes, and similarity to WC1 genes of cattle and among goat breeds. We found sequence for 17 complete WC1 genes and evidence for up to 30 SRCR a1 or d1 domains which represent distinct signature domains for individual genes. This suggests substantially more WC1 genes than in cattle. Moreover, goats had seven different WC1 gene structures of which 4 are unique to goats. Caprine WC1 genes also had multiple transcript splice variants of their intracytoplasmic domains that eliminated tyrosines shown previously to be important for signal transduction. The most distal WC1 SRCR a1 domains were highly conserved among goat breeds, but fewer were conserved between goats and cattle. Since goats have a greater number of WC1 genes and unique WC1 gene structures relative to cattle, goat WC1 molecules may have expanded functions. This finding may impact research on next-generation vaccines designed to stimulate γδ T cells.