Integrated analysis of lncRNA-mediated ceRNA network involved in immune regulation in the spleen of Meishan piglets

Comprehensive expression genome-wide association study of long non-coding RNAs in four porcine tissues

BACKGROUND

Long non-coding RNAs (lncRNAs), a type of non-coding RNA molecules, are known to play critical regulatory roles in various biological processes. However, the functions of the majority of lncRNAs remain largely unknown, and little is understood about the regulation of lncRNA expression. In this study, high-throughput DNA genotyping and RNA sequencing were applied to investigate genomic regions associated with lncRNA expression, commonly referred to as lncRNA expression quantitative trait loci (eQTLs). We analyzed the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred sows to identify lncRNA transcripts, explore genomic regions that might influence lncRNA expression, and identify potential regulators interacting with these regions.

RESULT

We identified 6380 lncRNA transcripts and 3733 lncRNA genes. Correlation tests between the expression of lncRNAs and protein-coding genes were performed. Subsequently, functional enrichment analyses were carried out on protein-coding genes highly correlated with lncRNAs. Our correlation results of these protein-coding genes uncovered terms that are related to tissue specific functions. Additionally, heatmaps of lncRNAs and protein-coding genes at different correlation levels revealed several distinct clusters. An expression genome-wide association study (eGWAS) was conducted using 535,896 genotypes and 1829, 1944, 2089, and 2074 expressed lncRNA genes for liver, spleen, lung, and muscle, respectively. This analysis identified 520,562 significant associations and 6654, 4525, 4842, and 7125 eQTLs for the respective tissues. Only a small portion of these eQTLs were classified as cis-eQTLs. Fifteen regions with the highest eQTL density were selected as eGWAS hotspots and potential mechanisms of lncRNA regulation in these hotspots were explored. However, we did not identify any interactions between the transcription factors or miRNAs in the hotspots and the lncRNAs, nor did we observe a significant enrichment of regulatory elements in these hotspots. While we could not pinpoint the key factors regulating lncRNA expression, our results suggest that the regulation of lncRNAs involves more complex mechanisms.

CONCLUSION

Our findings provide insights into several features and potential functions of lncRNAs in various tissues. However, the mechanisms by which lncRNA eQTLs regulate lncRNA expression remain unclear. Further research is needed to explore the regulation of lncRNA expression and the mechanisms underlying lncRNA interactions with small molecules and regulatory proteins.

Regulation of cell-mediated immune responses in dairy bulls via long non-coding RNAs from submandibular lymph nodes, peripheral blood, and the spleen

Cell-mediated immune responses (CMIRs) are critical to building a robust immune system and reducing disease susceptibility in cattle. Long non-coding RNAs (lncRNAs) regulate various biological processes. However, to the best of our knowledge, the characterization and functions of lncRNAs and their regulations on the bovine CMIR have not been investigated comprehensively. In this study, experimental bulls were immunized with heat-killed preparation of Candida albicans (HKCA) to induce delayed-type hypersensitivity (DTH). Three bulls were classified as high- CMIR responders and three were low-CMIR responders, based on their classical DTH skin reactions. LncRNAs were identified in the submandibular lymph nodes, peripheral blood, and spleen of high- and low-CMIR animals using strand-specific RNA sequencing. A total of 21,003 putative lncRNAs were identified across tissues, and 420, 468, and 599 lncRNAs were differentially expressed between the two groups in the submandibular lymph node, peripheral blood, and spleen tissues, respectively. Functional analysis of the differentially expressed lncRNA (DElncRNA) target genes showed that a number of immune-related Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched, including immune response, cell adhesion, nucleosome, DNA packaging, antigen processing and presentation, and complement and coagulation cascades. Tissue specificity analysis indicated that lncRNA transcripts have stronger tissue specificity than mRNA. Furthermore, an interaction network was constructed based on DElncRNAs and DEGs, and 11, 14, and 11 promising lncRNAs were identified as potential candidate genes influencing immune response regulation in submandibular lymph nodes, peripheral blood, and spleen tissues, respectively. These results provide a foundation for further research into the biological functions of lncRNAs associated with bovine CMIR and identify candidate lncRNA markers for cell-mediated immune responses.

Whole transcriptome sequencing reveals key genes and ceRNA regulatory networks associated with pimpled eggs in hens

Eggshell is one of the most important indicators of egg quality, and due to low shell strength, pimple eggs (PE) are more susceptible to breakage, thus causing huge economic losses to the egg industry. At the current time, the molecular mechanisms that regulate the formation of pimple eggs are poorly understood. In this study, uterine tissues of PE-laying hens (n = 8) and normal egg (NE) -laying hens (n = 8) were analyzed by whole transcriptome sequencing, and a total of 619 differentially expressed mRNAs (DE mRNAs), 122 differentially expressed lncRNAs (DE lncRNAs) and 21 differentially expressed miRNAs (DE miRNAs) were obtained. Based on the targeting relationship among DE mRNAs, DE lncRNAs and DE miRNAs, we constructed a competitive endogenous RNA (ceRNA) network including 12 DE miRNAs, 19 DE lncRNAs, and 128 DE mRNAs. Considering the large amount of information contained in the network, we constructed a smaller ceRNA network to better understand the complex mechanisms of pimple egg formation. The smaller ceRNA network network contains 7 DE lncRNAs (LOC107056551, LOC121109367, LOC121108909, LOC121108862, LOC112530033, LOC121113165, LOC107054145), 5 DE miRNAs (gga-miR-6568-3p, gga-miR-31-5p, gga-miR-18b-3p, gga-miR-1759-3p, gga-miR-12240-3p) and 7 DE mRNAs (CABP1, DNAJC5, HCN3, HPCA, IBSP, KCNT1, OTOP3), and these differentially expressed genes may play key regulatory roles in the formation of pimpled eggs in hens. This study provides the overall expression profiles of mRNAs, lncRNAs and miRNAs in the uterine tissues of hens, which provides a theoretical basis for further research on the molecular mechanisms of pimpled egg formation, and has potential applications in improving eggshell quality.