Integration Analysis of Hair Follicle Transcriptome and Proteome Reveals the Mechanisms Regulating Wool Fiber Diameter in Angora Rabbits.

Fiber diameter is an important characteristic that determines the quality and economic value of rabbit wool. This study aimed to investigate the genetic determinants of wool fiber diameter through an integration analysis using transcriptomic and proteomic datasets from hair follicles of coarse and fine wool from Angora rabbits. Using a 4D label-free technique, we identified 423 differentially expressed proteins (DEPs) in hair follicles of coarse and fine wool in Angora rabbits. Eighteen DEPs were examined using parallel reaction monitoring, which verified the reliability of our proteomic data. Functional enrichment analysis revealed that a set of biological processes and signaling pathways related to wool growth and hair diameter were strongly enriched by DEPs with fold changes greater than two, such as keratinocyte differentiation, skin development, epidermal and epithelial cell differentiation, epidermis and epithelium development, keratinization, and estrogen signaling pathway. Association analysis and protein-protein interaction network analysis further showed that the keratin (KRT) family members, including KRT77, KRT82, KRT72, KRT32, and KRT10, as well as CASP14 and CDSN, might be key factors contributing to differences in fiber diameter. Our results identified DEPs in hair follicles of coarse and fine wool and promoted understanding of the molecular mechanisms underlying wool fiber diameter variation among Angora rabbits.

Transcriptome analysis of divergent residual feed intake phenotypes in the M. longissimus thoracis et lumborum of Wannan Yellow rabbits.

Introduction: Feed efficiency is an important economic trait in rabbit meat production. The identification of molecular mechanisms and candidate genes for feed efficiency may improve the economic and environmental benefits of the rabbit meat industry. As an alternative to the conventional feed conversion ratio, residual feed intake (RFI) can be used as an accurate indicator of feed efficiency. Methods: RNA sequencing was used to identify the differentially expressed genes (DEGs) in the M. longissimus thoracis et lumborum of eight Wannan Yellow rabbits with excessively high or low RFIs (HRFI or LRFI, respectively). Thereafter, Gene Ontology (GO) analysis, enrichment using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, gene set enrichment analysis (GSEA), and protein-protein interaction (PPI) network analysis was conducted. Results: In total, 445 DEGs were identified in the M. longissimus thoracis et lumborum of rabbits with high and low RFIs. The significantly enriched GO terms identified in these two groups were primarily involved in energy and mitochondrial metabolism and oxidation-reduction processes. KEGG analysis identified 11 significantly enriched pathways, including oxidative phosphorylation, PI3K-Akt signaling, and extracellular matrix-receptor interaction pathways. According to GSEA, the expressions of genes and pathways related to mitochondrial function were upregulated in HRFI rabbits, whereas genes with upregulated expressions in LRFI rabbits were related to immune response and energy metabolism. Additionally, PPI network analysis revealed five potential candidate genetic markers. Conclusion: Comparative analysis of the M. longissimus thoracis et lumborum transcriptomes in HRFI and LRFI rabbits revealed FOS, MYC, PRKACB, ITGA2, and FN1 as potential candidate genes that affect feed efficiency in rabbits. In addition, key signaling pathways involved in oxidative phosphorylation and PI3K-Akt and ECM-receptor interaction signaling impact rabbit feed efficiency. These findings will aid in breeding programs to improve feed efficiency and optimize RFI selection of rabbits for meat production.

Semen quality and seminal plasma metabolites in male rabbits (Oryctolagus cuniculus) under heat stress.

Heat stress causes infertility in male rabbits in summer. This study was conducted to determine the effects of heat stress on semen quality and seminal plasma metabolites of male rabbits. To achieve these objectives, the temperature and humidity index (THI) was used to determine the stress state of male rabbits during different months, thereby the rabbits were divided into heat stress and no heat stress groups. The quality of the semen and the biochemical indices of seminal plasma were then analyzed. Next the plasma metabolites of rabbits in both groups were evaluated using the ultra-high performance liquid chromatography-mass spectroscopy (UPLC-MS)/MS technique. Our results showed that the THI value of the rabbit housing in May was 20.94 (no heat stress). The THI value of the housing in August was 29.10 (heat stress group, n = 10). Compared with the non-heat stress group, the sperm motility, density, and pH in the heat stress group (n = 10) were significantly decreased (P < 0.01); the semen volume decreased significantly (P < 0.05); and the sperm malformation rate increased significantly (P < 0.01). The number of grade A sperm significantly decreased, while the numbers of B and C grade sperm significantly increased (P < 0.01). The total sperm output (TSO), total motile sperm (TMS), and total functional sperm fraction (TFSF) decreased significantly (P < 0.01). Heat stress protein 70 (HSP70) and acid phosphatase (ACP) in the seminal plasma of rabbits in the heat stress group (n = 20) were significantly increased (P < 0.01). Seminal plasma testosterone (T), α-glucosidase (α-Glu), and fructose decreased significantly (P < 0.01). The concentrations of Mg2+ (P < 0.05), Na+ (P < 0.01), and K+ (P < 0.01) in metal ions were significantly decreased. These findings indicated that heat stress severely affected the quality of the male rabbit semen. Furthermore, UPLC-MS/MS technology was used to analyze the seminal plasma samples of rabbits in the heat stress group and non-heat stress group (n = 9 for each group). In total, 346 metabolites were identified, with variable importance in project (VIP) > 1.0, fold change (FC) > 1.5 or < 0.667, and P < 0.05 as the threshold. A total of 71 differential metabolites were matched, including stearic acid, betaine, arachidonic acid, L-malic acid, and indole. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differential metabolites revealed 51 metabolic pathways, including synthesis and degradation of ketones, serine and threonine metabolism, tryptophan metabolism, and the citric acid cycle. Our study has shown that the sperm motility, sperm pH value, and sperm density of male rabbits decreased significantly under heat stress, and the sperm malformation rate increased significantly. Furthermore, the quality of semen was shown to deteriorate and the energy metabolism pathway was disturbed. These findings provide a theoretical reference for alleviating the adaptive heat stress in male rabbits.

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.

MeDIP-seq and RNA-seq analysis during porcine testis development reveals functional DMR at the promoter of LDHC.

Testis development requires tight regulation of gene expression programmed by epigenetic modifiers. However, their mechanism remains to be elucidated. Here, we investigated the genome-wide DNA methylation landscape in the Duroc and Meishan boar testes using methylated DNA immunoprecipitation sequencing (MeDIP-seq). We identified over 1100 promoter differential methylation genes (DMGs) before and after puberty, most of which are associated with testis development. Furthermore, we discovered that the expression of lactate dehydrogenase C (LDHC) gene during testis development is regulated by DNA methylation. The promoter of LDHC in pre-pubertal testes is substantially methylated, whereas considerably demethylated in post-pubertal testes. Artificial demethylation with the demethylating agent 5-Aza-CdR induced LDHC expression in immature Sertoli cells (SCs). Mechanistically, we confirmed the transcription factor SP1 was recruited to bind in hypomethylated differentially methylated regions (DMRs) in LDHC promoter, which upregulated the expression of LDHC. Functionally, we demonstrated that LDHC was activated in mature SCs (mSCs) and its overexpression significantly increases lactate secretion in SCs. In conclusion, our results highlight the function and regulation of dynamic DNA methylation in testis development.

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.

Integrated analysis of miRNA and mRNA expression profiles in testes of Duroc and Meishan boars.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNAs playing vital roles in regulating posttranscriptional gene expression. Elucidating the expression regulation of miRNAs underlying pig testis development will contribute to a better understanding of boar fertility and spermatogenesis. RESULTS: In this study, miRNA expression profile was investigated in testes of Duroc and Meishan boars at 20, 75, and 270 days of age by high-throughput sequencing. Forty-five differentially expressed miRNAs were identified from testes of Duroc and Meishan boars before and after puberty. Integrated analysis of miRNA and mRNA profiles predicted many miRNA-mRNA pairs. Gene ontology and biological pathway analyses revealed that predicted target genes of ssc-mir-423-5p, ssc-mir-34c, ssc-mir-107, ssc-mir-196b-5p, ssc-mir-92a, ssc-mir-320, ssc-mir-10a-5p, and ssc-mir-181b were involved in sexual reproduction, male gamete generation, and spermatogenesis, and GnRH, Wnt, and MAPK signaling pathway. Four significantly differentially expressed miRNAs and their predicted target genes were validated by quantitative real-time polymerase chain reaction, and phospholipase C beta 1 (PLCß1) gene was verified to be a target of ssc-mir-423-5p. CONCLUSIONS: This study provides an insight into the functional roles of miRNAs in testis development and spermatogenesis and offers useful resources for understanding differences in sexual function development caused by the change in miRNAs expression between Duroc and Meishan boars.

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.

Expression and regulation of GnRHR2 gene and testosterone secretion mediated by GnRH2 and GnRHR2 within porcine testes.

Gonadotropin-releasing hormone 2 receptor (GnRHR2) together with its cognate ligand involves in regulating reproductive behavior. However, little is known concerning the effect of transcription factor steroidogenic factor1 (SF-1) regulation on porcine GnRHR2 gene expression and GnRH2 regulation mechanism in testosterone secretion through GnRHR2. Our study demonstrated that GnRHR2 transcription levels were high in porcine testis. Immunohistochemistry analyses showed that GnRHR2 immunoreactivity was strong in the Leydig cells in boar testes. Two SF-1 binding sites were predicted in GnRHR2 promoter and the second site (-159/-149) was considered to be important for GnRHR2 promoter activity through site-directed mutagenesis. The binding of SF-1 to GnRHR2 promoter was confirmed by electrophoretic mobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP). Overexpression and knockdown experiments revealed that SF-1 could up-regulate porcine GnRHR2 expression. DNA methylation of GnRHR2 promoter CpG island also specifically regulated GnRHR2 expression. Meanwhile, our study also demonstrated GnRH2 treatment promoted the expression of SF-1 and steroidogenic acute regulatory protein (StAR), and that this treatment stimulated cAMP responsive element binding protein (CREB) phosphorylation, regulated the expression of GnRHR2, especially that of GnRHR2-X1, and promoted testosterone secretion in porcine Leydig cells. We speculated that testosterone secretion mediated by GnRH2 and GnRHR2 (mainly GnRHR2-X1) was regulated by phosphorylated CREB interacting with SF-1 to control StAR expression. Taken together, the present study indicates that SF-1 and GnRH2 are the essential regulatory factors for GnRHR2 expression. This study also explores the regulation mechanism of testosterone secretion mediated by GnRH2 and GnRHR2 in porcine Leydig cells.

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

BACKGROUND: Hair fibre length is an important economic trait of rabbits in fur production. However, molecular mechanisms regulating rabbit hair growth have remained elusive. RESULTS: Here we aimed to characterise the skin traits and gene expression profiles of short-hair and long-hair rabbits by histological and transcriptome analyses. Haematoxylin-eosin staining was performed to observe the histological structure of the skin of short-hair and long-hair rabbits. Compared to that in short-hair rabbits, a significantly longer anagen phase was observed in long-hair rabbits. In addition, by RNA sequencing, we identified 951 genes that were expressed at significantly different levels in the skin of short-hair and long-hair rabbits. Nine significantly differentially expressed genes were validated by quantitative real-time polymerase chain reaction. A gene ontology analysis revealed that epidermis development, hair follicle development, and lipid metabolic process were significantly enriched. Further, we identified potential functional genes regulating follicle development, lipid metabolic, and apoptosis as well as important pathways including extracellular matrix-receptor interaction and basal cell carcinoma pathway. CONCLUSIONS: The present study provides transcriptome evidence for the differences in hair growth between short-hair and long-hair rabbits and reveals that lipid metabolism and apoptosis might constitute major factors contributing to hair length.

Histological and transcriptome analyses of testes from Duroc and Meishan boars.

Meishan boars are known for their early sexual maturity. However, they exhibit a significantly smaller testicular size and a reduced proportion of Sertoli cells and daily sperm production compared with Duroc boars. The testes of Duroc and Meishan boars at 20, 75 and 270 days of age were used for histological and transcriptome analyses. Haematoxylin-eosin staining was conducted to observe histological structure of the testes in Duroc and Meishan boars at different ages. Although spermatogenesis occurred prior to 75 days in Meishan boars, the number of spermatogonia and Sertoli cells in Meishan boars were less than in Duroc boars at adulthood. The diameters of the seminiferous tubules of the testes differed significantly during the initiation of development of the seminiferous tubules between the two breeds. We obtained differentially expressed functional genes and analysed seven pathways involved in male sexual maturity and spermatogenesis using RNA-seq. We also detected four main alternative splicing events and many single nucleotide polymorphisms from testes. Eight functionally important genes were validated by qPCR, and Neurotrophin 3 was subjected to quantification and cellular localization analysis. Our study provides the first transcriptome evidence for the differences in sexual function development between Meishan and Duroc boars.