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.

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.

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.

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.

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.

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.

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.

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/.

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.

Unraveling molecular mechanisms involved in the development of leptin resistance using the pig as a model

The increase in obesity worldwide underlines the need for research concerning its metabolic and genetic determinants. One of the most intriguing mechanisms regarding obesity involves leptin and its signaling cascade. Leptin is a key regulator contributing to the fine-tuned crosstalk between nutrient availability and appetite signaling in the central nervous system. Owing to ethical concerns, many human tissues are not readily available and pigs can serve as a good animal model owing to their comparable anatomy, metabolism and genetics. In the present study, we utilized the pig to investigate the possible impact of increased adiposity on the development of alterations within the leptin signaling pathway. Two divergent groups of pigs (High and Low) were defined based on a high and low amount of mesenteric fat. Cortex, cerebellum, hypothalamus, mesenteric, subcutaneous and retroperitoneal fat tissues were used to study changes in expression levels of 94 mRNA transcripts related to the leptin signaling pathway using the qPCR approach. No significant differences were found at the central nervous system, whereas the expression level of STAT1 was reduced in mesenteric fat and leptin (LEP) and interleukin 6 (IL6) were shown to be consistently increased in all analyzed fat compartments from pigs with a high amount of mesenteric fat. These results could imply the onset of leptin and pro-inflammatory cytokine overexpression at early stages of obesity in the analyzed pigs without affecting any key components in the central nervous system. Thus, these pigs showing a unique leptin deregulation in adipose tissues could be a useful translational resource for studies of obesity and leptin resistance phenotypes.

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.

Comparing the diversity of the casein genes in the Asian mouflon and domestic sheep

We aimed to determine whether casein variants that are currently segregating in ovine populations existed before the domestication of sheep or, to the contrary, if their emergence is much more recent. To this end, we have retrieved whole-genome sequences from Iranian and domestic sheep from Africa, Europe, South and East Asia and West Asia. Population structure analysis based on 55,352,935 SNPs revealed a clear separation between Iranian mouflons and domestic sheep. Moreover, we also observed a strong genetic differentiation between Iranian mouflons sampled in geographic areas close to Tehran and Tabriz. Based on sequence data, hundreds of SNPs mapping to the casein α_S1 (CSN1S1, 248 SNPs), casein α_S2 (CSN1S2, 268 SNPs), casein ß (CSN2, 146 SNPs) and casein κ (CSN3, 112 SNPs) genes were identified. Approximately 25-63.02% of the casein variation was shared between Iranian mouflons and domestic sheep, and the four domestic sheep populations also shared 44.2-57.4% of the casein polymorphic sites. These findings suggest that an important fraction of the casein variation present in domestic sheep was already segregating in the mouflon prior to its domestication. Genomic studies performed in horses and dogs are consistent with this view, suggesting that much of the diversity that we currently detect in domestic animals comes from standing variation already segregating in their wild ancestors.

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.

Discovery and annotation of novel microRNAs in the porcine genome by using a semi-supervised transductive learning approach

Despite the broad variety of available microRNA (miRNA) prediction tools, their application to the discovery and annotation of novel miRNA genes in domestic species is still limited. In this study we designed a comprehensive pipeline (eMIRNA) for miRNA identification in the yet poorly annotated porcine genome and demonstrated the usefulness of implementing a motif search positional refinement strategy for the accurate determination of precursor miRNA boundaries. The small RNA fraction from gluteus medius skeletal muscle of 48 Duroc gilts was sequenced and used for the prediction of novel miRNA loci. Additionally, we selected the human miRNA annotation for a homology-based search of porcine miRNAs with orthologous genes in the human genome. A total of 20 novel expressed miRNAs were identified in the porcine muscle transcriptome and 27 additional novel porcine miRNAs were also detected by homology-based search using the human miRNA annotation. The existence of three selected novel miRNAs (ssc-miR-483, ssc-miR484 and ssc-miR-200a) was further confirmed by reverse transcription quantitative real-time PCR analyses in the muscle and liver tissues of Göttingen minipigs. In summary, the eMIRNA pipeline presented in the current work allowed us to expand the catalogue of porcine miRNAs and showed better performance than other commonly used miRNA prediction approaches. More importantly, the flexibility of our pipeline makes possible its application in other yet poorly annotated non-model species.

Detection of homozygous genotypes for a putatively lethal recessive mutation in the porcine argininosuccinate synthase 1 (ASS1) gene

The sequencing of the pig genome revealed the existence of homozygous individuals for a nonsense mutation in the argininosuccinate synthase 1 (ASS1) gene (rs81212146, c.944T>A, L315X). Paradoxically, an AA homozygous genotype for this polymorphism is expected to abolish the function of the ASS1 enzyme that participates in the urea cycle, leading to citrullinemia, hyperammonemia, coma and death. Sequencing of five Duroc boars that sired a population of 350 Duroc barrows revealed the segregation of the c.944T>A polymorphism, so we aimed to investigate its phenotypic consequences. Genotyping of this mutation in the 350 Duroc barrows revealed the existence of seven individuals homozygous (AA) for the nonsense mutation. These AA pigs had a normal weight despite the fact that mild citrullinemia often involves impaired growth. Sequencing of the region surrounding the mutation in TT, TA and AA individuals revealed that the A substitution in the second position of the codon (c.944T>A) is in complete linkage disequilibrium with a C replacement (c.943T>C) in the first position of the codon. This second mutation would compensate for the potentially damaging effect of the c.944T>A replacement. In fact, this is the most probable reason why pigs with homozygous AA genotypes at the 944 site of the ASS1 coding region are alive. Our results illustrate the complexities of predicting the consequences of nonsense mutations on gene function and phenotypes, not only because of annotation issues but also owing to the existence of genetic mechanisms that sometimes limit the penetrance of highly harmful mutations.

An association analysis for 14 candidate genes mapping to meat quality quantitative trait loci in a Duroc pig population reveals that the ATP1A2 genotype is highly associated with muscle electric conductivity

In previous GWAS carried out in a Duroc commercial line (Lipgen population), we detected on pig chromosomes 3, 4 and 14 several QTL for gluteus medius muscle redness (GM a*), electric conductivity in the longissimus dorsi muscle (LD CE) and vaccenic acid content in the LD muscle (LD C18:1 n - 7), respectively. We have genotyped, in the Lipgen population, 19 SNPs mapping to 14 genes located within these QTL. Subsequently, association analyses have been performed. After correction for multiple testing, two SNPs in the TGFBRAP1 (rs321173745) and SELENOI (rs330820437) genes were associated with GM a*, whereas ACADSB (rs81449951) and GPR26 (rs343087568) genotypes displayed significant associations with LD vaccenic content. Moreover, the polymorphisms located at the ATP1A2 (rs344748241), ATP8B2 (rs81382410) and CREB3L4 (rs321278469 and rs330133789) genes showed significant associations with LD CE. We made a second round of association analyses including the SNPs mentioned above as well as other SNPs located in the chromosomes to which they map. After performing a correction for multiple testing, the only association that remained significant at the chromosome-wide level was that between the ATP1A2 genotype and LD CE. From a functional point of view, this association is meaningful because this locus encodes a subunit of the Na⁺ /K⁺ -ATPase responsible for maintaining an electrochemical gradient across the plasma membrane.

Integrating genome-wide co-association and gene expression to identify putative regulators and predictors of feed efficiency in pigs

BACKGROUND

Feed efficiency (FE) has a major impact on the economic sustainability of pig production. We used a systems-based approach that integrates single nucleotide polymorphism (SNP) co-association and gene-expression data to identify candidate genes, biological pathways, and potential predictors of FE in a Duroc pig population.

RESULTS

We applied an association weight matrix (AWM) approach to analyse the results from genome-wide association studies (GWAS) for nine FE associated and production traits using 31K SNPs by defining residual feed intake (RFI) as the target phenotype. The resulting co-association network was formed by 829 SNPs. Additive effects of this SNP panel explained 61% of the phenotypic variance of RFI, and the resulting phenotype prediction accuracy estimated by cross-validation was 0.65 (vs. 0.20 using pedigree-based best linear unbiased prediction and 0.12 using the 31K SNPs). Sixty-eight transcription factor (TF) genes were identified in the co-association network; based on the lossless approach, the putative main regulators were COPS5, GTF2H5, RUNX1, HDAC4, ESR1, USP16, SMARCA2 and GTF2F2. Furthermore, gene expression data of the gluteus medius muscle was explored through differential expression and multivariate analyses. A list of candidate genes showing functional and/or structural associations with FE was elaborated based on results from both AWM and gene expression analyses, and included the aforementioned TF genes and other ones that have key roles in metabolism, e.g. ESRRG, RXRG, PPARGC1A, TCF7L2, LHX4, MAML2, NFATC3, NFKBIZ, TCEA1, CDCA7L, LZTFL1 or CBFB. The most enriched biological pathways in this list were associated with behaviour, immunity, nervous system, and neurotransmitters, including melatonin, glutamate receptor, and gustation pathways. Finally, an expression GWAS allowed identifying 269 SNPs associated with the candidate genes' expression (eSNPs). Addition of these eSNPs to the AWM panel of 829 SNPs did not improve the accuracy of genomic predictions.

CONCLUSIONS

Candidate genes that have a direct or indirect effect on FE-related traits belong to various biological processes that are mainly related to immunity, behaviour, energy metabolism, and the nervous system. The pituitary gland, hypothalamus and thyroid axis, and estrogen signalling play fundamental roles in the regulation of FE in pigs. The 829 selected SNPs explained 61% of the phenotypic variance of RFI, which constitutes a promising perspective for applying genetic selection on FE relying on molecular-based prediction.

About the existence of common determinants of gene expression in the porcine liver and skeletal muscle

BACKGROUND

The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population.

RESULTS

The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals.

CONCLUSION

Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

Polymorphisms of the cryptochrome 2 and mitoguardin 2 genes are associated with the variation of lipid-related traits in Duroc pigs

The genetic factors determining the phenotypic variation of porcine fatness phenotypes are still largely unknown. We investigated whether the polymorphism of eight genes (MIGA2, CRY2, NPAS2, CIART, ARNTL2, PER1, PER2 and PCK1), which display differential expression in the skeletal muscle of fasted and fed sows, is associated with the variation of lipid and mRNA expression phenotypes in Duroc pigs. The performance of an association analysis with the GEMMA software demonstrated that the rs330779504 SNP in the MIGA2 gene is associated with LDL concentration at 190 days (LDL₂, corrected P-value = 0.057). Moreover, the rs320439526 SNP of the CRY2 gene displayed a significant association with stearic acid content in the longissimus dorsi muscle (LD C18:0, corrected P-value = 0.015). Both SNPs were also associated with the mRNA levels of the corresponding genes in the gluteus medius skeletal muscle. From a biological perspective these results are meaningful because MIGA2 protein plays an essential role in mitochondrial fusion, a process tightly connected with the energy status of the cell, while CRY2 is a fundamental component of the circadian clock. However, inclusion of these two SNPs in chromosome-wide association analyses demonstrated that they are not located at the peaks of significance for the two traits under study (LDL₂ for rs330779504 and LD C18:0 for rs320439526), thus implying that these two SNPs do not have causal effects.

Genomic analysis of the origins of extant casein variation in goats

The variation in the casein genes has a major impact on the milk composition of goats. Even though many casein polymorphisms have been identified so far, we do not know yet whether they are evolutionarily ancient (i.e., they existed before domestication) or young (i.e., they emerged after domestication). Herewith, we identified casein polymorphisms in a data set of 106 caprine whole-genome sequences corresponding to bezoars (Capra aegagrus, the ancestor of domestic goats) and 4 domestic goat (Capra hircus) populations from Europe, Africa, the Far East, and the Near East. Domestic and wild goat populations shared a substantial number of casein SNP, from 36.1% (CSN2) to 55.1% (CSN1S2). The comparison of casein variation among bezoars and the 4 domestic goat populations demonstrated that more than 50% of the casein SNP are shared by 2 or more populations, and 18 to 44% are shared by all populations. Moreover, the majority of casein alleles reported in domestic goats also segregate in the bezoar, including several alleles displaying significant associations with milk composition (e.g., the A/B alleles of the CSN1S1 and CSN3 genes, the A allele of the CSN2 gene). We conclude that much of the current diversity of the caprine casein genes comes from ancient standing variation segregating in the ancestor of modern domestic goats.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles

Background

Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles.

Results

A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110–114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq.

Conclusions

We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5557-9) contains supplementary material, which is available to authorized users.

Analysing the Expression of Eight Clock Genes in Five Tissues From Fasting and Fed Sows

In a previous study, we observed that circadian clock genes are differentially expressed in the skeletal muscle of fasting and fed sows. The goal of the current work was to investigate if these genes are also differentially expressed in tissues containing the central (hypothalamus) and peripheral (duodenum, dorsal fat, muscle, and liver) clocks. As animal material, we used 12 sows that fasted 12 h before slaughtering (T0) and 12 sows that were fed ad libitum 7 h prior slaughtering (T2). Tissue samples were collected immediately after slaughter and total RNA was subsequently extracted. The expression of the ARNTL, BHLHE40, CRY2, NPAS2, NR1D1, PER1, PER2, and SIK1 genes was measured by quantitative reverse transcription PCR. The numbers of clock genes showing differential expression before and after feeding varied depending on the tissue i.e., four in dorsal fat and duodenum, six in skeletal muscle, and seven in the liver. In contrast, none of the eight analysed genes displayed a significant differential expression in hypothalamus, the tissue where the central clock resides. This result supports that the differential expression of clock genes in the four tissues mentioned above is probably induced by nutrition and not by the central clock entrained by light. Moreover, we have observed that the NPAS2 and ARNTL genes display positive log2(FC) values in the five tissues under analysis, whilst the CRY2, PER1 (except dorsal fat) and PER2 (except hypothalamus) genes generally show negative log2(FC) values. Such result might be explained by the existence of a negative feedback loop between the ARNTL/NPAS2 and CRY/PER genes. Collectively, these results support that nutrition plays an important role in modulating the timing of porcine peripheral circadian clocks. Such regulation could be essential for coordinating the subsequent metabolic response to nutrient supply.

Nutrient supply affects the mRNA expression profile of the porcine skeletal muscle

Background

The genetic basis of muscle fat deposition in pigs is not well known. So far, we have only identified a limited number of genes involved in the absorption, transport, storage and catabolism of lipids. Such information is crucial to interpret, from a biological perspective, the results of genome-wide association analyses for intramuscular fat content and composition traits. Herewith, we have investigated how the ingestion of food changes gene expression in the gluteus medius muscle of Duroc pigs.

Results

By comparing the muscle mRNA expression of fasted pigs (T0) with that of pigs sampled 5 h (T1) and 7 h (T2) after food intake, we have detected differential expression (DE) for 148 (T0-T1), 520 (T0-T2) and 135 (T1-T2) genes (q-value <0.05 and a |FC| > of 1.5). Many of these DE genes were transcription factors, suggesting that we have detected the coordinated response of the skeletal muscle to nutrient supply. We also found DE genes with a dual role in oxidative stress and angiogenesis (THBS1, THBS2 and TXNIP), two biological processes that are probably activated in the post-prandial state. Finally, we have identified several loci playing a key role in the modulation of circadian rhythms (ARNTL, PER1, PER2, BHLHE40, NR1D1, SIK1, CIART and CRY2), a result that indicates that the porcine muscle circadian clock is modulated by nutrition.

Conclusion

We have shown that hundreds of genes change their expression in the porcine skeletal muscle in response to nutrient intake. Many of these loci do not have a known metabolic role, a result that suggests that our knowledge about the genetic basis of muscle energy homeostasis is still incomplete.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3986-x) contains supplementary material, which is available to authorized users.