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

A meta-analysis approach for annotation and identification of lncRNAs controlling perirenal fat deposition in suckling lambs

Long non-coding RNAs (lncRNAs) are being studied in farm animals due to their association with traits of economic interest, such as fat deposition. Based on the analysis of perirenal fat transcriptomes, this research explored the relevance of these regulatory elements to fat deposition in suckling lambs. To that end, meta-analysis techniques have been implemented to efficiently characterize and detect differentially expressed transcripts from two different RNA-seq datasets, one including samples of two sheep breeds that differ in fat deposition features, Churra and Assaf (n = 14), and one generated from Assaf suckling lambs with different fat deposition levels (n = 8). The joint analysis of the 22 perirenal fat RNA-seq samples with the FEELnc software allowed the detection of 3953 novel lncRNAs. After the meta-analysis, 251 differentially expressed genes were identified, 21 of which were novel lncRNAs. Additionally, a co-expression analysis revealed that, in suckling lambs, lncRNAs may play a role in controlling angiogenesis and thermogenesis, processes highlighted in relation to high and low fat deposition levels, respectively. Overall, while providing information that could be applied for the improvement of suckling lamb carcass traits, this study offers insights into the biology of perirenal fat deposition regulation in mammals.

Construction of a molecular regulatory network related to fat deposition by multi-tissue transcriptome sequencing of Jiaxian red cattle

Intramuscular fat (IMF) refers to the fat that accumulates between muscle bundles or within muscle cells, whose content significantly impacts the taste, tenderness, and flavor of meat products, making it a crucial economic characteristic in livestock production. However, the intricate mechanisms governing IMF deposition, involving non-coding RNAs (ncRNAs), genes, and complex regulatory networks, remain largely enigmatic. Identifying adipose tissue-specific genes and ncRNAs is paramount to unravel these molecular mysteries. This study, conducted on Jiaxian red cattle, harnessed whole transcriptome sequencing to unearth the nuances of circRNAs and miRNAs across seven distinct tissues. The interplay of these ncRNAs was assessed through differential expression analysis and network analysis. These findings are not only pivotal in unveiling the intricacies of fat deposition mechanisms but also lay a robust foundation for future research, setting the stage for enhancing IMF content in Jiaxian red cattle breeding.

Comparative epigenomics reveals the impact of ruminant-specific regulatory elements on complex traits

BACKGROUND

Insights into the genetic basis of complex traits and disease in both human and livestock species have been achieved over the past decade through detection of genetic variants in genome-wide association studies (GWAS). A majority of such variants were found located in noncoding genomic regions, and though the involvement of numerous regulatory elements (REs) has been predicted across multiple tissues in domesticated animals, their evolutionary conservation and effects on complex traits have not been fully elucidated, particularly in ruminants. Here, we systematically analyzed 137 epigenomic and transcriptomic datasets of six mammals, including cattle, sheep, goats, pigs, mice, and humans, and then integrated them with large-scale GWAS of complex traits.

RESULTS

Using 40 ChIP-seq datasets of H3K4me3 and H3K27ac, we detected 68,479, 58,562, 63,273, 97,244, 111,881, and 87,049 REs in the liver of cattle, sheep, goats, pigs, humans and mice, respectively. We then systematically characterized the dynamic functional landscapes of these REs by integrating multi-omics datasets, including gene expression, chromatin accessibility, and DNA methylation. We identified a core set (n = 6359) of ruminant-specific REs that are involved in liver development, metabolism, and immune processes. Genes with more complex cis-REs exhibited higher gene expression levels and stronger conservation across species. Furthermore, we integrated expression quantitative trait loci (eQTLs) and GWAS from 44 and 52 complex traits/diseases in cattle and humans, respectively. These results demonstrated that REs with different degrees of evolutionary conservation across species exhibited distinct enrichments for GWAS signals of complex traits.

CONCLUSIONS

We systematically annotated genome-wide functional REs in liver across six mammals and demonstrated the evolution of REs and their associations with transcriptional output and conservation. Detecting lineage-specific REs allows us to decipher the evolutionary and genetic basis of complex phenotypes in livestock and humans, which may benefit the discovery of potential biomedical models for functional variants and genes of specific human diseases.

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.

Identification of Potential Long Non-Coding RNA Candidates that Contribute to Triple-Negative Breast Cancer in Humans through Computational Approach

Breast cancer (BC) is the most frequent malignancy identified in adult females, resulting in enormous financial losses worldwide. Owing to the heterogeneity as well as various molecular subtypes, the molecular pathways underlying carcinogenesis in various forms of BC are distinct. Therefore, the advancement of alternative therapy is required to combat the ailment. Recent analyses propose that long non-coding RNAs (lncRNAs) perform an essential function in controlling immune response, and therefore, may provide essential information about the disorder. However, their function in patients with triple-negative BC (TNBC) has not been explored in detail. Here, we analyzed the changes in the genomic expression of messenger RNA (mRNA) and lncRNA in standard control in response to cancer metastasis using publicly available single-cell RNA-Seq data. We identified a total of 197 potentially novel lncRNAs in TNBC patients of which 86 were differentially upregulated and 111 were differentially downregulated. In addition, among the 909 candidate lncRNA transcripts, 19 were significantly differentially expressed (DE) of which three were upregulated and 16 were downregulated. On the other hand, 1901 mRNA transcripts were significantly DE of which 1110 were upregulated and 791 were downregulated by TNBCs subtypes. The Gene Ontology (GO) analyses showed that some of the host genes were enriched in various biological, molecular, and cellular functions. The Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis showed that some of the genes were involved in only one pathway of prostate cancer. The lncRNA-miRNA-gene network analysis showed that the lncRNAs TCONS_00076394 and TCONS_00051377 interacted with breast cancer-related micro RNAs (miRNAs) and the host genes of these lncRNAs were also functionally related to breast cancer. Thus, this study provides novel lncRNAs as potential biomarkers for the therapeutic intervention of this cancer subtype.

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.

Identification of sheep lncRNAs related to the immune response to vaccines and aluminium adjuvants

Background

Long non-coding RNAs (lncRNAs) are involved in several immune processes, including the immune response to vaccination, but most of them remain uncharacterised in livestock species. The mechanism of action of aluminium adjuvants as vaccine components is neither not fully understood.

Results

We built a transcriptome from sheep PBMCs RNA-seq data in order to identify unannotated lncRNAs and analysed their expression patterns along protein coding genes. We found 2284 novel lncRNAs and assessed their conservation in terms of sequence and synteny. Differential expression analysis performed between animals inoculated with commercial vaccines or aluminium adjuvant alone and the co-expression analysis revealed lncRNAs related to the immune response to vaccines and adjuvants. A group of co-expressed genes enriched in cytokine signalling and production highlighted the differences between different treatments. A number of differentially expressed lncRNAs were correlated with a divergently located protein-coding gene, such as the OSM cytokine. Other lncRNAs were predicted to act as sponges of miRNAs involved in immune response regulation.

Conclusions

This work enlarges the lncRNA catalogue in sheep and puts an accent on their involvement in the immune response to repetitive vaccination, providing a basis for further characterisation of the non-coding sheep transcriptome within different immune cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08086-z.

Identification of Regulatory Functions of LncRNAs Associated With T. circumcincta Infection in Adult Sheep

Several recent studies have demonstrated the role of long non-coding RNAs (lncRNAs) in regulating the defense mechanism against parasite infections, but no studies are available that investigated their relevance for immune response to nematode infection in sheep. Thus, the aim of the current study was to (i) detect putative lncRNAs that are expressed in the abomasal lymph node of adult sheep after an experimental infection with the gastrointestinal nematode (GIN) Teladorsagia circumcincta and (ii) to elucidate their potential functional role associated with the differential host immune response. We hypothesized that putative lncRNAs differentially expressed (DE) between samples from animals that differ in resistance to infection may play a significant regulatory role in response to nematode infection in adult sheep. To obtain further support for our hypothesis, we performed co-expression and functional gene enrichment analyses with the differentially expressed lncRNAs (DE lncRNAs). In a conservative approach, we included for this predictive analysis only those lncRNAs that are confirmed and supported by documentation of expression in gastrointestinal tissues in the current sheep gene atlas. We identified 9,105 putative lncRNA transcripts corresponding to 7,124 gene loci. Of these, 457 were differentially expressed lncRNA loci (DELs) with 683 lncRNA transcripts. Based on a gene co-expression analysis via weighted gene co-expression network analysis, 12 gene network modules (GNMs) were found significantly correlated with at least one of 10 selected target DE lncRNAs. Based on the principle of “guilt-by-association,” the DE genes from each of the three most significantly correlated GNMs were subjected to a gene enrichment analysis. The significant pathways associated with DE lncRNAs included ERK5 Signaling, SAPK/JNK Signaling, RhoGDI Signaling, EIF2 Signaling, Regulation of eIF4 and p70S6K Signaling and Oxidative Phosphorylation pathways. They belong to signaling pathway categories like Cellular Growth, Proliferation and Development, Cellular Stress and Injury, Intracellular and Second Messenger Signaling and Apoptosis. Overall, this lncRNA study conducted in adult sheep after GIN infection provided first insights into the potential functional role of lncRNAs in the differential host response to nematode infection.

Lantern: an integrative repository of functional annotations for lncRNAs in the human genome

BACKGROUND

With advancements in omics technologies, the range of biological processes where long non-coding RNAs (lncRNAs) are involved, is expanding extensively, thereby generating the need to develop lncRNA annotation resources. Although, there are a plethora of resources for annotating genes, despite the extensive corpus of lncRNA literature, the available resources with lncRNA ontology annotations are rare.

RESULTS

We present a lncRNA annotation extractor and repository (Lantern), developed using PubMed's abstract retrieval engine and NCBO's recommender annotation system. Lantern's annotations were benchmarked against lncRNAdb's manually curated free text. Benchmarking analysis suggested that Lantern has a recall of 0.62 against lncRNAdb for 182 lncRNAs and precision of 0.8. Additionally, we also annotated lncRNAs with multiple omics annotations, including predicted cis-regulatory TFs, interactions with RBPs, tissue-specific expression profiles, protein co-expression networks, coding potential, sub-cellular localization, and SNPs for ~ 11,000 lncRNAs in the human genome, providing a one-stop dynamic visualization platform.

CONCLUSIONS

Lantern integrates a novel, accurate semi-automatic ontology annotation engine derived annotations combined with a variety of multi-omics annotations for lncRNAs, to provide a central web resource for dissecting the functional dynamics of long non-coding RNAs and to facilitate future hypothesis-driven experiments. The annotation pipeline and a web resource with current annotations for human lncRNAs are freely available on sysbio.lab.iupui.edu/lantern.

lncRNA OGFRP1 promotes tumor progression by activating the AKT/mTOR pathway in human gastric cancer

As biomolecules of great clinical value, lncRNAs play a crucial role as regulators in the processes of tumor origin, metastasis, and recurrence. Thus, lncRNAs are urgently needed for research in gastric cancer. We elucidated the specific function of OGFRP1, both in vitro and in vivo. OGFRP1 was expressed at abnormally high levels in gastric cancer samples (n = 408) compared to normal samples (n = 211). Similar results were obtained in 30 clinical case samples. Interference of OGFRP1 markedly blocked cell proliferation and migration, and it induced cell cycle arrest and the apoptosis of gastric cancer cells in vitro. Phosphorylation of AKT was inhibited in cells transfected with OGFRP1 siRNA, as compared to their control cells. The in vivo results further confirmed the antitumor effects of OGFRP1 knockdown on gastric cancer. Decreases in tumor volume (104.23±62.27 mm³) and weight (0.1006±0.0488 g) in nude mice were observed during the OGFRP1 interference, as compared with the control group (418.96±211.96 mm³ and 0.2741±0.0769 g). OGFRP1 promotes tumor progression through activating the AKT/mTOR pathway. Our findings provide a new potential target for the clinical treatment of human gastric cancer.

Long non-coding RNA signatures as predictors of prognosis in thyroid cancer: a narrative review

Thyroid cancer (TC) is the most common endocrine malignancy, with high incidence rates in recent decades. Most TC cases have good prognoses, but a high risk of recurrence and metastases poses challenges, especially for patients with high-risk factors. Currently used prognostic markers for TC involve a combination of genetic factors and overexpressed proteins. Long non-coding RNAs (lncRNAs) regulate several integral biologic processes by playing key roles in the transcription of several downstream targets maintaining cellular behavior. Prior studies have revealed that lncRNAs promote tumor cell proliferation, invasion, metastasis, and angiogenesis, making them important targets for therapeutic intervention in cancer. While the exact molecular mechanisms underlying the role of lncRNAs in modulating TC progression and recurrence is still unclear, it is important to note that some lncRNAs are upregulated in certain cancers, while others are downregulated. In the present study, we review several key lncRNAs, their association with cancer progression, and the important roles they may play as tumor suppressors or tumor promoters in tumorigenesis. We discuss the potential mechanisms of lncRNA-mediated pathogenesis that can be targeted for the treatment of TC, the existing and potential benefits of using lncRNAs as diagnostic and prognostic measures for cancer detection, and tumor burden in patients.

Whole transcriptome approach to evaluate the effect of aluminium hydroxide in ovine encephalon

Aluminium hydroxide adjuvants are crucial for livestock and human vaccines. Few studies have analysed their effect on the central nervous system in vivo. In this work, lambs received three different treatments of parallel subcutaneous inoculations during 16 months with aluminium-containing commercial vaccines, an equivalent dose of aluminium hydroxide or mock injections. Brain samples were sequenced by RNA-seq and miRNA-seq for the expression analysis of mRNAs, long non-coding RNAs and microRNAs and three expression comparisons were made. Although few differentially expressed genes were identified, some dysregulated genes by aluminium hydroxide alone were linked to neurological functions, the lncRNA TUNA among them, or were enriched in mitochondrial energy metabolism related functions. In the same way, the miRNA expression was mainly disrupted by the adjuvant alone treatment. Some differentially expressed miRNAs had been previously linked to neurological diseases, oxidative stress and apoptosis. In brief, in this study aluminium hydroxide alone altered the transcriptome of the encephalon to a higher degree than commercial vaccines that present a milder effect. The expression changes in the animals inoculated with aluminium hydroxide suggest mitochondrial disfunction. Further research is needed to elucidate to which extent these changes could have pathological consequences.

Single-molecule real-time sequencing of the full-length transcriptome of loquat under low-temperature stress

In this study, third-generation full-length (FL) transcriptome sequencing was performed of loquat using single-molecule real-time(SMRT) sequencing from the pooled cDNA of embryos of young loquat fruit under different low temperatures (three biological replicates for treatments of 1°C, -1°C, and -3°C, for 12 h or 24 h) and the control group(three biological replicates for treatments of room temperature), Illumina sequencing was used to correct FL transcriptome sequences. A total of 3 PacBio Iso-Seq libraries (1–2 kb, 2–3 kb and 3–6 kb) and 21 Illumina transcriptome libraries were constructed, a total of 13.41 Gb of clean reads were generated, which included 215,636 reads of insert (ROIs) and 121,654 FL, non-chimaric (FLNC) reads. Transcript clustering analysis of the FLNC reads revealed 76,586 consensus isoforms, and a total of 12,520 high-quality transcript sequences corrected with non-FL sequences were used for subsequent analysis. After the redundant reads were removed, 38,435 transcripts were obtained. A total of 27,905 coding DNA sequences (CDSs) were identified, and 407 long non-coding RNAs (lncRNAs) were ultimately predicted. Additionally, 24,832 simple sequence repeats (SSRs) were identified, and a total of 1,295 alternative splicing (AS) events were predicted. Furthermore, 37,993 transcripts were annotated in eight functional databases. This is the first study to perform SMRT sequencing of the FL transcriptome of loquat. The obtained transcriptomic data are conducive for further exploration of the mechanism of loquat freezing injury and thus serve as an important theoretical basis for generating new loquat material and for identifying new ways to improve loquat cold resistance.

Species-Specificity of Transcriptional Regulation and the Response to Lipopolysaccharide in Mammalian Macrophages

Mammalian macrophages differ in their basal gene expression profiles and response to the toll-like receptor 4 (TLR4) agonist, lipopolysaccharide (LPS). In human macrophages, LPS elicits a temporal cascade of transient gene expression including feed forward activators and feedback regulators that limit the response. Here we present a transcriptional network analysis of the response of sheep bone marrow-derived macrophages (BMDM) to LPS based upon RNA-seq at 0, 2, 4, 7, and 24 h post-stimulation. The analysis reveals a conserved transcription factor network with humans, and rapid induction of feedback regulators that constrain the response at every level. The gene expression profiles of sheep BMDM at 0 and 7 h post LPS addition were compared to similar data obtained from goat, cow, water buffalo, horse, pig, mouse and rat BMDM. This comparison was based upon identification of 8,200 genes annotated in all species and detected at >10TPM in at least one sample. Analysis of expression of transcription factors revealed a conserved transcriptional millieu associated with macrophage differentiation and LPS response. The largest co-expression clusters, including genes encoding cell surface receptors, endosome–lysosome components and secretory activity, were also expressed in all species and the combined dataset defines a macrophage functional transcriptome. All of the large animals differed from rodents in lacking inducible expression of genes involved in arginine metabolism and nitric oxide production. Instead, they expressed inducible transporters and enzymes of tryptophan and kynurenine metabolism. BMDM from all species expressed high levels of transcripts encoding transporters and enzymes involved in glutamine metabolism suggesting that glutamine is a major metabolic fuel. We identify and discuss transcripts that were uniquely expressed or regulated in rodents compared to large animals including ACOD1, CXC and CC chemokines, CD163, CLEC4E, CPM, CSF1, CSF2, CTSK, MARCO, MMP9, SLC2A3, SLC7A7, and SUCNR1. Conversely, the data confirm the conserved regulation of multiple transcripts for which there is limited functional data from mouse models and knockouts. The data provide a resource for functional annotation and interpretation of loci involved in susceptibility to infectious and inflammatory disease in humans and large animal species.

Pan-tissue transcriptome analysis of long noncoding RNAs in the American beaver Castor canadensis

BACKGROUND

Long noncoding RNAs (lncRNAs) have roles in gene regulation, epigenetics, and molecular scaffolding and it is hypothesized that they underlie some mammalian evolutionary adaptations. However, for many mammalian species, the absence of a genome assembly precludes the comprehensive identification of lncRNAs. The genome of the American beaver (Castor canadensis) has recently been sequenced, setting the stage for the systematic identification of beaver lncRNAs and the characterization of their expression in various tissues. The objective of this study was to discover and profile polyadenylated lncRNAs in the beaver using high-throughput short-read sequencing of RNA from sixteen beaver tissues and to annotate the resulting lncRNAs based on their potential for orthology with known lncRNAs in other species.

RESULTS

Using de novo transcriptome assembly, we found 9528 potential lncRNA contigs and 187 high-confidence lncRNA contigs. Of the high-confidence lncRNA contigs, 147 have no known orthologs (and thus are putative novel lncRNAs) and 40 have mammalian orthologs. The novel lncRNAs mapped to the Oregon State University (OSU) reference beaver genome with greater than 90% sequence identity. While the novel lncRNAs were on average shorter than their annotated counterparts, they were similar to the annotated lncRNAs in terms of the relationships between contig length and minimum free energy (MFE) and between coverage and contig length. We identified beaver orthologs of known lncRNAs such as XIST, MEG3, TINCR, and NIPBL-DT. We profiled the expression of the 187 high-confidence lncRNAs across 16 beaver tissues (whole blood, brain, lung, liver, heart, stomach, intestine, skeletal muscle, kidney, spleen, ovary, placenta, castor gland, tail, toe-webbing, and tongue) and identified both tissue-specific and ubiquitous lncRNAs.

CONCLUSIONS

To our knowledge this is the first report of systematic identification of lncRNAs and their expression atlas in beaver. LncRNAs-both novel and those with known orthologs-are expressed in each of the beaver tissues that we analyzed. For some beaver lncRNAs with known orthologs, the tissue-specific expression patterns were phylogenetically conserved. The lncRNA sequence data files and raw sequence files are available via the web supplement and the NCBI Sequence Read Archive, respectively.

A Mini-Atlas of Gene Expression for the Domestic Goat (Capra hircus)

Goats (Capra hircus) are an economically important livestock species providing meat and milk across the globe. They are of particular importance in tropical agri-systems contributing to sustainable agriculture, alleviation of poverty, social cohesion, and utilisation of marginal grazing. There are excellent genetic and genomic resources available for goats, including a highly contiguous reference genome (ARS1). However, gene expression information is limited in comparison to other ruminants. To support functional annotation of the genome and comparative transcriptomics, we created a mini-atlas of gene expression for the domestic goat. RNA-Seq analysis of 17 transcriptionally rich tissues and 3 cell-types detected the majority (90%) of predicted protein-coding transcripts and assigned informative gene names to more than 1000 previously unannotated protein-coding genes in the current reference genome for goat (ARS1). Using network-based cluster analysis, we grouped genes according to their expression patterns and assigned those groups of coexpressed genes to specific cell populations or pathways. We describe clusters of genes expressed in the gastro-intestinal tract and provide the expression profiles across tissues of a subset of genes associated with functional traits. Comparative analysis of the goat atlas with the larger sheep gene expression atlas dataset revealed transcriptional similarities between macrophage associated signatures in the sheep and goats sampled in this study. The goat transcriptomic resource complements the large gene expression dataset we have generated for sheep and contributes to the available genomic resources for interpretation of the relationship between genotype and phenotype in small ruminants.

Microsyntenic Clusters Reveal Conservation of lncRNAs in Chordates Despite Absence of Sequence Conservation

Homologous long non-coding RNAs (lncRNAs) are elusive to identify by sequence similarity due to their fast-evolutionary rate. Here we develop LincOFinder, a pipeline that finds conserved intergenic lncRNAs (lincRNAs) between distant related species by means of microsynteny analyses. Using this tool, we have identified 16 bona fide homologous lincRNAs between the amphioxus and human genomes. We characterized and compared in amphioxus and Xenopus the expression domain of one of them, Hotairm1, located in the anterior part of the Hox cluster. In addition, we analyzed the function of this lincRNA in Xenopus, showing that its disruption produces a severe headless phenotype, most probably by interfering with the regulation of the Hox cluster. Our results strongly suggest that this lincRNA has probably been regulating the Hox cluster since the early origin of chordates. Our work pioneers the use of syntenic searches to identify non-coding genes over long evolutionary distances and helps to further understand lncRNA evolution.

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.

Screening and evaluating of long non-coding RNAs in prenatal and postnatal pituitary gland of sheep

Long non-coding RNAs are transcribed into RNA molecules that are >200 nucleotides in length. However, the expression and function analysis of lncRNAs in the sheep pituitary gland are still lacking. In this study, we identified 1755 lncRNAs (545 annotated lncRNAs and 1210 novel lncRNAs) from RNA-seq data in the pituitary gland of embryonic and adult sheep. A total of 235 lncRNAs were differentially expressed between embryonic and adult group. We verified the presence of some lncRNAs using RT-PCR and DNA sequencing, and identified some differentially expressed lncRNAs using qPCR. We also investigated the role of cis-acting lncRNAs on target genes. GO and KEGG enrichment analysis revealed that the target genes of lncRNAs were involved in the regulation of hormones secretion and some signaling pathways in the sheep pituitary gland. Our study provides comprehensive expression profiles of lncRNAs and valuable resource for understanding their function in the pituitary gland.

Genome-wide identification and characterization of long non-coding RNAs expressed during sheep fetal and postnatal hair follicle development

Long non-coding RNAs (lncRNAs), >200 nt in length, are transcribed from mammalian genomes. They play important regulatory roles in various biological processes; However, the function and expression profile of lncRNAs involved in the development of hair follicles in the fetus, have been relatively under-explored area. To investigate the specific role of lncRNAs and mRNAs that regulate hair follicle development, we herein performed a comprehensive study on the lncRNA and mRNA expression profiles of sheep at multiple embryonic days (E65, E85, E105, and E135) and six lambs aged one week (D7) and one month (D30) using RNA-seq technology. The number of genes (471 lncRNAs and 12,812 mRNAs) differentially expressed and potential targets of differentially expressed lncRNAs were predicted. Differentially expressed lncRNAs were grouped into 10 clusters based on their expression pattern by K-means clustering. Moreover, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that some differentially expressed mRNAs, such as DKK1, DSG4, FOXE1, Hoxc13, SFRP1, SFRP2, and Wnt10A overlapped with lncRNAs targets, and enriched in important hair follicle developmental pathways, including Wnt, TNF, and MAPK signaling pathways. In addition, 9 differentially expressed lncRNAs and 4 differentially expressed mRNAs were validated using quantitative real-time PCR (qRT-PCR). This study helps enrich the Ovis lncRNA databases and provides a comprehensive lncRNA transcriptome profile of fetal and postnatal skin of sheep. Additionally, it provides a foundation for further experiments on the role of lncRNAs in the regulation of hair growth in sheep.

DNA methylation and functional characterization of the XIST gene during in vitro early embryo development in cattle

XIST, in association with the shorter ncRNA RepA, are essential for the initiation of X chromosome inactivation (XCI) in mice. The molecular mechanisms controlling XIST and RepA expression are well characterized in that specie. However, little is known in livestock. We aimed to characterize the DNA methylation status along the 5' portion of XIST and to characterize its transcriptional profile during early development in cattle. Three genomic regions of XIST named here as promoter, RepA and DMR1 had their DNA methylation status characterized in gametes and embryos. Expression profile of XIST was evaluated, including sense and antisense transcription. Oocytes showed higher levels of methylation than spermatozoa that was demethylated. DMR1 was hypermethylated throughout oogenesis. At the 8-16-cell embryo stage DMR1 was completed demethylated. Interestingly, RepA gain methylation during oocyte maturation and was demethylated at the blastocyst stage, later than DMR1. These results suggest that DMR1 and RepA are transient differentially methylated regions in cattle. XIST RNA was detected in matured oocytes and in single cells from the 2-cell to the morula stage, confirming the presence of maternal and embryonic transcripts. Sense and antisense transcripts were detected along the XIST in blastocyst. In silico analysis identified 63 novel transcript candidates at bovine XIST locus from both the plus and minus strands. Taking together these results improve our understanding of the molecular mechanisms involved in XCI initiation in cattle. This information may be useful for the improvement of assisted reproductive technologies in livestock considering that in vitro conditions may impair epigenetic reprogramming.

Identification and Expression Analysis of Long Noncoding RNAs in Fat-Tail of Sheep Breeds

Emerging evidence suggests that long non-coding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. However, most studies have been focused on a few established model organisms and little is known about lncRNAs in fat-tail development in sheep. Here, the first profile of lncRNA in sheep fat-tail along with their possible roles in fat deposition were investigated, based on a comparative transcriptome analysis between fat-tailed (Lori-Bakhtiari) and thin-tailed (Zel) Iranian sheep breeds. Among all identified lncRNAs candidates, 358 and 66 transcripts were considered novel intergenic (lincRNAs) and novel intronic (ilncRNAs) corresponding to 302 and 58 gene loci, respectively. Our results indicated that a low percentage of the novel lncRNAs were conserved. Also, synteny analysis identified 168 novel lincRNAs with the same syntenic region in human, bovine and chicken. Only seven lncRNAs were identified as differentially expressed genes between fat and thin tailed breeds. Q-RT-PCR results were consistent with the RNA-Seq data and validated the findings. Target prediction analysis revealed that the novel lncRNAs may act in cis or trans and regulate the expression of genes that are involved in the lipid metabolism. A gene regulatory network including lncRNA-mRNA interactions were constructed and three significant modules were found, with genes relevant to lipid metabolism, insulin and calcium signaling pathway. Moreover, integrated analysis with AnimalQTLdb database further suggested six lincRNAs and one ilncRNAs as candidates of sheep fat-tail development. Our results highlighted the putative contributions of lncRNAs in regulating expression of genes associated with fat-tail development in sheep.

Comprehensive Transcriptional Profiling of the Gastrointestinal Tract of Ruminants from Birth to Adulthood Reveals Strong Developmental Stage Specific Gene Expression

One of the most significant physiological challenges to neonatal and juvenile ruminants is the development and establishment of the rumen. Using a subset of RNA-Seq data from our high-resolution atlas of gene expression in sheep (Ovis aries) we have provided the first comprehensive characterization of transcription of the entire gastrointestinal (GI) tract during the transition from pre-ruminant to ruminant. The dataset comprises 164 tissue samples from sheep at four different time points (birth, one week, 8 weeks and adult). Using network cluster analysis we illustrate how the complexity of the GI tract is reflected in tissue- and developmental stage-specific differences in gene expression. The most significant transcriptional differences between neonatal and adult sheep were observed in the rumen complex. Comparative analysis of gene expression in three GI tract tissues from age-matched sheep and goats revealed species-specific differences in genes involved in immunity and metabolism. This study improves our understanding of the transcriptomic mechanisms involved in the transition from pre-ruminant to ruminant by identifying key genes involved in immunity, microbe recognition and metabolism. The results form a basis for future studies linking gene expression with microbial colonization of the developing GI tract and provide a foundation to improve ruminant efficiency and productivity through identifying potential targets for novel therapeutics and gene editing.

Analysis of long non-coding RNA and mRNA expression in bovine macrophages brings up novel aspects of Mycobacterium avium subspecies paratuberculosis infections

Paratuberculosis is a major disease in cattle that severely affects animal welfare and causes huge economic losses worldwide. Development of alternative diagnostic methods is of urgent need to control the disease. Recent studies suggest that long non-coding RNAs (lncRNAs) play a crucial role in regulating immune function and may confer valuable information about the disease. However, their role has not yet been investigated in cattle with respect to infection towards Paratuberculosis. Therefore, we investigated the alteration in genomic expression profiles of mRNA and lncRNA in bovine macrophages in response to Paratuberculosis infection using RNA-Seq. We identified 397 potentially novel lncRNA candidates in macrophages of which 38 were differentially regulated by the infection. A total of 820 coding genes were also significantly altered by the infection. Co-expression analysis of lncRNAs and their neighbouring coding genes suggest regulatory functions of lncRNAs in pathways related to immune response. For example, this included protein coding genes such as TNIP3, TNFAIP3 and NF-κB2 that play a role in NF-κB2 signalling, a pathway associated with immune response. This study advances our understanding of lncRNA roles during Paratuberculosis infection.

The Structure, Evolution, and Gene Expression Within the Caprine Leukocyte Receptor Complex

The leukocyte receptor complex (LRC) encodes a large number of immunoglobulin (Ig)-like receptors involved in the immune response, particularly in modulating natural killer (NK) cell function. The killer cell Ig-like receptors (KIR), the leukocyte Ig-like receptors (LILR), and a recently described novel Ig-like receptor family are highly variable between species, which is consistent with rapid evolution driven by selection pressure from pathogens. Among the species studied to date, only simians (such as humans) and bovids (such as cattle and goats) have an expanded complement of KIR genes and represent an interesting model to study KIR evolution. Using recently improved genome assemblies and an assembly of bacterial artificial chromosomes, we describe the structure of the LRC, and the KIR region in particular, in goats and compare this to sheep as the assemblies allow. These species diverged from a common ancestor ~10 million years ago and from cattle ~25 million years ago. We identified conserved KIR genes common to both goats and sheep and confirm a partial sheep haplotype shared between the Rambouillet and Texel breeds. Goats and sheep have independently expanded two novel KIR subgroups, and unlike cattle or any other mammal, they do not appear to possess a functional 3DL-lineage KIR gene. Investigation of LRC gene expression using available transcriptomic data for various sheep and goat tissues largely confirmed putative gene annotation and revealed that a relatively conserved caprinae-specific KIR subgroup is expressed in macrophages. The LILR and novel Ig-like receptors were also highly expressed across a diverse range of tissues. This further step toward our understanding of the LRC receptor repertoire will help inform future studies investigating immune response variation in these species.

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.