Comprehensive Atlas of Alternative Splicing Reveals NSRP1 Promoting Adipogenesis through CCDC18.

Alternative splicing (AS) plays a crucial role in regulating gene expression, function, and diversity. However, limited reports exist on the identification and comparison of AS in Eastern and Western pigs. Here, we analyzed 243 transcriptome data from eight tissues, integrating information on transcription factors (TFs), selection signals, splicing factors (SFs), and quantitative trait loci (QTL) to comprehensively study alternative splicing events (ASEs) in pigs. Five ASE types were identified, with Mutually Exclusive Exon (MXE) and Skipped Exon (SE) ASEs being the most prevalent. A significant portion of genes with ASEs (ASGs) showed conservation across all eight tissues (63.21-76.13% per tissue). Differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) exhibited tissue specificity, with blood and adipose tissues having more DASGs. Functional enrichment analysis revealed coDASG_DEGs in adipose were enriched in pathways associated with adipose deposition and immune inflammation, while coDASG_DEGs in blood were enriched in pathways related to immune inflammation and metabolism. Adipose deposition in Eastern pigs might be linked to the down-regulation of immune-inflammation-related pathways and reduced insulin resistance. The TFs, selection signals, and SFs appeared to regulate ASEs. Notably, ARID4A (TF), NSRP1 (SF), ANKRD12, IFT74, KIAA2026, CCDC18, NEXN, PPIG, and ROCK1 genes in adipose tissue showed potential regulatory effects on adipose-deposition traits. NSRP1 could promote adipogenesis by regulating alternative splicing and expression of CCDC18. Conducting an in-depth investigation into AS, this study has successfully identified key marker genes essential for pig genetic breeding and the enhancement of meat quality, which will play important roles in promoting the diversity of pork quality and meeting market demand.

Transcriptional analyses provide novel insights into the transgenerational effects of Poly (I:C) on chickens.

Pathogenic microorganisms that are ubiquitous in the environment threaten human health and food safety. Polyinosinic:polycytidylic acid (Poly (I:C)) is a macromolecule with a double-stranded RNA structure, which is often used to simulate viruses. Our previous study found that Poly (I:C) maternal stimulation could affect the reproduction of laying hens and their offspring, but the underlying mechanism needed to be explored. In the present study, splenic transcriptomes were sequenced and analyzed from two groups (Poly (I:C) treatment as the challenged group and saline treatment as the control) and in three generations (maternal stimulated F0 hens, unchallenged F1 and F2 generations). The results showed that Poly (I:C) maternal stimulation affected gene expression patterns in laying hens and their offspring. A total of 27 differentially expressed genes (DEGs) with the same regulating trend were discovered in the F0 and F1 generations, indicating an influence of the intergenerational transmission effect. Functional enrichment analysis of Gene Ontology (GO) showed that lymphocyte differentiation, positive regulation of leukocyte differentiation, positive regulation of MAPK cascade, and T cell differentiation were the common biological processes between F0 and F1 generations, revealing Poly (I:C) could affect the immunity of the treated F0 hens and the unchallenged subsequent generations. Further study showed that pathways associated with growth, development, biosynthesis, and metabolism of F2 chicks were also affected by Poly (I:C) maternal stimulation. Correlation analysis between DEGs and reproductive traits revealed that PHLDA2 (pleckstrin homology-like domain family A member 2) and PODN (podocan) with inheritable effect were highly correlated with egg-laying rate and egg weight in F1 hens, suggesting their potential long-term role in regulating reproductive traits. ARHGAP40, FGB, HRH4, PHLDA2, PODN, NTSR1, and NMU were supposed to play important roles in regulating chickens' immunity and reproductive traits. This study reveals the far-reaching effect on transcriptome induced by Poly (I:C), reflecting the influence of the mother's living environment on the offspring. It is an important reference for future research into the multi-generational transmission of maternal stimulation and harmful environmental factors.

Research Note: Association of single nucleotide polymorphism of AKT3 with egg production traits in White Muscovy ducks (Cairina moschata).

Prior studies on transcriptomes of hypothalamus and ovary revealed that AKT3 is one of the candidate genes that might affect egg production in White Muscovy ducks. The role of AKT3 in the uterus during reproductive processes cannot be overemphasized. However, functional role of this gene in the tissues and on egg production traits of Muscovy ducks remains unknown. To identify the relationship between AKT3 and egg production traits in ducks, relative expression profile was first examined prior to identifying the variants within AKT3 that may underscore egg production traits [age at first egg (AFE), number of eggs at 300 d (N300D), and number of eggs at 59 wk (N59W)] in 549 ducks. The mRNA expression of AKT3 gene in high producing (HP) ducks was significantly higher than low producing (LP) ducks in the ovary, oviduct, and hypothalamus (P < 0.05 or 0.001). Three variants in AKT3 (C-3631A, C-3766T, and C-3953T) and high linkage block between C-3766T and C-3953T which are significantly (P < 0.05) associated with N300D and N59W were discovered. This study elucidates novel knowledge on the molecular mechanism of AKT3 that might be regulating egg production traits in Muscovy ducks.

MiR-743a-5p regulates differentiation of myoblast by targeting Mob1b in skeletal muscle development and regeneration.

The microRNAs (miRNAs) play an important role in regulating myogenesis by targeting mRNA. However, the understanding of miRNAs in skeletal muscle development and diseases is unclear. In this study, we firstly performed the transcriptome profiling in differentiating C2C12 myoblast cells. Totally, we identified 187 miRNAs and 4260 mRNAs significantly differentially expressed that were involved in myoblast differentiation. We carried out validation of microarray data based on 5 mRNAs and 5 miRNAs differentially expressed and got a consistent result. Then we constructed and validated the significantly up- and down-regulated mRNA-miRNA interaction networks. Four interaction pairs (miR-145a-5p-Fscn1, miR-200c-5p-Tmigd1, miR-27a-5p-Sln and miR-743a-5p-Mob1b) with targeted relationships in differentiated myoblast cells were demonstrated. They are all closely related to myoblast development. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated cell cycle signals important for exploring skeletal muscle development and disease. Functionally, we discovered that miR-743a targeting gene Mps One Binder Kinase Activator-Like 1B (Mob1b) gene in differentiated C2C12. The up-regulated miR-743a can promote the differentiation of C2C12 myoblast. While the down-regulated Mob1b plays a negative role in differentiation. In addition, the expression profile of miR-743a and Mob1b are consistent with skeletal muscle recovery after Cardiotoxin (CTX) injury. Our study revealed that miR-743a-5p regulates myoblast differentiation by targeting Mob1b involved in skeletal muscle development and regeneration. Our findings made a further exploration for mechanisms in myogenesis and might provide potential possible miRNA-based target therapies for skeletal muscle regeneration and disease in the near future.

Long-lasting effects of lipopolysaccharide on the reproduction and splenic transcriptome of hens and their offspring.

Lipopolysaccharide (LPS) is ubiquitous in the environment and is released after the death of gram-negative bacteria, which may be related to inflammation and immunosuppression. However, its impact on the reproduction of animals and their offspring, especially the underlying mechanism need further elucidation. Here, we used laying hens as a model organism to investigate the effects of maternal exposure to LPS (LPS maternal stimulation) on animal and their offspring's immunity and reproductive performance, as well as the regulatory role of the transcriptome. We found that the LPS maternal stimulation could reduce the egg-laying rate of hens and their offspring, especially during the early and late laying stages. The transcriptome study of the spleen in F0, F1 and F2 generations showed that the maternal stimulation of the LPS affects the patterns of gene expression in laying hens, and this change has a long-lasting effect. Further analysis of DEGs and their enrichment pathways found that the LPS maternal stimulation mainly affects the reproduction and immunity of laying hens and their offspring. The DEGs such as AVD, HPS5, CATHL2, S100A12, EXFABP, RSFR, LY86, PKD4, XCL1, FOS, TREM2 and MST1 may play an essential role in the regulation of the immunity and egg-laying rate of hens. Furthermore, the MMR1L3, C3, F13A1, LY86 and GDPD2 genes with heritable effects are highly correlated with the egg-laying rate, may have an important reference value for further research. Our study reveals the profound implications of LPS exposure on immunity and reproduction of offspring, elaborating the impact of immune alteration on the egg-laying rate, emphasizing the regulatory role of intergenerational transmission of the transcriptome, implying that the environment parents being exposed to has an important impact on offspring.

circRNAome profiling reveals circFgfr2 regulates myogenesis and muscle regeneration via a feedback loop.

BACKGROUND: Circular RNAs (circRNAs) represent a novel class of non-coding RNAs formed by a covalently closed loop and play crucial roles in many biological processes. Several circRNAs associated with myogenesis have been reported. However, the dynamic expression, function, and mechanism of circRNAs during myogenesis and skeletal muscle development are largely unknown. METHODS: Strand-specific RNA-sequencing (RNA-seq) and microarray datasets were used to profile the dynamic circRNAome landscape during skeletal muscle development and myogenic differentiation. Bioinformatics analyses were used to characterize the circRNAome and identify candidate circRNAs associated with myogenesis. Bulk and single-cell RNA-seq were performed to identify the downstream genes and pathways of circFgfr2. The primary myoblast cells, C2C12 cells, and animal model were used to assess the function and mechanism of circFgfr2 in myogenesis and muscle regeneration in vitro or in vivo by RT-qPCR, western blotting, dual-luciferase activity assay, RNA immunoprecipitation, RNA fluorescence in situ hybridization, and chromatin immunoprecipitation. RESULTS: We profiled the dynamic circRNAome in pig skeletal muscle across 27 developmental stages and detected 52 918 high-confidence circRNAs. A total of 2916 of these circRNAs are conserved across human, mouse, and pig, including four circRNAs (circFgfr2, circQrich1, circMettl9, and circCamta1) that were differentially expressed (|log2 fold change| > 1 and adjusted P value < 0.05) in various myogenesis systems. We further focused on a conserved circRNA produced from the fibroblast growth factor receptor 2 (Fgfr2) gene, termed circFgfr2, which was found to inhibit myoblast proliferation and promote differentiation and skeletal muscle regeneration. Mechanistically, circFgfr2 acted as a sponge for miR-133 to regulate the mitogen-activated protein kinase kinase kinase 20 (Map3k20) gene and JNK/MAPK pathway. Importantly, transcription factor Kruppel like factor 4 (Klf4), the downstream target of the JNK/MAPK pathway, directly bound to the promoter of circFgfr2 and affected its expression via an miR-133/Map3k20/JNK/Klf4 auto-regulatory feedback loop. RNA binding protein G3BP stress granule assembly factor 1 (G3bp1) inhibited the biogenesis of circFgfr2. CONCLUSIONS: The present study provides a comprehensive circRNA resource for skeletal muscle study. The functional and mechanistic analysis of circFgfr2 uncovered a circRNA-mediated auto-regulatory feedback loop regulating myogenesis and muscle regeneration, which provides new insight to further understand the regulatory mechanism of circRNAs.

A genome-wide single nucleotide polymorphism scan reveals genetic markers associated with fertility rate in Chinese Jing Hong chicken.

The function of the sperm storage tubules is directly correlated with the fertility of laying hens. However, little is known about the molecular mechanisms regulating the fertility traits in chicken. To identify genetic markers associated with reproductive traits, we calculated fertility rate at 61 to 69 wk (51 D) of Jing Hong chickens parent generation as the phenotype and the genotype were detected by the chicken 600K Affymetrix Axiom High Density single nucleotide polymorphisms (SNP)-array. The genome-wide association study using 190 Jing Hong hens showed that the 20 SNP in chromosomes 3 and 13 were significantly associated with fertility rate. To verify these results, a total of 1900 Jing Hong laying hens from 2 populations (P1 and P2) were further genotyped by polymerase chain reaction-restriction fragments length polymorphisms method. The association analysis results revealed that 12 polymorphisms (AX-75769978, AX-76582632, AX-75730546, AX-75730496, AX-75730588, AX-76530282, AX-76530329, AX-76529310, AX-75769906, AX-75755394, AX-80813697 and AX-76582809) out of 20 showed highly significant effects (P < 0.0001) on fertility rate in P1, P2 and P1+P2. Six haplotypes (TTAA, TTGG, TTAG, CTAA, CTGG, and CTAG) were inferred based on significant loci (AX-75730546 and AX-76530282) also showed significant association with fertility rate, where haplotype CTAG was shown to be markedly associated with the significantly highest (P < 0.0001) fertility rate (in P1, 86.42 ± 0.59; P2, 85.98 ± 0.59 and P1+P2, 86.16 ± 0.42) followed by other haplotypes for the irrespective of population studied. Collectively, we report for the first time that 12 SNP in the chromosomes 3 and 13 were significantly associated with fertility rate during the later stage of egg production, which could be used as the potential genetic markers that would be able to facilitate in the selection and improvement of fertility rate through chicken breeding.

RAI14 in the blood feather regulates chicken pigmentation.

A genome-wide association study (GWAS) was performed on a resource family consisting of white and colored chickens for identification of genes related to plumage coloration using the Fixed and random model Circulating Probability Unification (FarmCPU) package. GWAS identified three chromosomal single-nucleotide polymorphisms (SNPs), demonstrating the polygenic basis of plumage phenotypes. Herein, retinoic acid-induced protein 14 (RAI14), a developmentally regulated gene that encodes a protein containing many ankyrin repeats, was identified as a candidate gene involved in plumage color. In this study, mRNA expression profiles of chicken RAI14 were determined, indel (insertion-deletion) variants were identified, and their association was analyzed in white and colored chickens. RA114 mRNA was expressed in all tissues tested (brain, spleen, liver, heart, oviduct, kidney, lung, pituitary gland, ovary, muscle, feather bulb, and skin). A relatively high RAI14 expression in white feather bulb compared to colored feather bulb ( P < 0.01 ) indicated a potential association with plumage color. Additionally, statistical analysis revealed that a 4 bp indel genetic variation in RAI14 was associated with plumage phenotypes ( P < 0.01 ). Together, our analysis of the identification of the RAI14 gene will enable us to understand the genetic mechanisms behind chicken pigmentation.

Association of single nucleotide polymorphism at long non-coding RNA 8138.1 with duration of fertility in egg-laying hens.

A previous genome-wide transcriptional analysis identified long non-coding RNA 8138.1 (lncRNA8138.1) as a candidate gene related to hen duration of the fertility (DF) trait. LncRNA8138.1 gene response to growth factor and reproductive system development suggests it has a vital role in reproduction. In this study, we investigated the lncRNA8138.1 gene sequence in a population of egg-laying hens. The sequence analysis of the lncRNA8138.1 gene containing about 1.6 k nucleotides (nt) was observed with four single nucleotide polymorphisms (SNPs) and 7 nt indel including r.4937159A > G; r.4937219T > C; r.4937258G > C; r.4937318C > G and g.4937319_4937325delinsTGTGTGG. Next, the genomic DNAs from laying hen populations were subjected to polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) to detect a region of 457 bp carrying lncRNA8138.1 r.4937159A > G substitution. Further inspection of the region containing r.4937159A > G mutation revealed three genotypes viz., AA, AG, and GG were observed with respective frequencies of 0.106, 0.607, and 0.287 in laying hen population 1 (P1) (n = 1, 042) and respective frequencies of 0.176, 0.708, and 0.116 in laying hen population 2 (P2) (n = 826). Moreover, to further examining the frequencies of r.4937159A > G genotypes in P1 and P2, and their additive and dominance effects; r.4937159A > G locus was significantly associated with DF-trait in both P1 and P2 (EN: the number of eggs, FN: the number of fertile eggs after a single AI), and DN (the number of days post-insemination until last fertile egg). In testing for additive and dominance effects, additive effect was significant (P < 0.05) in both P1 and P2 for DF-trait, and the dominance effect was significant (P < 0.05) for EN and FN traits, suggesting that r.4937159A > G polymorphism is a potential biomarker for DF-trait. However, the identified novel r.4937159A > G mutation and others require further investigation to confirm phenotypic causality and potential genetic relationships with reproductive traits. Overall, our findings suggest the significance of genetic variation in long non-coding RNAs may assist in future breeding programs to improve selection for prolonged DF-trait.

Gene expression profiling of porcine skeletal muscle satellite cells after poly(I:C) stimulation.

Porcine satellite cells (PSCs) play a vital role in the construction, development and self-renewal of skeletal muscle. In this study, PSCs were exposed to poly(I:C) stimulation to mimic viral infection during the proliferation and differentiation phases at 0, 12, 24 and 48 hours (h) of the stimulation. The untreated and treated PSCs were analyzed by the RNA-Seq technology. There were 88, 119, 104 and 95 genes being differentially expressed in 0 h vs 12 h treated, 12 h vs 24 h treated, 0 h vs 24 h treated and 24 h vs 48 h untreated comparison libraries, respectively. The GO terms analysis results showed that during the proliferation phase of treated PSCs, the up-regulated genes related to the immune system were highly expressed. In addition, the gene expressions associated with muscle structure development in response to growth factor emerged during the differentiation phase of untreated PSCs. The biological pathways associated with Influenza A, Toll-like receptor and chemokine signaling were revealed in PSCs following poly(I:C) stimulation. The differentially expressed genes were confirmed by quantitative real-time PCR. These findings expanded our understanding of gene expressions and signaling pathways about the infiltrated mechanism of the virus into PSCs.

Transcriptional Profiling of Leucocyte Count Variation from Porcine Peripheral Blood Reveals Differential Gene Expression.

Leucocytes have tremendous health-check importance related to the individual antiviral capacity of pigs and other mammals. However, the molecular mechanism of the immune response of blood leucocytes in pigs is not completely known. This study investigated the leucocyte-count variation before and after poly I:C stimulation in a Duroc-Erhualian F2 population. Pigs with increased and decreased differences in leucocyte counts were coded as increased responder (IR) and decreased responder (DR), respectively. Then, we used microarray technology to compare the gene-expression profiles of both groups of pigs. Transcriptomic analysis identified 129 differentially expressed genes (DEGs) in IR pigs and 136 DEGs in DR pigs. Forty-one common DEGs showed that both groups had similar expression patterns of immune responses. These results illustrated a differential expression in both groups. Furthermore, qPCR experiment was performed to verify the differential-expression profile. Functional annotation of the DEGs indicated that both IR and DR pigs were similar in several biological processes, including innate immune response, and also exhibited distinct differences in biological processes, molecular function, and pathways. These results provided insights into the mechanism underlying the antiviral capacity of pigs. Trial registration number is CAS Registry Number 24939-03-5.

Transcriptome sequencing reveals key potential long non-coding RNAs related to duration of fertility trait in the uterovaginal junction of egg-laying hens.

Duration of fertility, (DF) is an important functional trait in poultry production and lncRNAs have emerged as important regulators of various process including fertility. In this study we applied a genome-guided strategy to reconstruct the uterovaginal junction (UVJ) transcriptome of 14 egg-laying birds with long- and short-DF (n = 7); and sought to uncover key lncRNAs related to duration of fertility traits by RNA-sequencing technology. Examination of RNA-seq data revealed a total of 9977 lncRNAs including 2576 novel lncRNAs. Differential expression (DE) analysis of lncRNA identified 223 lncRNAs differentially expressed between the two groups. DE-lncRNA target genes prediction uncovered over 200 lncRNA target genes and functional enrichment tests predict a potential function of DE-lncRNAs. Gene ontology classification and pathway analysis revealed 8 DE-lncRNAs, with the majority of their target genes enriched in biological functions such as reproductive structure development, developmental process involved in reproduction, response to cytokine, carbohydrate binding, chromatin organization, and immune pathways. Differential expression of lncRNAs and target genes were confirmed by qPCR. Together, these results significantly expand the utility of the UVJ transcriptome and our analysis identification of key lncRNAs and their target genes regulating DF will form the baseline for understanding the molecular functions of lncRNAs regulating DF.

Exploring the Genetic Resistance to Gastrointestinal Nematodes Infection in Goat Using RNA-Sequencing.

Gastrointestinal nematodes (GINs) are one of the most economically important parasites of small ruminants and a major animal health concern in many regions of the world. However, the molecular mechanisms of the host response to GIN infections in goat are still little known. In this study, two genetically distinct goat populations, one relatively resistant and the other susceptible to GIN infections, were identified in Yichang goat and then four individuals in each group were chosen to compare mRNA expression profiles using RNA-seq. Field experiment showed lower worm burden, delayed and reduced egg production in the relatively resistant group than the susceptible group. The analysis of RNA-seq showed that 2369 genes, 1407 of which were up-regulated and 962 down-regulated, were significantly (p < 0.001) differentially expressed between these two groups. Functional annotation of the 298 genes more highly expressed in the resistant group yielded a total of 46 significant (p < 0.05) functional annotation clusters including 31 genes (9 in innate immunity, 13 in immunity, and 9 in innate immune response) related to immune biosynthetic process as well as transforming growth factor (TGF)-ß, mitogen-activated protein kinase (MAPK), and cell adhesion molecules (CAMs) pathways. Our findings provide insights that are immediately relevant for the improvement of host resistance to GIN infections and which will make it possible to know the mechanisms underlying the resistance of goats to GIN infections.