Expression profiling of preadipocyte microRNAs by deep sequencing on chicken lines divergently selected for abdominal fatness

Resequencing of a Pekin duck breeding population provides insights into the genomic response to short-term artificial selection

BACKGROUND

Short-term, intense artificial selection drives fast phenotypic changes in domestic animals and leaves imprints on their genomes. However, the genetic basis of this selection response is poorly understood. To better address this, we employed the Pekin duck Z2 pure line, in which the breast muscle weight was increased nearly 3-fold after 10 generations of breeding. We denovo assembled a high-quality reference genome of a female Pekin duck of this line (GCA_003850225.1) and identified 8.60 million genetic variants in 119 individuals among 10 generations of the breeding population.

RESULTS

We identified 53 selected regions between the first and tenth generations, and 93.8% of the identified variations were enriched in regulatory and noncoding regions. Integrating the selection signatures and genome-wide association approach, we found that 2 regions covering 0.36 Mb containing UTP25 and FBRSL1 were most likely to contribute to breast muscle weight improvement. The major allele frequencies of these 2 loci increased gradually with each generation following the same trend. Additionally, we found that a copy number variation region containing the entire EXOC4 gene could explain 1.9% of the variance in breast muscle weight, indicating that the nervous system may play a role in economic trait improvement.

CONCLUSIONS

Our study not only provides insights into genomic dynamics under intense artificial selection but also provides resources for genomics-enabled improvements in duck breeding.

The landscape of the long non-coding RNAs and circular RNAs of the abdominal fat tissues in the chicken lines divergently selected for fatness

BACKGROUND

Excessive deposition of abdominal fat poses serious problems in broilers owing to rapid growth. Recently, the evolution of the existing knowledge on long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) have established their indispensable roles in multiple physiological metabolic processes, including adipogenesis and fat deposition. However, not much has been explored on their profiles in the abdominal fat tissues of broilers to date. In the study, we aimed to characterize the vital candidates of lncRNAs and circRNAs and their underlying regulations for abdominal fat deposition in broilers.

RESULTS

The present study sequenced the lncRNAs and circRNAs expression profiles in the abdominal fat tissues isolated from 7-week-old broilers, who were divergently selected for their fatness. It identified a total of 3359 lncRNAs and 176 circRNAs, demonstrating differential expressed (DE) 30 lncRNAs and 17 circRNAs between the fat- and lean-line broilers (|log2FC| ≥ 1, P < 0.05). Subsequently, the 20 cis-targets and 48 trans-targets of the candidate DE lncRNAs were identified for depositing abdominal fat by adjacent gene analysis and co-expression analysis, respectively. In addition, the functional enrichment analysis showed the DE lncRNAs targets and DE circRNAs host genes to be mainly involved in the cellular processes, amino/fatty acid metabolism, and immune inflammation-related pathways and GO terms. Finally, the vital 16 DE lncRNAs located in cytoplasm and specifically expressed in fat/lean line and their targets were used to construct the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) regulatory network, comprising 7 DE lncRNAs, 28 miRNAs, 11 DE mRNAs. Notably, three lncRNAs including XR_001468036.2, XR_003077610.1 and XR_001466431.2 with the most connected degrees might play hub regulatory roles in abdominal fat deposition of broilers.

CONCLUSIONS

This study characterized the whole expression difference of lncRNAs and circRNAs between the two lines broilers with divergently ability of abdominal fat. The vital candidate DE lncRNAs/circRNAs and ceRNA regulations were identified related to the deposition of abdominal fat in chicken. These results might further improve our understanding of regulating the non-coding RNAs in obesity.

Whole Transcriptome Analysis Reveals a Potential Regulatory Mechanism of LncRNA-FNIP2/miR-24-3p/FNIP2 Axis in Chicken Adipogenesis

Lipid biosynthesis is a complex process, which is regulated by multiple factors including lncRNA. However, the role of lncRNA in chicken abdominal fat accumulation is still unclear. In this research, we collected liver tissues from six high abdominal fat rate Sanhuang broilers and six low abdominal fat rate Sanhuang broilers to perform lncRNA sequencing and small RNA sequencing. A total of 2,265 lncRNAs, 245 miRNAs, and 5,315 mRNAs were differently expressed. Among of them, 1,136 differently expressed genes were enriched in the metabolic process. A total of 36 differently expressed genes, which were considered as differently expressed lncRNAs' targets, were enriched in the metabolic process. In addition, we also found out that eight differently expressed miRNAs could target 19 differently expressed genes. FNIP2 and PEX5L were shared in a cis-regulatory network and a differently expressed miRNA target relationship network. LncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as a potential candidate that may participate in lipid synthesis. Experimentally, the objective reality of lncRNA-FNIP2/miR-24-3p/FNIP2 axis was clarified and the regulation effect of lncRNA-FNIP2/miR-24-3p/FNIP2 axis on synthesis was validated. In brief, our study reveals a potential novel regulatory mechanism that lncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as being involved in lipid synthesis during chicken adipogenesis in liver.

Molecular Regulation of Lipogenesis, Adipogenesis and Fat Deposition in Chicken

In the poultry industry, excessive fat deposition is considered an undesirable factor, affecting feed efficiency, meat production cost, meat quality, and consumer’s health. Efforts to reduce fat deposition in economically important animals, such as chicken, can be made through different strategies; including genetic selection, feeding strategies, housing, and environmental strategies, as well as hormone supplementation. Recent investigations at the molecular level have revealed the significant role of the transcriptional and post-transcriptional regulatory networks and their interaction on modulating fat metabolism in chickens. At the transcriptional level, different transcription factors are known to regulate the expression of lipogenic and adipogenic genes through various signaling pathways, affecting chicken fat metabolism. Alternatively, at the post-transcriptional level, the regulatory mechanism of microRNAs (miRNAs) on lipid metabolism and deposition has added a promising dimension to understand the structural and functional regulatory mechanism of lipid metabolism in chicken. Therefore, this review focuses on the progress made in unraveling the molecular function of genes, transcription factors, and more notably significant miRNAs responsible for regulating adipogenesis, lipogenesis, and fat deposition in chicken. Moreover, a better understanding of the molecular regulation of lipid metabolism will give researchers novel insights to use functional molecular markers, such as miRNAs, for selection against excessive fat deposition to improve chicken production efficiency and meat quality.

Effects of miR-125b-5p on Preadipocyte Proliferation and Differentiation in Chicken

Our previous studies have shown that miR-125b-5p was highly expressed and significantly upregulated during abdominal fat deposition in chickens. However, the role of miR-125b in the regulation of adipogenesis is not clear in chickens. Therefore, we evaluated the effects of miR-125b-5p on preadipocyte proliferation and differentiation and the interaction between miR-125b-5p and the acyl-CoA synthetase bubblegum family member 2 (ACSBG2) gene in adipogenesis in chicken abdominal adipose tissue. Here, transfection tests of miR-125b-5p mimic/inhibitor were performed in preadipocytes, and the effects of miR-125b-5p on preadipocytes proliferation and differentiation were analyzed. The target site of miR-125b-5p in the 3'UTR (untranslated region) of ACSBG2 were verified by a luciferase reporter assay. Our results showed that miR-125b-5p overexpression inhibited proliferation and reduced the number of cells in S phase and G2/M phase in preadipocytes; conversely, miR-125b-5p inhibition promoted the proliferation and increased the number of cells in S phase and G2/M phase. In adipocytes after induction, miR-125b-5p overexpression led to a notable increase in the accumulation of lipid droplets as well as in the concentration of triglycerides, while miR-125b-5p inhibition had the opposite effect. Furthermore, miR-125b-5p could directly bind to the 3'UTR of ACSBG2, and its overexpression could significantly repress the mRNA and protein expression of ACSBG2. These results indicate that miR-125b-5p can inhibit preadipocyte proliferation and can promote preadipocyte differentiation to affect adipogenesis in chicken abdominal adipose tissues, at least partially by downregulating ACSBG2.

Functional microRNA screening for dietary vitamin E regulation of abdominal fat deposition in broilers

1. Functional microRNA (miRNA) screening for abdominal fat tissue with different dietary vitamin E (VE) levels was performed to reveal miRNAs, genes and metabolic pathways involved in abdominal fat deposition in broilers. 2. A total of 240, one-day-old healthy female chicks were randomly allocated into five dietary treatments containing either 0, 20, 50, 75 or 100 IU DL-α-tocopherol acetate. The sequencing of miRNAs from abdominal fat tissues was performed. The target genes of miRNAs were predicted and enrichment analysis for these genes was performed. Diets supplemented with 50 IU VE significantly diminished abdominal fat deposition in broilers at day 35 of age. 3. A total of 29 miRNAs were differentially expressed between control and 50 IU VE treatment. Ten of the 23 target genes were enriched in four signalling pathways: tight junction, SNARE interactions in vesicular transport, regulation of autophagy and proteasome. 4. This study identified miRNA, target genes and pathways in dietary VE treatment for broilers, providing new insights into the miRNA regulation of abdominal fat deposition in broilers.

GATA Binding Protein 3 Is a Direct Target of Kruppel-Like Transcription Factor 7 and Inhibits Chicken Adipogenesis

Kruppel-like transcription factor 7 (KLF7) is a negative regulator of adipogenesis, however, its precise mechanism is poorly understood. Our previous KLF7 ChIP-seq analysis showed that one of the KLF7 binding peaks was present upstream of GATA binding protein 3 (GATA3) in chicken preadipocytes. In the present study, we identified GATA3 as a target of KLF7. Overexpression analysis showed KLF7 markedly enhanced the endogenous expression of GATA3 in the immortalized chicken preadipcyte cell line (ICP2), and the luciferase reporter assay showed that KLF7 overexpression increased the reporter gene activity of the cloned upstream region (-5285/-4336 relative to the translation initiation codon ATG) of GATA3 in ICP2 and DF1 cells, and mutation of the putative KLF7 binding site abolished the promotive effect of KLF7 overexpression on the reporter gene activity of the cloned GATA3 upstream region. ChIP-qPCR further demonstrated that KLF7 directly bound to the GATA3 upstream region. Gene expression analysis showed that GATA3 mRNA expression in abdominal adipose tissue was significantly higher in lean chicken line than in the fat line at 2, 3, and 6 weeks of age. In addition, GATA3 mRNA expression markedly decreased during the preadipocyte differentiation. Furthermore, a functional study showed that GATA3 overexpression inhibited the differentiation of the ICP2 cells. Taken together, our results demonstrated that KLF7 inhibits chicken adipogenesis, at least in part through direct upregulation of GATA3.

MicroRNAs and their regulatory networks in Chinese Gushi chicken abdominal adipose tissue during postnatal late development

Background

Abdominal fat is the major adipose tissue in chickens. The growth status of abdominal fat during postnatal late development ultimately affects meat yield and quality in chickens. MicroRNAs (miRNAs) are endogenous small noncoding RNAs that regulate gene expression at the post-transcriptional level. Studies have shown that miRNAs play an important role in the biological processes involved in adipose tissue development. However, few studies have investigated miRNA expression profiles and their interaction networks associated with the postnatal late development of abdominal adipose tissue in chickens.

Results

We constructed four small RNA libraries from abdominal adipose tissue obtained from Chinese domestic Gushi chickens at 6, 14, 22, and 30 weeks. A total of 507 known miRNAs and 53 novel miRNAs were identified based on the four small RNA libraries. Fifty-one significant differentially expressed (SDE) miRNAs were identified from six combinations by comparative analysis, and the expression patterns of these SDE miRNAs were divided into six subclusters by cluster analysis. Gene ontology enrichment analysis showed that the SDE miRNAs were primarily involved in the regulation of fat cell differentiation, regulation of lipid metabolism, regulation of fatty acid metabolism, and unsaturated fatty acid metabolism in the lipid metabolism- or deposition-related biological process categories. In addition, we constructed differentially expressed miRNA–mRNA interaction networks related to abdominal adipose development. The results showed that miRNA families, such as mir-30, mir-34, mir-199, mir-8, and mir-146, may have key roles in lipid metabolism, adipocyte proliferation and differentiation, and cell junctions during abdominal adipose tissue development in chickens.

Conclusions

This study determined the dynamic miRNA transcriptome and characterized the miRNA–mRNA interaction networks in Gushi chicken abdominal adipose tissue for the first time. The results expanded the number of known miRNAs in abdominal adipose tissue and provide novel insights and a valuable resource to elucidate post-transcriptional regulation mechanisms during postnatal late development of abdominal adipose tissue in chicken.

MicroRNA-15a Regulates the Differentiation of Intramuscular Preadipocytes by Targeting ACAA1, ACOX1 and SCP2 in Chickens

Our previous studies showed that microRNA-15a (miR-15a) was closely related to intramuscular fat (IMF) deposition in chickens; however, its regulatory mechanism remains unclear. Here, we evaluated the expression characteristics of miR-15a and its relationship with the expression of acetyl-CoA acyltransferase 1 (ACAA1), acyl-CoA oxidase 1 (ACOX1) and sterol carrier protein 2 (SCP2) by qPCR analysis in Gushi chicken breast muscle at 6, 14, 22, and 30 weeks old, where we performed transfection tests of miR-15a mimics in intramuscular preadipocytes and verified the target gene of miR-15a in chicken fibroblasts (DF1). The miR-15a expression level at 30 weeks increased 13.5, 4.5, and 2.7-fold compared with the expression levels at 6, 14, and 22 weeks, respectively. After 6 days of induction, miR-15a over-expression significantly promoted intramuscular adipogenic differentiation and increased cholesterol and triglyceride accumulation in adipocytes. Meanwhile, 48 h after transfection with miR-15a mimics, the expression levels of ACAA1, ACOX1 and SCP2 genes decreased by 56.52%, 31.18% and 37.14% at the mRNA level in intramuscular preadipocytes. In addition, the co-transfection of miR-15a mimics and ACAA1, ACOX1 and SCP2 3′UTR (untranslated region) dual-luciferase vector significantly inhibited dual-luciferase activity in DF1 cells. Taken together, our data demonstrate that miR-15a can reduce fatty acid oxidation by targeting ACAA1, ACOX1, and SCP2, which subsequently indirectly promotes the differentiation of chicken intramuscular preadipocytes.

Hyperglycemia Affects miRNAs Expression Pattern during Adipogenesis of Human Visceral Adipocytes-Is Memorization Involved?

microRNAs are increasingly analyzed in adipogenesis, whose deregulation, especially visceral, contributes to the development of diabetes. Hyperglycemia is known to affect cells while occurring acutely and chronically. Therefore, we aimed to evaluate the effect of hyperglycemia on human visceral pre/adipocytes from the perspective of microRNAs. The relative expression of 78 microRNAs was determined by TaqMan Low Density Arrays at three stages of HPA-v adipogenesis conducted under normoglycemia, chronic, and intermittent hyperglycemia (30 mM). Hierarchical clustering/Pearson correlation revealed the relationship between various microRNAs' expression profiles, while functional analysis identified the genes and signaling pathways regulated by differentially expressed microRNAs. Hyperglycemia affected microRNAs' expression patterns during adipogenesis, and at the stage of pre-adipocytes, differentiated and matured adipocytes compared to normoglycemia. Interestingly, the changes that were evoked upon hyperglycemic exposure during one adipogenesis stage resembled those observed upon chronic hyperglycemia. At least 15 microRNAs were modulated during normoglycemic and/or hyperglycemic adipogenesis and/or upon intermittent/chronic hyperglycemia. Bioinformatics analysis revealed the involvement of these microRNAs in cell cycles, lipid metabolism, ECM⁻receptor interaction, oxidative stress, signaling of insulin, MAPK, TGF-β, p53, and more. The obtained data suggests that visceral pre/adipocytes exposed to chronic/intermittent hyperglycemia develop a microRNAs' expression pattern, which may contribute to further visceral dysfunction, the progression of diabetic phenotype, and diabetic complications possibly involving "epi"-memory.

Genomic Insights Into the Multiple Factors Controlling Abdominal Fat Deposition in a Chicken Model

Genetic selection for an increased growth rate in meat-type chickens has been accompanied by excessive fat accumulation particularly in abdominal cavity. These progressed to indirect and often unhealthy effects on meat quality properties and increased feed cost. Advances in genomics technology over recent years have led to the surprising discoveries that the genome is more complex than previously thought. Studies have identified multiple-genetic factors associated with abdominal fat deposition. Meanwhile, the obesity epidemic has focused attention on adipose tissue and the development of adipocytes. The aim of this review is to summarize the current understanding of genetic/epigenetic factors associated with abdominal fat deposition, or as it relates to the proliferation and differentiation of preadipocytes in chicken. The results discussed here have been identified by different genomic approaches, such as QTL-based studies, the candidate gene approach, epistatic interaction, copy number variation, single-nucleotide polymorphism screening, selection signature analysis, genome-wide association studies, RNA sequencing, and bisulfite sequencing. The studies mentioned in this review have described multiple-genetic factors involved in an abdominal fat deposition. Therefore, it is inevitable to further study the multiple-genetic factors in-depth to develop novel molecular markers or potential targets, which will provide promising applications for reducing abdominal fat deposition in meat-type chicken.

The differentiation of preadipocytes and gene expression related to adipogenesis in ducks (Anas platyrhynchos)

Meat quality is closely related to adipose tissues in ducks, and adipogenesis is controlled by a complex network of transcription factors tightly acting at different stages of differentiation especially in ducks. The aim of this study was to establish the preadipocyte in vitro culture system and understand the biological characteristics of expansion of duck adipocyte tissue at the cellular and molecular level. We isolated pre-adipocytes from the subcutaneous fat of three breeds of duck and differentiated them into mature adipocytes using a mixture of insulin, rosiglitazone, dexamethasone, 3-isobutyl-1-methylxanthine, and oleic acid over 0,2, 4, 6, and 8 days. Successful differentiation was confirmed from the development of lipid droplets and their response to Oil Red O, and increasing numbers of lipid droplets were stained red over time. The expression of key marker genes, including peroxisome proliferator activated receptor γ (PPARγ), CCAAT/enhancer binding protein-α (C/EBPα), adipocyte fatty acid binding protein 4 (FABP4), and fatty acid synthetase (FAS), gradually increased during pre-adipocyte differentiation. Furthermore, it was verified by interference experiments that the knockdown of PPARγ directly reduced lipid production. Meanwhile we analyzed the role of unsaturated fatty acids in the production of poultry fat using different concentrations of oleic acid and found that lipid droplet deposition was highest when the concentration of oleic acid was 300 μM. We also compared the level of differentiated pre-adipocytes that were isolated from Jianchang ducks (fatty-meat duck), Cherry Valley ducks (lean-meat duck) and White-crested ducks (egg-producing duck). The proliferation and differentiation rate of pre-adipocytes derived from Jianchang ducks was higher than that of White-crested ducks. These results provide the foundation for further research into waterfowl adipogenesis.

MiR-17-5p and miR-20a promote chicken cell proliferation at least in part by upregulation of c-Myc via MAP3K2 targeting

The miR-17-92 cluster has been well studied in mammals but less extensively studied in birds. Here, we demonstrated that miR-17-92 cluster overexpression promoted the proliferation of DF1 cells and immortalized chicken preadipocytes (ICPA-1), and miR-17-5p and miR-20a, members of the miR-17-92 cluster, targeted MAP3K2. Further analysis showed that MAP3K2 overexpression reduced the proliferation of DF1 and ICPA-1 cells and attenuated the promotive effect of the miR-17-92 cluster on cell proliferation. Downstream gene expression analysis of the MAPK signalling pathway showed that MAP3K2 overexpression decreased c-Myc expression; in contrast, MAP3K2 knockdown using RNA interference and miR-17-92 cluster overexpression increased c-Myc expression. Furthermore, c-Myc overexpression promoted miR-17-92 cluster expression and DF1 cell proliferation. Taken together, these data indicated that miR-17-92 promotes chicken cell proliferation at least in part by the upregulation of c-Myc via targeting MAP3K2, and the miR-17-92 cluster, c-Myc and E2F1 form a complex regulatory network in chicken cell proliferation.

MicroRNA-101-2-5p targets the ApoB gene in the liver of chicken (Gallus Gallus)

Apolipoprotein B (ApoB) is a major protein component of plasma lipoproteins. It is involved in many important biological processes such as lipid transportation, enzyme activity regulation, and receptor recognition. Extensive studies have shown that the expression of ApoB is regulated at multiple levels. However, the regulation of ApoB expression by microRNAs (miRNAs) still remains unknown. In the present study, identified are miRNAs that are predicted to interact with ApoB in chicken. The predicted relationship between the identified miRNAs and ApoB was verified through dual luciferase reporter assay in chicken DF1 cells, and the effect of miRNAs on ApoB expression was analyzed in chicken embryo hepatocytes stimulated by 17β-estradiol. The results show that miR-101-2-5p was predicted to interact with ApoB. Dual luciferase reporter assay together with the miR-101-2-5p mimics study demostrate that ApoB is the target of miR-101-2-5p, which suppresses the expression of ApoB through binding with the 3′UTR of ApoB. Our experiments suggest that miR-101-2-5p might be involved in lipid metabolism through binding to the 3′UTR of ApoB in the liver of egg-laying chickens.

Genome-Wide Analysis of lncRNA and mRNA Expression During Differentiation of Abdominal Preadipocytes in the Chicken

Long noncoding RNAs (lncRNAs) regulate adipogenesis and other processes associated with metabolic tissue development and function. However, little is known about the function and profile of lncRNAs during preadipocyte differentiation in the chicken (Gallus gallus). Herein, lncRNA and mRNA expression in preadipocytes at different stages of differentiation were analyzed using RNA sequencing. A total of 1,300,074,528 clean reads and 27,023 novel lncRNAs were obtained from 12 samples. The number of genes (1336 lncRNAs and 1759 mRNAs; 3095 in total) differentially expressed across various stages declined as differentiation progressed. Differentially expressed genes were found to be involved in several pathways related to preadipocyte differentiation that have been extensively studied, including glycerolipid metabolism, and the mammalian target of rapamycin, peroxisome proliferator-activated receptor, and mitogen-activated protein kinase signaling pathways. To our knowledge, some pathways are being reported for the first time, including the propanoate metabolism, fatty acid metabolism, and oxidative phosphorylation pathways. Furthermore, 3095 differentially expressed genes were clustered into eight clusters, and their expression patterns were determined through K-means clustering. Genes involved in the K2 cluster likely play important roles in preadipocyte differentiation. Six stage-specific modules related to A0 (day 0), A2 (day 2), and A6 (day 6) stages were identified, using weighted coexpression network analysis. Nine central, highly connected .genes in stage-specific modules were subsequently identified, including XLOC_068731, XLOC_022661, XLOC_045161, XLOC_070302, CHD6, LLGL1, NEURL1B, KLHL38, and ACTR6. This study provides a valuable resource for further study of chicken lncRNA and facilitates a better understanding of preadipocyte differentiation in the chicken

miR-429 Inhibits Differentiation and Promotes Proliferation in Porcine Preadipocytes

MicroRNAs (miRNAs) are crucial regulatory molecules for adipogenesis. They contribute to the controlling of proliferation and differentiation of preadipocytes. Previous studies revealed an important role of miR-429 in cell invasion, migration, and apoptosis. Our previous work has shown that the expression of miR-429 in subcutaneous fat can be observed in newly born (3-day-old) Rongchang piglets rather than their adult counterparts (180-day-old). This expression pattern suggests that miR-429 might be functionally related to postnatal adipogenesis. However, we currently lack a mechanistic understanding of miR-429 within the context of preadipocyte differentiation. In this study, we investigated the function of miR-429 in porcine subcutaneous and intramuscular preadipocyte proliferation and differentiation. In our porcine preadipocyte differentiation model, miR-429 expression decreased remarkably upon adipogenic induction. Overexpression of miR-429 notably down-regulated the expression of adipogenic marker genes: PPARγ, aP2, FAS and impaired the triglyceride accumulation, while the expression of lipolytic gene ATGL was not affected. In addition, we observed that miR-429 significantly promoted the proliferation of porcine preadipocytes. We also found that miR-429 could directly bind to the 3′-UTRs of KLF9 and p27, which have been well documented to promote preadipocyte differentiation and repress cell cycle progression. Taken together, our data support a novel role of miR-429 in regulating porcine preadipocyte differentiation and proliferation, and KLF9 and p27 are potent targets of miR-429 during these processes.

Differential fecal microbiota are retained in broiler chicken lines divergently selected for fatness traits

Our study combined 16S rRNA-pyrosequencing and whole genome sequencing to analyze the fecal metagenomes of the divergently selected lean (LL) and fat (FL) line chickens. Significant structural differences existed in both the phylogenic and functional metagenomes between the two chicken lines. At phylum level, the FL group had significantly less Bacteroidetes. At genus level, fourteen genera of different relative abundance were identified, with some known short-chain fatty acid producers (including Subdoligranulum, Butyricicoccus, Eubacterium, Bacteroides, Blautia) and a potentially pathogenic genus (Enterococcus). Redundancy analysis identified 190 key responsive operational taxonomic units (OTUs) that accounted for the structural differences between the phylogenic metagenome of the two groups. Four Cluster of Orthologous Group (COG) categories (Amino acid transport and metabolism, E; Nucleotide transport and metabolism, F; Coenzyme transport and metabolism, H; and Lipid transport and metabolism, I) were overrepresented in LL samples. Fifteen differential metabolic pathways (Biosynthesis of amino acids, Pyruvate metabolism, Nitrotoluene degradation, Lipopolysaccharide biosynthesis, Peptidoglycan biosynthesis, Pantothenate and CoA biosynthesis, Glycosaminoglycan degradation, Thiamine metabolism, Phosphotransferase system, Two-component system, Bacterial secretion system, Flagellar assembly, Bacterial chemotaxis, Ribosome, Sulfur relay system) were identified. Our data highlighted interesting variations between the gut metagenomes of these two chicken lines.

Gga-miR-21 inhibits chicken pre-adipocyte proliferation in part by down-regulating Kruppel-like factor 5

Gga-miR-21 is abundantly expressed in chicken pre-adipocytes, but its role is unclear. The present study investigated the role of gga-miR-21 in chicken pre-adipocyte proliferation. Cell proliferation assay and gene expression analysis of proliferating cell nuclear antigen (PCNA) showed that the gga-miR-21 mimic inhibited pre-adipocyte proliferation. In contrast, the gga-miR-21 inhibitor enhanced pre-adipocyte proliferation. The subsequent investigation identified Kruppel-like factor 5 (KLF5) mRNA as a target of gga-miR-21. The gga-miR-21 mimic inhibited KLF5 3^'UTR reporter activity and decreased endogenous KLF5 expression in primary pre-adipocytes. KLF5 knockdown using RNAi had a similar effect to that of the gga-miR-21 mimic on cell proliferation. The promoting effect of the gga-miR-21 inhibitor on pre-adipocyte proliferation was partially attenuated by KLF5 knockdown. Taken together, our results demonstrated that miR-21 inhibits chicken pre-adipocyte proliferation, at least in part, by targeting KLF5.

Chicken gga-miR-130a targets HOXA3 and MDFIC and inhibits Marek's disease lymphoma cell proliferation and migration

Marek's disease (MD) is an infectious disease of chickens caused by MD virus (MDV), which is a herpesvirus that initiates tumor formation. Studies have indicated that microRNAs (miRNAs) are linked with the development of cancers or tumors. Previously, gga-miR-130a was discovered downregulated in MDV-infected tissues. Here, we aimed to explore the further function of gga-miR-130a in MD. The expression of gga-miR-130a in MDV-infected and uninfected spleens was detected by quantitative real-time PCR (qRT-PCR). Subsequently, proliferation and migration assays of MDV-transformed lymphoid cells (MSB1) were carried out by transfecting gga-miR-130a. The target genes of gga-miR-130a were predicted using TargetScan and miRDB and clustered through Gene Ontology analysis. The target genes were validated by western blot, qRT-PCR, and a dual luciferase reporter assay. Our results show that the expression of gga-miR-130a was reduced in MDV-infected spleens. Gga-miR-130a showed an inhibitory effect on MSB1 cell proliferation and migration. Two target genes, homeobox A3 (HOXA3) and MyoD family inhibitor domain containing (MDFIC), were predicted and clustered to cell proliferation. Results indicate that gga-miR-130a regulates HOXA3 and MDFIC at the protein level but not at the mRNA level. Moreover, the gga-miR-130a binding sites of two target genes have been confirmed. We conclude that gga-miR-130a can arrest MSB1 cell proliferation and migration, and target HOXA3 and MDFIC, which are both involved in the regulation of cell proliferation. Collectively, gga-miR-130a plays a critical role in the tumorigenesis associated with chicken MD.

Chicken gga-miR-103-3p Targets CCNE1 and TFDP2 and Inhibits MDCC-MSB1 Cell Migration

Marek’s disease (MD) is a highly contagious viral neoplastic disease caused by Marek’s disease virus (MDV), which can lead to huge economic losses in the poultry industry. Recently, microRNAs (miRNAs) have been found in various cancers and tumors. In recent years, 994 mature miRNAs have been identified through deep sequencing in chickens, but only a few miRNAs have been investigated further in terms of their function. Previously, gga-miR-103-3p was found downregulated in MDV-infected samples by using Solexa deep sequencing. In this study, we further verified the expression of gga-miR-103-3p among MDV-infected spleen, MD lymphoma from liver, noninfected spleen, and noninfected liver, by qPCR. The results showed that the expression of gga-miR-103-3p was decreased in MDV-infected tissues, which was consistent with our previous study. Furthermore, two target genes of gga-miR-103-3p, cyclin E1 (CCNE1) and transcription factor Dp-2 (E2F dimerization partner 2) (TFDP2), were predicted and validated by luciferase reporter assay, qPCR, and western blot analysis. The results suggested that CCNE1 and TFDP2 are direct targets of gga-miR-103-3p in chickens. Subsequent cell proliferation and migration assay showed that gga-miR-103-3p suppressed MDCC-MSB1 migration, but did not obviously modulate MDCC-MSB1 cell proliferation. In conclusion, gga-miR-103-3p targets the CCNE1 and TFDP2 genes, and suppresses cell migration, which indicates that it might play an important role in MD tumor transformation.