Systematic Identification and Molecular Characteristics of Long Noncoding RNAs in Pig Tissues

T-2 toxin exposure induces ovarian damage in sows: lncRNA CUFF.253988.1 promotes cell apoptosis by inhibiting the SIRT3/PGC1α pathway

T-2 toxin, a mycotoxin found in foods and feeds, poses a threat to female reproductive health in both humans and animals. LncRNA CUFF.253988.1 (CUFF.253988.1), highly expressed in pigs, has an undisclosed regulatory role. This study aimed to establish a model of T-2 toxin-induced ovarian injury in sows, both in vivo and in vitro, and to explore the regulatory role and potential mechanisms of CUFF.253988.1. The results showed that feeding T-2 toxin-contaminated feed (1 mg/kg) induced ovarian follicle atresia and mitochondrial structural damage, accompanied by a significant upregulation of CUFF.253988.1 expression in the ovaries. Additionally, T-2 toxin inhibited the SIRT3/PGC1-α pathway associated with mitochondrial function. Moreover, T-2 toxin induced cell apoptosis by upregulating the expression of Cyt c, Bax, cleaved-caspase-9, and cleaved-caspase-3 proteins. In T-2 toxin-induced injury to the ovarian granulosa AVG-16 cells at concentrations of 10, 40 and 160 nM, not only were the previously mentioned effects observed, but also a decrease in mitochondrial membrane potential, ATP content, and an elevation in ROS levels. However, downregulating CUFF.253988.1 reversed T-2 toxin's inhibition of the SIRT3/PGC1-α pathway, alleviating mitochondrial dysfunction and reducing cell apoptosis. Notably, this may be attributed to the inhibition of T-2 toxin-induced enrichment of CUFF.253988.1 in mitochondria. In conclusion, CUFF.253988.1 plays a pivotal role in T-2 toxin-induced ovarian damage, operating through the inhibition of the SIRT3/PGC1-α pathway and promotion of cell apoptosis.

Whole genome sequencing of mouse lines divergently selected for fatness (FLI) and leanness (FHI) revealed several genetic variants as candidates for novel obesity genes

BACKGROUND

Analysing genomes of animal model organisms is widely used for understanding the genetic basis of complex traits and diseases, such as obesity, for which only a few mouse models exist, however, without their lean counterparts.

OBJECTIVE

To analyse genetic differences in the unique mouse models of polygenic obesity (Fat line) and leanness (Lean line) originating from the same base population and established by divergent selection over more than 60 generations.

METHODS

Genetic variability was analysed using WGS. Variants were identified with GATK and annotated with Ensembl VEP. g.Profiler, WebGestalt, and KEGG were used for GO and pathway enrichment analysis. miRNA seed regions were obtained with miRPathDB 2.0, LncRRIsearch was used to predict targets of identified lncRNAs, and genes influencing adipose tissue amount were searched using the IMPC database.

RESULTS

WGS analysis revealed 6.3 million SNPs, 1.3 million were new. Thousands of potentially impactful SNPs were identified, including within 24 genes related to adipose tissue amount. SNP density was highest in pseudogenes and regulatory RNAs. The Lean line carries SNP rs248726381 in the seed region of mmu-miR-3086-3p, which may affect fatty acid metabolism. KEGG analysis showed deleterious missense variants in immune response and diabetes genes, with food perception pathways being most enriched. Gene prioritisation considering SNP GERP scores, variant consequences, and allele comparison with other mouse lines identified seven novel obesity candidate genes: 4930441H08Rik, Aff3, Fam237b, Gm36633, Pced1a, Tecrl, and Zfp536.

CONCLUSION

WGS revealed many genetic differences between the lines that accumulated over the selection period, including variants with potential negative impacts on gene function. Given the increasing availability of mouse strains and genetic polymorphism catalogues, the study is a valuable resource for researchers to study obesity.

MALAT1: A Long Non-Coding RNA with Multiple Functions and Its Role in Processes Associated with Fat Deposition

Metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) belongs to the lncRNA molecules, which are involved in transcriptional and epigenetic regulation and the control of gene expression, including the mechanism of chromatin remodeling. MALAT1 was first discovered during carcinogenesis in lung adenocarcinoma, hence its name. In humans, 66 of its isoforms have been identified, and in pigs, only 2 are predicted, for which information is available in Ensembl databases (Ensembl Release 111). MALAT1 is expressed in numerous tissues, including adipose, adrenal gland, heart, kidney, liver, ovary, pancreas, sigmoid colon, small intestine, spleen, and testis. MALAT1, as an lncRNA, shows a wide range of functions. It is involved in the regulation of the cell cycle, where it has pro-proliferative effects and high cellular levels during the G1/S and mitotic (M) phases. Moreover, it is involved in invasion, metastasis, and angiogenesis, and it has a crucial function in alternative splicing during carcinogenesis. In addition, MALAT1 plays a significant role in the processes of fat deposition and adipogenesis. The human adipose tissue stem cells, during differentiation into adipocytes, secrete MALAT1 as one the most abundant lncRNAs in the exosomes. MALAT1 expression in fat tissue is positively correlated with adipogenic FABP4 and LPL. This lncRNA is involved in the regulation of PPARγ at the transcription stage, fatty acid metabolism, and insulin signaling. The wide range of MALAT1 functions makes it an interesting target in studies searching for drugs to prevent obesity development in humans. In turn, in farm animals, it can be a source of selection markers to control the fat tissue content.

Identification of a Novel Long Non-Coding RNA G8110 That Modulates Porcine Adipogenic Differentiation and Inflammatory Responses

Long non-coding RNAs (lncRNAs) have been extensively studied, and their crucial roles in adipogenesis, lipid metabolism, and gene expression have been revealed. However, the exact regulatory or other mechanisms by which lncRNAs influence the functioning of mesenteric adipose tissue (MAT) remain largely unknown. In this paper, we report the identification of a new lncRNA, named G8110, from the MAT of Bama pigs. The coordinated expression levels of lncRNA G8110 and NFE2L1 were significantly decreased in the MAT of obese Bama pigs compared with those in the MAT of lean pigs. Using a bone mesenchymal stem cell adipogenic differentiation model, we found that lncRNA G8110 played a role in adipocyte differentiation by positively regulating NFE2L1. We also found that lncRNA G8110 inhibited the formation of intracellular lipid synthesis, promoted lipid metabolism, and inhibited the expression of inflammatory cytokines. Our findings regarding lipid synthesis may further promote the role of lncRNAs in driving adipose tissue remodeling and maintaining metabolic health.

Genomic distribution and polymorphism of G-quadruplex motifs occupying ovine promoters and enhancers

G-quadruplexes are non-canonical DNA structures that are formed in regions with short runs of guanines. During the last decade they have gained considerable attention due to their involvement in basic cellular processes, linking them to several physiological processes and pathological conditions. Regulation of gene transcription is among the crucial roles that G-quadruplexes play in the cells. Several ways in which these structures affect transcription have been described, both negatively and positively. Recently, G-quadruplexes have been shown to be implicated in the three-dimensional rearrangement of the chromosomes that enables the interaction of enhancers and gene promoters during regulation of transcription. Sheep is a species for which almost no G-quadruplex-related studies have been conducted and thus research on this species is kept out from the progress that has been made in the G-quadruplex field. In this context, we investigated the DNA sequences with potential to form G-quadruplexes (G4-motifs) in the ovine enhancers and promoters. We describe the distribution of G4-motifs within the regulatory regions which is shown to be enriched in G4-motifs in a way similar to other mammals. Furthermore, our data suggest that G4-motifs promote promoter-enhancer interactions in sheep. The single nucleotide polymorphisms colocalizing with promoter- and enhancer-associated ovine G4-motifs constitute a considerable pool of polymorphism and given the crucial role of these specific G4-motifs on regulation of transcription, we suggest this polymorphism as an interesting target for ovine genetic studies.

Flnc: Machine Learning Improves the Identification of Novel Long Noncoding RNAs from Stand-Alone RNA-Seq Data

Long noncoding RNAs (lncRNAs) play critical regulatory roles in human development and disease. Although there are over 100,000 samples with available RNA sequencing (RNA-seq) data, many lncRNAs have yet to be annotated. The conventional approach to identifying novel lncRNAs from RNA-seq data is to find transcripts without coding potential but this approach has a false discovery rate of 30-75%. Other existing methods either identify only multi-exon lncRNAs, missing single-exon lncRNAs, or require transcriptional initiation profiling data (such as H3K4me3 ChIP-seq data), which is unavailable for many samples with RNA-seq data. Because of these limitations, current methods cannot accurately identify novel lncRNAs from existing RNA-seq data. To address this problem, we have developed software, Flnc, to accurately identify both novel and annotated full-length lncRNAs, including single-exon lncRNAs, directly from RNA-seq data without requiring transcriptional initiation profiles. Flnc integrates machine learning models built by incorporating four types of features: transcript length, promoter signature, multiple exons, and genomic location. Flnc achieves state-of-the-art prediction power with an AUROC score over 0.92. Flnc significantly improves the prediction accuracy from less than 50% using the conventional approach to over 85%. Flnc is available via GitHub platform.

Integrative analysis of transcriptome complexity in pig granulosa cells by long-read isoform sequencing

Background

In intensive and large-scale farms, abnormal estradiol levels in sows can cause reproductive disorders. The high incidence rate of reproductive disturbance will induce the elimination of productive sows in large quantities, and the poor management will bring great losses to the pig farms. The change in estradiol level has an important effect on follicular development and estrus of sows. To solve this practical problem and improve the productive capacity of sows, it is significant to further clarify the regulatory mechanism of estradiol synthesis in porcine granulosa cells (GCs). The most important function of granulosa cells is to synthesize estradiol. Thus, the studies about the complex transcriptome in porcine GCs are significant. As for precursor-messenger RNAs (pre-mRNAs), their post-transcriptional modification, such as alternative polyadenylation (APA) and alternative splicing (AS), together with long non-coding RNAs (lncRNAs), may regulate the functions of granulosa cells. However, the above modification events and their function are unclear within pig granulosa cells.

Methods

Combined PacBio long-read isoform sequencing (Iso-Seq) was conducted in this work for generating porcine granulosa cells' transcriptomic data. We discovered new transcripts and possible gene loci via comparison against reference genome. Later, combined Iso-Seq data were adopted to uncover those post-transcriptional modifications such as APA or AS, together with lncRNA within porcine granulosa cells. For confirming that the Iso-Seq data were reliable, we chose four AS genes and analyzed them through RT-PCR.

Results

The present article illustrated that pig GCs had a complex transcriptome, which gave rise to 8,793 APA, 3,465 AS events, 703 candidate new gene loci, as well as 92 lncRNAs. The results of this study revealed the complex transcriptome in pig GCs. It provided a basis for the interpretation of the molecular mechanism in GCs.

Transcriptome analysis reveals that long noncoding RNAs contribute to developmental differences between medium-sized ovarian follicles of Meishan and Duroc sows

Ovulation rate is an extremely important factor affecting litter size in sows. It differs greatly among pig breeds with different genetic backgrounds. Long non-coding RNAs (lncRNAs) can regulate follicle development, granulosa cell growth, and hormone secretion, which in turn can affect sow litter size. In this study, we identified 3554 lncRNAs and 25,491 mRNAs in M2 follicles of Meishan and Duroc sows. The lncRNA sequence and open reading frame lengths were shorter than mRNAs, and lncRNAs had fewer exons, were less abundant, and more conserved than protein-coding RNAs. Furthermore, 201 lncRNAs were differentially expressed (DE) between breeds, and quantitative trait loci analysis of DE lncRNAs were performed. A total of 127 DE lncRNAs were identified in 119 reproduction trait-related loci. In addition, the potential target genes of lncRNAs in cis or trans configurations were predicted. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that some potential target genes were involved in follicular development and hormone secretion-related biological processes or pathways, such as progesterone biosynthetic process, estrogen metabolic process, ovarian steroidogenesis, and PI3K-Akt signaling pathway. Furthermore, we also screened 19 differentially expressed lncRNAs in the PI3K-Akt signaling pathway as candidates. This study provides new insights into the roles of lncRNAs in follicular growth and development in pigs.

The genome variation and developmental transcriptome maps reveal genetic differentiation of skeletal muscle in pigs

Natural and artificial directional selections have resulted in significantly genetic and phenotypic differences across breeds in domestic animals. However, the molecular regulation of skeletal muscle diversity remains largely unknown. Here, we conducted transcriptome profiling of skeletal muscle across 27 time points, and performed whole-genome re-sequencing in Landrace (lean-type) and Tongcheng (obese-type) pigs. The transcription activity decreased with development, and the high-resolution transcriptome precisely captured the characterizations of skeletal muscle with distinct biological events in four developmental phases: Embryonic, Fetal, Neonatal, and Adult. A divergence in the developmental timing and asynchronous development between the two breeds was observed; Landrace showed a developmental lag and stronger abilities of myoblast proliferation and cell migration, whereas Tongcheng had higher ATP synthase activity in postnatal periods. The miR-24-3p driven network targeting insulin signaling pathway regulated glucose metabolism. Notably, integrated analysis suggested SATB2 and XLOC_036765 contributed to skeletal muscle diversity via regulating the myoblast migration and proliferation, respectively. Overall, our results provide insights into the molecular regulation of skeletal muscle development and diversity in mammals.

Analysis of mRNA and Long Non-Coding RNA Expression Profiles in Developing Yorkshire Pig Spleens

Simple Summary

Epidemic disease is a prominent problem in intensive pig production. The spleen is a blood bank and the largest immune organ, and most of the diseases in pig farms will be reflected as spleen abnormality. The results showed how the mRNA–lncRNA expression profiles in Yorkshire spleens varied with age (seven, 90, and 180 days after birth). Our study shows that 90 days after birth the gene-expression profile of pig spleen no longer changes significantly. The results are helpful for a better understanding of the transcriptome and functional genomics of spleen tissue in farm animals and could provide reference for precision pig disease research and prevention and control in pig farms.

Abstract

Epidemic diseases cause great economic loss in pig farms each year; some of these diseases are characterized mainly in the spleen, but mRNA and lncRNA (long non-coding RNA) expression networks in developing Yorkshire pig spleens remain obscure. Here, we profiled the systematic characters of mRNA and lncRNA repertoires in three groups of spleens from nine Yorkshire pigs, each three aged at seven days, 90 days, and 180 days. By using a precise mRNA and lncRNA identification pipeline, we identified 19,647 genes and 219 known and 3219 putative lncRNA transcripts; 1729 genes and 64 lncRNAs therein were found to express differentially. The gene expression characteristics of genes and lncRNAs were found to be basically fixed before 90 days after birth. Three large gene expression modules were detected. The enrichment analyses of differentially expressed genes and the potential target genes of differentially expressed lncRNAs both displayed the crucial roles of up-regulation in immune activation and hematopoiesis, and down-regulation in cell replication and division in 90 days and 180 days compared to seven days. ENSSSCT00000001325 was the only lncRNA transcript that existed in the three groups. CDK1, PCNA, and PLK were detected to be node genes that varied with age. This study contributes to a further understanding of mRNA and lncRNA expression in different developmental pig spleens.

Review: Long non-coding RNA in livestock

Less than 2% of mammalian genomes code for proteins, but 'the majority of its bases can be found in primary transcripts' - a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species - Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals' genomes is critical in bridging the gap between genotype and phenotype.

Dynamic Transcriptome Analysis Reveals Potential Long Non-coding RNAs Governing Postnatal Pineal Development in Pig

Postnatal development and maturation of pineal gland is a highly dynamic period of tissue remodeling and phenotype maintenance, which is genetically controlled by programmed gene expression regulations. However, limited molecular characterization, particularly regarding long noncoding RNAs (lncRNA), is available for postnatal pineal at a whole transcriptome level. The present study first characterized the comprehensive pineal transcriptome profiles using strand-specific RNA-seq to illustrate the dynamic mRNA/lncRNA expression at three developmental stages (infancy, puberty, and adulthood). The results showed that 21,448 mRNAs and 8,166 novel lncRNAs were expressed in pig postnatal pineal gland. Among these genes, 3,573 mRNAs and 851 lncRNAs, including the 5-hydroxytryptamine receptors, exhibited significant dynamic regulation along maturation process, while the expression of homeobox genes didn't show significant differences. Gene Ontology analysis revealed that the differentially expressed genes (DEGs) were significantly enriched in ion transport and synaptic transmission, highlighting the critical role of calcium signaling in postnatal pineal development. Additionally, co-expression analysis revealed the DEGs could be grouped into 12 clusters with distinct expression patterns. Many differential lncRNAs were functionally enriched in co-expressed clusters of genes related to ion transport, transcription regulation, DNA binding, and visual perception. Our study first provided an overview of postnatal pineal transcriptome dynamics in pig and demonstrated that dynamic lncRNA regulation of developmental transitions impact pineal physiology.

Transcriptome analysis to identify long non coding RNA (lncRNA) and characterize their functional role in back fat tissue of pig

Long non coding RNAs (lncRNA) have been previously found to be involved in important cellular activities like epigenetics, implantation, cell growth etc. in pigs. However, comprehensive analysis of lncRNA in back fat tissues at different developmental stages in pigs is still lacking. In this study we conducted transcriptome analysis in the back fat tissue of a F1 crossbred Korean Native Pig (KNP) × Yorkshire Pig to identify lncRNA. We investigated their role in 16 pigs at two different growth stages; stage 1 (10 weeks, n = 8) and stage 2 (26 weeks, n = 8). After quality assessment of sequencing reads, we got a total of 1,641,165 assembled transcripts out of eight paired end read from each stage. Among them, 6808 lncRNA transcripts were identified by filtering on the basis of multiple parameters like read length ≥ 200 nucleotides, exon numbers ≥2, FPKM ≥0.5, coding potential score < 0 etc. PFAM and RFAM were used to filter out all possible protein coding genes and housekeeping RNAs respectively. A total of 103 lncRNAs and 1057 mRNAs were found to be differentially expressed (DE) between the two stages (|log2FC| > 2, q < 0.05). We also identified 306 genes located around 100 kb upstream and 234 genes downstream around these DE lncRNA transcripts. The expression of top eleven DE lncRNAs (COL4A6, LY7S, MYH2, OXCT1, SMPDL3A, TMEM182, TTC36, RFOOOO4, RFOOO15, RFOOO45, CADM2) had been validating by qRT-PCR. Pathway and GO terms analysis showed that, positive regulation of biosynthetic process, Wnt signaling pathway, cellular protein modification process, and positive regulation of nitrogen compound were differentially enriched. Our results suggested that, KEGG pathways such as protein digestion and absorption, Arrhythmogenic right ventricular cardiomyopathy (ARVC) to be significantly enriched in both DE lncRNAs as well as DE mRNAs and involved in back fat tissues development. It also suggests that, identified lncRNAs are involved in regulation of important adipose tissues development pathways.

Effects of melatonin on the synthesis of estradiol and gene expression in pig granulosa cells

The interaction of granulosa cells (GCs) with oocytes is important to regulate follicle development. The exogenous melatonin promoting the maturation of oocytes by GCs has been approved in pig, however, the transcriptome profile and the functions of the genes regulated by melatonin in GCs have not yet to be fully characterized. In this study, we found melatonin could stimulate the synthesis of estradiol in pig GCs. The RNA-seq was used to explore the effects of melatonin on gene expression, a total of 89 differentially expressed genes (DEGs) were identified. Gene ontology analysis showed DEGs which associated with regulation of cell proliferation, cell cycle, and anti-apoptosis were significantly enriched. The functions of two DEGs, NOTCH2 and FILIP1L, were studied in pig GCs. The results showed that NOTCH2 inhibited the synthesis of estradiol, but FILIP1L promoted the synthesis of estradiol. Furthermore, inhibiting NOTCH2 in granulosa cells cocultured with cumulus-oocyte-complexes had no obvious effect on the maturation of pig oocyte, but could upregulate the cleavage rate of oocyte. We proved that FILIP1L had no effect on the maturation and cleavage of pig oocytes. Our work deepens the understanding of melatonin's effects on GCs and oocyte. The DEGs we found will be beneficial to reveal mechanisms of melatonin acting on GCs and oocytes and design the pharmacological interventions.

Genome-wide epigenetic landscape of pig lincRNAs and their evolution during porcine domestication

AIM

We aimed to identify previously unreported long intergenic noncoding RNAs (lincRNAs) in the porcine liver, an important metabolic tissue, and further illustrate the epigenomic landscapes and the evolution of lincRNAs.

MATERIALS & METHODS

We used porcine omics data and comprehensively analyzed and identified lincRNAs and their methylation, expression and evolutionary patterns during pig domestication.

RESULTS

LincRNAs exhibit highly methylated promoter and downstream regions, as well as lower expression levels and higher tissue specificity than protein-coding genes. We identified a batch of lincRNAs with selection signals that are associated with pig domestication, which are more highly expressed in the liver than in other tissues (19:10/8/6/3/2/1/1). Interestingly, the lincRNA linc-sscg1779 and its target gene C6, which is crucial in liver metabolism, are differentially expressed during pig domestication.

CONCLUSION

Although they may originate from noisy transcripts, lincRNAs may be subjected to artificial selection. This phenomenon implies the functional importance of lincRNAs in pig domestication.

Identification of long noncoding RNAs in Schistosoma mansoni and Schistosoma japonicum

Schistosomiasis is a major parasitic disease caused by 3 principal species of schistosome. Studies of schistosome transcriptomes have focused on protein-coding transcripts and although miRNAs are attracting increased attention, few reports have concerned the long noncoding RNAs (lncRNAs). These have been shown to play key roles in the regulation of gene expression through interactions with mRNAs, proteins and miRNAs. In this study, we first identified lncRNAs from RNA-seq data in Schistosoma mansoni and Schistosoma japonicum: 3247 and 3033 potential lncRNAs were found in these two species respectively. ChIP-seq analysis to determine H3K4me3 profiles along the gene regions corresponding to lncRNAs showed that in 12% of cases this mark was enriched in regions proximal to the transcription start sites, supporting their validity as actively transcribed genes. Besides, the sequence conservation of lncRNAs between schistosome species was much lower than that of mRNAs, but higher than that of the randomly selected genomic sequences, which is consistent with that in mammals. Our results demonstrate that lncRNAs form a significant part of the schistosome transcriptome and suggest that they play an important role in the biology of the parasite.

Identification and Functional Analysis of Long Intergenic Non-coding RNAs Underlying Intramuscular Fat Content in Pigs

Intramuscular fat (IMF) content is an important trait that can affect pork quality. Previous studies have identified many genes that can regulate IMF. Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in various biological processes. However, lincRNAs related to IMF in pig are largely unknown, and the mechanisms by which they regulate IMF are yet to be elucidated. Here we reconstructed 105,687 transcripts and identified 1,032 lincRNAs in pig longissimus dorsi muscle (LDM) of four stages with different IMF contents based on published RNA-seq. These lincRNAs show typical characteristics such as shorter length and lower expression compared with protein-coding genes. Combined with methylation data, we found that both the promoter and genebody methylation of lincRNAs can negatively regulate lincRNA expression. We found that lincRNAs exhibit high correlation with their protein-coding neighbors in expression. Co-expression network analysis resulted in eight stage-specific modules, gene ontology and pathway analysis of them suggested that some lincRNAs were involved in IMF-related processes, such as fatty acid metabolism and peroxisome proliferator-activated receptor signaling pathway. Furthermore, we identified hub lincRNAs and found six of them may play important roles in IMF development. This work detailed some lincRNAs which may affect of IMF development in pig, and facilitated future research on these lincRNAs and molecular assisted breeding for pig.