Malat1 affects transcription and splicing through distinct pathways in mouse embryonic stem cells.

Malat1 is a long-noncoding RNA with critical roles in gene regulation and cancer metastasis, however its functional role in stem cells is largely unexplored. We here perform a nuclear knockdown of Malat1 in mouse embryonic stem cells, causing the de-regulation of 320 genes and aberrant splicing of 90 transcripts, some of which potentially affecting the translated protein sequence. We find evidence that Malat1 directly interacts with gene bodies and aberrantly spliced transcripts, and that it locates upstream of down-regulated genes at their putative enhancer regions, in agreement with functional genomics data. Consistent with this, we find these genes affected at both exon and intron levels, suggesting that they are transcriptionally regulated by Malat1. Besides, the down-regulated genes are regulated by specific transcription factors and bear both activating and repressive chromatin marks, suggesting that some of them might be regulated by bivalent promoters. We propose a model in which Malat1 facilitates the transcription of genes involved in chromatid dynamics and mitosis in one pathway, and affects the splicing of transcripts that are themselves involved in RNA processing in a distinct pathway. Lastly, we compare our findings with Malat1 perturbation studies performed in other cell systems and in vivo.

The landscape of long noncoding RNA expression in the goat brain.

The brain regulates multiple metabolic processes, such as food intake, energy expenditure, insulin secretion, hepatic glucose production, and glucose and fatty acid metabolism in adipose tissue, which are fundamental for the maintenance of energy and glucose homeostasis during lactation and pregnancy. In addition, brain expression has a fundamental impact on the development of maternal behavior. Although brain functions are partly regulated by long noncoding RNAs (lncRNAs), their expression profiles have not been characterized in depth in any ruminant species. We have sequenced the transcriptome of 12 brain tissues from 3 goats that were 1 mo pregnant and 4 nonpregnant goats to investigate their lncRNA expression patterns. Between 4,363 (adenohypophysis) and 4,604 (olfactory bulb) lncRNAs were expressed in brain tissues, leading us to establish a set of 794 already annotated lncRNAs and 5,098 novel lncRNA candidates. The detected lncRNAs shared features with those of other mammals, and tissue-specific lncRNAs were enriched in brain development-related terms. Differential expression analyses between goats that were 1 mo pregnant and nonpregnant goats showed that the lncRNA expression profiles of certain brain regions experience substantial changes associated with early pregnancy (238 lncRNAs are differentially expressed in the olfactory bulb), but others do not. Enrichment analysis showed that differentially expressed lncRNAs from the olfactory bulb are co-expressed with genes previously linked to behavioral changes related to pregnancy. These findings provide a first characterization of the landscape of lncRNA expression in the goat brain and provides valuable clues to understand the molecular events triggered by early pregnancy in the central nervous system.

Detection of transcriptome-wide microRNA-target interactions in single cells with agoTRIBE.

MicroRNAs (miRNAs) exert their gene regulatory effects on numerous biological processes based on their selection of target transcripts. Current experimental methods available to identify miRNA targets are laborious and require millions of cells. Here we have overcome these limitations by fusing the miRNA effector protein Argonaute2 to the RNA editing domain of ADAR2, allowing the detection of miRNA targets transcriptome-wide in single cells. miRNAs guide the fusion protein to their natural target transcripts, causing them to undergo A>I editing, which can be detected by sensitive single-cell RNA sequencing. We show that agoTRIBE identifies functional miRNA targets, which are supported by evolutionary sequence conservation. In one application of the method we study microRNA interactions in single cells and identify substantial differential targeting across the cell cycle. AgoTRIBE also provides transcriptome-wide measurements of RNA abundance and allows the deconvolution of miRNA targeting in complex tissues at the single-cell level.

Exon-intron split analysis reveals posttranscriptional regulatory signals induced by high and low n-6/n-3 polyunsaturated fatty acid ratio diets in piglets.

Polyunsaturated fatty acids (PUFA), such as omega-6 (n-6) and omega-3 (n-3), play a vital role in nutrient metabolism, inflammatory response, and gene regulation. microRNAs (miRNA), which can potentially degrade targeted messenger RNAs (mRNA) and/or inhibit their translation, might play a relevant role in PUFA-related changes in gene expression. Although differential expression analyses can provide a comprehensive picture of gene expression variation, they are unable to disentangle when in the mRNA life cycle the regulation of expression is taking place, including any putative functional miRNA-driven repression. To capture this, we used an exon-intron split analysis (EISA) approach to account for posttranscriptional changes in response to extreme values of n-6/n-3 PUFA ratio. Longissimus dorsi muscle samples of male and female piglets from sows fed with n-6/n-3 PUFA ratio of 13:1 (SOY) or 4:1 (LIN), were analyzed in a bidirectional contrast (LIN vs. SOY, SOY vs. LIN). Our results allowed the identification of genes showing strong posttranscriptional downregulation signals putatively targeted by significantly upregulated miRNA. Moreover, we identified genes primarily involved in the regulation of lipid-related metabolism and immune response, which may be associated with the pro- and anti-inflammatory functions of the n-6 and n-3 PUFA, respectively. EISA allowed us to uncover regulatory networks complementing canonical differential expression analyses, thus providing a more comprehensive view of muscle metabolic changes in response to PUFA concentration.

The relationship between dietary lipids, such as omega-6 and omega-3 polyunsaturated fatty acids (PUFA), and gene expression regulation was explored in piglet muscle. While these PUFA can influence nutrient metabolism and inflammatory response, small regulatory molecules called microRNAs (miRNA) can also influence the activity of genes. In this experiment, we used a computational approach dubbed exon­intron split analysis (EISA) to fully understand the role of miRNA in this context and how the genes and miRNA respond to changes in PUFA levels. Our findings demonstrated that some genes involved in lipid metabolism and immune response were affected by different PUFA concentrations and that EISA provides a more comprehensive view of how genes are regulated throughout their life cycle.

Historical RNA expression profiles from the extinct Tasmanian tiger.

Paleogenomics continues to yield valuable insights into the evolution, population dynamics, and ecology of our ancestors and other extinct species. However, DNA sequencing cannot reveal tissue-specific gene expression, cellular identity, or gene regulation, which are only attainable at the transcriptional level. Pioneering studies have shown that useful RNA can be extracted from ancient specimens preserved in permafrost and historical skins from extant canids, but no attempts have been made so far on extinct species. We extract, sequence, and analyze historical RNA from muscle and skin tissue of a ∼130-year-old Tasmanian tiger (Thylacinus cynocephalus) preserved in desiccation at room temperature in a museum collection. The transcriptional profiles closely resemble those of extant species, revealing specific anatomical features such as slow muscle fibers or blood infiltration. Metatranscriptomic analysis, RNA damage, tissue-specific RNA profiles, and expression hotspots genome-wide further confirm the thylacine origin of the sequences. RNA sequences are used to improve protein-coding and noncoding annotations, evidencing missing exonic loci and the location of ribosomal RNA genes while increasing the number of annotated thylacine microRNAs from 62 to 325. We discover a thylacine-specific microRNA isoform that could not have been confirmed without RNA evidence. Finally, we detect traces of RNA viruses, suggesting the possibility of profiling viral evolution. Our results represent the first successful attempt to obtain transcriptional profiles from an extinct animal species, providing thought-to-be-lost information on gene expression dynamics. These findings hold promising implications for the study of RNA molecules across the vast collections of natural history museums and from well-preserved permafrost remains.

A protein-coding gene expression atlas from the brain of pregnant and non-pregnant goats.

Background: The brain is an extraordinarily complex organ with multiple anatomical structures involved in highly specialized functions related with behavior and physiological homeostasis. Our goal was to build an atlas of protein-coding gene expression in the goat brain by sequencing the transcriptomes of 12 brain regions in seven female Murciano-Granadina goats, from which three of them were 1-month pregnant. Results: Between 14,889 (cerebellar hemisphere) and 15,592 (pineal gland) protein-coding genes were expressed in goat brain regions, and most of them displayed ubiquitous or broad patterns of expression across tissues. Principal component analysis and hierarchical clustering based on the patterns of mRNA expression revealed that samples from certain brain regions tend to group according to their position in the anterior-posterior axis of the neural tube, i.e., hindbrain (pons and medulla oblongata), midbrain (rostral colliculus) and forebrain (frontal neocortex, olfactory bulb, hypothalamus, and hippocampus). Exceptions to this observation were cerebellum and glandular tissues (pineal gland and hypophysis), which showed highly divergent mRNA expression profiles. Differential expression analysis between pregnant and non-pregnant goats revealed moderate changes of mRNA expression in the frontal neocortex, hippocampus, adenohypophysis and pons, and very dramatic changes in the olfactory bulb. Many genes showing differential expression in this organ are related to olfactory function and behavior in humans. Conclusion: With the exception of cerebellum and glandular tissues, there is a relationship between the cellular origin of sampled regions along the anterior-posterior axis of the neural tube and their mRNA expression patterns in the goat adult brain. Gestation induces substantial changes in the mRNA expression of the olfactory bulb, a finding consistent with the key role of this anatomical structure on the development of maternal behavior.

Identifying miRNA-mRNA regulatory networks on extreme n-6/n-3 polyunsaturated fatty acid ratio expression profiles in porcine skeletal muscle.

Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acids (PUFAs) are essential fatty acids with antagonistic inflammatory functions that play vital roles in metabolic health and immune response. Current commercial swine diets tend to over-supplement with n-6 PUFAs, which may increase the likelihood of developing inflammatory diseases and affect the overall well-being of the animals. However, it is still poorly understood how n-6/n-3 PUFA ratios affect the porcine transcriptome expression and how messenger RNAs (mRNAs) and microRNAs (miRNAs) might regulate biological processes related to PUFA metabolism. On account of this, we selected a total of 20 Iberian × Duroc crossbred pigs with extreme values for n-6/n-3 FA ratio (10 high vs 10 low), and longissimus dorsi muscle samples were used to identify differentially expressed mRNAs and miRNAs. The observed differentially expressed mRNAs were associated to biological pathways related to muscle growth and immunomodulation, while the differentially expressed microRNAs (ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-15b and ssc-miR-7142-3p) were correlated to adipogenesis and immunity. Relevant miRNA-to-mRNA regulatory networks were also predicted (i.e., mir15b to ARRDC3; mir-7142-3p to METTL21C), and linked to lipolysis, obesity, myogenesis, and protein degradation. The n-6/n-3 PUFA ratio differences in pig skeletal muscle revealed genes, miRNAs and enriched pathways involved in lipid metabolism, cell proliferation and inflammation.

MicroRNA profiling of cerebrospinal fluid from dogs with steroid responsive meningitis-arteritis and meningoencephalitis of unknown origin.

Introduction: Non-infectious inflammatory diseases of the central nervous system in dogs, such as steroid responsive meningitis-arteritis (SRMA) and meningoencephalitis of unknown origin (MUO), represent a common clinical challenge that needs extensive and multimodal work-up to reach a presumptive diagnosis. Both diseases are presumably caused by dysregulations of the immune system, but further research is needed in order to understand the molecular mechanisms behind each disease and to optimize treatment. Methods: By next-generation sequencing and subsequent quantitative real-time PCR (qPCR) verification, we designed a prospective case-control pilot study to analyze the small RNA profiles of cerebrospinal fluid from dogs suffering from MUO (N = 5), dogs suffering from SRMA (N = 8), and healthy dogs (N = 5) presented for elective euthanasia used as the Control group. Results: Our results showed an overall enrichment in Y-RNA fragments across all samples, followed by microRNAs (miRNAs) and ribosomal RNAs as the major findings. Additional traces of short RNA reads mapped to long non-coding RNAs and protein-coding genes were also found. From the detected canine miRNAs, miR-21, miR-486, miR-148a, miR-99a, miR-191 and miR-92a were among the most abundant. Dogs with SRMA showed higher differences in miRNA abundance than dogs with MUO when compared to healthy dogs, and miR-142-3p was consistently detected as differentially upregulated in both diseases, although at a low concentration. Moreover, miR-405-5p and miR-503-5p showed different profiles between SRMA and MUO dogs. Subsequent qPCR analyses confirmed miR-142-5p, miR-191-5p and miR-92a-3p as significantly upregulated miRNAs in dogs with SRMA and/or MUO. Discussion: Cerebrospinal fluid is a challenging biological material to use for profiling miRNAs due to the low content of circulating RNAs. Despite this, we could confirm several miRNAs being differentially abundant when comparing healthy dogs and dogs with MUO and SRMA, respectively. The results of this study indicate a potential role of miRNAs in the underlying molecular mechanisms of these diseases and establish the basis for further studies.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs.

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Urinary microRNAome in healthy cats and cats with pyelonephritis or other urological conditions.

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression at the post-transcriptional level. miRNAs have been found in urine and have shown diagnostic potential in human nephropathies. Here, we aimed to characterize, for the first time, the feline urinary miRNAome and explore the use of urinary miRNA profiles as non-invasive biomarkers for feline pyelonephritis (PN). Thirty-eight cats were included in a prospective case-control study and classified in five groups: healthy Control cats (n = 11), cats with PN (n = 10), cats with subclinical bacteriuria or cystitis (SB/C, n = 5), cats with ureteral obstruction (n = 7) and cats with chronic kidney disease (n = 5). By small RNA sequencing we identified 212 miRNAs in cat urine, including annotated (n = 137) and putative novel (n = 75) miRNAs. The 15 most highly abundant urinary miRNAs accounted for nearly 71% of all detected miRNAs, most of which were previously identified in feline kidney. Ninety-nine differentially abundant (DA) miRNAs were identified when comparing Control cats to cats with urological conditions and 102 DA miRNAs when comparing PN to other urological conditions. Tissue clustering analysis revealed that the majority of urine samples clustered close to kidney, which confirm the likely cellular origin of the secreted urinary miRNAs. Relevant DA miRNAs were verified by quantitative real-time PCR (qPCR). Eighteen miRNAs discriminated Control cats from cats with a urological condition. Of those, seven miRNAs were DA by both RNAseq and qPCR methods between Control and PN cats (miR-125b-5p, miR-27a-3p, miR-21-5p, miR-27b-3p, miR-125a-5p, miR-17-5p and miR-23a-3p) or DA between Control and SB/C cats (miR-125b-5p). Six additional miRNAs (miR-30b-5p, miR-30c, miR-30e-5p, miR-27a-3p, miR-27b-39 and miR-222) relevant for discriminating PN from other urological conditions were identified by qPCR alone (n = 4) or by both methods (n = 2) (P<0.05). This panel of 13 miRNAs has potential as non-invasive urinary biomarkers for diagnostic of PN and other urological conditions in cats.

MirGeneDB 2.1: toward a complete sampling of all major animal phyla.

We describe an update of MirGeneDB, the manually curated microRNA gene database. Adhering to uniform and consistent criteria for microRNA annotation and nomenclature, we substantially expanded MirGeneDB with 30 additional species representing previously missing metazoan phyla such as sponges, jellyfish, rotifers and flatworms. MirGeneDB 2.1 now consists of 75 species spanning over ∼800 million years of animal evolution, and contains a total number of 16 670 microRNAs from 1549 families. Over 6000 microRNAs were added in this update using ∼550 datasets with ∼7.5 billion sequencing reads. By adding new phylogenetically important species, especially those relevant for the study of whole genome duplication events, and through updating evolutionary nodes of origin for many families and genes, we were able to substantially refine our nomenclature system. All changes are traceable in the specifically developed MirGeneDB version tracker. The performance of read-pages is improved and microRNA expression matrices for all tissues and species are now also downloadable. Altogether, this update represents a significant step toward a complete sampling of all major metazoan phyla, and a widely needed foundation for comparative microRNA genomics and transcriptomics studies. MirGeneDB 2.1 is part of RNAcentral and Elixir Norway, publicly and freely available at http://www.mirgenedb.org/.

Assessing the levels of intraspecific admixture and interspecific hybridization in Iberian wild goats (Capra pyrenaica).

Iberian wild goats (Capra pyrenaica, also known as Iberian ibex, Spanish ibex, and Spanish wild goat) underwent strong genetic bottlenecks during the 19th and 20th centuries due to overhunting and habitat destruction. From the 1970s to 1990s, augmentation translocations were frequently carried out to restock Iberian wild goat populations (very often with hunting purposes), but they were not systematically planned or recorded. On the other hand, recent data suggest the occurrence of hybridization events between Iberian wild goats and domestic goats (Capra hircus). Augmentation translocations and interspecific hybridization might have contributed to increase the diversity of Iberian wild goats. With the aim of investigating this issue, we have genotyped 118 Iberian wild goats from Tortosa-Beceite, Sierra Nevada, Muela de Cortes, Gredos, Batuecas, and Ordesa and Monte Perdido by using the Goat SNP50 BeadChip (Illumina). The analysis of genotypic data indicated that Iberian wild goat populations are strongly differentiated and display low diversity. Only three Iberian wild goats out from 118 show genomic signatures of mixed ancestry, a result consistent with a scenario in which past augmentation translocations have had a limited impact on the diversity of Iberian wild goats. Besides, we have detected eight Iberian wild goats from Tortosa-Beceite with signs of domestic goat introgression. Although rare, hybridization with domestic goats could become a potential threat to the genetic integrity of Iberian wild goats; hence, measures should be taken to avoid the presence of uncontrolled herds of domestic or feral goats in mountainous areas inhabited by this iconic wild ungulate.

Integrating multi-taxon palaeogenomes and sedimentary ancient DNA to study past ecosystem dynamics.

Ancient DNA (aDNA) has played a major role in our understanding of the past. Important advances in the sequencing and analysis of aDNA from a range of organisms have enabled a detailed understanding of processes such as past demography, introgression, domestication, adaptation and speciation. However, to date and with the notable exception of microbiomes and sediments, most aDNA studies have focused on single taxa or taxonomic groups, making the study of changes at the community level challenging. This is rather surprising because current sequencing and analytical approaches allow us to obtain and analyse aDNA from multiple source materials. When combined, these data can enable the simultaneous study of multiple taxa through space and time, and could thus provide a more comprehensive understanding of ecosystem-wide changes. It is therefore timely to develop an integrative approach to aDNA studies by combining data from multiple taxa and substrates. In this review, we discuss the various applications, associated challenges and future prospects of such an approach.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes.

BACKGROUND: Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs. RESULTS: By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits. CONCLUSIONS: Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

A gene co-association network regulating gut microbial communities in a Duroc pig population.

BACKGROUND: Analyses of gut microbiome composition in livestock species have shown its potential to contribute to the regulation of complex phenotypes. However, little is known about the host genetic control over the gut microbial communities. In pigs, previous studies are based on classical "single-gene-single-trait" approaches and have evaluated the role of host genome controlling gut prokaryote and eukaryote communities separately. RESULTS: In order to determine the ability of the host genome to control the diversity and composition of microbial communities in healthy pigs, we undertook genome-wide association studies (GWAS) for 39 microbial phenotypes that included 2 diversity indexes, and the relative abundance of 31 bacterial and six commensal protist genera in 390 pigs genotyped for 70 K SNPs. The GWAS results were processed through a 3-step analytical pipeline comprised of (1) association weight matrix; (2) regulatory impact factor; and (3) partial correlation and information theory. The inferred gene regulatory network comprised 3561 genes (within a 5 kb distance from a relevant SNP-P < 0.05) and 738,913 connections (SNP-to-SNP co-associations). Our findings highlight the complexity and polygenic nature of the pig gut microbial ecosystem. Prominent within the network were 5 regulators, PRDM15, STAT1, ssc-mir-371, SOX9 and RUNX2 which gathered 942, 607, 588, 284 and 273 connections, respectively. PRDM15 modulates the transcription of upstream regulators of WNT and MAPK-ERK signaling to safeguard naive pluripotency and regulates the production of Th1- and Th2-type immune response. The signal transducer STAT1 has long been associated with immune processes and was recently identified as a potential regulator of vaccine response to porcine reproductive and respiratory syndrome. The list of regulators was enriched for immune-related pathways, and the list of predicted targets includes candidate genes previously reported as associated with microbiota profile in pigs, mice and human, such as SLIT3, SLC39A8, NOS1, IL1R2, DAB1, TOX3, SPP1, THSD7B, ELF2, PIANP, A2ML1, and IFNAR1. Moreover, we show the existence of host-genetic variants jointly associated with the relative abundance of butyrate producer bacteria and host performance. CONCLUSIONS: Taken together, our results identified regulators, candidate genes, and mechanisms linked with microbiome modulation by the host. They further highlight the value of the proposed analytical pipeline to exploit pleiotropy and the crosstalk between bacteria and protists as significant contributors to host-microbiome interactions and identify genetic markers and candidate genes that can be incorporated in breeding program to improve host-performance and microbial traits. Video Abstract.

Ancient microRNA profiles of a 14,300-year-old canid samples confirm taxonomic origin and give glimpses into tissue-specific gene regulation from the Pleistocene.

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue-specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 years, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by a) the presence of canid / Caniformia-specific sequences that never evolved outside of this clade, b) tissue-specific expression patterns (cartilage, liver and muscle) that resemble those of modern dogs and c) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing. .

A genome-wide association analysis for body, udder, and leg conformation traits recorded in Murciano-Granadina goats.

Morphological traits are of great importance to dairy goat production given their effect on phenotypes of economic interest. However, their underlying genomic architecture has not yet been extensively characterized. Herein, we aimed to identify genomic regions associated with body, udder, and leg conformation traits recorded in 825 Murciano-Granadina goats. We genotyped this resource population using the GoatSNP50 BeadChip (Illumina Inc., San Diego, CA) and performed genome-wide association analyses using the GEMMA software. We found 2 genome-wide significant associations between markers rs268273468 [Capra hircus (CHI) 16:69617700] and rs268249346 (CHI 28:18321523) and medial suspensory ligament. In contrast, we did not detect any genome-wide significant associations for body and leg traits. Moreover, we found 12, 19, and 7 chromosome-wide significant associations for udder, body, and leg traits, respectively. Comparison of our data with previous studies revealed a low level of positional concordance between regions associated with morphological traits. In addition to technical factors, this lack of concordance could be due to a substantial level of genetic heterogeneity among breeds or to the strong polygenic background of morphological traits, which makes it difficult to detect genetic factors that have small phenotypic effects.

Analyzing the genomic and transcriptomic architecture of milk traits in Murciano-Granadina goats.

BACKGROUND: In this study, we aimed to investigate the molecular basis of lactation as well as to identify the genetic factors that influence milk yield and composition in goats. To achieve these two goals, we have analyzed how the mRNA profile of the mammary gland changes in seven Murciano-Granadina goats at each of three different time points, i.e. 78 d (T1, early lactation), 216 d (T2, late lactation) and 285 d (T3, dry period) after parturition. Moreover, we have performed a genome-wide association study (GWAS) for seven dairy traits recorded in the 1st lactation of 822 Murciano-Granadina goats. RESULTS: The expression profiles of the mammary gland in the early (T1) and late (T2) lactation were quite similar (42 differentially expressed genes), while strong transcriptomic differences (more than one thousand differentially expressed genes) were observed between the lactating (T1/T2) and non-lactating (T3) mammary glands. A large number of differentially expressed genes were involved in pathways related with the biosynthesis of amino acids, cholesterol, triglycerides and steroids as well as with glycerophospholipid metabolism, adipocytokine signaling, lipid binding, regulation of ion transmembrane transport, calcium ion binding, metalloendopeptidase activity and complement and coagulation cascades. With regard to the second goal of the study, the performance of the GWAS allowed us to detect 24 quantitative trait loci (QTLs), including three genome-wide significant associations: QTL1 (chromosome 2, 130.72-131.01 Mb) for lactose percentage, QTL6 (chromosome 6, 78.90-93.48 Mb) for protein percentage and QTL17 (chromosome 17, 11.20 Mb) for both protein and dry matter percentages. Interestingly, QTL6 shows positional coincidence with the casein genes, which encode 80% of milk proteins. CONCLUSIONS: The abrogation of lactation involves dramatic changes in the expression of genes participating in a broad array of physiological processes such as protein, lipid and carbohydrate metabolism, calcium homeostasis, cell death and tissue remodeling, as well as immunity. We also conclude that genetic variation at the casein genes has a major impact on the milk protein content of Murciano-Granadina goats.

Co-expression network analysis predicts a key role of microRNAs in the adaptation of the porcine skeletal muscle to nutrient supply.

BACKGROUND: The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake. RESULTS: We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation. CONCLUSIONS: As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.