Comprehensive analysis of differentially expressed mRNA profiles in chicken jejunum and cecum following Eimeria maxima infection.

Coccidiosis, a protozoan disease that substantially impacts poultry production, is characterized by an intracellular parasite. The study utilized 48 one-day-old Horro chickens, randomly divided into the infected (I) and control (C) groups. The challenge group of chickens were administered Eimeria maxima oocysts via oral gavage at 21-days-old, and each chicken received 2 mL containing 7×104 sporulated oocysts. The total RNAs of chicken jejunum and cecum tissues were isolated from three samples, each from I and C groups. Our study aimed to understand the host immune-parasite interactions and compare immune response mRNA profiles in chicken jejunum and cecum tissues at 4 and 7 days postinfection with Eimeria maxima. The results showed that 823 up- and 737 down-regulated differentially expressed mRNAs (DEmRNAs) in jejunum at 4 d infection and control (J4I vs. J4C), and 710 up- and 368 down-regulated DEmRNAs in jejunum at 7 days infection and control (J7I vs. J7C) were identified. In addition, DEmRNAs in cecum tissue, 1424 up- and 1930 down-regulated genes in cecum at 4 days infection and control (C4I vs. C4C), and 77 up- and 191 down-regulated genes in cecum at 7 days infection and control (C7I vs. C7C) were detected. The crucial DEmRNAs, including SLC7A5, IL1R2, GLDC, ITGB6, ADAMTS4, IL1RAP, TNFRSF11B, IMPG2, WNT9A, and FOXF1, played pivotal roles in the immune response during Eimeria maxima infection of chicken jejunum. In addition, the potential detection of FSTL3, RBP7, CCL20, DPP4, PRKG2, TFPI2, and CDKN1A in the cecum during the host immune response against Eimeria maxima infection is particularly noteworthy. Furthermore, our functional enrichment analysis revealed the primary involvement of DEmRNAs in small molecule metabolic process, immune response function, inflammatory response, and toll-like receptor 10 signaling pathway in the jejunum at 4 and 7 days postinfection. Similarly, in the cecum, DEmRNAs at 4 and 7 days postinfection were enriched in processes related to oxidative stress response and immune responses. Our findings provide new insights and contribute significantly to the field of poultry production and parasitology.

MicroRNA expression profile of chicken jejunum in different time points Eimeria maxima infection.

Coccidiosis stands as a protozoan disease of notable economic impact, characterized by an intracellular parasite that exerts substantial influence over poultry production. This invasion disrupts the integrity of the enteric mucosa, leading to the emergence of severe lesions and diminishing the efficiency of feed utilization in chickens. MicroRNA (miRNA) are short, non-coding RNA molecules with approximately 21-24 nucleotides long in size that play essential roles in various infectious diseases and inflammatory responses. However, the miRNA's expression patterns and roles in the context of Eimeria maxima infection of chicken intestines remain unclear. miRNA sequencing was employed to assess the miRNA expression profile in chicken jejunum during E. maxima infection. In this study, we analyzed miRNA expression profiles related to the host's immune response in the chicken jejunum during E. maxima infection. At 4 days infection and control (J4I versus J4C), 21 differentially expressed miRNAs in the jejunum were identified, comprising 9 upregulated and 12 downregulated miRNAs. Furthermore, in the jejunum, at 7 days infection and control (J7I versus J7C) groups, a total of 35 significantly differentially expressed miRNAs were observed, with 13 upregulated and 22 downregulated miRNAs. The regulatory networks were constructed between differentially expressed miRNA and mRNAs to offer insight into the interaction mechanisms between chickens and E. maxima coccidian infection. Furthermore, within the comparison group, we obtained 946, 897, and 281 GO terms that exhibited significant enrichment associated with host immunity in the following scenarios, J4I vs. J4C, J7I vs. J7C, and J4I vs. J7I, respectively. The KEGG pathway analysis indicated notable enrichment of differentially expressed miRNAs in the jejunum, particularly in J4I vs. J4C; enriched pathways included metabolic pathways, endocytosis, MAPK signaling pathway, regulation of actin cytoskeleton, and cytokine-cytokine receptor interaction. Moreover, in J7I vs. J7C, the KEGG pathway was significantly enriched, including metabolic pathways, protein processing in the endoplasmic reticulum, ubiquitin-mediated proteolysis, and FoxO signaling pathway. A comprehensive understanding of the host genetic basis of resistance with a combination of time-dependent infection to the Eimeria parasite is crucial for pinpointing resistance biomarkers for poultry production.

Transcriptome analysis of differentially expressed circRNAs miRNAs and mRNAs during the challenge of coccidiosis.

Avian coccidiosis is a common enzootic disease caused by infection of Eimeria species parasites. It causes huge economic losses in the global poultry industry. Current control using anticoccidial drugs or vaccination is limited due to drug resistance and the relatively high cost of vaccines. Improving host genetic resistance to Eimeria species is considered an effective strategy for improved control of coccidiosis. Circular RNAs (circRNAs) have been found to function as biomarkers or diagnoses of various kinds of diseases. The molecular biological functions of circRNAs, miRNAs, and mRNAs related to Sasso chicken have not yet been described during Eimeria species challenge. In this study, RNA-seq was used to profile the expression pattern of circRNAs, miRNAs, and mRNAs in spleens from Eimeria tenella-infected and non-infected commercial dual-purpose Sasso T445 breed chickens. Results showed a total of 40 differentially expressed circRNAs (DEcircRNAs), 31 differentially expressed miRNAs (DEmiRNAs), and 820 differentially expressed genes (DEmRNAs) between infected and non-infected chickens. Regulatory networks were constructed between differentially expressed circRNAs, miRNAs, and mRNAs to offer insights into the interaction mechanisms between chickens and Eimeria spp. Functional validation of a significantly differentially expressed circRNA, circMGAT5, revealed that circMGAT5 could sponge miR-132c-5p to promote the expression of the miR-132c-5p target gene monocyte to macrophage differentiation-associated (MMD) during the infection of E. tenella sporozoites or LPS stimulation. Pathologically, knockdown of circMGAT5 significantly upregulated the expression of macrophage surface markers and the macrophage activation marker, F4/80 and MHC-II, which indicated that circMGAT5 might inhibit the activation of macrophage. miR-132c-5p markedly facilitated the expression of F4/80 and MHC-II while circMGAT5 could attenuate the increase of F4/80 and MHC-II induced by miR-132c-5p, indicating that circMGAT5 exhibited function through the circMGAT5-miR-132c-5p-MMD axis. Together, our results indicate that circRNAs exhibit their resistance or susceptive roles during E. tenella infection. Among these, circMGAT5 may inhibit the activation of macrophages through the circMGAT5-miR-132c-5p-MMD axis to participate in the immune response induced by Eimeria infection.

Influence of Eimeria maxima coccidia infection on gut microbiome diversity and composition of the jejunum and cecum of indigenous chicken.

Coccidiosis is an economically significant protozoan disease and an intracellular parasite that significantly impacts poultry production. The gastrointestinal tract microbiota plays a central role in host health and metabolism, and these microbes enhance chickens' immune systems and nutrient absorption. In this study, we analyzed the abundance and diversity of microbiota of the jejunum and cecum of a dual-purpose indigenous Horro chicken following Eimeria maxima infection. We compared microbial abundance, composition, and diversity at the 4- and 7- days post-infection using 16S rRNA gene sequencing. We obtained, on average, 147,742 and 132,986 high-quality sequences per sample for jejunum and cecum content, respectively. Firmicutes, Proteobacteria, Campilobacterota and Bacteroidota were the major microbial phylum detected in the jejunum content. Firmicutes were the dominant phylum for 4- and 7-days jejunum control groups accounting for (>60% of the sequences). In the infected group Campilobacterota was the dominant phylum in the jejunum (> 24% of sequences) at 4-and 7-days post-infection groups, while Proteobacteria was predominant at 4- and 7-days post-infection of the cecum (> 40% of the sequences). The microbial genus Lactobacillus and Helicobacter were found in the jejunum, while Alistipes, Barnesiella and Faecalibacterium were detected in the cecum. In the jejunum, Helicobacter was dominant at 4 -and-7 days post-infection (≥24%), and Lactobacillus was dominant at 4 -and 7- days in the control group (> 50%). In 4- and 7-days post-infection, Alistipes genus was the more prevalent (> 38%) in the cecum. Thus, clear differences were observed in the bacterial microbiota distribution and abundance between the jejunum and cecum, as well as between infected and control groups for both tissues. The results indicate that chicken intestinal microbial imbalance (dysbiosis) is associated with Eimeria parasite infection and will likely affect the host-microbial non-pathogenic and pathogenic molecular interactions.

Long noncoding RNA SMUL suppresses SMURF2 production-mediated muscle atrophy via nonsense-mediated mRNA decay.

As the world population grows, muscle atrophy leading to muscle wasting could become a bigger risk. Long noncoding RNAs (lncRNAs) are known to play important roles in muscle growth and muscle atrophy. Meanwhile, it has recently come to light that many putative small open reading frames (sORFs) are hidden in lncRNAs; however, their translational capabilities and functions remain unclear. In this study, we uncovered 104 myogenic-associated lncRNAs translated, in at least a small peptide, by integrated transcriptome and proteomic analyses. Furthermore, an upstream ORF (uORF) regulatory network was constructed, and a novel muscle atrophy-associated lncRNA named SMUL (Smad ubiquitin regulatory factor 2 [SMURF2] upstream lncRNA) was identified. SMUL was highly expressed in skeletal muscle, and its expression level was downregulated during myoblast differentiation. SMUL promoted myoblast proliferation and suppressed differentiation in vitro. In vivo, SMUL induced skeletal muscle atrophy and promoted a switch from slow-twitch to fast-twitch fibers. In the meantime, translation of the SMUL sORF disrupted the stability of SMURF2 mRNA. Mechanistically, SMUL restrained SMURF2 production via nonsense-mediated mRNA decay (NMD), participating in the regulation of the transforming growth factor ß (TGF-ß)/SMAD pathway and further regulating myogenesis and muscle atrophy. Taken together, these results suggest that SMUL could be a novel therapeutic target for muscle atrophy.

Cellular function of chicken FOXO3 and its associations with chicken growth.

FOXO3 belongs to the Forkhead O transcription factor family and it is an important gene in multiple biological processes, such as cell cycle control, cell proliferation, cell apoptosis, human longevity, and oxidative stress. Previous studies have shown that FOXO3 is associated with skeletal muscle growth and adipose development in mammals. However, the sequence of chicken FOXO3 is still incomplete and the cellular functions of FOXO3 in chickens are poorly understood. Thus, we obtained the full-length sequence of chicken FOXO3 by 5' rapid amplification of cDNA ends (5' RACE) and the phylogenetic tree showed that the chicken FOXO3 sequence was homologous with those in other species. Flow cytometry analysis and 5-ethynyl-2'-deoxyuridine assays showed that FOXO3 repressed cellular proliferation and induced apoptosis in a chicken hepatocellular carcinoma cell line (LMH). Mutations were screened in the second exon of FOXO3 and 13 synonymous single nucleotide polymorphisms were found in the test population. Further analysis showed that rs317670452 and rs15379317 were associated with many growth and carcass traits, such as the body weight at different ages and breast muscle weight. Our results indicate that chicken FOXO3 has similar cellular functions to those found in mammals and it is significantly associated with chicken growth.

A Novel Dnmt3a1 Transcript Inhibits Adipogenesis.

DNA (cytosine-5)-methyltransferase 3a (Dnmt3a) is an enzyme that catalyzes the transfer of methyl groups to specific CpG forms in DNA. In mammals, two variant transcripts of Dnmt3a have been successfully identified. To the best of our knowledge, no Dnmt3a transcripts in an avian have been successfully identified. This study was performed to detect different transcripts of Dnmt3a in chickens and to examine whether a novel Dnmt3a transcript named Dnmt3a1 may regulate adipogenesis. In addition to cloning, sequencing, transcript detection, and expression studies, a novel Dnmt3a1 transcript overexpression and knockdown were conducted to explore the potential role of Dnmt3a1 in preadipocyte proliferation and the early stage of adipocyte differentiation. In chicken abdominal fat tissue, we detected a novel Dnmt3a1 transcript that differs from Dnmt3a by lacking 23 amino acids at the exon-1/exon-2 border. Dnmt3a1 mRNA was ubiquitously expressed in a variety of tissues or cells and highly expressed in chicken adipose tissue/cells. The expression of Dnmt3a1 was regulated under different physiological conditions including aging, fasting, and high-fat diet. In addition, overexpression of Dnmt3a1 significantly decreased preadipocyte proliferation and induced cell-cycle arrest while its inhibition increased cell proliferation and S-phase cells. Furthermore, the overexpression of Dnmt3a1 significantly upregulated the mRNA level of cell-cycle-related genes, such as CDKN1A, CDKN1B, CCNB3, CCND2, CCNG2, CDKN2B, and CDK9, or the protein level of CDKN1A, CDKN1B, and CCNG2. Conversely, the knockdown of Dnmt3a1 by siRNA had the opposite effects. Moreover, during early adipocyte differentiation, the overexpression of Dnmt3a1 significantly decreased the mRNA and the protein levels of PPAR-γ, C/EBP-α, ADIPOR1, and STAT3, and the mRNA levels of FAS, LEPR, LPL, PRKAB2, and ATGL. In contrast, their expression was significantly increased after the knockdown of Dnmt3a1. Taken together, we identified a novel transcript of Dnmt3a, and it played a potential role in adipogenesis.

Characterization of miRNA and their target gene during chicken embryo skeletal muscle development.

MicroRNAs (miRNAs) are non-coding RNAs that regulate mRNA expression by degradation or translational inhibition. We investigated the underlying molecular mechanisms of skeletal muscle development based on differentially expressed genes and miRNAs. We compared mRNA and miRNA from chicken skeletal muscle at embryonic day E11, E16 and one day post-hatch (P1). The interaction networks were constructed, according to target prediction results and integration analysis of up-regulated genes with down regulated miRNAs or down-regulated genes with up-regulated miRNAs with |log2fold change| ≥ 1.75, P < 0.005. The miRNA-mRNA integration analysis showed high number of mRNAs regulated by a few number of miRNAs. In the E11_VS_E16, comparison group we identified biological processes including muscle maintenance, myoblast proliferation and muscle thin filament formation. The E11_VS_P1 group comparison included negative regulation of axon extension, sarcomere organization, and cell redox homeostasis and kinase inhibitor activity. The E16_VS_P1 comparison group contained genes for the negative regulation of anti-apoptosis and axon extension as well as glomerular basement membrane development. Functional in vitro assays indicated that over expression of miR-222a and miR-126-5p in DF-1 cells significantly reduced the mRNA levels of the target genes CPEB3 and FGFR3, respectively. These integrated analyses provide several candidates for future studies concerning miRNAs-target function on regulation of embryonic muscle development and growth.

Circular RNAs are abundant and dynamically expressed during embryonic muscle development in chickens.

The growth and development of skeletal muscle is regulated by proteins as well as non-coding RNAs. Circular RNAs (circRNAs) are universally expressed in various tissues and cell types, and regulate gene expression in eukaryotes. To identify the circRNAs during chicken embryonic skeletal muscle development, leg muscles of female Xinghua (XH) chicken at three developmental time points 11 embryo age (E11), 16 embryo age (E16) and 1 day post hatch (P1) were performed RNA sequencing. We identified 13,377 circRNAs with 3,036 abundantly expressed and most were derived from coding exons. A total of 462 differentially expressed circRNAs were identified (fold change > 2; q-value < 0.05). Parental genes of differentially expressed circRNAs were related to muscle biological processes. There were 946 exonic circRNAs have been found that harbored one or more miRNA-binding site for 150 known miRNAs. We validated that circRBFOX2s promoted cell proliferation through interacted with miR-206. These data collectively indicate that circRNAs are abundant and dynamically expressed during embryonic muscle development and could play key roles through sequestering miRNAs as well as other functions.

Molecular characterization, expression profile of the FSHRgene and its association with egg production traits in muscovy duck.

Follicle-stimulating hormone (FSH) and its receptor play a key role in the follicular development and regulation of steroidogenesis in the ovary and spermatogenesis in the testis. The purpose of this study was to characterize themuscovy duck FSHR gene, identify SNPs and their association with egg production traits in muscovy ducks. Here, we cloned the complementary DNA (cDNA) sequence of FSHR, and examined the expression patterns of FSHR gene in adult female muscovy duck tissues. The cloned cDNA of the muscovy duck FSHR gene shared high similarity to those of pekin duck (Anas platyrhynchos) (95.7%) and chicken (93.2%). Three different muscovy duck FSHR transcripts were identified. Quantitative real-time PCR (RT-qPCR) results showed that the FSHR gene was expressed in all the 14 tested tissues, and the highest expression level was seen in the ovary. A total of 16 SNPs were identified, among which, four SNPs were located in the coding region of FSHR. The SNP C320T is significantly associated with egg production at 59 weeks of age (P < 0.05), whereas the SNP A227G is significantly associated with age at first egg stage (P < 0.05). These results suggest that the two SNPs (A227G and C320T) of FSHR gene are associated with egg production traits and could be potential markers that can be used for marker-assisted selection programmes to increase egg production in muscovy duck.

Identification, expression and variation of the GNPDA2 gene, and its association with body weight and fatness traits in chicken.

Background. The GNPDA2 (glucosamine-6-phosphate deaminase 2) gene is a member of Glucosamine-6-phosphate (GlcN6P) deaminase subfamily, which encoded an allosteric enzyme of GlcN6P. Genome-wide association studies (GWAS) have shown that variations of human GNPDA2 are associated with body mass index and obesity risk, but its function and metabolic implications remain to be elucidated.The object of this study was to characterize the gene structure, expression, and biological functions of GNPDA2 in chickens. Methods. Variant transcripts of chicken GNPDA2 and their expression were investigated using rapid amplification of cDNA ends (RACE) system and real-time quantitative PCR technology. We detected the GNPDA2 expression in hypothalamic, adipose, and liver tissue of Xinghua chickens with fasting and high-glucose-fat diet treatments, and performed association analysis of variations of GNPDA2 with productive traits in chicken. The function of GNPDA2 was further studied by overexpression and small interfering RNA (siRNA) methods in chicken preadipocytes. Results.Four chicken GNPDA2 transcripts (cGNPDA2-a∼cGNPDA2-d) were identified in this study. The complete transcript GNPDA2-a was predominantly expressed in adipose tissue (subcutaneous fat and abdominal fat), hypothalamus, and duodenum. In fasting chickens, the mRNA level of GNPDA2 was decreased by 58.8% (P < 0.05) in hypothalamus, and returned to normal level after refeeding. Chicken fed a high-glucose-fat diet increased GNPDA2 gene expression about 2-fold higher in adipose tissue (P < 0.05) than that in the control (fed a basal diet), but decreased its expression in hypothalamus. Two single-nucleotide polymorphisms of the GNPDA2 gene were significantly associated with body weight and a number of fatness traits in chicken (P < 0.05). Conclusion. Our findings indicated that the GNPDA2 gene has a potential role in the regulation of body weight, fat and energy metabolism in chickens.

Effects of abhydrolase domain containing 5 gene (ABHD5) expression and variations on chicken fat metabolism.

Abhydrolase domain containing 5 gene (ABHD5), also known as comparative gene identification 58 (CGI-58), is a member of the α/ß-hydrolase family as a protein cofactor of ATGL stimulating its triacylglycerol hydrolase activity. In this study, we aim to characterize the expression and variations of ABHD5 and to study their functions in chicken fat metabolism. We compared the ABHD5 expression level in various tissues and under different nutrition conditions, identified the variations of ABHD5, and associated them with production traits in an F2 resource population of chickens. Overexpression analysis with two different genotypes and siRNA interfering analysis of ABHD5 were performed in chicken preadipocytes. Chicken ABDH5 was expressed widely and most predominantly in adipose tissue. Five SNPs of the ABHD5 gene were identified and genotyped in the F2 resource population. The c.490C > T SNP was associated with subcutaneous fat thickness (P < 0.01), carcass weight (P < 0.05), body weight (P < 0.05), shank diameter (P < 0.05), and shank length (P < 0.05). The c.423T > C SNP was also associated with chicken body weight (P < 0.05) and shank diameter (P < 0.05). In chicken preadipocytes, overexpression of wild type ABDH5 did not affect the mRNA level of ATGL (adipose triglyceride lipase) but markedly decreased (P < 0.05) the TG (triglyceride) content of the cell, whereas overexpression of mutation type ABHD5 did not affect either ATGL expression or TG content of the cell. The expression of ATGL and TG content of the cell were decreased (P < 0.05) after ABHD5 knockdown in preadipocytes. The mRNA level of ABHD5 was regulated by both feeding and fasting, and by consumption of a high fat diet. It was increased greatly by fasting (P < 0.05) and was returned to control levels after re-feeding in the adipose tissues, and down-regulated in abdominal fat (P < 0.05) and the liver (P < 0.01) of chickens with a high fat diet. These results suggest that expression and variations of ABHD5 may affect fat metabolism through regulating the activity of ATGL in chickens.