Coloration in Equine: Overview of Candidate Genes Associated with Coat Color Phenotypes.

Variation in coat color among equids has attracted significant interest in genetics and breeding research. The range of colors is primarily determined by the type, concentration, and distribution of melanin pigments, with the balance between eumelanin and pheomelanin influenced by numerous genetic factors. Advances in genomic and sequencing technologies have enabled the identification of several candidate genes that influence coat color, thereby clarifying the genetic basis of these diverse phenotypes. In this review, we concisely categorize coat coloration in horses and donkeys, focusing on the biosynthesis and types of melanin involved in pigmentation. Moreover, we highlight the regulatory roles of some key candidate genes, such as MC1R, TYR, MITF, ASIP, and KIT, in coat color variation. Moreover, the review explores how coat color relates to selective breeding and specific equine diseases, offering valuable insights for developing breeding strategies that enhance both the esthetic and health aspects of equine species.

Elucidating the Role of Transcriptomic Networks and DNA Methylation in Collagen Deposition of Dezhou Donkey Skin.

DNA methylation represents a predominant epigenetic modification with broad implications in various biological functions. Its role is particularly significant in the process of collagen deposition, a fundamental aspect of dermal development in donkeys. Despite its critical involvement, the mechanistic insights into how DNA methylation influences collagen deposition in donkey skin remain limited. In this study, we employed whole genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) to investigate the epigenetic landscape and gene expression profiles in the dorsal skin tissues of Dezhou donkeys across three developmental stages: embryonic (YD), juvenile (2-year-old, MD), and mature (8-year-old, OD). Our analysis identified numerous differentially methylated genes that play pivotal roles in skin collagen deposition and overall skin maturation, including but not limited to COL1A1, COL1A2, COL3A1, COL4A1, COL4A2, GLUL, SFRP2, FOSL1, SERPINE1, MMP1, MMP2, MMP9, and MMP13. Notably, we observed an inverse relationship between gene expression and DNA methylation proximal to transcription start sites (TSSs), whereas a direct correlation was detected in regions close to transcription termination sites (TTSs). Detailed bisulfite sequencing analyses of the COL1A1 promoter region revealed a low methylation status during the embryonic stage, correlating with elevated transcriptional activity and gene expression levels. Collectively, our findings elucidate key genetic markers associated with collagen deposition in the skin of Dezhou donkeys, underscoring the significant regulatory role of DNA methylation. This research work contributes to the foundational knowledge necessary for the genetic improvement and selective breeding of Dezhou donkeys, aiming to enhance skin quality attributes.

Comprehensive transcriptomic analysis unveils the interplay of mRNA and LncRNA expression in shaping collagen organization and skin development in Dezhou donkeys.

The primary focus of donkey hide gelatin processing lies in the dermal layer of donkey hide due to its abundant collagen content. However, the molecular mechanism involved in collagen organization and skin development in donkey skin tissue across various developmental stages remains incomplete. The current study aims to investigate the transcriptomic screening of lncRNAs and mRNA associated with skin development and collagen organization across different ages in Dezhou donkeys' skin. In the pursuit of this objective, we used nine skin tissue samples obtained from Dezhou donkeys at various ages including 8-month fetal stage, followed by 2 and 8 years. RNA-seq analysis was performed for the transcriptomic profiling of differentially expressed genes (DEGs) and lncRNAs associated with skin development in different age groups. Our investigation revealed the presence of 6,582, 6,455, and 405 differentially expressed genes and 654, 789, and 29 differentially expressed LncRNAs within the skin tissues of Dezhou donkeys when comparing young donkeys (YD) vs. middle-aged donkeys (MD), YD vs. old donkeys (OD), and MD vs. OD, respectively. Furthermore, we identified Collagen Type I Alpha 1 Chain (COL1A1), Collagen Type III Alpha 1 Chain (COL3A1), and Collagen Type VI Alpha 5 Chain (COL6A5) as key genes involved in collagen synthesis, with COL1A1 being subject to cis-regulation by several differentially expressed LncRNAs, including ENSEAST00005041187, ENSEAST00005038497, and MSTRG.17248.1, among others. Interestingly, collagen organizational and skin development linked pathways including Protein digestion and absorption, metabolic pathways, Phosphatidylinositol 3-Kinase-Protein Kinase B signaling pathway (PI3K-Akt signaling pathway), Extracellular Matrix-Receptor Interaction (ECM-receptor interaction), and Relaxin signaling were also reported across different age groups in Dezhou donkey skin. These findings enhance our comprehension of the molecular mechanisms underlying Dezhou donkey skin development and collagen biosynthesis and organization, thus furnishing a solid theoretical foundation for future research endeavors in this domain.

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits.

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology.

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.

[Activity and transcriptional regulatory elements of the promoter in Arctic fox (Vulpes lagopus) ß-defensin103 gene].

The aim of this study was to screen the active regions and transcription factor binding sites in the promoter of the CBD103 gene related to Arctic fox coat color, and to provide a basis for revealing the molecular genetic mechanism of CBD103 gene regulating the coat color formation. The 5'-flanking region fragment 2 123 bp of Arctic fox CBD103 gene was cloned, and 4 truncated promoter reporter vectors of different lengths were constructed. The promoter activity was detected by the dual-luciferase reporter assay system. Point mutations were performed on the 3 predicted specificity protein 1 (Sp1) transcription factor binding sites in the highest promoter active region, and 3 mutant vectors were constructed. The activity was then detected by the dual-luciferase reporter assay system. The results showed that the region 1 656 (-1 604/+51) had the highest activity in the 4 truncated promoters of different lengths, and the promoter activity of the three mutant vectors constructed in this region were significantly lower than that of the wild type (fragment 1 656). The region of -1 604 /+51 was the core promoter region of CBD103 gene in Arctic fox and -1 552/-1 564, -1 439/-1 454 and -329/-339 regions were positive regulatory regions. This study successfully obtained the core promoter region and positive regulation regions of the Arctic fox CBD103 gene, which laid a foundation for further study on the molecular genetic mechanism of this gene regulating Arctic fox coat color.

Correction to: Hair follicles transcriptome profiles in Bashang long-tailed chickens with different plumage colors.

The words 'hair follicles' was replaced with 'feather follicles' in the title and the main text.

Feather follicles transcriptome profiles in Bashang long-tailed chickens with different plumage colors.

Despite the rich variety in plumage color found in nature, genetic studies on how feather follicles affect pigmentation are often limited to animals that have black and white pigment. To test how gene expression influences plumage color, transcriptomes of chicken feather follicles with white, black, hemp, reed catkins, silvery grey, and landscape plumage colors were generated using Illumina sequencing. We generated six RNA-Seq libraries with over 25 million paired-end clean reads per library with percentage of paired-end clean reads ranging from 96.73 to 96.98%. 78% of the reads mapped to the chicken genome, and approximately 70% of the reads were mapped to exons and 6% mapped to introns. Transcriptomes of feather follicles producing hemp and land plumage were similar, but these two showed moderate differences compared with gray and reed colored plumage. The black and white follicle transcriptomes were most divergent from the other colors. We identified several candidate genes, including GPNMB, PMEL, TYRP1, GPR143, OCA2, SOX10, SLC45A2, KRT75, and TYR. All of these genes are known to induce pigment formation in mice. White feathers result from the lack of pigment formation, and our results suggest that the white chickens due to the recessive insertion mutation of TYR. The formation of black area size and color depth may be due to the expression levels of GPNMB, PMEL, TYRP1, GPR143, OCA2, SOX10, SLC45A2, KRT75, and TYR. The GO analysis of the differentially expressed genes (DEGs) revealed that DEGs in our transcriptome analysis were enriched in cytoskeleton and cell structure related pathways. The black plumage transcriptome showed significant differences in melanogenesis, tyrosine metabolism, and riboflavin metabolism compared with transcriptomes of other plumage colors. The transcriptome profiles of the different chicken plumage colors provide a valuable resource to understand how gene expression influences plumage color, and will be an important resource for identifying candidate genes in breeding programs.

Genome-wide differential expression of long noncoding RNAs and mRNAs in ovarian follicles of two different chicken breeds.

Bashang long-tail chickens are an indigenous breed with dual purpose in China (meat and eggs) but have low egg laying performance. To improve the low egg laying performance, a genome-wide analysis of mRNAs and long noncoding RNAs (lncRNAs) from Bashang long-tail chickens and Hy-Line brown layers was performed. A total of 16,354 mRNAs and 8691 lncRNAs were obtained from ovarian follicles. Between the breeds, 160 mRNAs and 550 lncRNAs were found to be significantly differentially expressed. Integrated network analysis suggested some differentially expressed genes were involved in ovarian follicular development through oocyte meiosis, progesterone-mediated oocyte maturation, and cell cycle. The impact of lncRNAs on cis and trans target genes, indicating some lncRNAs may play important roles in ovarian follicular development. The current results provided a catalog of chicken ovarian follicular lncRNAs and genes for further study to understand their roles in regulation of egg laying performance.

[Identification of the core promoter of the pmel gene of Bashang long-tail chickens].

To explore the activity of the pmel core promoter of Bashang long-tail chickens, we constructed dual-luciferase expression vectors and transiently transfected into DF1 cells with Lipofectamine 2000. We measured the luciferase activity with the dual-luciferase detection kit. The 1 268 bp fragment in 5-flanking region of the pmel gene in Bashang long-tail chickens was cloned. The region from -1 200 bp to +68 bp included 2 CpG islands and multiple transcription factor binding sites. We constructed 9 expression vectors with different promoter regions and a mutant vector of the core promoter region of the pmel gene of Bashang long-tail chickens. The core promoter region from -840 bp to +68 bp was identified in the pmel gene. The region from -590 to -525 bp negatively regulated the pmel gene during the transcription process. The -840--590 bp and -525--266 bp regions were positive regulatory regions. The polymorphic sites (-456, -435, -410, -374 and -341) had a significant effect on the promoter activity of the pmel gene.

The Emerging Role of Zfp217 in Adipogenesis.

Zinc finger protein 217 (Zfp217), a member of the krüppel-type zinc finger protein family, plays diverse roles in cell differentiation and development of mammals. Despite extensive research on the functions of Zfp217 in cancer, pluripotency and reprogramming, its physiological roles in adipogenesis remain unknown. Our previous RNA sequencing data suggest the involvement of Zfp217 in adipogenesis. In this study, the potential function of Zfp217 in adipogenesis was investigated through bioinformatics analysis and a series of experiments. The expression of Zfp217 was found to be gradually upregulated during the adipogenic differentiation in C3H10T1/2 cells, which was consistent with that of the adipogenic marker gene Pparg2. Furthermore, there was a positive, significant relationship between Zfp217 expression and adipocyte differentiation. It was also observed that Zfp217 could not only trigger proliferative defect in C3H10T1/2 cells, but also interact with Ezh2 and suppress the downstream target genes of Ezh2. Besides, three microRNAs (miR-503-5p, miR-135a-5p and miR-19a-3p) which target Zfp217 were found to suppress the process of adipogenesis. This is the first report showing that Zfp217 has the capacity to regulate adipogenesis.

MicroRNA-215 impairs adipocyte differentiation and co-represses FNDC3B and CTNNBIP1.

MicroRNAs (miRNAs) are small ∼22 nucleotide regulatory RNAs that regulate the stability and translation of cognate mRNAs. MiRNAs participate in the regulation of adipogenesis, and identification of the full repertoire of miRNAs expressed in adipose tisse is likely to improve our understanding of adipose tissue growth and development significantly. In the present study, miR-215-5p was found to inhibit adipocyte differentiation of 3T3-L1 cells. Moreover, fibronectin type III domain containing 3B (FNDC3B) and catenin, beta interacting protein 1 (CTNNBIP1) were found to be direct targets of miR-215-5p. Further studies in mouse 3T3-L1 cell-line suggests that miR-215-5p is a negative regulator of adipocyte differentiation through post-transcriptional regulation of FNDC3B and CTNNBIP1 during early adipogenesis.

STAT5a promotes the transcription of mature mmu-miR-135a in 3T3-L1 cells by binding to both miR-135a-1 and miR-135a-2 promoter elements.

Despite extensive research on the role of miR-135a in biological processes, very little attention has been paid to the regulation of its transcription. We have previously reported that miR-135a suppresses 3T3-L1 preadipocyte differentiation and adipogenesis by directly targeting the adenomatous polyposis coli (APC) gene and activating the canonical Wnt/ß-catenin signaling pathway, but the regulatory elements that regulate the expression of the two isoforms of miR-135a (miR-135a-1 and miR-135a-2) remain poorly understood. Here, by using deletion analysis, we predicted two binding sites (-874/-856 and -2020/-2002) for the transcription factor Signal Transducers and Activators of Transcription 5a (STAT5a) within the core promoters of miR-135a-1 and miR-135a-2 (-1128/-556 and -2264/-1773), and the subsequent site-directed mutagenesis indicated that the two STAT5a binding sites regulated the activity of the miR-135a-1 and miR-135a-2 promoters. The binding of STAT5a to the miR-135a-1/2 core promoters in vitro and in cell culture was identified by electrophoretic mobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP) assays. Overexpression and RNAi knockdown of STAT5a showed that the transcription factor regulated the endogenous miR-135a expression. Additionally, The expression time frame of STAT5a and APC indicated a potential negative feedback between them. In sum, the overall results from this study indicate that STAT5a regulates miR-135a transcription by binding to both miR-135a-1 and miR135a-2 promoter elements and the findings provide novel insights into the molecular regulatory mechanisms of miR-135a during adipogenesis.

miR-135a-5p inhibits 3T3-L1 adipogenesis through activation of canonical Wnt/ß-catenin signaling.

MicroRNAs are endogenous, conserved, and non-coding small RNAs that function as post-transcriptional regulators of fat development and adipogenesis. Adipogenic marker genes, such as CCAAT/enhancer binding protein α (Cebpa), peroxisome proliferator-activated receptor γ (Pparg), adipocyte fatty acid binding protein (Ap2), and fatty acid synthase (Fas), are regarded as the essential transcriptional regulators of preadipocyte differentiation and lipid storage in mature adipocytes. Canonical Wnt/ß-catenin signaling is recognized as a negative molecular switch during adipogenesis. In the present work we found that miR-135a-5p is markedly downregulated during the process of 3T3-L1 preadipocyte differentiation. Overexpression of miR-135a-5p impairs the expressions of adipogenic marker genes as well as lipid droplet accumulation and triglyceride content, indicating the importance of miR-135a-5p for adipogenic differentiation and adipogenesis. Further studies show that miR-135a-5p directly targets adenomatous polyposis coli (Apc), contributes to the translocation of ß-catenin from cytoplasm to nucleus, and then activates the expressions of cyclin D1 (Ccnd1) and Cmyc, indicating the induction of canonical Wnt/ß-catenin signaling. In addition, inhibition of APC with siRNA exhibits the same effects as overexpression of miR-135a-5p. Our findings demonstrate that miR-135a-5p suppresses 3T3-L1 preadipocyte differentiation and adipogenesis through the activation of canonical Wnt/ß-catenin signaling by directly targeting Apc. Taken together, these results offer profound insights into the adipogenesis mechanism and the development of adipose tissue.

MicroRNAs: emerging roles in adipogenesis and obesity.

Obesity is a serious health problem worldwide associated with an increased risk of life-threatening diseases such as type 2 diabetes, atherosclerosis, and certain types of cancer. Understanding the molecular basis of adipogenesis and fat cell development in obesity is essential to identify new biomarkers and therapeutic targets for the development of anti-obesity drugs. Recent computational and experimental studies have shown that microRNAs (miRNAs) appear to play regulatory roles in many biological processes associated with obesity, including adipocyte differentiation and lipid metabolism. In addition, many miRNAs are dysregulated in metabolic tissues from obese animals and humans, which potentially contributes to the pathogenesis of obesity-associated complications. The discovery of circulating miRNAs has highlighted their potential as both endocrine signaling molecules and disease markers. The potential of miRNA based therapeutics targeting obesity is highlighted as well as recommendations for future research which could lead to a breakthrough in the treatment of obesity.

MiR-224 impairs adipocyte early differentiation and regulates fatty acid metabolism.

MicroRNAs (miRNAs) are small ~22 nucleotide regulatory RNAs that regulate the stability and translation of cognate messenger RNAs (mRNAs). MicroRNAs participate in the regulation of adipogenesis and identification of the full repertoire of MicroRNAs expressed in adipose tissue is likely to improve our understanding of adipose tissue growth and development significantly. In the present study, it is found that miR-224-5p abundance decreases first and then increases during adipogenesis of 3T3-L1 cells. And early growth response 2 (EGR2) and Acyl-CoA synthetase long-chain family member 4 (ACSL4) are direct targets of miR-224-5p. Further studies in mouse 3T3-L1 cell-line shows that miR-224-5p is a novel negative regulator of adipocyte differentiation through post-transcriptional regulation of early growth response 2 during early adipogenesis. Furthermore, miR-224-5p could regulate fatty acid metabolism through Acyl-CoA synthetase long-chain family member 4 at terminal differentiation. It indicates that miR-224 plays different roles on different stages of adipogenesis.

Molecular cloning, tissue expression, and analysis with genome DNA methylation of porcine LSD1 gene.

Lysine-specific demethylase 1 (LSD1) functioned as a demethyl methylase gene, underlying a wide range of biological processes, including cancer, cell apoptosis, differentiation, and development. To further understand the functions of the porcine LSD1 gene, we first obtained cDNA sequence of porcine LSD1 gene, using in silico cloning method. We further found that the porcine LSD1 gene has two transcripts, in which cDNA sequences are 2,716 and 2,656 bp, ORF are 2,622 and 2,562 bp, respectively. Then, RT-PCR analysis showed that the LSD1 gene is expressed in various tissues and relatively higher in the tissues of ovary, kidney, and spleen. Besides, the LSD1 gene was expressed higher in the growth nonage and peaked at 3 days in muscle tissue. Meanwhile, the expression of two transcript variants of LSD1 gene presented the same change trend. Besides, the level of DNA methylation was approximately fourfold higher in a 3-day muscle than in an old pig (180 days), significantly positive related to the gene expression of LSD1 (R = 0.9362, P < 0.05), and declined with growing age. Cloning, expression pattern, and analysis with genome DNA methylation of porcine LSD1 gene laid a foundation to clarify the molecular mechanisms of porcine growth and development and also for further work on animal breeding.