Analyses of histological and transcriptome differences in the skin of short-hair and long-hair rabbits

Molecular insights into Pashmina fiber production: comparative skin transcriptomic analysis of Changthangi goats and sheep

Ladakh, one of the highest inhabited regions globally, hosts the unique Changthangi goat, renowned for producing Pashmina, the world's most luxurious natural fiber. In comparison, the fiber derived from Changthangi sheep is considered next only to Pashmina. This research endeavors to compare the skin transcriptome profiles of Changthangi goats and Changthangi sheep, aiming to discern the molecular determinants behind the recognition of Changthangi goats as the source of Pashmina. Drawing upon previously conducted studies, a collective of 225 genes correlated with fiber characteristics were extracted from the differentially expressed genes noticed between the two species (p-value of ≤ 0.05 and a log2 fold change of ≥ 1.5). These genes were analyzed using DAVID software to understand their biological functions and to identify enriched KEGG and Reactome pathways. The protein-protein interaction networks were constructed using Cytoscape, cytoHubba, and STRING to focus on key genes and infer their biological significance. Comparative transcriptome analysis revealed significantly higher expression of genes involved in signaling pathways like Wnt, MAPK, PI3K-Akt, Hedgehog, associated with fiber development and quality in Changthangi goats. These pathways play crucial roles in hair follicle (HF) formation, maintenance of epidermal stem cells, and fiber characteristics. Findings also highlight the enrichment of cell adhesion molecules and ECM-receptor interaction, emphasizing their roles in HF structure, growth, and signaling. This investigation offers an in-depth understanding of the molecular intricacies governing Pashmina production in Changthangi goats, providing valuable insights into their unique genetic makeup and underlying mechanisms influencing the exceptional quality of Pashmina fibers.

Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differential isoform in skin tissues of different hair length Yak

BACKGROUND

The hair follicle development process is regulated by sophisticated genes and signaling networks, and the hair grows from the hair follicle. The Tianzhu white yak population exhibits differences in hair length, especially on the forehead and shoulder region. However, the genetic mechanism is still unclear. Isoform sequencing (Iso-seq) technology with advantages in long reads sequencing. Hence, we combined the Iso-seq and RNA-seq methods to investigate the transcript complexity and difference between long-haired yak (LHY) and normal-haired yak (NHY).

RESULTS

The hair length measurement result showed a significant difference between LHY and NHY on the forehead and the shoulder (P-value < 0.001). The skin samples from the forehead and the shoulder of LHY and NHY were pooled for isoform sequencing (Iso-seq). We obtained numerous long transcripts, including novel isoforms, long non-coding RNA, alternative splicing events, and alternative polyadenylation events. Combined with RNA-seq data, we performed differential isoforms (DEIs) analysis between LHY and NHY. We found that some hair follicle and skin development-related DEIs, like BMP4, KRT2, IGF2R, and COL1A2 in the forehead skin; BMP1, KRT1, FGF5, COL2A1, and IGFBP5 in the shoulder skin. Enrichment analysis revealed that DEIs in both two comparable groups significantly participated in skin and hair follicle development-related pathways, like ECM-receptor interaction, focal adhesion, and PI3K-Akt signaling pathways. The results indicated that the hair follicle development of Tianzhu white yak may influence the hair length difference. Besides, the protein-protein interaction (PPI) network of DEIs showed COL2A1 and COL3A1 exhibited a high degree of centrality, and these two genes were suggested as potential candidates for the hair length growth of Tianzhu white yak.

CONCLUSIONS

The results provided a comprehensive analysis of the transcriptome complexity and identified differential transcripts that enhance our understanding of the molecular mechanisms underlying the variation in hair length growth in Tianzhu white yak.

Single nucleotide polymorphisms in the KRT82 promoter region modulate irregular thickening and patchiness in the dorsal skin of New Zealand rabbits

BACKGROUND

While rabbits are used as models in skin irritation tests, the presence of irregular patches and thickening on the dorsal skin can affect precise evaluation. In this study, genes associated with patchiness or non-patchiness on the dorsal skin of New Zealand rabbits were investigated to identify potential regulators of the patchiness phenotype.

RESULTS

The results showed that parameters associated with hair follicles (HFs), such as HF density, skin thickness, and HF depth, were augmented in rabbits with the patchiness phenotype relative to the non-patchiness phenotype. A total of 592 differentially expressed genes (DEGs) were identified between the two groups using RNA-sequencing. These included KRT72, KRT82, KRT85, FUT8, SOX9, and WNT5B. The functions of the DEGs were investigated by GO and KEGG enrichment analyses. A candidate gene, KRT82, was selected for further molecular function verification. There was a significant positive correlation between KRT82 expression and HF-related parameters, and KRT82 overexpression and knockdown experiments with rabbit dermal papilla cells (DPCs) showed that it regulated genes related to skin and HF growth and development. Investigation of single nucleotide polymorphisms (SNPs) in the exons and promoter region of KRT82 identified four SNPs in the promoter region but none in the exons. The G.-631G > T, T.-696T > C, G.-770G > T and A.-873 A > C alleles conformed to the Hardy - Weinberg equilibrium, and three identified haplotypes showed linkage disequilibrium. Luciferase reporter assays showed that the core promoter region of KRT82 was located in the - 600 to - 1200 segment, in which the four SNPs were located.

CONCLUSIONS

The morphological characteristics of the patchiness phenotype were analyzed in New Zealand rabbits and DEGs associated with this phenotype were identified by RNA-sequencing. The biological functions of the gene KRT82 associated with this phenotype were analyzed, and four SNPs were identified in the promoter region of the gene. These findings suggest that KRT82 may be a potential biomarker for the breeding of experimental New Zealand rabbits.

Identification of N6-Methyladenosine-Related Factors and the Prediction of the Regulatory Mechanism of Hair Follicle Development in Rex and Hycole Rabbits

Hair follicle development directly affects the development of the rabbit fur industry. The growth and development of a hair follicle is modified and regulated by many genes and mechanisms. M6A is an important RNA modification. However, there are few studies on the effects of the regulation of m6A on hair follicle growth and development. In this study, hematoxylin-eosin (HE) staining was used to explore the difference in hair follicle development between Rex rabbits and Hycole rabbits, and we performed m6A sequencing to identify the key genes with m6A modification in hair follicle growth. The results showed that the hair length, coarse hair percentage, primary hair follicle ratio, and skin thickness of Hycole rabbits were significantly higher than those of Rex rabbits. However, the proportion of secondary hair follicles in Hycole rabbits was significantly lower than that in Rex rabbits. In addition, we found five differential methylases, 20 differential genes, and 24 differential signaling pathways related to hair growth and development. The results of the Sankey diagram showed that 12 genes were related to 13 signal pathways. Finally, we found that five methylases regulated the development of hair follicles through differential genes/signal pathways. These findings laid a molecular foundation for the function of m6A modification in hair development.

A homozygous missense mutation in the fibroblast growth factor 5 gene is associated with the long-hair trait in Angora rabbits

BACKGROUND

Rabbits are well-domesticated animals. As a crucial economic animal, rabbit has been successfully bred into wool-use, meat-use and fur-use breeds. Hair length is one of the most economically important traits affecting profitability in wool rabbits. In this study, to identify selection signatures with the long-hair trait, whole-genomic resequencing of long-haired rabbits (Angora rabbits) and short-haired rabbits (Rex and New Zealand rabbits) was performed.

RESULTS

By genome-wide selective sweeping analysis based on population comparison, we identified a total of 5.85 Mb regions (containing 174 candidate genes) with strong selection signals. Six of these genes (Dusp1, Ihh, Fam134a, Map3k1, Spata16, and Fgf5) were enriched in the MAPK signalling and Hedgehog signalling pathways, both of which are closely associated with hair growth regulation. Among these genes, Fgf5 encodes the FGF5 protein, which is a well-established regulator of hair growth. There was a nonsynonymous nucleotide substitution (T19234C) in the Fgf5 gene. At this locus, the C allele was present in all of the tested Angora rabbits, while the T allele was dominant in New Zealand and Rex rabbits. We further confirmed that the C allele was conserved in Angora rabbits by screening an additional 135 rabbits. Moreover, the results of functional predictions and co-immunoprecipitation revealed that the T19234C mutation impaired the binding capacity of FGF5 to its receptor FGFR1.

CONCLUSIONS

We discovered that the homozygous missense mutation T19234C within Fgf5 might contribute to the long-hair trait of Angora rabbits by reducing its receptor binding capacity. This finding will provide new insights into the genetic basis underlying the genetic improvement of Angora rabbits and benefit the improvement of rabbit breeding in the future.

Characterization of Growth Secondary Hair in Min Pig Activated by Follicle Stem Cell Stimulated by Wnt and BMP Signaling Pathway

In China, the national-level protected pig, the Min pig, is characterized by the development of secondary hairs and hair follicles in winter. Factors that dominate the genotype in the growth of secondary hairs are not clear through the concrete cell signaling pathways. This study compared hair phenotypes based on morphological structure, transcriptomics, and potential targeting molecules in the breeds of Min, Berkshire, and Yorkshire pigs. The results indicated that Min pigs have specific characteristics for the growth of secondary hairs compared with the Berkshire and Yorkshire pigs. The transcriptome analyses and quantitative reverse transcription-polymerase chain reaction results revealed that secondary hair growth was activated by follicle stem cells. The specific inhibitors of Wnt and BMP were studied using respective signals. The density of follicles, activity of follicle stem cells, and relative gene expression results have shown that Wnt and BMP stimulate the activity of follicle stem cells, and the Wnt signaling molecule has a significantly better effect than the BMP signaling molecule on stem cells. Wnt and BMP can promote the growth of local secondary hair and gene expression. Therefore, this study was conducted to verify the development mechanisms of secondary hairs, which have potential applications in laboratory animals and comparative medicine.

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.

Assessment of 2 Hair Removal Methods in New Zealand White Rabbits (Oryctolagus cuniculus)

Rabbits are frequently used as surgical models in research. However, studies assessing the effects of various hair removal methods on wound healing and surgical site infection (SSI) in rabbits are sparse. Here we evaluated the effects of 2 hair removal methods-clipping with electric clippers and using a commercial depilatory agent-on wound healing and SSI as assessed via wound scoring and histology. Incisions were assigned ASEPSIS scores on days 3 and 7. To assess whether the hair removal methods influenced aseptic preparation, swabs for bacterial culture were obtained just after hair removal on day 0, after aseptic skin preparation on day 0, and on day 1. For histopathologic assessment, full-thickness punch biopsies were obtained on days 0, 1, 3, 7, and 21. Histopathology revealed significant differences between the 2 methods, with the depilatory method having consistently higher scores (that is, more abnormalities). We conclude that for a surgical preparation regimen, clipping is safer, more efficacious, and less traumatic to tissues in rabbits.

Analysis of histology and long noncoding RNAs involved in the rabbit hair follicle density using RNA sequencing

BACKGROUND

Hair follicle density influences wool fibre production, which is one of the most important traits of the Wan Strain Angora rabbit. However, molecular mechanisms regulating hair follicle density have remained elusive.

RESULTS

In this study, hair follicle density at different body sites of Wan Strain Angora rabbits with high and low wool production (HWP and LWP) was investigated by histological analysis. Haematoxylin-eosin staining showed a higher hair follicle density in the skin of the HWP rabbits. The long noncoding RNA (lncRNA) profile was investigated by RNA sequencing, and 50 and 38 differentially expressed (DE) lncRNAs and genes, respectively, were screened between the HWP and LWP groups. A gene ontology analysis revealed that phospholipid, lipid metabolic, apoptotic, lipid biosynthetic, and lipid and fatty acid transport processes were significantly enriched. Potential functional lncRNAs that regulate lipid metabolism, amino acid synthesis, as well as the Janus kinase (JAK)-signal transducer and activator of transcription (STAT) and hedgehog signalling pathways, were identified. Consequently, five lncRNAs (LNC_002171, LNC_000797, LNC_005567, LNC_013595, and LNC_020367) were considered to be potential regulators of hair follicle density and development. Three DE lncRNAs and genes were validated by quantitative real-time polymerase chain reaction (q-PCR).

CONCLUSIONS

LncRNA profiles provide information on lncRNA expression to improve the understanding of molecular mechanisms involved in the regulation of hair follicle density.

Rabbits - their domestication and molecular genetics of hair coat development and quality

The European rabbit (Oryctolagus cuniculus) is the only representative of its genus living in present-day Europe and North Africa, and all domestic rabbits are descendants of this one species, which is native to the Iberian Peninsula. There are over 300 breeds of rabbits that differ in size, coat color, length of ears and type of fur. Rabbits are bred for various reasons, such as for laboratory animals and a source of meat, wool and fur, as well as for pets and exhibition animals. The hair coat is a important economic trait of rabbits. Its development and quality are influenced by various factors, both environmental and genetic. The genetic mechanisms underlying its development have not been thoroughly researched. The aim of this review is to discuss the domestication of rabbits and the different aspects of rabbit genetics. A brief review of the properties of rabbit hair coat, hair coat development and hair cycle will be provided, followed by discussion of the factors regulating hair coat development, molecular control of hair coat development and the role of non-coding RNAs in the regulation of gene expression in the hair follicles of rabbits. Information about genetic regulation of pathways could provide useful tools for improving hair coat quality and be of practical use in rabbit breeding.

Transcriptome-based selection and validation of optimal house-keeping genes for skin research in goats (Capra hircus)

BACKGROUND

In quantitative real-time polymerase chain reaction (qRT-PCR) experiments, accurate and reliable target gene expression results are dependent on optimal amplification of house-keeping genes (HKGs). RNA-seq technology offers a novel approach to detect new HKGs with improved stability. Goat (Capra hircus) is an economically important livestock species and plays an indispensable role in the world animal fiber and meat industry. Unfortunately, uniform and reliable HKGs for skin research have not been identified in goat. Therefore, this study seeks to identify a set of stable HKGs for the skin tissue of C. hircus using high-throughput sequencing technology.

RESULTS

Based on the transcriptome dataset of 39 goat skin tissue samples, 8 genes (SRP68, NCBP3, RRAGA, EIF4H, CTBP2, PTPRA, CNBP, and EEF2) with relatively stable expression levels were identified and selected as new candidate HKGs. Commonly used HKGs including SDHA and YWHAZ from a previous study, and 2 conventional genes (ACTB and GAPDH) were also examined. Four different experimental variables: (1) different development stages, (2) hair follicle cycle stages, (3) breeds, and (4) sampling sites were used for determination and validation. Four algorithms (geNorm, NormFinder, BestKeeper, and ΔCt method) and a comprehensive algorithm (ComprFinder, developed in-house) were used to assess the stability of each HKG. It was shown that NCBP3 + SDHA + PTPRA were more stably expressed than previously used genes in all conditions analysis, and that this combination was effective at normalizing target gene expression. Moreover, a new algorithm for comprehensive analysis, ComprFinder, was developed and released.

CONCLUSION

This study presents the first list of candidate HKGs for C. hircus skin tissues based on an RNA-seq dataset. We propose that the NCBP3 + SDHA + PTPRA combination could be regarded as a triplet set of HKGs in skin molecular biology experiments in C. hircus and other closely related species. In addition, we also encourage researchers who perform candidate HKG evaluations and who require comprehensive analysis to adopt our new algorithm, ComprFinder.

Skin transcriptome profiling of Changthangi goats highlights the relevance of genes involved in Pashmina production

Pashmina, the world's finest natural fiber is derived from secondary hair follicles of Changthangi goats which are domesticated in Ladakh region of Jammu and Kashmir by nomadic pastoralists. Complex epithelial-mesenchymal interactions involving numerous signal molecules and signaling pathways govern hair follicle morphogenesis and mitosis across different species. The present study involved transcriptome profiling of skin from fiber type Changthangi goats and meat type Barbari goats to unravel gene networks and metabolic pathways that might contribute to Pashmina development. In Changthangi goats, 525 genes were expressed at significantly higher levels and 54 at significantly lower levels with fold change >2 (p_adj < 0.05). Functional annotation and enrichment analysis identified significantly enriched pathways to be formation of the cornified envelope, keratinization and developmental biology. Expression of genes for keratins (KRTs) and keratin-associated proteins (KRTAPs) was observed to be much higher in Changthangi goats. A host of transcriptional regulator genes for hair follicle keratin synthesis such as GPRC5D, PADI3, HOXC13, FOXN1, LEF1 and ELF5 showed higher transcript abundance in Pashmina producing goats. Positive regulation of Wnt signaling pathway and negative regulation of Oncostatin M signaling pathway may be speculated to be important contributors to hair follicle development and hair shaft differentiation in Changthangi goats.

CRISPR/Cas9-mediated Disruption of Fibroblast Growth Factor 5 in Rabbits Results in a Systemic Long Hair Phenotype by Prolonging Anagen

Hair growth and morphology are generally regulated by the hair cycle in mammals. Fibroblast Growth Factor 5 (FGF5), which is a hair cycle regulator, has a role in regulating the hair cycle during the transition from the anagen phase to the catagen phase, and a hereditary long hair phenotype has been widely reported when FGF5 is mutated in humans and other species. However, there has been no such report in rabbits. Thus, the first exon of rabbit FGF5 was disrupted by the CRISPR/Cas9 system, and the phenotype of FGF5^-/- rabbits was characterized while using hematoxylin and eosin (H&E) staining, immunohistochemistry, quantitative PCR, scanning electron microscopy, and western blotting. The results showed a significant and systemic long hair phenotype in the FGF5^-/- rabbits, which indicated that FGF5 is a negative regulator of hair growth. In addition, a decreased diameter of the fiber and a higher area proportion of hair follicle clusters were determined in FGF5^-/- rabbits as compared with the WT rabbits. Further investigation verified that prolonging the anagen phase in rabbits, with decreased BMP2/4 pathway signaling and increased VERSICAN pathway signaling, caused the systemic long hair phenotype. Taken together, these results indicate a systemic long hair phenotype by prolonging anagen in FGF5^-/- rabbits, which could be widely used for Fur production and an ideal model for studying the mechanism of long hair in the future.

Analysis of histological and microRNA profiles changes in rabbit skin development

The periodic regrowth of rabbit fur is economically important. Here, we aimed to characterise the histological traits and microRNA (miRNA) expression profiles in the skin tissue of Wan Strain Angora rabbits at different weeks after plucking. Haematoxylin-eosin staining showed that hair follicles were in the telogen phase in the first week, while they were in the anagen phase from the fourth to twenty-fourth weeks. In addition, two small RNA libraries derived from skin samples of Wan Strain Angora rabbits at telogen and anagen stages yielded over 24 million high-quality reads. Specifically, 185 miRNAs were differentially expressed between the telogen and anagen phases. The function of the differentially expressed miRNAs was explored by comparing them with known mammalian miRNAs and by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis of their predicted targets. Five new functional miRNAs were validated using quantitative real-time PCR. Moreover, the fibroblast growth factor 5 (FGF5) gene was verified to be a target of conservative_NC_013672.1_9290 and conservative_NC_013675.1_10734. We investigated differential miRNA profiles between the telogen and anagen phases of the hair cycle and our findings provide a basis for future studies focusing on the mechanisms of miRNA-mediated regulation of rabbit hair follicle cycling.