Identification of long non-coding RNA in the horse transcriptome

Evaluation of Muscle Long Non-Coding RNA Profile during Rearing and Finishing Phase of Bulls Subjected to Different Prenatal Nutritional Strategies

Maternal nutrition has the ability of influence critical processes in fetal life, including muscle development. Also, in this period, epigenetic sensitivity to external stimuli is higher and produces long-lasting effects. Thus, the aim of this study was to investigate epigenetic mechanisms, including the identification and characterization of long non-coding RNA (lncRNA) from animals that had undergone different strategies of prenatal supplementation. A group of Nellore cows (n = 126) were separated into three nutritional plans: NP (control)-Not Programmed, without protein-energy supplementation; PP-Partially Programmed, protein-energy supplementation in the final third of pregnancy; and CP-Complete Programming, protein-energy supplementation during the full period of gestation. A total of 63 male offspring were used in this study, of which 15 (5 per treatment) had Longissimus thoracis muscle at 15 (biopsy) and 22 months (slaughter). Biopsy samples were subjected to RNA extraction and sequencing. Differential expression (DE) of remodeling factors and chromatin-modifying enzyme genes were performed. For the identification and characterization of lncRNA, a series of size filters and protein coding potential tests were performed. The lncRNAs identified had their differential expression and regulatory potential tested. Regarding DE of epigenetic mechanisms, no differentially expressed gene was found (p > 0.1). Identification of potential lncRNA was successful, identifying 1823 transcripts at 15 months and 1533 at 22 months. Among these, four were considered differentially expressed between treatments at 15 months and 6 were differentially expressed at 22 months. Yet, when testing regulatory potential, 13 lncRNAs were considered key regulators in the PP group, and 17 in the CP group. PP group lncRNAs possibly regulate fat-cell differentiation, in utero embryonic development, and transforming growth factor beta receptor, whereas lncRNA in the CP group regulates in utero embryonic development, fat-cell differentiation and vasculogenesis. Maternal nutrition had no effect on differential expression of epigenetic mechanisms; however, it seems to impair lncRNA regulation of epigenetics.

Genome-based prediction of cross-protective, HLA-DR-presented epitopes as putative vaccine antigens for multiple Bordetella species

IMPORTANCE

Pertussis, caused by Bordetella pertussis, can cause debilitating respiratory symptoms, so whole-cell pertussis vaccines (wPVs) were introduced in the 1940s. However, reactogenicity of wPV necessitated the development of acellular pertussis vaccines (aPVs) that were introduced in the 1990s. Since then, until the COVID-19 pandemic began, reported pertussis incidence was increasing, suggesting that aPVs do not induce long-lasting immunity and may not effectively prevent transmission. Additionally, aPVs do not provide protection against other Bordetella species that are observed during outbreaks. The significance of this work is in determining potential new vaccine antigens for multiple Bordetella species that are predicted to elicit long-term immune responses. Genome-based approaches have aided the development of novel vaccines; here, these methods identified Bordetella vaccine candidates that may be cross-protective and predicted to induce strong memory responses. These targets can lead to an improved vaccine with a strong safety profile while also strengthening the longevity of the immune response.

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.

Mutation of a major CG methylase alters genome-wide lncRNA expression in rice

Plant long non-coding RNAs (lncRNAs) function in diverse biological processes, and lncRNA expression is under epigenetic regulation, including by cytosine DNA methylation. However, it remains unclear whether 5-methylcytosine (5mC) plays a similar role in different sequence contexts (CG, CHG, and CHH). In this study, we characterized and compared the profiles of genome-wide lncRNA profiles (including long intergenic non-coding RNAs [lincRNAs] and long noncoding natural antisense transcripts [lncNATs]) of a null mutant of the rice DNA methyltransferase 1, OsMET1-2 (designated OsMET1-2-/-) and its isogenic wild type (OsMET1-2+/+). The En/Spm transposable element (TE) family, which was heavily methylated in OsMET1-2+/+, was transcriptionally de-repressed in OsMET1-2-/- due to genome-wide erasure of CG methylation, and this led to abundant production of specific lncRNAs. In addition, RdDM-mediated CHH hypermethylation was increased in the 5'-upstream genomic regions of lncRNAs in OsMET1-2-/-. The positive correlation between the expression of lincRNAs and that of their proximal protein-coding genes was also analyzed. Our study shows that CG methylation negatively regulates the TE-related expression of lncRNA and demonstrates that CHH methylation is also involved in the regulation of lncRNA expression.

Derivation of Intermediate Pluripotent Stem Cells Amenable to Primordial Germ Cell Specification

Dynamic pluripotent stem cell (PSC) states are in vitro adaptations of pluripotency continuum in vivo. Previous studies have generated a number of PSCs with distinct properties. To date, however, no known PSCs have demonstrated dual competency for chimera formation and direct responsiveness to primordial germ cell (PGC) specification, a unique functional feature of formative pluripotency. Here, by modulating fibroblast growth factor (FGF), transforming growth factor β (TGF-β), and WNT pathways, we derived PSCs from mice, horses, and humans (designated as XPSCs) that are permissive for direct PGC-like cell induction in vitro and are capable of contributing to intra- or inter-species chimeras in vivo. XPSCs represent a pluripotency state between naive and primed pluripotency and harbor molecular, cellular, and phenotypic features characteristic of formative pluripotency. XPSCs open new avenues for studying mammalian pluripotency and dissecting the molecular mechanisms governing PGC specification. Our method may be broadly applicable for the derivation of analogous stem cells from other mammalian species.

An integrative atlas of chicken long non-coding genes and their annotations across 25 tissues

Long non-coding RNAs (LNC) regulate numerous biological processes. In contrast to human, the identification of LNC in farm species, like chicken, is still lacunar. We propose a catalogue of 52,075 chicken genes enriched in LNC (http://www.fragencode.org/), built from the Ensembl reference extended using novel LNC modelled here from 364 RNA-seq and LNC from four public databases. The Ensembl reference grew from 4,643 to 30,084 LNC, of which 59% and 41% with expression ≥ 0.5 and ≥ 1 TPM respectively. Characterization of these LNC relatively to the closest protein coding genes (PCG) revealed that 79% of LNC are in intergenic regions, as in other species. Expression analysis across 25 tissues revealed an enrichment of co-expressed LNC:PCG pairs, suggesting co-regulation and/or co-function. As expected LNC were more tissue-specific than PCG (25% vs. 10%). Similarly to human, 16% of chicken LNC hosted one or more miRNA. We highlighted a new chicken LNC, hosting miR155, conserved in human, highly expressed in immune tissues like miR155, and correlated with immunity-related PCG in both species. Among LNC:PCG pairs tissue-specific in the same tissue, we revealed an enrichment of divergent pairs with the PCG coding transcription factors, as for example LHX5, HXD3 and TBX4, in both human and chicken.

Adding exogenous biglycan or decorin improves tendon formation for equine peritenon and tendon proper cells in vitro

Background

Tendon injuries amount to one of the leading causes of career-ending injuries in horses due to the inability for tendon to completely repair and the high reinjury potential. As a result, novel therapeutics are necessary to improve repair with the goal of decreasing leg lameness and potential reinjury. Small leucine-rich repeat proteoglycans (SLRPs), a class of regulatory molecules responsible for collagen organization and maturation, may be one such therapeutic to improve tendon repair. Before SLRP supplementation can occur in vivo, proper evaluation of the effect of these molecules in vitro needs to be assessed. The objective of this study was to evaluate the effectiveness of purified bovine biglycan or decorin on tendon proper and peritenon cell populations in three-dimensional tendon constructs.

Methods

Equine tendon proper or peritenon cell seeded fibrin three-dimensional constructs were supplemented with biglycan or decorin at two concentrations (5 nM or 25 nM). The functionality and ultrastructural morphology of the constructs were assessed using biomechanics, collagen content analysis, transmission electron microscopy (TEM), and gene expression by real time – quantitative polymerase chain reaction (RT-qPCR).

Results

SLRP supplementation affected both tendon proper and peritenon cells-seeded constructs. With additional SLRPs, material and tensile properties of constructs strengthened, though ultrastructural analyses indicated production of similar-sized or smaller fibrils. Overall expression of tendon markers was bolstered more in peritenon cells supplemented with either SLRP, while supplementation of SLRPs to TP cell-derived constructs demonstrated fewer changes in tendon and extracellular matrix markers. Moreover, relative to non-supplemented tendon proper cell-seeded constructs, SLRP supplementation of the peritenon cells showed increases in mechanical strength, material properties, and collagen content.

Conclusions

The SLRP-supplemented peritenon cells produced constructs with greater mechanical and material properties than tendon proper seeded constructs, as well as increased expression of matrix assembly molecules. These findings provide evidence that SLRPs should be further investigated for their potential to improve tendon formation in engineered grafts or post-injury.

Comparison of Poly-A+ Selection and rRNA Depletion in Detection of lncRNA in Two Equine Tissues Using RNA-seq

Long non-coding RNAs (lncRNAs) are untranslated regulatory transcripts longer than 200 nucleotides that can play a role in transcriptional, post-translational, and epigenetic regulation. Traditionally, RNA-sequencing (RNA-seq) libraries have been created by isolating transcriptomic RNA via poly-A+ selection. In the past 10 years, methods to perform ribosomal RNA (rRNA) depletion of total RNA have been developed as an alternative, aiming for better coverage of whole transcriptomic RNA, both polyadenylated and non-polyadenylated transcripts. The purpose of this study was to determine which library preparation method is optimal for lncRNA investigations in the horse. Using liver and cerebral parietal lobe tissues from two healthy Thoroughbred mares, RNA-seq libraries were prepared using standard poly-A+ selection and rRNA-depletion methods. Averaging the two biologic replicates, poly-A+ selection yielded 327 and 773 more unique lncRNA transcripts for liver and parietal lobe, respectively. More lncRNA were found to be unique to poly-A+ selected libraries, and rRNA-depletion identified small nucleolar RNA (snoRNA) to have a higher relative expression than in the poly-A+ selected libraries. Overall, poly-A+ selection provides a more thorough identification of total lncRNA in equine tissues while rRNA-depletion may allow for easier detection of snoRNAs.

Identification and Comparative Analysis of Long Non-Coding RNA in the Skeletal Muscle of Two Dezhou Donkey Strains

Long non-coding RNA (lncRNA) has been extensively studied in many livestock. However, compared with other livestock breeds, there is less research regarding donkey lncRNA function. It has been reported that lncRNA plays an important role in the timing control of development, aging, and death of livestock. Therefore, the study of donkey skeletal muscle lncRNA is of great significance for exploring donkey meat production performance. In this study, RNA-Seq was used to perform high-throughput sequencing of skeletal muscle (longissimus dorsi and gluteus) of two Dezhou donkey strains (SanFen and WuTou). The differentially expressed lncRNAs were screened between different strains and tissues. Then candidate genes for conjoint analysis were screened based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Finally, the accuracy of the sequencing data was verified by real-time quantitative polymerase chain reaction (RT-qPCR). Herein, we identified 3869 novel lncRNAs and 73 differentially expressed lncRNAs. Through the comparison between groups, the specific expression of lncRNAs were found in different strains and muscle tissues. Importantly, we constructed the lncRNA-miRNA-mRNA interaction network and found three important candidate lncRNAs (MSTRG.9787.1, MSTRG.3144.1, and MSTRG.9886.1) and four candidate genes (ACTN1, CDON, FMOD, and BMPR1B). Analysis of the KEGG pathway indicates that the TGF-β signaling pathway plays a pivotal role in the growth and development of skeletal muscle in Dezhou donkey strains.

miRNAs and lncRNAs in Echinococcus and Echinococcosis

Echinococcosis are considered to be potentially lethal zoonotic diseases that cause serious damage to hosts. The metacestode of Echinococcus multilocularis and E. granulosus can result in causing the alveolar and cystic echinococcoses, respectively. Recent studies have shown that non-coding RNAs are widely expressed in Echinococcus spp. and hosts. In this review, the two main types of non-coding RNAs—long non-coding RNAs (lncRNAs) and microRNAs (miRNAs)—and the wide-scale involvement of these molecules in these parasites and their hosts were discussed. The expression pattern of miRNAs in Echinococcus spp. is species- and developmental stage-specific. Furthermore, common miRNAs were detected in three Echinococcus spp. and their intermediate hosts. Here, we primarily focus on recent insights from transcriptome studies, the expression patterns of miRNAs and lncRNAs, and miRNA-related databases and techniques that are used to investigate miRNAs in Echinococcus and echinococcosis. This review provides new avenues for screening therapeutic and diagnostic markers.

Identification of Novel lncRNAs Differentially Expressed in Placentas of Chinese Ningqiang Pony and Yili Horse Breeds

As a nutrient sensor, the placenta plays a key role in regulating fetus growth and development. Long non-coding RNAs (lncRNAs) have been shown to regulate growth-related traits. However, the biological function of lncRNAs in horse placentas remains unclear. To compare the expression patterns of lncRNAs in the placentas of the Chinese Ningqiang (NQ) and Yili (YL) breeds, we performed a transcriptome analysis using RNA sequencing (RNA-seq) technology. NQ is a pony breed with an average adult height at the withers of less than 106 cm, whereas that of YL is around 148 cm. Based on 813 million high-quality reads and stringent quality control procedures, 3011 transcripts coding for 1464 placental lncRNAs were identified and mapped to the horse reference genome. We found 107 differentially expressed lncRNAs (DELs) between NQ and YL, including 68 up-regulated and 39 down-regulated DELs in YL. Six (TBX3, CACNA1F, EDN3, KAT5, ZNF281, TMED2, and TGFB1) out of the 233 genes targeted by DELs were identified as being involved in limb development, skeletal myoblast differentiation, and embryo development. Two DELs were predicted to target the TBX3 gene, which was found to be under strong selection and associated with small body size in the Chinese Debao pony breed. This finding suggests the potential functional significance of placental lncRNAs in regulating horse body size.

A Gene Expression Atlas of the Domestic Water Buffalo (Bubalus bubalis)

The domestic water buffalo (Bubalus bubalis) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.

Distribution, Diversity, and Evolution of Endogenous Retroviruses in Perissodactyl Genomes

The evolution of mammalian genomes has been shaped by interactions with endogenous retroviruses (ERVs). In this study, we investigated the distribution and diversity of ERVs in the mammalian order Perissodactyla, with a view to understanding their impact on the evolution of modern equids (family Equidae). We characterize the major ERV lineages in the horse genome in terms of their genomic distribution, ancestral genome organization, and time of activity. Our results show that subsequent to their ancestral divergence from rhinoceroses and tapirs, equids acquired four novel ERV lineages. We show that two of these ERV lineages proliferated extensively in the lineage leading to modern horses, and one contains loci that are actively transcribed in specific tissues. In addition, we show that the white rhinoceros has resisted germ line colonization by retroviruses for more than 54 million years-longer than any other extant mammalian species. The map of equine ERVs that we provide here will be of great utility to future studies aiming to investigate the potential functional roles of equine ERVs and their impact on equine evolution.IMPORTANCE ERVs in the host genome are highly informative about the long-term interactions of retroviruses and hosts. They are also interesting because they have influenced the evolution of mammalian genomes in various ways. In this study, we derive a calibrated timeline describing the process through which ERV diversity has been generated in the equine germ line. We determined the distribution and diversity of perissodactyl ERV lineages and inferred their retrotranspositional activity during evolution, thereby gaining insight into the long-term coevolutionary history of retroviruses and mammals. Our study provides a platform for future investigations to identify equine ERV loci involved in physiological processes and/or pathological conditions.

Functional Annotation of Animal Genomes (FAANG): Current Achievements and Roadmap

Functional annotation of genomes is a prerequisite for contemporary basic and applied genomic research, yet farmed animal genomics is deficient in such annotation. To address this, the FAANG (Functional Annotation of Animal Genomes) Consortium is producing genome-wide data sets on RNA expression, DNA methylation, and chromatin modification, as well as chromatin accessibility and interactions. In addition to informing our understanding of genome function, including comparative approaches to elucidate constrained sequence or epigenetic elements, these annotation maps will improve the precision and sensitivity of genomic selection for animal improvement. A scientific community-driven effort has already created a coordinated data collection and analysis enterprise crucial for the success of this global effort. Although it is early in this continuing process, functional data have already been produced and application to genetic improvement reported. The functional annotation delivered by the FAANG initiative will add value and utility to the greatly improved genome sequences being established for domesticated animal species.

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

Background

Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs.

Results

Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen.

Conclusions

LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5037-7) contains supplementary material, which is available to authorized users.

Cross-species inference of long non-coding RNAs greatly expands the ruminant transcriptome

Background

mRNA-like long non-coding RNAs (lncRNAs) are a significant component of mammalian transcriptomes, although most are expressed only at low levels, with high tissue-specificity and/or at specific developmental stages. Thus, in many cases lncRNA detection by RNA-sequencing (RNA-seq) is compromised by stochastic sampling. To account for this and create a catalogue of ruminant lncRNAs, we compared de novo assembled lncRNAs derived from large RNA-seq datasets in transcriptional atlas projects for sheep and goats with previous lncRNAs assembled in cattle and human. We then combined the novel lncRNAs with the sheep transcriptional atlas to identify co-regulated sets of protein-coding and non-coding loci.

Results

Few lncRNAs could be reproducibly assembled from a single dataset, even with deep sequencing of the same tissues from multiple animals. Furthermore, there was little sequence overlap between lncRNAs that were assembled from pooled RNA-seq data. We combined positional conservation (synteny) with cross-species mapping of candidate lncRNAs to identify a consensus set of ruminant lncRNAs and then used the RNA-seq data to demonstrate detectable and reproducible expression in each species. In sheep, 20 to 30% of lncRNAs were located close to protein-coding genes with which they are strongly co-expressed, which is consistent with the evolutionary origin of some ncRNAs in enhancer sequences. Nevertheless, most of the lncRNAs are not co-expressed with neighbouring protein-coding genes.

Conclusions

Alongside substantially expanding the ruminant lncRNA repertoire, the outcomes of our analysis demonstrate that stochastic sampling can be partly overcome by combining RNA-seq datasets from related species. This has practical implications for the future discovery of lncRNAs in other species.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0391-0) contains supplementary material, which is available to authorized users.