Exploring right ovary degeneration in duck and goose embryos by histology and transcriptome dynamics analysis.

BACKGROUND: The development of asymmetric chick gonads involves separate developmental programs in the left and right gonads. In contrast to the left ovary developing into a fully functional reproductive organ, the right ovary undergoes gradual degeneration. However, the molecular mechanisms underlying the the degeneration of the right ovary remain incompletely understood. In the present study, we investigated the histomorphological and transcriptomic changes in the right ovary of ducks and geese during the the embryonic stage up to post-hatching day 1. RESULT: Hematoxylin-eosin stainings revealed that the right ovary developed until embryonic day 20 in ducks (DE20) or embryonic day 22 in geese (GE22), after which it started to regress. Further RNA-seq analyses revealed that both the differentially expressed genes (DEGs) in ducks and geese right ovary developmental stage were significantly enriched in cell adhesion-related pathway (ECM-receptor interaction, Focal adhesion pathway) and Cellular senescence pathway. Then during the degeneration stage, the DEGs were primarily enriched in pathways associated with inflammation, including Herpes simplex virus 1 infection, Influenza A, and Toll-like receptor signaling pathway. Moreover, duck-specific DEGs showed enrichment in Steroid hormone biosynthesis, Base excision repair, and the Wnt signaling pathway, while geese-specifically DEGs were found to be enriched in apoptosis and inflammation-related pathways, such as Ferroptosis, Necroptosis, RIG-I-like receptor signaling pathway, and NOD-like receptor signaling pathway. These findings suggest that the degeneration process of the right ovary in ducks occurs at a slower pace compared to that in geese. Additionally, the observation of the left ovary of the geese varying degeneration rates in the right ovary after hatching indicated that the development of the left ovary may be influenced by the degeneration of the right ovary. CONCLUSION: The data presented in this study provide valuable insights into the dynamic changes in histological structure and transcriptome during the degeneration of the right ovary in ducks and geese. In addition, through the analysis of shared characteristics in the degeneration process of the right ovary in both ducks and geese, we have uncovered the patterns of degradation and elucidated the molecular mechanisms involved in the regression of the right ovary in poultry. Furthermore, we have also made initial discoveries regarding the relationship between the degeneration of the right ovary and the development of the left ovary.

Identification of the genetic basis of the duck growth rate in multiple growth stages using genome-wide association analysis.

BACKGROUND: The genetic locus responsible for duck body size has been fully explained before, but the growth trait-related genetic basis is still waiting to be explored. For example, the genetic site related to growth rate, an important economic trait affecting marketing weight and feeding cost, is still unclear. Here, we performed genome wide association study (GWAS) to identify growth rate-associated genes and mutations. RESULT: In the current study, the body weight data of 358 ducks were recorded every 10 days from hatching to 120 days of age. According to the growth curve, we evaluated the relative and absolute growth rates (RGR and AGR) of 5 stages during the early rapid growth period. GWAS results for RGRs identified 31 significant SNPs on autosomes, and these SNPs were annotated by 24 protein-coding genes. Fourteen autosomal SNPs were significantly associated with AGRs. In addition, 4 shared significant SNPs were identified as having an association with both AGR and RGR, which were Chr2: 11483045 C>T, Chr2: 13750217 G>A, Chr2: 42508231 G>A and Chr2: 43644612 C>T. Among them, Chr2: 11483045 C>T, Chr2: 42508231 G>A, and Chr2: 43644612 C>T were annotated by ASAP1, LYN and CABYR, respectively. ASAP1 and LYN have already been proven to play roles in the growth and development of other species. In addition, we genotyped every duck using the most significant SNP (Chr2: 42508231 G>A) and compared the growth rate difference among each genotype population. The results showed that the growth rates of individuals carrying the Chr2: 42508231 A allele were significantly lower than those without this allele. Moreover, the results of the Mendelian randomization (MR) analysis supported the idea that the growth rate and birth weight had a causal effect on the adult body weight, with the growth rate having a greater effect size. CONCLUSION: In this study, 41 SNPs significantly related to growth rate were identified. In addition, we considered that the ASAP1 and LYN genes are essential candidate genes affecting the duck growth rate. The growth rate also showed the potential to be used as a reliable predictor of adult weight, providing a theoretical reference for preselection.

Comparison of the Effects of Recombinant and Native Prolactin on the Proliferation and Apoptosis of Goose Granulosa Cells.

In poultry, prolactin (PRL) plays a key role in the regulation of incubation behavior, hormone secretion, and reproductive activities. However, previous in vitro studies have focused on the actions of PRL in ovarian follicles of poultry, relying on the use of exogenous or recombinant PRL, and the true role of PRL in regulating ovarian granulosa cell (GC) functions in poultry awaits a further investigation using endogenous native PRL. Therefore, in this study, we first isolated and purified recombinant goose PRL protein (rPRL) and native goose PRL protein (nPRL) using Ni-affinity chromatography and rabbit anti-rPRL antibodies-filled immunoaffinity chromatography, respectively. Then, we analyzed and compared the effects of rPRL and nPRL at different concentrations (0, 3, 30, or 300 ng/mL) on the proliferation and apoptosis of both GCs isolated from goose ovarian pre-hierarchical follicles (phGCs) and from hierarchical follicles (hGCs). Our results show that rPRL at lower concentrations increased the viability and proliferation of both phGCs and hGCs, while it exerted anti-apoptotic effects in phGCs by upregulating the expression of Bcl-2. On the other hand, nPRL increased the apoptosis of phGCs in a concentration-dependent manner by upregulating the expressions of caspase-3 and Fas and downregulating the expressions of Bcl-2 and Becn-1. In conclusion, this study not only obtained a highly pure nPRL for the first time, but also suggested a dual role of PRL in regulating the proliferation and apoptosis of goose GCs, depending on its concentration and the stage of follicle development. The data presented here can be helpful in purifying native proteins of poultry and enabling a better understanding of the roles of PRL during the ovarian follicle development in poultry.

MiR-202-5p Regulates Geese Follicular Selection by Targeting BTBD10 to Regulate Granulosa Cell Proliferation and Apoptosis.

The regulation of granulosa cells (GCs) proliferation and apoptosis is the key step in follicular selection which determines the egg production performance of poultry. miR-202-5p has been reported to be involved in regulating the proliferation and apoptosis of mammalian ovarian GCs. However, its role in regulating the proliferation and apoptosis of goose GCs is still unknown. In the present study, the GCs of pre-hierarchical follicles (phGCs, 8-10 mm) and those of hierarchical follicles (hGCs, F2-F4) were used to investigate the role of miR-202-5p in cell proliferation and apoptosis during follicle selection. In phGCs and hGCs cultured in vitro, miR-202-5p was found to negatively regulate cell proliferation and positively regulate cell apoptosis. The results of RNA-seq showed that BTB Domain Containing 10 (BTBD10) is predicted to be a key target gene for miR-202-5p to regulate the proliferation and apoptosis of GCs. Furthermore, it is confirmed that miR-202-5p can inhibit BTBD10 expression by targeting its 3'UTR region, and BTBD10 was revealed to promote the proliferation and inhibit the apoptosis of phGCs and hGCs. Additionally, co-transfection with BTBD10 effectively prevented miR-202-5p mimic-induced cell apoptosis and the inhibition of cell proliferation. Meanwhile, miR-202-5p also remarkably inhibited the expression of Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta (PIK3CB) and AKT Serine/Threonine Kinase 1 (AKT1), while it was significantly restored by BTBD10. Overall, miR-202-5p suppresses the proliferation and promotes the apoptosis of GCs through the downregulation of PIK3CB/AKT1 signaling by targeting BTBD10 during follicular selection. Our study provides a theoretical reference for understanding the molecular mechanism of goose follicular selection, as well as a candidate gene for molecular marker-assisted breeding to improve the geese' egg production performance.

Comparative transcriptome analysis identifies crucial candidate genes and pathways in the hypothalamic-pituitary-gonadal axis during external genitalia development of male geese.

BACKGROUND: All birds reproduce via internal fertilization, but only ~3% of male birds possess the external genitalia that allows for intromission. Waterfowl (e.g., duck and goose) are representatives of them, and the external genitalia development of male geese is directly related to mating ability. Notably, some male geese show abnormal external genitalia development during ontogenesis. However, until now little is known about the molecular mechanisms of the external genitalia development in goose. In the present study, comparative transcriptomic analyses were performed on the hypothalamus, pituitary gland, testis, and external genitalia isolated from the 245-day-old male Tianfu meat geese showing normal (NEGG, n = 3) and abnormal (AEGG, n = 3) external genitals in order to provide a better understanding of the mechanisms controlling the development of the external genitalia in aquatic bird species. RESULTS: There were 107, 284, 2192, and 1005 differentially expressed genes (DEGs) identified in the hypothalamus, pituitary gland, testis and external genitalia between NEGG and AEGG. Functional enrichment analysis indicated that the DEGs identified in the hypothalamus were mainly enriched in the ECM-receptor interaction pathway. The ECM-receptor interaction, focal adhesion, and neuroactive ligand-receptor interaction pathways were significantly enriched by the DEGs in the pituitary gland. In the testis, the DEGs were enriched in the neuroactive ligand-receptor interaction, cell cycle, oocyte meiosis, and purine metabolism. In the external genitalia, the DEGs were enriched in the metabolic, neuroactive ligand-receptor interaction, and WNT signaling pathways. Furthermore, through integrated analysis of protein-protein interaction (PPI) network and co-expression network, fifteen genes involved in the neuroactive ligand-receptor interaction and WNT signaling pathways were identified, including KNG1, LPAR2, LPAR3, NPY, PLCB1, AVPR1B, GHSR, GRM3, HTR5A, FSHB, FSHR, WNT11, WNT5A, WIF1, and WNT7B, which could play crucial roles in the development of goose external genitalia. CONCLUSIONS: This study is the first systematically comparing the hypothalamus, pituitary gland, testis, and external genitalia transcriptomes of male geese exhibiting normal and abnormal external genitals. Both bioinformatic analysis and validation experiments indicated that the neuroactive ligand-receptor interaction pathway could regulate the WNT signaling pathway through PLCB1 to control male goose external genitalia development.

Effects of floor- and net-rearing systems on intestinal growth and microbial diversity in the ceca of ducks.

BACKGROUND: Rearing systems can affect livestock production directly, but whether they have effects on intestinal growth states and ceca microorganisms in ducks is largely unclear. The current study used Nonghua ducks to estimate the effects of rearing systems on the intestines by evaluating differences in intestinal growth indices and cecal microorganisms between ducks in the floor-rearing system (FRS) and net-rearing system (NRS). RESULTS: The values of relative weight (RW), relative length (RL) and RW/RL of the duodenum, jejunum, ileum and ceca in the FRS were significantly higher than those in the NRS during weeks 4, 8 and 13 (p < 0.05). A total of 157 genera were identified from ducks under the two systems, and the dominant microorganisms in both treatments were Firmicutes, Bacteroidetes, Actinobacteria and Proteobacteria at the phylum level. The distribution of microorganisms in the ceca of the two treatments showed significant separation during the three time periods, and the value of the Simpson index in the FRS was significantly higher than that in the NRS at 13 weeks (p < 0.05). Five differential microorganisms and 25 differential metabolic pathways were found in the ceca at week 4, seven differential microorganisms and 25 differential metabolic pathways were found in the ceca at week 8, and four differential microorganisms and two differential metabolic pathways were found in the ceca at week 13. CONCLUSIONS: The rearing system influences duck intestinal development and microorganisms. The FRS group had higher intestinal RL, RW and RW/RL and obviously separated ceca microorganisms compared to those of the NRS group. The differential metabolic pathways of cecal microorganisms decreased with increasing age, and the abundance of translation pathways was higher in the NRS group at week 13, while cofactor and vitamin metabolism were more abundant in the FRS group.

Association of leptin mRNA expression with meat quality trait in Tianfu black rabbits.

Leptin is a hormone synthesized and secreted primarily in adipocyte which can help to regulate energy balance. In this experiment, three tissue samples of Tianfu black rabbits at four growing periods were selected. The expression levels of leptin gene in different tissues were detected by real-time fluorescence quantitative PCR. The correlation analysis showed that the correlation coefficient between the expression levels of leptin gene in perirenal fat and intramuscular fat content in 84-day-old male rabbits was 0.73 (p < 0.05); the correlation coefficients between the expression levels of leptin gene in left biceps femoris and intramuscular fat and 24-hour pH in 84-day-old male rabbits were 0.95 (p < 0.01) and 0.85 (p < 0.05), respectively. Besides, the correlation coefficient between the expression levels of leptin gene in male left biceps femoris and cooked meat rate was 0.83 (p < 0.05). According to the analysis results, we inferred the expression levels of leptin gene in Tianfu black rabbits can influence meat quality and the meat quality of high expression levels of leptin gene in Tianfu black rabbits is better. These results revealed the leptin gene may be one of the important candidate genes for meat quality traits of molecular markers.

Integrated Transcriptome and Metabolome Analysis Reveals the Regulatory Mechanisms of FASN in Geese Granulosa Cells.

FASN plays a critical role in lipid metabolism, which is involved in regulating ovarian follicular development. However, the molecular mechanisms of how FASN regulate the function of ovarian follicular cells still remain elusive. In this study, by overexpression or interference of FASN in pre-hierarchical follicle granulosa cells (phGCs) and hierarchical follicle granulosa cells (hGCs), we analyzed their effects on the granulosa cell transcriptome and metabolome profiles using RNA-Seq and LC-MS/MS, respectively. The results showed that overexpression of FASN promoted proinflammatory factors expression by activating TLR3/IRF7 and TLR3/NF-κB pathways in phGCs, but only by activating TLR3/IRF7 pathways in hGCs. Then, necroptosis and apoptosis were triggered through the JAK/STAT1 pathway (induced by inflammatory factors) and BAK/caspase-7 pathway, respectively. The combined analysis of the metabolome and transcriptome revealed that FASN affected the demand of GCs for 5-hydroxytryptamine (5-HT) by activating the neuroactive ligand-receptor interaction pathway in two categorized GCs and only altering the metabolic pathway of tryptophan in phGCs, and ultimately participated in regulating the physiological function of geese GCs. Taken together, this study showed that the mechanisms of FASN regulating the physiological function of geese phGCs and hGCs were similar, but they also had some different characteristics.

Oestrogen promotes lipids transportation through oestrogen receptor α in hepatic steatosis of geese in vitro.

Evidence has shown that oestrogen suppresses lipids deposition in the liver of mammals. However, the molecular mechanism of oestrogen action in hepatic steatosis of geese liver has yet to be determined. This study aimed to investigate the effect of oestrogen on lipid homeostasis at different states of geese hepatocytes in vitro. The results showed that an in vitro model of hepatic steatosis was induced by 1.5 mM sodium oleate via detecting the viability of hepatocytes and content of lipids. When the normal hepatocytes were administrated with different concentrations of oestrogen (E2 ), the expression levels of diacylglycerol acyltransferase 2 (DGAT2), microsomal triglyceride transfer protein (MTTP) and oestrogen receptors (ERs, alpha and beta) were up-regulated only at high concentrations of E2 , whereas the lipid content was not a significant difference. In goose hepatocytes of hepatic steatosis, however, the expression levels of MTTP, apolipoprotein B (apoB) and ERα/ß significantly increased at 10-7 or 10-6  M E2 . Meanwhile, the lipids content significantly increased at 10-9 and 10-8  M E2 and decreased at 80 µM E2 . Further heatmap analysis showed that ERα was clustered with apoB and MTTP in either normal hepatocytes or that of hepatic steatosis. Taken together, E2  might bind to ERα to up-regulate the expression levels of apoB and MTTP, promoting the transportation of lipids and alleviating lipids overload in hepatic steatosis of geese in vitro.

Integrative analysis of histomorphology, transcriptome and whole genome resequencing identified DIO2 gene as a crucial gene for the protuberant knob located on forehead in geese.

BACKGROUND: During domestication, remarkable changes in behavior, morphology, physiology and production performance have taken place in farm animals. As one of the most economically important poultry, goose owns a unique appearance characteristic called knob, which is located at the base of the upper bill. However, neither the histomorphology nor the genetic mechanism of the knob phenotype has been revealed in geese. RESULTS: In the present study, integrated radiographic, histological, transcriptomic and genomic analyses revealed the histomorphological characteristics and genetic mechanism of goose knob. The knob skin was developed, and radiographic results demonstrated that the knob bone was obviously protuberant and pneumatized. Histologically, there were major differences in structures in both the knob skin and bone between geese owing knob (namely knob-geese) and those devoid of knob (namely non-knob geese). Through transcriptome analysis, 592 and 952 genes differentially expressed in knob skin and bone, and significantly enriched in PPAR and Calcium pathways in knob skin and bone, respectively, which revealed the molecular mechanisms of histomorphological differences of the knob between knob- and non-knob geese. Furthermore, integrated transcriptomic and genomic analysis contributed to the identification of 17 and 21 candidate genes associated with the knob formation in the skin and bone, respectively. Of them, DIO2 gene could play a pivotal role in determining the knob phenotype in geese. Because a non-synonymous mutation (c.642,923 G > A, P265L) changed DIO2 protein secondary structure in knob geese, and Sanger sequencing further showed that the AA genotype was identified in the population of knob geese, and was prevalent in a crossing population which was artificially selected for 10 generations. CONCLUSIONS: This study was the first to uncover the knob histomorphological characteristics and genetic mechanism in geese, and DIO2 was identified as the crucial gene associated with the knob phenotype. These data not only expand and enrich our knowledge on the molecular mechanisms underlying the formation of head appendages in both mammalian and avian species, but also have important theoretical and practical significance for goose breeding.

Molecular Profiling of DNA Methylation and Alternative Splicing of Genes in Skeletal Muscle of Obese Rabbits.

DNA methylation and the alternative splicing of precursor messenger RNAs (pre-mRNAs) are two important genetic modification mechanisms. However, both are currently uncharacterized in the muscle metabolism of rabbits. Thus, we constructed the Tianfu black rabbit obesity model (obese rabbits fed with a 10% high-fat diet and control rabbits from 35 days to 70 days) and collected the skeletal muscle samples from the two groups for Genome methylation sequencing and RNA sequencing. DNA methylation data showed that the promoter regions of 599 genes and gene body region of 2522 genes had significantly differential methylation rates between the two groups, of which 288 genes had differential methylation rates in promoter and gene body regions. Analysis of alternative splicing showed 555 genes involved in exon skipping (ES) patterns, and 15 genes existed in differential methylation regions. Network analysis showed that 20 hub genes were associated with ubiquitinated protein degradation, muscle development pathways, and skeletal muscle energy metabolism. Our findings suggest that the two types of genetic modification have potential regulatory effects on skeletal muscle development and provide a basis for further mechanistic studies in the rabbit.

miR-194-5p negatively regulates the proliferation and differentiation of rabbit skeletal muscle satellite cells.

Skeletal muscle satellite cells (SMSCs), also known as a multipotential stem cell population, play a crucial role during muscle growth and regeneration. In recent years, numerous miRNAs have been associated with the proliferation and differentiation of SMSCs in a number of mammalian species; however, the regulatory mechanisms of miR-194-5p in rabbit SMSCs still remain scarce. In this study, miR-194-5p was first observed to be highly expressed in the rabbit leg muscle. Furthermore, both the mimics and inhibitor of miR-194-5p were used to explore its role in the proliferation and differentiation of rabbit SMSCs cultured in vitro. Results from both EdU and CCK8 assays showed that miR-194-5p inhibited the proliferation of SMSCs. Meanwhile, Mef2c was identified as a target gene of miR-194-5p based on the dual-luciferase reporter assay results. In addition, upregulation of miR-194-5p decreased the expression levels of Mef2c and MyoG during rabbit SMSCs differentiation on Days 3 and 7 of in vitro culture. Taken together, these data demonstrated that miR-194-5p negatively regulates the proliferation and differentiation of rabbit SMSCs by targeting Mef2c.

Co-culture model reveals the characteristics of theca cells and the effect of granulosa cells on theca cells at different stages of follicular development.

Theca cells (TCs) play an important role in follicular development, which cannot be separated from granulosa cells (GCs). However, compared with mammals, the TCs and the effects of GCs on TCs at different follicular development stages (FDSs) have specific characteristics in avian species, but none of them have been clearly defined. In this study, we established an in vitro co-culture (with GC at the corresponding stage) model of goose TCs at different FDSs (pre-hierarchical, hierarchical and F1) by using a transwell system. The properties of TCs in co-culture at the three FDSs, including cell morphology, activity and intracellular lipid content, as well as the expression of key genes involved in de novo lipogenesis, steroidogenesis, proliferation and apoptosis, were examined and defined. We further compared the mono-culture and co-culture groups. After co-culture, the activity of TCs showed significant (p < .01) increases in all stages; moreover, in pre-hierarchical TCs, the expression levels of FAS, SREBP, 3ß-HSD and CCND1 were promoted, and PPARγ, CYP19, BCL2 and CAS3 were inhibited (p < .05); in the hierarchical TCs, the expression levels of PPARγ, FAS, CYP19, CCND1 and BCL2 were promoted, and SREBP, STAR, 3ß-HSD and CAS3 were inhibited (p < .05), whereas in the F1 TCs, the expression levels of PPARγ, FAS, 3ß-HSD, CYP19 and CCND1 were promoted, and STAR and CAS3 were inhibited (p < .05). These results suggested that GCs at the three FDSs have dynamic and complex influences on the physiological characteristics of TCs, and the influences on TCs at the three FDSs were varied.

Genome-wide identification and characterization of long non-coding RNAs during differentiation of visceral preadipocytes in rabbit.

Emerging evidence suggests that long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes, including adipogenesis. Despite being considered an ideal animal model for studying adipogenesis, little is known about the roles of lncRNAs in the regulation of rabbit preadipocyte differentiation. In the present study, visceral preadipocytes isolated from newborn rabbits were cultured in vitro and induced for differentiation, and global lncRNA expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. A total of 2066 lncRNAs were identified from nine RNA-seq libraries. Compared to protein-coding transcripts, lncRNA transcripts exhibited characteristics of a longer length and lower expression level. Furthermore, 486 and 357 differentially expressed (DE) lncRNAs were identified when comparing day 3 vs. day 0 and day 9 vs. day 3, respectively. Target genes of DE lncRNAs were predicted by the cis-regulating approach. Prediction of functions revealed that DE lncRNAs when comparing day 3 vs. day 0 were involved in gene ontology (GO) terms of developmental growth, growth, developmental cell growth, and stem cell proliferation, and involved in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of PI3K-Akt signaling pathway, fatty acid biosynthesis, and the insulin signaling pathway. The DE lncRNAs when comparing day 9 vs. day 3 were involved in GO terms that associated with epigenetic modification and were involved in the KEGG pathway of cAMP signaling pathway. This study provides further insight into the regulatory function of lncRNAs in rabbit visceral adipose and facilitates a better understanding of different stages of preadipocyte differentiation.

Paternal weight of ducks may have an influence on offspring' small intestinal function and cecal microorganisms.

BACKGROUND: In animals, many factors affect the small intestinal function and cecal microorganisms, including body weight and genetic background. However, whether paternal weight impacts the small intestinal function and cecal microorganisms remains unknown to date. The current study used Nonghua sheldrake to estimate the effect of paternal weight on the intestine of the offspring by evaluating differences in small intestinal morphology, digestive enzyme activity, and cecal microorganisms between the offspring of male parents with high body weight (group H) and that of male parents with low body weight (group L). RESULTS: The results of the analysis of small intestinal morphology showed that the villus height of the jejunum of group H ducks was higher than that of group L ducks, and the difference was significant for ducks at 10 weeks of age. Moreover, the villus height/crypt depth of the duodenum in group H significantly exceeded that of group L at a duck age of 2 weeks. The amylase activity in the jejunum content of group H exceeded that of group L at 5 and 10 weeks of age. Furthermore, the proportion of the Firmicutes to Bacteroidetes was significantly higher in group H (duck age of 2 weeks). Among the genera with a relative abundance exceeding 1%, the relative abundances of genera Desulfovibrio, Megamonas, Alistipes, Faecalibacterium, and Streptococcus observed in group H were significantly different between group H and group L. CONCLUSIONS: For the first time, this study identifies the effect of paternal weight on offspring small intestinal function and cecal microorganisms. Consequently, this lays a foundation for further research on the relationship between male parents and offspring intestinal function.

MiR-9-5p promotes rabbit preadipocyte differentiation by suppressing leptin gene expression.

BACKGROUND: MicroRNAs (miRNAs) are a class of small non-coding RNAs, which participate in the regulation of cell differentiation. Previous studies have demonstrated that miR-9-5p plays a key role in cancer cell development, but the mechanisms by which miR-9-5p regulates adipogenesis remain poorly understood. The present study intended to investigate its significance in producing rabbits with high-quality meat by observing the regulatory effect of miR-9-5p in preadipocytes and finding the related targets. METHODS: In this study, a dual-luciferase reporter assay was employed to validate the targeting relationship between miR-9-5p and leptin gene. We also utilized quantitative reverse transcription PCR (qRT-PCR), western blot, oil red-O staining assay, and determination of triglyceride content to analyze the regulation of miR-9-5p and leptin gene during adipocyte differentiation. RESULTS: The analysis demonstrated that during preadipocyte differentiation, miR-9-5p was up-regulated and the fat formation related biomarkers, i.e., fatty acid-binding protein 4 (FABP4), CCAAT-enhancer binding protein α (C/EBPα), and peroxisome proliferator activated receptor γ (PPARγ) were also up-regulated. Meanwhile, the oil red-O staining assay revealed that the accumulation of lipid droplets increased. We also explored the expression pattern and role of miR-9-5p in adipogenesis using white pre-adipocytes. The results showed that miR-9-5p was up-regulated during preadipocyte differentiation, and overexpression of miR-9-5p enhanced lipid accumulation. Furthermore, we found that the overexpression of miR-9-5p significantly up- regulated the expression of marker genes, PPARγ, C/EBPα and FABP4, and increased the protein levels of PPARγ and triglyceride content. The results suggest that miR-9-5p might be involved in the regulation of rabbit preadipocyte differentiation. We predicted that leptin is the target gene of miR-9-5p, by using bioinformatics tools and the conclusion was validated by a luciferase reporter assay. Finally, we verified that the knock-down of leptin by si-leptin promoted preadipocyte differentiation in rabbits. CONCLUSION: The results of the present study indicate that miR-9-5p regulates white preadipocyte differentiation in rabbits by targeting the leptin gene.

Screening and identification of MicroRNAs expressed in perirenal adipose tissue during rabbit growth.

BACKGROUND: MicroRNAs (miRNAs) regulate adipose tissue development, which are closely related to subcutaneous and intramuscular fat deposition and adipocyte differentiation. As an important economic and agricultural animal, rabbits have low adipose tissue deposition and are an ideal model to study adipose regulation. However, the miRNAs related to fat deposition during the growth and development of rabbits are poorly defined. METHODS: In this study, miRNA-sequencing and bioinformatics analyses were used to profile the miRNAs in rabbit perirenal adipose tissue at 35, 85 and 120 days post-birth. Differentially expressed (DE) miRNAs between different stages were identified by DEseq in R. Target genes of DE miRNAs were predicted by TargetScan and miRanda. To explore the functions of identified miRNAs, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. RESULTS: Approximately 1.6 GB of data was obtained by miRNA-seq. A total of 987 miRNAs (780 known and 207 newly predicted) and 174 DE miRNAs were identified. The miRNAs ranged from 18 nt to 26 nt. GO enrichment and KEGG pathway analyses revealed that the target genes of the DE miRNAs were mainly involved in zinc ion binding, regulation of cell growth, MAPK signaling pathway, and other adipose hypertrophy-related pathways. Six DE miRNAs were randomly selected, and their expression profiles were validated by q-PCR. CONCLUSIONS: This is the first report of the miRNA profiles of adipose tissue during different growth stages of rabbits. Our data provide a theoretical reference for subsequent studies on rabbit genetics, breeding and the regulatory mechanisms of adipose development.

Effect of thermal manipulation during embryogenesis on the promoter methylation and expression of myogenesis-related genes in duck skeletal muscle.

Avian embryos are an ideal system to investigate the effect of incubation temperature on embryonic development, but the characteristics and mechanisms of temperature effects on poultry embryonic myogenesis are unclear. In this study, we investigated the effect of increasing the incubation temperature by 1 °C on the expression of nine myogenesis-related genes in ducks and then explored the correlation between the alteration of promoter methylation and the expression of two of the nine genes under thermal manipulation (TM). The qRT-PCR results showed that TM during embryonic days (ED) 1-10 promoted (P < 0.05) the expression of genes in breast muscle (PAX3, PAX7, MYOG, MCK, SIX1, TNNC1) and leg muscle (MYOD, MYOG, MYF5, MCK, AKIRIN2, TNNC1). TM during ED10-20 promoted the expression of PAX3, MYF5 and MCK and inhibited AKIRIN2 expression in breast muscle (P < 0.05); however, it inhibited the expression of PAX3, PAX7, MYOD, MYOG, MYF5, SIX1, AKIRIN2 and TNNC1 and promoted MCK expression in leg muscle (P < 0.05). TM during ED20-27 inhibited the expression of genes in breast muscle (PAX7) and leg muscle (MYOD, MYOG, MYF5, TNNC1) and promoted MCK expression in breast and leg muscle (P < 0.05). Furthermore, with the Sequenom MassARRAY platform, it was observed that the average methylation level of AKIRIN2 (ED10) and TNNC1 (ED20) in leg muscle decreased (P < 0.05) after TM. Notably, we found significant (P < 0.05) inverse correlations between the methylation and mRNA levels of AKIRIN2 under TM during ED1-10 (r = - 0.969) and ED10-20 (r = - 0.805). Taken together, TM during ED1-10 was more favorable for improving duck myogenesis-related gene expression than TM during ED10-20 and ED20-27. TM during duck embryogenesis seemed to have a greater effect on the development of leg muscle than breast muscle and might alter AKIRIN2 expression by changing its promoter methylation status. These findings may be helpful to understand temperature effects on the muscle development of avian embryos and to explore the role of epigenetic regulation during this process.

Molecular cloning, polymorphism, and expression analysis of the LKB1/STK11 gene and its association with non-specific digestive disorder in rabbits.

Liver kinase B1 (LKB1, also called STK11) encodes a serine/threonine kinase mutated in Peutz-Jeghers cancer syndrome characterized by gastrointestinal polyposis. Although LKB1 plays an important role in regulating energy homeostasis, cell growth, and metabolism via activation of adenosine monophosphate (AMP)-activated protein kinase (AMPK), nothing is known about its molecular characteristics and possible involvement in non-specific digestive disorder (NSDD) of rabbits. In the present study, we first cloned the coding sequence (CDS) of rabbit LKB1, which consisted of 1317 bp encoding 438 amino acids (AAs) and contained a highly conserved S_TKc kinase domain. Its deduced AA sequence showed 87.93-91.10% similarities with that of other species. In order to determine its involvement in NSDD, a NSDD rabbit model was built by a dietary fiber deficiency. The polymorphic site of LKB1 was then investigated in both healthy and NSDD groups using directing sequencing. Our results suggested that a synonymous variant site (840 c. G > C, CCC→CCG) existed in its S_TKc domain, which was associated with susceptibility to NSDD. Furthermore, qPCR was utilized to examine the mRNA levels of LKB1 and its downstream targets (i.e., PRKAA2, mTOR and NF-kß) in several intestinal-related tissues from both healthy and NSDD groups. Significant changes in their expression levels between two groups indicated that impaired LKB1 signaling contributed to the intestinal abnormality in NSDD rabbits. Taken together, it could be concluded that LKB1 might be a potential candidate gene affecting the occurrence of rabbit NSDD. This information may serve as a basis for further investigations on rabbit digestive diseases.

Genome-wide identification and characterization of long non-coding RNAs during postnatal development of rabbit adipose tissue.

BACKGROUND: The rabbit is widely used as an important experimental model for biomedical research, and shows low adipose tissue deposition during growth. Long non-coding RNAs (lncRNAs) are associated with adipose growth, but little is known about the function of lncRNAs in the rabbit adipose tissue. METHODS: Deep RNA-sequencing and comprehensive bioinformatics analyses were used to characterize the lncRNAs of rabbit visceral adipose tissue (VAT) at 35, 85 and 120 days after birth. Differentially expressed (DE) lncRNAs were identified at the three growth stages by DESeq. The cis and trans prediction ways predicted the target genes of the DE lncRNAs. To explore the function of lncRNAs, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed on the candidate genes. RESULTS: A total of 991,157,544 clean reads were generated after RNA-Seq of the three growth stages, of which, 30,353 and 107 differentially expressed (DE) lncRNAs were identified. Compared to the protein-coding transcripts, the rabbit lncRNAs shared some characteristics such as shorter length and fewer exons. Cis and trans target gene prediction revealed, 43 and 64 DE lncRNAs respectively, corresponding to 72 and 20 protein-coding genes. GO enrichment and KEGG pathway analyses revealed that the candidate DE lncRNA target genes were involved in oxidative phosphorylation, glyoxylate and dicarboxylate metabolism, and other adipose growth-related pathways. Six DE lncRNAs were randomly selected and validated by q-PCR. CONCLUSIONS: This study is the first to profile the potentially functional lncRNAs in the adipose tissue growth in rabbits, and contributes to our understanding of mammalian adipogenesis.