Genome-wide identification of tissue-specific long non-coding RNA in three farm animal species

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.

An mRNA expression atlas for the duck with public RNA-seq datasets

BACKGROUND

Ducks are globally important poultry species and a major source of farm animal products, including meat, eggs, and feathers. A thorough understanding of the functional genomic and transcriptomic sequences is crucial for improving production efficiency.

RESULT

This study constructed the largest duck mRNA expression atlas among all waterfowl species to date. The atlas encompasses 1,257 tissue samples across 30 tissue types, representing all major organ systems. Using advanced clustering analysis, we established co-expression network clusters to describe the transcriptional features in the duck mRNA expression atlas and, when feasible, assign these features to unique tissue types or pathways. Additionally, we identified 27 low-variance, highly expressed housekeeping genes suitable for gene expression experiments. Furthermore, in-depth analysis revealed potential sex-biased gene expression patterns within tissues and specific gene expression profiles in meat-type and egg-type ducks, providing valuable resources to understand the genetic basis of sex differences and particular phenotypes. This research elucidates the biological processes affecting duck productivity.

CONCLUSION

This study presents the most extensive gene expression atlas for any waterfowl species to date. These findings are of significant value for advancing duck biological research and industrial applications.

Unveiling the role of long non-coding RNAs in chicken immune response to highly pathogenic avian influenza H5N1 infection

Avian influenza viruses (AIVs) pose a significant threat to global poultry production, necessitating effective control strategies to mitigate economic losses and ensure animal welfare. Long non-coding RNAs (lncRNAs) have emerged as crucial regulators of immune responses, yet their roles in AIV-infected chickens remain poorly understood. This study aimed to investigate the expression profiles of lncRNAs and their targets in Vietnamese Ri chickens infected with the highly pathogenic AIV (HPAIV) H5N1. Through RNA sequencing, we identified novel lncRNAs and analyzed differentially expressed (DE) transcripts at 1 and 3 days post-infection (dpi) in chicken lung tissue. Our results revealed a higher number of DE lncRNAs and mRNAs at 1 dpi and 3 dpi, respectively, compared to control, with resistant chickens exhibiting a notably stronger immune response than susceptible chickens at 3 dpi. Functional analysis implicated these lncRNAs in immune-related pathways crucial for host responses to H5N1 viral infection. Furthermore, we identified lncRNA-mRNA interactions associated with antiviral responses and immune function. Notably, several genes involved in antiviral resistance and immune responses showed higher expression in resistant chickens, confirming their stronger antiviral response. Overall, our study provides insights into the role of lncRNAs in the host's response to HPAIV H5N1 infection in chickens and highlights potential candidates for further investigation into host-pathogen interactions. These findings could drive the development of novel control strategies for AIVs, significantly enhancing poultry health and biosecurity.

Transcriptomic Profiling Reveals Altered Expression of Genes Involved in Metabolic and Immune Processes in NDV-Infected Chicken Embryos

OBJECTIVE

The poultry industry is significantly impacted by viral infections, particularly Newcastle Disease Virus (NDV), which leads to substantial economic losses. It is essential to comprehend how the sequence of development affects biological pathways and how early exposure to infections might affect immune responses.

METHODS

This study employed transcriptome analysis to investigate host-pathogen interactions by analyzing gene expression changes in NDV-infected chicken embryos' lungs.

RESULT

RNA-Seq reads were aligned with the chicken reference genome (Galgal7), revealing 594 differentially expressed genes: 264 upregulated and 330 downregulated. The most overexpressed genes, with logFC between 8.15 and 8.75, included C8A, FGG, PIT54, FETUB, APOC3, and FGA. Notably, downregulated genes included BPIFB3 (-4.46 logFC) and TRIM39.1 (-4.26 logFC). The analysis also identified 29 novel transcripts and 20 lncRNAs that were upregulated. Gene Ontology and KEGG pathways' analyses revealed significant alterations in gene expression related to immune function, metabolism, cell cycle, nucleic acid processes, and mitochondrial activity due to NDV infection. Key metabolic genes, such as ALDOB (3.27 logFC), PRPS2 (2.66 logFC), and XDH (2.15 logFC), exhibited altered expression patterns, while DCK2 (-1.99 logFC) and TK1 (-2.11 logFC) were also affected. Several immune-related genes showed significant upregulation in infected lung samples, including ALB (6.15 logFC), TLR4 (1.86 logFC), TLR2 (2.79 logFC), and interleukin receptors, such as IL1R2 (3.15 logFC) and IL22RA2 (1.37 logFC). Conversely, genes such as CXCR4 (-1.49 logFC), CXCL14 (-2.57 logFC), GATA3 (-1.51 logFC), and IL17REL (-2.93 logFC) were downregulated. The higher expression of HSP genes underscores their vital role in immune responses.

CONCLUSION

Comprehension of these genes' interactions is essential for regulating viral replication and immune responses during infections, potentially aiding in the identification of candidate genes for poultry breed improvement amidst NDV challenges.

In search of epigenetic hallmarks of different tissues: an integrative omics study of horse liver, lung, and heart

DNA methylation and microRNA (miRNA) expression are epigenetic mechanisms essential for regulating tissue-specific gene expression and metabolic processes. However, high-resolution transcriptome, methylome, or miRNAome data is only available for a few model organisms and selected tissues. Up to date, only a few studies have reported on gene expression, DNA methylation, or miRNA expression in adult equine tissues at the genome-wide level. In the present study, we used RNA-Seq, miRNA-seq, and reduced representation bisulfite sequencing (RRBS) data from the heart, lung, and liver tissues of healthy cold-blooded horses to identify differentially expressed genes (DEGs), differentially expressed miRNA (DE miRNA) and differentially methylated sites (DMSs) between three types of horse tissues. Additionally, based on integrative omics analysis, we described the observed interactions of epigenetic mechanisms with tissue-specific gene expression alterations. The obtained data allowed identification from 4067 to 6143 DMSs, 9733 to 11,263 mRNAs, and 155 to 185 microRNAs, differentially expressed between various tissues. We pointed out specific genes whose expression level displayed a negative correlation with the level of CpG methylation and miRNA expression and revealed biological processes that they enrich. Furthermore, we confirmed and validated the accuracy of the Next-Generation Sequencing (NGS) results with bisulfite sequencing PCR (BSP) and quantitative PCR (qPCR). This comprehensive analysis forms a strong foundation for exploring the epigenetic mechanisms involved in tissue differentiation, especially the growth and development of the equine heart, lungs, and liver.

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.

Long non-coding RNA variability in porcine skeletal muscle

Recently, numerous studies including various tissues have been carried out on long non-coding RNAs (lncRNAs), but still, its variability has not yet been fully understood. In this study, we characterised the inter-individual variability of lncRNAs in pigs, in the context of number, length and expression. Transcriptomes collected from muscle tissue belonging to six Polish Landrace boars (PL1-PL6), including half-brothers (PL1-PL3), were investigated using bioinformatics (lncRNA identification and functional analysis) and statistical (lncRNA variability) methods. The number of lncRNA ranged from 1289 to 3500 per animal, and the total number of common lncRNAs among all boars was 232. The number, length and expression of lncRNAs significantly varied between individuals, and no consistent pattern has been found between pairs of half-brothers. In detail, PL5 exhibits lower expression than the others, while PL4 has significantly higher expression than PL2-PL3 and PL5-PL6. Noteworthy, comparing the inter-individual variability of lncRNA and mRNA expression, they exhibited concordant patterns. The enrichment analysis for common lncRNA target genes determined a variety of biological processes that play fundamental roles in cell biology, and they were mostly related to whole-body homeostasis maintenance, energy and protein synthesis as well as dynamics of multiple nucleoprotein complexes. The high variability of lncRNA landscape in the porcine genome has been revealed in this study. The inter-individual differences have been found in the context of three aspects: the number, length and expression of lncRNAs, which contribute to a better understanding of its complex nature.

GEGA (Gallus Enriched Gene Annotation): an online tool providing genomics and functional information across 47 tissues for a chicken gene-enriched atlas gathering Ensembl and Refseq genome annotations

GEGA is a user-friendly tool designed to navigate through various genomic and functional information related to an enriched gene atlas in chicken that integrates the gene catalogues from the two reference databases, NCBI-RefSeq and EMBL-Ensembl/GENCODE, along with four additional rich resources such as FAANG and NONCODE. Using the latest GRCg7b genome assembly, GEGA encompasses a total of 78 323 genes, including 24 102 protein-coding genes (PCGs) and 44 428 long non-coding RNAs (lncRNAs), significantly increasing the number of genes provided by each resource independently. However, GEGA is more than just a gene database. It offers a range of features that allow us to go deeper into the functional aspects of these genes. Users can explore gene expression and co-expression profiles across 47 tissues from 36 datasets and 1400 samples, discover tissue-specific variations and their expression as a function of sex or age and extract orthologous genes or their genomic configuration relative to the closest gene. For the communities interested in a specific gene, a list of genes or a quantitative trait locus region in chicken, GEGA's user-friendly interface facilitates efficient gene analysis, easy downloading of results and a multitude of graphical representations, from genomic information to detailed visualization of expression levels.

TransTEx: novel tissue-specificity scoring method for grouping human transcriptome into different expression groups

MOTIVATION

Although human tissues carry out common molecular processes, gene expression patterns can distinguish different tissues. Traditional informatics methods, primarily at the gene level, overlook the complexity of alternative transcript variants and protein isoforms produced by most genes, changes in which are linked to disease prognosis and drug resistance.

RESULTS

We developed TransTEx (Transcript-level Tissue Expression), a novel tissue-specificity scoring method, for grouping transcripts into four expression groups. TransTEx applies sequential cut-offs to tissue-wise transcript probability estimates, subsampling-based P-values and fold-change estimates. Application of TransTEx on GTEx mRNA-seq data divided 199 166 human transcripts into different groups as 17 999 tissue-specific (TSp), 7436 tissue-enhanced, 36 783 widely expressed (Wide), 79 191 lowly expressed (Low), and 57 757 no expression (Null) transcripts. Testis has the most (13 466) TSp isoforms followed by liver (890), brain (701), pituitary (435), and muscle (420). We found that the tissue specificity of alternative transcripts of a gene is predominantly influenced by alternate promoter usage. By overlapping brain-specific transcripts with the cell-type gene-markers in scBrainMap database, we found that 63% of the brain-specific transcripts were enriched in nonneuronal cell types, predominantly astrocytes followed by endothelial cells and oligodendrocytes. In addition, we found 61 brain cell-type marker genes encoding a total of 176 alternative transcripts as brain-specific and 22 alternative transcripts as testis-specific, highlighting the complex TSp and cell-type specific gene regulation and expression at isoform-level. TransTEx can be adopted to the analysis of bulk RNA-seq or scRNA-seq datasets to find tissue- and/or cell-type specific isoform-level gene markers.

AVAILABILITY AND IMPLEMENTATION

TransTEx database: https://bmi.cewit.stonybrook.edu/transtexdb/ and the R package is available via GitHub: https://github.com/pallavisurana1/TransTEx.

A Synergistic, Cultivator Model of De Novo Gene Origination

The origin and fixation of evolutionarily young genes is a fundamental question in evolutionary biology. However, understanding the origins of newly evolved genes arising de novo from noncoding genomic sequences is challenging. This is partly due to the low likelihood that several neutral or nearly neutral mutations fix prior to the appearance of an important novel molecular function. This issue is particularly exacerbated in large effective population sizes where the effect of drift is small. To address this problem, we propose a regulation-focused, cultivator model for de novo gene evolution. This cultivator-focused model posits that each step in a novel variant's evolutionary trajectory is driven by well-defined, selectively advantageous functions for the cultivator genes, rather than solely by the de novo genes, emphasizing the critical role of genome organization in the evolution of new genes.

A first complete catalog of highly expressed genes in eight chicken tissues reveals uncharacterized gene families specific for the chicken testis

Based on next-generation sequencing, we established a repertoire of differentially overexpressed genes (DoEGs) in eight adult chicken tissues: the testis, brain, lung, liver, kidney, muscle, heart, and intestine. With 4,499 DoEGs, the testis had the highest number and proportion of DoEGs compared with the seven somatic tissues. The testis DoEG set included the highest proportion of long noncoding RNAs (lncRNAs; 1,851, representing 32% of the lncRNA genes in the whole genome) and the highest proportion of protein-coding genes (2,648, representing 14.7% of the protein-coding genes in the whole genome). The main significantly enriched Gene Ontology terms related to the protein-coding genes were "reproductive process," "tubulin binding," and "microtubule cytoskeleton." Using real-time quantitative reverse transcription-polymerase chain reaction, we confirmed the overexpression of genes that encode proteins already described in chicken sperm [such as calcium binding tyrosine phosphorylation regulated (CABYR), spermatogenesis associated 18 (SPATA18), and CDK5 regulatory subunit associated protein (CDK5RAP2)] but whose testis origin had not been previously confirmed. Moreover, we demonstrated the overexpression of vertebrate orthologs of testis genes not yet described in the adult chicken testis [such as NIMA related kinase 2 (NEK2), adenylate kinase 7 (AK7), and CCNE2]. Using clustering according to primary sequence homology, we found that 1,737 of the 2,648 (67%) testis protein-coding genes were unique genes. This proportion was significantly higher than the somatic tissues except muscle. We clustered the other 911 testis protein-coding genes into 495 families, from which 47 had all paralogs overexpressed in the testis. Among these 47 testis-specific families, eight contained uncharacterized duplicated paralogs without orthologs in other metazoans except birds: these families are thus specific for chickens/birds.NEW & NOTEWORTHY Comparative next-generation sequencing analysis of eight chicken tissues showed that the testis has highest proportion of long noncoding RNA and protein-coding genes of the whole genome. We identified new genes in the chicken testis, including orthologs of known mammalian testicular genes. We also identified 47 gene families in which all the members were overexpressed, if not exclusive, in the testis. Eight families, organized in duplication clusters, were unknown, without orthologs in metazoans except birds, and are thus specific for chickens/birds.

Plant-LncPipe: a computational pipeline providing significant improvement in plant lncRNA identification

Long non-coding RNAs (lncRNAs) play essential roles in various biological processes, such as chromatin remodeling, post-transcriptional regulation, and epigenetic modifications. Despite their critical functions in regulating plant growth, root development, and seed dormancy, the identification of plant lncRNAs remains a challenge due to the scarcity of specific and extensively tested identification methods. Most mainstream machine learning-based methods used for plant lncRNA identification were initially developed using human or other animal datasets, and their accuracy and effectiveness in predicting plant lncRNAs have not been fully evaluated or exploited. To overcome this limitation, we retrained several models, including CPAT, PLEK, and LncFinder, using plant datasets and compared their performance with mainstream lncRNA prediction tools such as CPC2, CNCI, RNAplonc, and LncADeep. Retraining these models significantly improved their performance, and two of the retrained models, LncFinder-plant and CPAT-plant, alongside their ensemble, emerged as the most suitable tools for plant lncRNA identification. This underscores the importance of model retraining in tackling the challenges associated with plant lncRNA identification. Finally, we developed a pipeline (Plant-LncPipe) that incorporates an ensemble of the two best-performing models and covers the entire data analysis process, including reads mapping, transcript assembly, lncRNA identification, classification, and origin, for the efficient identification of lncRNAs in plants. The pipeline, Plant-LncPipe, is available at: https://github.com/xuechantian/Plant-LncRNA-pipline.

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

Deep Conservation and Unexpected Evolutionary History of Neighboring lncRNAs MALAT1 and NEAT1

Long non-coding RNAs (lncRNAs) have begun to receive overdue attention for their regulatory roles in gene expression and other cellular processes. Although most lncRNAs are lowly expressed and tissue-specific, notable exceptions include MALAT1 and its genomic neighbor NEAT1, two highly and ubiquitously expressed oncogenes with roles in transcriptional regulation and RNA splicing. Previous studies have suggested that NEAT1 is found only in mammals, while MALAT1 is present in all gnathostomes (jawed vertebrates) except birds. Here we show that these assertions are incomplete, likely due to the challenges associated with properly identifying these two lncRNAs. Using phylogenetic analysis and structure-aware annotation of publicly available genomic and RNA-seq coverage data, we show that NEAT1 is a common feature of tetrapod genomes except birds and squamates. Conversely, we identify MALAT1 in representative species of all major gnathostome clades, including birds. Our in-depth examination of MALAT1, NEAT1, and their genomic context in a wide range of vertebrate species allows us to reconstruct the series of events that led to the formation of the locus containing these genes in taxa from cartilaginous fish to mammals. This evolutionary history includes the independent loss of NEAT1 in birds and squamates, since NEAT1 is found in the closest living relatives of both clades (crocodilians and tuataras, respectively). These data clarify the origins and relationships of MALAT1 and NEAT1 and highlight an opportunity to study the change and continuity in lncRNA structure and function over deep evolutionary time.

Long non-coding RNAs mediate fish gene expression in response to ocean acidification

The majority of the transcribed genome does not have coding potential but these non-coding transcripts play crucial roles in transcriptional and post-transcriptional regulation of protein-coding genes. Regulation of gene expression is important in shaping an organism's response to environmental changes, ultimately impacting their survival and persistence as population or species face global change. However, the roles of long non-coding RNAs (lncRNAs), when confronted with environmental changes, remain largely unclear. To explore the potential role of lncRNAs in fish exposed to ocean acidification (OA), we analyzed publicly available brain RNA-seq data from a coral reef fish Acanthochromis polyacanthus. We annotated the lncRNAs in its genome and examined the expression changes of intergenic lncRNAs (lincRNAs) between A. polyacanthus samples from a natural CO2 seep and a nearby control site. We identified 4728 lncRNAs, including 3272 lincRNAs in this species. Remarkably, 93.03% of these lincRNAs were species-specific. Among the 125 highly expressed lincRNAs and 403 differentially expressed lincRNAs in response to elevated CO2, we observed that lincRNAs were either neighboring or potentially trans-regulating differentially expressed coding genes associated with pH regulation, neural signal transduction, and ion transport, which are known to be important in the response to OA in fish. In summary, lncRNAs may facilitate fish acclimation and mediate the responses of fish to OA by modulating the expression of crucial coding genes, which offers insight into the regulatory mechanisms underlying fish responses to environmental changes.

Combined analysis of lncRNA and mRNA emphasizes the potential role of tryptophan-mediated regulation of muscle development in weaned piglets by lncRNA

Pork is an important high-value protein source that fulfills the nutritional requirements for normal growth development, repair, and metabolism. Tryptophan (Trp), a crucial amino acid for piglet growth performance and muscle development, has an essential yet unclear regulatory mechanism. To investigate the biological basis of Trp regulation of piglet muscle development and identify the related regulatory pathways, we studied 20 weaned piglets. The piglets were divided into control (CON, 0.14% Trp) and high Trp (HT, 0.35% Trp) groups. They were fed with different Trp concentrations for 28 d, after which we collected the longissimus dorsi (LD) muscle for histomorphometric analysis and RNA extraction. Our results showed that the HT diet significantly increased the average daily weight gain, myocyte number, and muscle fiber density in weaned piglets. We then analyzed the differentially expressed (DE) genes in the LD muscle through RNA sequencing (RNA-seq). We identified 253 lncRNAs and 1,055 mRNAs mainly involved in myoblast proliferation and myofiber formation, particularly through the FoxO and AMPK signaling pathways and metabolism. Further analysis of the DE lncRNA targeting relationship and construction of a protein-protein interaction network resulted in the discovery of a novel lncRNA, XLOC_021675, or FRPMD, and elucidated its role in regulating piglet muscle development. Finally, we confirmed the RNA-seq results by reverse transcription polymerase chain reaction (RT-PCR). This study provides valuable insights into the regulatory mechanism of lncRNA-mediated Trp regulation of muscle development in weaned piglets offering a theoretical basis for optimizing piglet dietary ratios and enhancing pork production.

Transcriptome analysis of ovarian tissues highlights genes controlling energy homeostasis and oxidative stress as potential drivers of heterosis for egg number and clutch size in crossbred laying hens

Heterosis is the major benefit of crossbreeding and has been exploited in laying hens breeding for a long time. This genetic phenomenon has been linked to various modes of nonadditive gene action. However, the molecular mechanism of heterosis for egg production in laying hens has not been fully elucidated. To fill this research gap, we sequenced mRNAs and lncRNAs of the ovary stroma containing prehierarchical follicles in White Leghorn, Rhode Island Red chickens as well as their reciprocal crossbreds that demonstrated heterosis for egg number and clutch size. We further delineated the modes of mRNAs and lncRNAs expression to identify their potential functions in the observed heterosis. Results showed that dominance was the principal mode of nonadditive expression exhibited by mRNAs and lncRNAs in the prehierarchical follicles of crossbred hens. Specifically, low-parent dominance was the main mode of mRNA expression, while high-parent dominance was the predominant mode of lncRNA expression. Important pathways enriched by genes that showed higher expression in crossbreds compared to either one or both parental lines were cell adhesion molecules, tyrosine and purine metabolism. In contrast, ECM-receptor interaction, focal adhesion, PPAR signaling, and ferroptosis were enriched in genes with lower expression in the crossbred. Protein network interaction identified nonadditively expressed genes including apolipoprotein B (APOB), transferrin, acyl-CoA synthetase medium-chain family member (APOBEC) 3, APOBEC1 complementation factor, and cathepsin S as hub genes. Among these potential hub genes, APOB was the only gene with underdominance expression common to the 2 reciprocal crossbred lines, and has been linked to oxidative stress. LncRNAs with nonadditive expression in the crossbred hens targeted natriuretic peptide receptor 1, epidermal differentiation protein beta, spermatogenesis-associated gene 22, sperm-associated antigen 16, melanocortin 2 receptor, dolichol kinase, glycine amiinotransferase, and prolactin releasing hormone receptor. In conclusion, genes with nonadditive expression in the crossbred may play crucial roles in follicle growth and atresia by improving follicle competence and increasing oxidative stress, respectively. These 2 phenomena could underpin heterosis for egg production in crossbred laying hens.

Investigation of lncRNA in Bos taurus Mammary Tissue during Dry and Lactation Periods

This study aims to collect RNA-Seq data from Bos taurus samples representing dry and lactating mammary tissue, identify lncRNA transcripts, and analyze findings for their features and functional annotation. This allows for connections to be drawn between lncRNA and the lactation process. RNA-Seq data from 103 samples of Bos taurus mammary tissue were gathered from publicly available databases (60 dry, 43 lactating). The samples were filtered to reveal 214 dry mammary lncRNA transcripts and 517 lactating mammary lncRNA transcripts. The lncRNAs met common lncRNA characteristics such as shorter length, fewer exons, lower expression levels, and less sequence conservation when compared to the genome. Interestingly, several lncRNAs showed sequence similarity to genes associated with strong hair keratin intermediate filaments. Human breast cancer research has associated strong hair keratin filaments with mammary tissue cellular resilience. The lncRNAs were also associated with several genes/proteins that linked to pregnancy using expression correlation and gene ontology. Such findings indicate that there are crucial relationships between the lncRNAs found in mammary tissue and the development of the tissue, to meet both the animal's needs and our own production needs; these relationships should be further investigated to ensure that we continue to breed the most resilient, efficient dairy cattle.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Long Non-Coding RNAs and Their "Discrete" Contribution to IBD and Johne's Disease-What Stands out in the Current Picture? A Comprehensive Review

Non-coding RNAs (ncRNA) have paved the way to new perspectives on the regulation of gene expression, not only in biology and medicine, but also in associated fields and technologies, ensuring advances in diagnostic means and therapeutic modalities. Critical in this multistep approach are the associations of long non-coding RNA (lncRNA) with diseases and their causal genes in their networks of interactions, gene enrichment and expression analysis, associated pathways, the monitoring of the involved genes and their functional roles during disease progression from one stage to another. Studies have shown that Johne's Disease (JD), caused by Mycobacterium avium subspecies partuberculosis (MAP), shares common lncRNAs, clinical findings, and other molecular entities with Crohn's Disease (CD). This has been a subject of vigorous investigation owing to the zoonotic nature of this condition, although results are still inconclusive. In this review, on one hand, the current knowledge of lncRNAs in cells is presented, focusing on the pathogenesis of gastrointestinal-related pathologies and MAP-related infections and, on the other hand, we attempt to dissect the associated genes and pathways involved. Furthermore, the recently characterized and novel lncRNAs share common pathologies with IBD and JD, including the expression, molecular networks, and dataset analysis results. These are also presented in an attempt to identify potential biomarkers pertinent to cattle and human disease phenotypes.

LncDLSM: Identification of Long Non-Coding RNAs With Deep Learning-Based Sequence Model

Long non-coding RNAs (LncRNAs) serve a vital role in regulating gene expressions and other biological processes. Differentiation of lncRNAs from protein-coding transcripts helps researchers dig into the mechanism of lncRNA formation and its downstream regulations related to various diseases. Previous works have been proposed to identify lncRNAs, including traditional bio-sequencing and machine learning approaches. Considering the tedious work of biological characteristic-based feature extraction procedures and inevitable artifacts during bio-sequencing processes, those lncRNA detection methods are not always satisfactory. Hence, in this work, we presented lncDLSM, a deep learning-based framework differentiating lncRNA from other protein-coding transcripts without dependencies on prior biological knowledge. lncDLSM is a helpful tool for identifying lncRNAs compared with other biological feature-based machine learning methods and can be applied to other species by transfer learning achieving satisfactory results. Further experiments showed that different species display distinct boundaries among distributions corresponding to the homology and the specificity among species, respectively.

Transcriptome annotation of 17 porcine tissues using nanopore sequencing technology

The annotation of animal genomes plays an important role in elucidating molecular mechanisms behind the genetic control of economically important traits. Here, we employed long-read sequencing technology, Oxford Nanopore Technology, to annotate the pig transcriptome across 17 tissues from two Yorkshire littermate pigs. More than 9.8 million reads were obtained from a single flow cell, and 69 781 unique transcripts at 50 108 loci were identified. Of these transcripts, 16 255 were found to be novel isoforms, and 22 344 were found at loci that were novel and unannotated in the Ensembl (release 102) and NCBI (release 106) annotations. Novel transcripts were mostly expressed in cerebellum, followed by lung, liver, spleen, and hypothalamus. By comparing the unannotated transcripts to existing databases, there were 21 285 (95.3%) transcripts matched to the NT database (v5) and 13 676 (61.2%) matched to the NR database (v5). Moreover, there were 4324 (19.4%) transcripts matched to the SwissProt database (v5), corresponding to 11 356 proteins. Tissue-specific gene expression analyses showed that 9749 transcripts were highly tissue-specific, and cerebellum contained the most tissue-specific transcripts. As the same samples were used for the annotation of cis-regulatory elements in the pig genome, the transcriptome annotation generated by this study provides an additional and complementary annotation resource for the Functional Annotation of Animal Genomes effort to comprehensively annotate the pig genome.

The potential regulatory role of the lncRNA-miRNA-mRNA axis in teleost fish

Research over the past two decades has confirmed that noncoding RNAs (ncRNAs), which are abundant in cells from yeast to vertebrates, are no longer "junk" transcripts but functional regulators that can mediate various cellular and physiological processes. The dysregulation of ncRNAs is closely related to the imbalance of cellular homeostasis and the occurrence and development of various diseases. In mammals, ncRNAs, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), have been shown to serve as biomarkers and intervention targets in growth, development, immunity, and disease progression. The regulatory functions of lncRNAs on gene expression are usually mediated by crosstalk with miRNAs. The most predominant mode of lncRNA-miRNA crosstalk is the lncRNA-miRNA-mRNA axis, in which lncRNAs act as competing endogenous RNAs (ceRNAs). Compared to mammals, little attention has been given to the role and mechanism of the lncRNA-miRNA-mRNA axis in teleost species. In this review, we provide current knowledge about the teleost lncRNA-miRNA-mRNA axis, focusing on its physiological and pathological regulation in growth and development, reproduction, skeletal muscle, immunity to bacterial and viral infections, and other stress-related immune responses. Herein, we also explored the potential application of the lncRNA-miRNA-mRNA axis in the aquaculture industry. These findings contribute to an enhanced understanding of ncRNA and ncRNA-ncRNA crosstalk in fish biology to improve aquaculture productivity, fish health and quality.

Investigation of rumen long noncoding RNA before and after weaning in cattle

Background

This study aimed to identify long non-coding RNA (lncRNA) from the rumen tissue in dairy cattle, explore their features including expression and conservation levels, and reveal potential links between lncRNA and complex traits that may indicate important functional impacts of rumen lncRNA during the transition to the weaning period.

Results

A total of six cattle rumen samples were taken with three replicates from before and after weaning periods, respectively. Total RNAs were extracted and sequenced with lncRNA discovered based on size, coding potential, sequence homology, and known protein domains. As a result, 404 and 234 rumen lncRNAs were identified before and after weaning, respectively. However, only nine of them were shared under two conditions, with 395 lncRNAs found only in pre-weaning tissues and 225 only in post-weaning samples. Interestingly, none of the nine common lncRNAs were differentially expressed between the two weaning conditions. LncRNA averaged shorter length, lower expression, and lower conservation scores than the genome overall, which is consistent with general lncRNA characteristics. By integrating rumen lncRNA before and after weaning with large-scale GWAS results in cattle, we reported significant enrichment of both pre- and after-weaning lncRNA with traits of economic importance including production, reproduction, health, and body conformation phenotypes.

Conclusions

The majority of rumen lncRNAs are uniquely expressed in one of the two weaning conditions, indicating a functional role of lncRNA in rumen development and transition of weaning. Notably, both pre- and post-weaning lncRNA showed significant enrichment with a variety of complex traits in dairy cattle, suggesting the importance of rumen lncRNA for cattle performance in the adult stage. These relationships should be further investigated to better understand the specific roles lncRNAs are playing in rumen development and cow performance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08758-4.

Transcriptomic meta-analysis reveals unannotated long non-coding RNAs related to the immune response in sheep

Long non-coding RNAs (lncRNAs) are involved in several biological processes, including the immune system response to pathogens and vaccines. The annotation and functional characterization of lncRNAs is more advanced in humans than in livestock species. Here, we take advantage of the increasing number of high-throughput functional experiments deposited in public databases in order to uniformly analyse, profile unannotated lncRNAs and integrate 422 ovine RNA-seq samples from the ovine immune system. We identified 12302 unannotated lncRNA genes with support from independent CAGE-seq and histone modification ChIP-seq assays. Unannotated lncRNAs showed low expression levels and sequence conservation across other mammal species. There were differences in expression levels depending on the genomic location-based lncRNA classification. Differential expression analyses between unstimulated and samples stimulated with pathogen infection or vaccination resulted in hundreds of lncRNAs with changed expression. Gene co-expression analyses revealed immune gene-enriched clusters associated with immune system activation and related to interferon signalling, antiviral response or endoplasmic reticulum stress. Besides, differential co-expression networks were constructed in order to find condition-specific relationships between coding genes and lncRNAs. Overall, using a diverse set of immune system samples and bioinformatic approaches we identify several ovine lncRNAs associated with the response to an external stimulus. These findings help in the improvement of the ovine lncRNA catalogue and provide sheep-specific evidence for the implication in the general immune response for several lncRNAs.

Prediction of transcript isoforms in 19 chicken tissues by Oxford Nanopore long-read sequencing

To identify and annotate transcript isoforms in the chicken genome, we generated Nanopore long-read sequencing data from 68 samples that encompassed 19 diverse tissues collected from experimental adult male and female White Leghorn chickens. More than 23.8 million reads with mean read length of 790 bases and average quality of 18.2 were generated. The annotation and subsequent filtering resulted in the identification of 55,382 transcripts at 40,547 loci with mean length of 1,700 bases. We predicted 30,967 coding transcripts at 19,461 loci, and 16,495 lncRNA transcripts at 15,512 loci. Compared to existing reference annotations, we found ∼52% of annotated transcripts could be partially or fully matched while ∼47% were novel. Seventy percent of novel transcripts were potentially transcribed from lncRNA loci. Based on our annotation, we quantified transcript expression across tissues and found two brain tissues (i.e., cerebellum and cortex) expressed the highest number of transcripts and loci. Furthermore, ∼22% of the transcripts displayed tissue specificity with the reproductive tissues (i.e., testis and ovary) exhibiting the most tissue-specific transcripts. Despite our wide sampling, ∼20% of Ensembl reference loci were not detected. This suggests that deeper sequencing and additional samples that include different breeds, cell types, developmental stages, and physiological conditions, are needed to fully annotate the chicken genome. The application of Nanopore sequencing in this study demonstrates the usefulness of long-read data in discovering additional novel loci (e.g., lncRNA loci) and resolving complex transcripts (e.g., the longest transcript for the TTN locus).

LncRNA affects epigenetic reprogramming of porcine embryo development by regulating global epigenetic modification and the downstream gene SIN3A

The study of preimplantation development is of great significance to reproductive biology and regenerative medicine. With the development of high-throughput deep sequencing technology, it has been found that lncRNAs play a very important role in the regulation of embryonic development. In this study, key lncRNAs that regulate embryonic development were screened by analyzing the expression pattern of lncRNAs in porcine in vivo fertilization (IVV) embryos. By knocking down lncRNA expression in in vitro fertilization (IVF) embryos, we investigated its function and mechanism of regulating embryonic development. The results showed that the expression pattern of lncRNA was consistent with the time of gene activation. The lncRNAs were highly expressed in the 4-cell to blastocyst stage but barely expressed in the oocytes and 2-cell stage. So we speculated this part of lncRNAs may regulate gene expression. The lncRNA LOC102165808 (named lncT because the gene near this lncRNA is TFAP2C) was one of them. The knockdown (KD) of lncT inhibited embryonic development, resulting in decreased H3K4me3, H3K4me2, and H3K9me3, and increased DNA methylation. Meanwhile, RNAseq showed SIN3A was the top decreased gene in lncT-KD embryos. There was a severe blastocyst formation defect in SIN3A-KD embryos. Both lncT and SIN3A could affect NANOG and induce more cell apoptosis. In conclusion, the knockdown of lncT inhibits embryonic development by regulating H3K4me3, H3K4me2, DNA methylation, pluripotency gene, and apoptosis, and SIN3A is one of the downstream genes of lncT in regulating embryonic development.

Transcriptome-guided annotation and functional classification of long non-coding RNAs in Arabidopsis thaliana

Long non-coding RNAs (lncRNAs) are a prominent class of eukaryotic regulatory genes. Despite the numerous available transcriptomic datasets, the annotation of plant lncRNAs remains based on dated annotations that have been historically carried over. We present a substantially improved annotation of Arabidopsis thaliana lncRNAs, generated by integrating 224 transcriptomes in multiple tissues, conditions, and developmental stages. We annotate 6764 lncRNA genes, including 3772 that are novel. We characterize their tissue expression patterns and find 1425 lncRNAs are co-expressed with coding genes, with enriched functional categories such as chloroplast organization, photosynthesis, RNA regulation, transcription, and root development. This improved transcription-guided annotation constitutes a valuable resource for studying lncRNAs and the biological processes they may regulate.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle

BACKGROUND

A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle.

RESULTS

Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits.

CONCLUSIONS

We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits.

Mammary epithelial cell transcriptome reveals potential roles of lncRNAs in regulating milk synthesis pathways in Jersey and Kashmiri cattle

Background

Long noncoding RNAs (lncRNAs) are now proven as essential regulatory elements, playing diverse roles in many biological processes including mammary gland development. However, little is known about their roles in the bovine lactation process.

Results

To identify and characterize the roles of lncRNAs in bovine lactation, high throughput RNA sequencing data from Jersey (high milk yield producer), and Kashmiri cattle (low milk yield producer) were utilized. Transcriptome data from three Kashmiri and three Jersey cattle throughout their lactation stages were utilized for differential expression analysis. At each stage (early, mid and late) three samples were taken from each breed. A total of 45 differentially expressed lncRNAs were identified between the three stages of lactation. The differentially expressed lncRNAs were found co-expressed with genes involved in the milk synthesis processes such as GPAM, LPL, and ABCG2 indicating their potential regulatory effects on milk quality genes. KEGG pathways analysis of potential cis and trans target genes of differentially expressed lncRNAs indicated that 27 and 48 pathways were significantly enriched between the three stages of lactation in Kashmiri and Jersey respectively, including mTOR signaling, PI3K-Akt signaling, and RAP1 signaling pathways. These pathways are known to play key roles in lactation biology and mammary gland development.

Conclusions

Expression profiles of lncRNAs across different lactation stages in Jersey and Kashmiri cattle provide a valuable resource for the study of the regulatory mechanisms involved in the lactation process as well as facilitate understanding of the role of lncRNAs in bovine lactation biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08406-x.

Long non-coding RNAs and microRNAs as crucial regulators in cardio-oncology

Cancer is one of the leading causes of morbidity and mortality worldwide. Significant improvements in the modern era of anticancer therapeutic strategies have increased the survival rate of cancer patients. Unfortunately, cancer survivors have an increased risk of cardiovascular diseases, which is believed to result from anticancer therapies. The emergence of cardiovascular diseases among cancer survivors has served as the basis for establishing a novel field termed cardio-oncology. Cardio-oncology primarily focuses on investigating the underlying molecular mechanisms by which anticancer treatments lead to cardiovascular dysfunction and the development of novel cardioprotective strategies to counteract cardiotoxic effects of cancer therapies. Advances in genome biology have revealed that most of the genome is transcribed into non-coding RNAs (ncRNAs), which are recognized as being instrumental in cancer, cardiovascular health, and disease. Emerging studies have demonstrated that alterations of these ncRNAs have pathophysiological roles in multiple diseases in humans. As it relates to cardio-oncology, though, there is limited knowledge of the role of ncRNAs. In the present review, we summarize the up-to-date knowledge regarding the roles of long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) in cancer therapy-induced cardiotoxicities. Moreover, we also discuss prospective therapeutic strategies and the translational relevance of these ncRNAs.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

Expression of lncRNAs in response to bacterial infections of goat mammary epithelial cells reveals insights into mammary gland diseases

Mastitis or inflammation of the mammary gland is a highly economic and deadly alarming disease for the dairy sector as well as policymakers caused by microbial infection. Transcriptomic and proteomic approaches have been widely employed to identify the underlying molecular mechanisms of bacterial infections in the mammary gland. Numerous differentially expressed mRNAs, miRNAs, and proteins together with their associated signaling pathways have been identified during bacterial infection, paving the way for analysis of their biological functions. Long noncoding RNAs (lncRNAs) are important regulators of multiple biological processes. However, little is known regarding their role in bacterial infection in mammary epithelial cells. Hence, RNA-sequencing was performed by infecting primary mammary epithelial cells (pMECs) with both gram-negative (E. coli) and gram-positive bacteria (S. aureus). Using stringent pipeline, a set of 1957 known and 1175 novel lncRNAs were identified, among which, 112 lncRNAs were found differentially expressed in bacteria challenged PMECs compared with the control. Additionally, potential targets of the lncRNAs were predicted in cis- and trans-configuration. KEGG analysis revealed that DE lncRNAs were associated with at least 15 immune-related pathways. Therefore, our study revealed that bacterial challenge triggers the expression of lncRNAs associated with immune response and defense mechanisms in goat mammary epithelial cells.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

SARS-CoV-2, Cardiovascular Diseases, and Noncoding RNAs: A Connected Triad

Coronavirus Disease 2019 (COVID-19), caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), is characterized by important respiratory impairments frequently associated with severe cardiovascular damages. Moreover, patients with pre-existing comorbidity for cardiovascular diseases (CVD) often present a dramatic increase in inflammatory cytokines release, which increases the severity and adverse outcomes of the infection and, finally, mortality risk. Despite this evident association at the clinical level, the mechanisms linking CVD and COVID-19 are still blurry and unresolved. Noncoding RNAs (ncRNAs) are functional RNA molecules transcribed from DNA but usually not translated into proteins. They play an important role in the regulation of gene expression, either in relatively stable conditions or as a response to different stimuli, including viral infection, and are therefore considered a possible important target in the design of specific drugs. In this review, we introduce known associations and interactions between COVID-19 and CVD, discussing the role of ncRNAs within SARS-CoV-2 infection from the perspective of the development of efficient pharmacological tools to treat COVID-19 patients and taking into account the equally dramatic associated consequences, such as those affecting the cardiovascular system.

Transcriptome Analysis Revealed Potential Mechanisms of Resistance to Trichomoniasis gallinae Infection in Pigeon (Columba livia)

Trichomoniasis gallinae (T. gallinae) is one of the most pathogenic parasites in pigeon, particularly in squabs. Oral cavity is the main site for the host-parasite interaction. Herein, we used RNA-sequencing technology to characterize lncRNA and mRNA profiles and compared transcriptomic dynamics of squabs, including four susceptible birds (S) from infected group, four tolerant birds (T) without parasites after T. gallinae infection, and three birds from uninfected group (N), to understand molecular mechanisms underlying host resistance to this parasite. We identified 29,809 putative lncRNAs and characterized their genomic features subsequently. Differentially expressed (DE) genes, DE-lncRNAs and cis/trans target genes of DE-lncRNAs were further compared among the three groups. The KEGG analysis indicated that specific intergroup DEGs were involved in carbon metabolism (S vs. T), metabolic pathways (N vs. T) and focal adhesion pathway (N vs. S), respectively. Whereas, the cis/trans genes of DE-lncRNAs were enriched in cytokine-cytokine receptor interaction, toll-like receptor signaling pathway, p53 signaling pathway and insulin signaling pathway, which play crucial roles in immune system of the host animal. This suggests T. gallinae invasion in pigeon mouth may modulate lncRNAs expression and their target genes. Moreover, co-expression analysis identified crucial lncRNA-mRNA interaction networks. Several DE-lncRNAs including MSTRG.82272.3, MSTRG.114849.42, MSTRG.39405.36, MSTRG.3338.5, and MSTRG.105872.2 targeted methylation and immune-related genes, such as JCHAIN, IL18BP, ANGPT1, TMRT10C, SAMD9L, and SOCS3. This implied that DE-lncRNAs exert critical influence on T. gallinae infections. The quantitative exploration of host transcriptome changes induced by T. gallinae infection broaden both transcriptomic and epigenetic insights into T. gallinae resistance and its pathological mechanism.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

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.

PacBio Single-Molecule Long-Read Sequencing Provides New Light on the Complexity of Full-Length Transcripts in Cattle

Cattle (Bos taurus) is one of the most widely distributed livestock species in the world, and provides us with high-quality milk and meat which have a huge impact on the quality of human life. Therefore, accurate and complete transcriptome and genome annotation are of great value to the research of cattle breeding. In this study, we used error-corrected PacBio single-molecule real-time (SMRT) data to perform whole-transcriptome profiling in cattle. Then, 22.5 Gb of subreads was generated, including 381,423 circular consensus sequences (CCSs), among which 276,295 full-length non-chimeric (FLNC) sequences were identified. After correction by Illumina short reads, we obtained 22,353 error-corrected isoforms. A total of 305 alternative splicing (AS) events and 3,795 alternative polyadenylation (APA) sites were detected by transcriptome structural analysis. Furthermore, we identified 457 novel genes, 120 putative transcription factors (TFs), and 569 novel long non-coding RNAs (lncRNAs). Taken together, this research improves our understanding and provides new insights into the complexity of full-length transcripts in cattle.

Integrated Analysis of lncRNA and mRNA in Subcutaneous Adipose Tissue of Ningxiang Pig

Simple Summary

This study shows the transcription profiles and the functional network in lncRNA and mRNA in the subcutaneous adipose tissue of Ningxiang piglets in four stages of development (piglets, nursery pigs, early fattening, and late fattening). A total of 2872 novel lncRNAs have now been determined. A total of 10,084 DEmRNAs and 931 DElncRNAs were determined. Interestingly, most DEmRNAs were up-regulated in the piglet stage and, in contrast, DElncRNAs were up-regulated in the late fattening stage. A complicated interaction between mRNAs and lncRNAs was determined via STEM and WGCNA, demonstrating that lncRNAs are an essential regulatory component in mRNAs. Modules 2 and 5 shows a similar mode of transcriptions for both mRNA and lncRNA, which are mainly involved in steroid biosynthesis, glycosphingolipid biosynthesis, metabolic pathways, and glycerolipid metabolism. The transcription levels of mRNAs and lncRNAs for both modules were higher in the early and late fattening stage. This may be explained by the active fatty acids, sterols, steroids, and lipid-related metabolic activity in the subcutaneous adipose tissue during the early and late fattening stage.

Abstract

Ningxiang pigs, a Chinese bred pig known for its tender meat and high quality unsaturated fatty acids. This study discovers the transcription profiles and functional networks in long non-coding RNA (lncRNA) and messenger RNA (mRNA) in subcutaneous adipose tissue. Subcutaneous adipose tissue was collected from piglet, nursery pig, early fattening, and late fattening stage of Ningxiang piglets, and lncRNA and mRNA transcription of each stage was profiled. A total of 339,204,926 (piglet), 315,609,246 (nursery), 266,798,202 (early fattening), and 343,740,308 (late fattening) clean reads were generated, and 2872 novel lncRNAs were identified. Additionally, 10,084 differential mRNAs (DEmRNAs) and 931 differential lncRNAs were determined. Most DEmRNAs were up-regulated in the piglet stage, while they were down-regulated in late fattening stage. A complicated interaction between mRNAs and lncRNAs was identified via STEM and WGCNA, demonstrated that lncRNAs are a significant regulatory component in mRNAs. The findings showed that modules 2 and 5 have a similar mode of transcription for both mRNA and lncRNA, and were mainly participated in steroid biosynthesis, glycosphingolipid biosynthesis, metabolic pathways, and glycerolipid metabolism. The mRNAs and lncRNAs transcription levels of both modules was higher in the early and late fattening stage, which may be due to the active activity of the metabolism in relation to fatty acids, sterols, steroids, and lipids in the subcutaneous adipose tissue during the early and late fattening stage. These findings could be expected to result in further research of the functional properties of lncRNA from subcutaneous adipose tissue at different stages of development in Ningxiang pigs.

Identification and expression analysis of lncRNA in seven organs of Rhinopithecus roxellana

Long non-coding RNA (lncRNA) represents a new direction to identify expression profiles and regulatory mechanisms in various organisms. Here, we report the first dataset of lncRNAs of the golden snub-nosed monkey (GSM), including 12,557 putative lncRNAs identified from seven organs. Compared with mRNA, GSM lncRNA had fewer exons and isoforms, and longer length. LncRNA showed more obvious tissue-specific expression than mRNA. However, for the top ten most abundant genes in each organ, mRNAs expression was more tissue-specific than lncRNAs. By identification of specifically expressed lncRNAs and mRNAs in each organ, it indicates that the expression of SEG-lncRNA (specifically expressed lncRNA) and SEG-mRNA (specifically expressed mRNA) had high correlation. In particular, combined our lncRNA and mRNA data, we identified 92 heart SEG-lncRNAs targeted ten mRNA genes in the oxidative phosphorylation pathway and upregulated the expression of these target genes such as ND4, ATP6, and ATP8. These may contribute to GSM adaption to its high-elevation environment. We also identified 171 liver SEG-lncRNAs, which targeted 27 genes associated with the metabolism of xenobiotics and leaded to high expression of these target genes in liver. These lncRNAs may play important roles in GSM adaptation to a folivory diet.

Tissue Resources for the Functional Annotation of Animal Genomes

In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

Profiling and integrated analysis of whole-transcriptome changes in uterine caruncles of pregnant and non-pregnant buffaloes

Improved reproductive performance in buffaloes can be achieved by understanding the basic mechanism governing the embryonic attachment and feto-maternal communication. Considering this, trascriptomic profiling and integrative analysis of long intergenic non-coding RNAs were carried out in the uterine caruncles of pregnant and non-pregnant buffaloes. Transcriptome data of pregnant and non-pregnant uterine caruncles after quality control was used to perform the analysis. Total of 86 novel lincRNAs expressed in uterine caruncular tissues were identified and characterized. Differential expression analysis revealed that 447 mRNAs and 185 mRNAs were up- and down- regulated, respectively. The number of up- and down- regulated lincRNAs were 114 and 13, respectively. Of the identified 86 novel lincRNAs, six novel lincRNAs were up-regulated in the pregnant uterine caruncles. GO terms (biological process) and PANTHER pathways associated with reproduction and embryogenesis were over-represented in differentially expressed genes. Through miRNA interaction analysis, interactions of 16 differentially expressed lincRNAs with mi-RNAs involved in reproduction were identified. This study has provided a catalogue of differentially expressed genes and novel regions previously unknown to play a significant role in buffalo reproduction. The results from the current study extends the buffalo uterine lncRNAs database and provides candidate regulators for future molecular genetic studies on buffalo uterine physiology to improve the embryo implantation and successful completion of pregnancy.

Transcriptome Analysis Reveals the Profile of Long Non-coding RNAs During Chicken Muscle Development

The developmental complexity of muscle arises from elaborate gene regulation. Long non-coding RNAs (lncRNAs) play critical roles in muscle development through the regulation of transcription and post-transcriptional gene expression. In chickens, previous studies have focused on the lncRNA profile during the embryonic periods, but there are no studies that explore the profile from the embryonic to post-hatching period. Here, we reconstructed 14,793 lncRNA transcripts and identified 2,858 differentially expressed lncRNA transcripts and 4,282 mRNAs from 12-day embryos (E12), 17-day embryos (E17), 1-day post-hatch chicks (D1), 14-day post-hatch chicks (D14), 56-day post-hatch chicks (D56), and 98-day post-hatch chicks (D98), based on our published RNA-seq datasets. We performed co-expression analysis for the differentially expressed lncRNAs and mRNAs, using STEM, and identified two profiles with opposite expression trends: profile 4 with a downregulated pattern and profile 21 with an upregulated pattern. The cis- and trans-regulatory interactions between the lncRNAs and mRNAs were predicted within each profile. Functional analysis of the lncRNA targets showed that lncRNAs in profile 4 contributed to the cell proliferation process, while lncRNAs in profile 21 were mainly involved in metabolism. Our work highlights the lncRNA profiles involved in the development of chicken breast muscle and provides a foundation for further experiments on the role of lncRNAs in the regulation of muscle development.

Identification of Long Non-coding RNA Isolated From Naturally Infected Macrophages and Associated With Bovine Johne's Disease in Canadian Holstein Using a Combination of Neural Networks and Logistic Regression

Mycobacterium avium ssp. paratuberculosis (MAP) causes chronic enteritis in most ruminants. The pathogen MAP causes Johne's disease (JD), a chronic, incurable, wasting disease. Weight loss, diarrhea, and a gradual drop in milk production characterize the disease's clinical phase, culminating in death. Several studies have characterized long non-coding RNA (lncRNA) in bovine tissues, and a previous study characterizes (lncRNA) in macrophages infected with MAP in vitro. In this study, we aim to characterize the lncRNA in macrophages from cows naturally infected with MAP. From 15 herds, feces and blood samples were collected for each cow older than 24 months, twice yearly over 3–5 years. Paired samples were analyzed by fecal PCR and blood ELISA. We used RNA-seq data to study lncRNA in macrophages from 33 JD(+) and 33 JD(–) dairy cows. We performed RNA-seq analysis using the “new Tuxedo” suite. We characterized lncRNA using logistic regression and multilayered neural networks and used DESeq2 for differential expression analysis and Panther and Reactome classification systems for gene ontology (GO) analysis. The study identified 13,301 lncRNA, 605 of which were novel lncRNA. We found seven genes close to differentially expressed lncRNA, including CCDC174, ERI1, FZD1, TWSG1, ZBTB38, ZNF814, and ZSCAN4. None of the genes associated with susceptibility to JD have been cited in the literature. LncRNA target genes were significantly enriched for biological process GO terms involved in immunity and nucleic acid regulation. These include the MyD88 pathway (TLR5), GO:0043312 (neutrophil degranulation), GO:0002446 (neutrophil-mediated immunity), and GO:0042119 (neutrophil activation). These results identified lncRNA with potential roles in host immunity and potential candidate genes and pathways through which lncRNA might function in response to MAP infection.

Parsing the microRNA genetics basis regulating skeletal muscle fiber types and meat quality traits in pigs

Muscle fibers are closely related to human diseases and livestock meat quality. However, the genetics basis of microRNAs (miRNAs) in regulating muscle fibers is not completely understood. In this study, we constructed the whole genome-wide miRNA expression profiles of porcine fast-twitch muscle [biceps femoris (Bf)] and slow-twitch muscle [soleus (Sol)], and identified hundreds of miRNAs, including four skeletal muscle-highly expressed miRNAs, ssc-miR-378, ssc-let-7f, ssc-miR-26a, and ssc-miR-27b-3p. Moreover, we identified 63 differentially expressed (DE) miRNAs between biceps femoris vs. soleus, which are the key candidate miRNAs regulating the skeletal muscle fiber types. In addition, we found that the expression of DE ssc-miR-499-5p was significantly correlated to the expression of Myoglobin (r = 0.6872, P < 0.0001) and Myosin heavy chain 7 (MYH7; r = 0.5408, P = 0.0020), and pH_45 min (r = 0.3806, P = 0.0380) and glucose content (r = -0.4382, P = 0.0154); while the expression of DE ssc-miR-499-3p was significantly correlated to the expression of Myoglobin (r = 0.5340, P = 0.0024) and pH_45 min (r = 0.4857, P = 0.0065). Taken together, our data established a sound foundation for further studies on the regulatory mechanisms of miRNAs in skeletal muscle fiber conversion and meat quality traits in livestock, and could provide a genetic explanation of the role of miRNAs in human muscular diseases.

Insight into Liver lncRNA and mRNA Profiling at Four Developmental Stages in Ningxiang Pig

Ningxiang pigs, a fat-type pig, are native to Ningxiang County in Hunan Province, with thousands of years of breeding history. This study aims to explore the expression profiles and functional networks on messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) in the liver. Liver tissue of Ningxiang piglets was collected at 30, 90, 150, and 210 days after birth (four development stages), and the mRNA and lncRNA expression was profiled. Compared to mRNA and lncRNA expression profiles, most differentially expressed mRNAs (DEmRNAs) were upregulated at 30 days; however, most DElncRNAs were downregulated at 210 days. Via Short Time-series Expression Miner (STEM) analysis and weighted gene co-expression network analysis (WGCNA), a complex interaction between mRNAs and lncRNAs was identified, indicating that lncRNAs may be a critical regulatory element for mRNAs. One module of genes in particular (module profile 4) was related to fibril organization, vasculogenesis, GTPase activator activity, and regulation of kinase activity. The mRNAs and lncRNAs in module profile 4 had a similar pattern of expression, indicating that they have functional and regulatory relationships. Only CAV1, PACSIN2, and CDC42 in the particular mRNA profile 4 were the target genes of lncRNAs in that profile, which shows the possible regulatory relationship between lncRNAs and mRNAs. The expression of these genes and lncRNAs in profile 4 was the highest at 30 days, and it is believed that these RNAs may play a critical role during the suckling period in order to meet the dietary requirements of piglets. In the lncRNA-mRNA co-expression network, the identified gene hubs and associated lncRNAs were shown to be involved in saccharide, lipid, and glucose metabolism, which may play an important role in the development and health of the liver. This result will lead to further investigation of liver lncRNA functions at various stages of development in Ningxiang pigs.