Differences in caecal microbiota composition and Salmonella carriage between experimentally infected inbred lines of chickens.

BACKGROUND: Salmonella Enteritidis (SE) is one of the major causes of human foodborne intoxication resulting from consumption of contaminated poultry products. Genetic selection of animals that are more resistant to Salmonella carriage and modulation of the gut microbiota are two promising ways to decrease individual Salmonella carriage. The aims of this study were to identify the main genetic and microbial factors that control the level of Salmonella carriage in chickens (Gallus gallus) under controlled experimental conditions. Two-hundred and forty animals from the White Leghorn inbred lines N and 61 were infected by SE at 7 days of age. After infection, animals were kept in isolators to reduce recontamination of birds by Salmonella. Caecal contents were sampled at 12 days post-infection and used for DNA extraction. Microbiota DNA was used to measure individual counts of SE by digital PCR and to determine the bacterial taxonomic composition, using a 16S rRNA gene high-throughput sequencing approach. RESULTS: Our results confirmed that the N line is more resistant to Salmonella carriage than the 61 line, and that intra-line variability is higher for the 61 line. Furthermore, the 16S analysis showed strong significant differences in microbiota taxonomic composition between the two lines. Among the 617 operational taxonomic units (OTU) observed, more than 390 were differentially abundant between the two lines. Furthermore, within the 61 line, we found a difference in the microbiota taxonomic composition between the high and low Salmonella carriers, with 39 differentially abundant OTU. Using metagenome functional prediction based on 16S data, several metabolic pathways that are potentially associated to microbiota taxonomic differences (e.g. short chain fatty acids pathways) were identified between high and low carriers. CONCLUSIONS: Overall, our findings demonstrate that the caecal microbiota composition differs between genetic lines of chickens. This could be one of the reasons why the investigated lines differed in Salmonella carriage levels under experimental infection conditions.

Influence of genetics and the pre-vaccination blood transcriptome on the variability of antibody levels after vaccination against Mycoplasma hyopneumoniae in pigs.

BACKGROUND: The impact of individual genetic and genomic variations on immune responses is an emerging lever investigated in vaccination strategies. In our study, we used genetic and pre-vaccination blood transcriptomic data to study vaccine effectiveness in pigs. RESULTS: A cohort of 182 Large White pigs was vaccinated against Mycoplasma hyopneumoniae (M. hyo) at weaning (28 days of age), with a booster 21 days later. Vaccine response was assessed by measuring seric M. hyo antibodies (Ab) at 0 (vaccination day), 21 (booster day), 28, 35, and 118 days post-vaccination (dpv). Inter-individual variability of M. hyo Ab levels was observed at all time points and the corresponding heritabilities ranged from 0.46 to 0.57. Ab persistence was higher in females than in males. Genome-wide association studies with a 658 K SNP panel revealed two genomic regions associated with variations of M. hyo Ab levels at 21 dpv at positions where immunity-related genes have been mapped, DAB2IP on chromosome 1, and ASAP1, CYRIB and GSDMC on chromosome 4. We studied covariations of Ab responses with the pre-vaccination blood transcriptome obtained by RNA-Seq for a subset of 82 pigs. Weighted gene correlation network and differential expression analyses between pigs that differed in Ab responses highlighted biological functions that were enriched in heme biosynthesis and platelet activation for low response at 21 dpv, innate antiviral immunity and dendritic cells for high response at 28 and 35 dpv, and cell adhesion and extracellular matrix for high response at 118 dpv. Sparse partial least squares discriminant analysis identified 101 genes that efficiently predicted divergent responders at all time points. We found weak negative correlations of M. hyo Ab levels with body weight traits, which revealed a trade-off that needs to be further explored. CONCLUSIONS: We confirmed the influence of the host genetics on vaccine effectiveness to M. hyo and provided evidence that the pre-vaccination blood transcriptome co-varies with the Ab response. Our results highlight that both genetic markers and blood biomarkers could be used as potential predictors of vaccine response levels and more studies are required to assess whether they can be exploited in breeding programs.

Identification of rumen microbial biomarkers linked to methane emission in Holstein dairy cows.

Mitigation of greenhouse gas emissions is relevant for reducing the environmental impact of ruminant production. In this study, the rumen microbiome from Holstein cows was characterized through a combination of 16S rRNA gene and shotgun metagenomic sequencing. Methane production (CH4 ) and dry matter intake (DMI) were individually measured over 4-6 weeks to calculate the CH4 yield (CH4 y = CH4 /DMI) per cow. We implemented a combination of clustering, multivariate and mixed model analyses to identify a set of operational taxonomic unit (OTU) jointly associated with CH4 y and the structure of ruminal microbial communities. Three ruminotype clusters (R1, R2 and R3) were identified, and R2 was associated with higher CH4 y. The taxonomic composition on R2 had lower abundance of Succinivibrionaceae and Methanosphaera, and higher abundance of Ruminococcaceae, Christensenellaceae and Lachnospiraceae. Metagenomic data confirmed the lower abundance of Succinivibrionaceae and Methanosphaera in R2 and identified genera (Fibrobacter and unclassified Bacteroidales) not highlighted by metataxonomic analysis. In addition, the functional metagenomic analysis revealed that samples classified in cluster R2 were overrepresented by genes coding for KEGG modules associated with methanogenesis, including a significant relative abundance of the methyl-coenzyme M reductase enzyme. Based on the cluster assignment, we applied a sparse partial least-squares discriminant analysis at the taxonomic and functional levels. In addition, we implemented a sPLS regression model using the phenotypic variation of CH4 y. By combining these two approaches, we identified 86 discriminant bacterial OTUs, notably including families linked to CH4 emission such as Succinivibrionaceae, Ruminococcaceae, Christensenellaceae, Lachnospiraceae and Rikenellaceae. These selected OTUs explained 24% of the CH4 y phenotypic variance, whereas the host genome contribution was ~14%. In summary, we identified rumen microbial biomarkers associated with the methane production of dairy cows; these biomarkers could be used for targeted methane-reduction selection programmes in the dairy cattle industry provided they are heritable.

Host genotype and amoxicillin administration affect the incidence of diarrhoea and faecal microbiota of weaned piglets during a natural multiresistant ETEC infection.

Enterotoxigenic Escherichia coli (ETEC) is the aetiological agent of postweaning diarrhoea (PWD) in piglets. The SNPs located on the Mucine 4 (MUC4) and Fucosyltransferase 1 (FUT1) genes have been associated with the susceptibility to ETEC F4 and ETEC F18, respectively. The interplay between the MUC4 and FUT1 genotypes to ETEC infection and the use of amoxicillin in modifying the intestinal microbiota during a natural infection by multiresistant ETEC strains have never been investigated. The aim of this study was to evaluate the effects of the MUC4 and FUT1 genotypes and the administration of amoxicillin through different routes on the presence of diarrhoea and the faecal microbiota composition in piglets naturally infected with ETEC. Seventy-one piglets were divided into three groups: two groups differing by amoxicillin administration routes-parenteral (P) or oral (O) and a control group without antibiotics (C). Faecal scores, body weight, presence of ETEC F4 and F18 were investigated 4 days after the arrival in the facility (T0), at the end of the amoxicillin administration (T1) and after the withdrawal period (T2). The faecal bacteria composition was assessed by sequencing the 16S rRNA gene. We described that MUC4 and FUT1 genotypes were associated with the presence of ETEC F4 and ETEC F18. The faecal microbiota was influenced by the MUC4 genotypes at T0. We found the oral administration to be associated with the presence of diarrhoea at T1 and T2. Furthermore, the exposure to amoxicillin resulted in significant alterations of the faecal microbiota. Overall, MUC4 and FUT1 were confirmed as genetic markers for the susceptibility to ETEC infections in pigs. Moreover, our data highlight that group amoxicillin treatment may produce adverse outcomes on pig health in course of multiresistant ETEC infection. Therefore, alternative control measures able to maintain a healthy faecal microbiota in weaners are recommended.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing.

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

Genome-wide immunity studies in the rabbit: transcriptome variations in peripheral blood mononuclear cells after in vitro stimulation by LPS or PMA-Ionomycin.

BACKGROUND: Our purpose was to obtain genome-wide expression data for the rabbit species on the responses of peripheral blood mononuclear cells (PBMCs) after in vitro stimulation by lipopolysaccharide (LPS) or phorbol myristate acetate (PMA) and ionomycin. This transcriptome profiling was carried out using microarrays enriched with immunity-related genes, and annotated with the most recent data available for the rabbit genome. RESULTS: The LPS affected 15 to 20 times fewer genes than PMA-Ionomycin after both 4 hours (T4) and 24 hours (T24), of in vitro stimulation, in comparison with mock-stimulated PBMCs. LPS induced an inflammatory response as shown by a significant up-regulation of IL12A and CXCL11 at T4, followed by an increased transcription of IL6, IL1B, IL1A, IL36, IL37, TNF, and CCL4 at T24. Surprisingly, we could not find an up-regulation of IL8 either at T4 or at T24, and detected a down-regulation of DEFB1 and BPI at T24. A concerted up-regulation of SAA1, S100A12 and F3 was found upon stimulation by LPS. PMA-Ionomycin induced a very early expression of Th1, Th2, Treg, and Th17 responses by PBMCs at T4. The Th1 response increased at T24 as shown by the increase of the transcription of IFNG and by contrast to other cytokines which significantly decreased from T4 to T24 (IL2, IL4, IL10, IL13, IL17A, CD69) by comparison to mock-stimulation. The granulocyte-macrophage colony-stimulating factor (CSF2) was by far the most over-expressed gene at both T4 and T24 by comparison to mock-stimulated cells, confirming a major impact of PMA-Ionomycin on cell growth and proliferation. A significant down-regulation of IL16 was observed at T4 and T24, in agreement with a role of IL16 in PBMC apoptosis. CONCLUSIONS: We report new data on the responses of PBMCs to LPS and PMA-Ionomycin in the rabbit species, thus enlarging the set of mammalian species for which such reports exist. The availability of the rabbit genome assembly together with high throughput genomic tools should pave the way for more intense genomic studies for this species, which is known to be a very relevant biomedical model in immunology and physiology.

A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis.

During the last 3 years, a number of approaches for the normalization of RNA sequencing data have emerged in the literature, differing both in the type of bias adjustment and in the statistical strategy adopted. However, as data continue to accumulate, there has been no clear consensus on the appropriate normalization method to be used or the impact of a chosen method on the downstream analysis. In this work, we focus on a comprehensive comparison of seven recently proposed normalization methods for the differential analysis of RNA-seq data, with an emphasis on the use of varied real and simulated datasets involving different species and experimental designs to represent data characteristics commonly observed in practice. Based on this comparison study, we propose practical recommendations on the appropriate normalization method to be used and its impact on the differential analysis of RNA-seq data.

Epigenetic regulation of the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs.

BACKGROUND: In previous studies on an Iberian x Landrace cross, we have provided evidence that supported the porcine ELOVL6 gene as the major causative gene of the QTL on pig chromosome 8 for palmitic and palmitoleic acid contents in muscle and backfat. The single nucleotide polymorphism (SNP) ELOVL6:c.-533C > T located in the promoter region of ELOVL6 was found to be highly associated with ELOVL6 expression and, accordingly, with the percentages of palmitic and palmitoleic acids in longissimus dorsi and adipose tissue. The main goal of the current work was to further study the role of ELOVL6 on these traits by analyzing the regulation of the expression of ELOVL6 and the implication of ELOVL6 polymorphisms on meat quality traits in pigs. RESULTS: High-throughput sequencing of BAC clones that contain the porcine ELOVL6 gene coupled to RNAseq data re-analysis showed that two isoforms of this gene are expressed in liver and adipose tissue and that they differ in number of exons and 3'UTR length. Although several SNPs in the 3'UTR of ELOVL6 were associated with palmitic and palmitoleic acid contents, this association was lower than that previously observed with SNP ELOVL6:c.-533C > T. This SNP is in full linkage disequilibrium with SNP ELOVL6:c.-394G > A that was identified in the binding site for estrogen receptor alpha (ERα). Interestingly, the ELOVL6:c.-394G allele is associated with an increase in methylation levels of the ELOVL6 promoter and with a decrease of ELOVL6 expression. Therefore, ERα is clearly a good candidate to explain the regulation of ELOVL6 expression through dynamic epigenetic changes in the binding site of known regulators of ELOVL6 gene, such as SREBF1 and SP1. CONCLUSIONS: Our results strongly suggest the ELOVL6:c.-394G > A polymorphism as the causal mutation for the QTL on pig chromosome 8 that affects fatty acid composition in pigs.

Interrelation between gut microbiota, SCFA, and fatty acid composition in pigs.

The gut microbiota is a key player in the host metabolism. Some bacteria are able to ferment non-digestible compounds and produce short-chain fatty acids that the host can later transform and accumulate in tissue. In this study, we aimed to better understand the relationships between the microorganisms and the short-chain fatty acid composition of the rectal content, including the possible linkage with the fatty acid composition in backfat and muscle of the pig. We studied a Duroc × Iberian crossbred population, and we found significant correlations between different bacterial and archaeal genera and the fatty acid profile. The abundance of n-butyric acid in the rectal content was positively associated with Prevotella spp. and negatively associated with Akkermansia spp., while conversely, the abundance of acetic acid was negatively and positively associated with the levels of Prevotella spp. and Akkermansia spp., respectively. The most abundant genus, Rikenellaceae RC9 gut group, had a positive correlation with palmitic acid in muscle and negative correlations with stearic acid in backfat and oleic acid in muscle. These results suggest the possible role of Prevotella spp. and Akkermansia spp. as biomarkers for acetic and n-butyric acids, and the relationship of Rikenellaceae RC9 gut group with the lipid metabolism, building up the potential, although indirect, role of the microbiota in the modification of the backfat and muscle fatty acid composition of the host.IMPORTANCEThe vital role of the gut microbiota on its host metabolism makes it essential to know how its modulation is mirrored on the fatty acid composition of the host. Our findings suggest Prevotella spp. and Akkermansia spp. as potential biomarkers for the levels of beneficial short-chain fatty acids and the possible influence of Rikenellaceae RC9 gut group in the backfat and muscle fatty acid composition of the pig.

Malignant features of minipig melanomas prior to spontaneous regression.

In MeLiM minipigs, melanomas develop around birth, can metastasize, and have histopathologic characteristics similar to humans. Interestingly, MeLiM melanomas eventually regress. This favorable outcome raises the question of their malignancy, which we investigated. We clinically followed tens of tumors from onset to first signs of regression. Transcriptome analysis revealed an enrichment of all cancer hallmarks in melanomas, although no activating or suppressing somatic mutation were found in common driver genes. Analysis of tumor cell genomes revealed high mutation rates without UV signature. Canonical proliferative, survival and angiogenic pathways were detected in MeLiM tumor cells all along progression stages. Functionally, we show that MeLiM melanoma cells are capable to grow in immunocompromised mice, with serial passages and for a longer time than in MeLiM pigs. Pigs set in place an immune response during progression with dense infiltration by myeloid cells while melanoma cells are deficient in B2M expression. To conclude, our data on MeLiM melanomas reveal several malignancy characteristics. The combination of these features with the successful spontaneous regression of these tumors make it an outstanding model to study an efficient anti-tumor immune response.

The peripheral blood transcriptome reflects variations in immunity traits in swine: towards the identification of biomarkers.

BACKGROUND: Immune traits (ITs) are potentially relevant criteria to characterize an individual's immune response. Our aim was to investigate whether the peripheral blood transcriptome can provide a significant and comprehensive view of IT variations in pig. RESULTS: Sixty-day-old Large White pigs classified as extreme for in vitro production of IL2, IL10, IFNγ and TNFα, phagocytosis activity, in vivo CD4⁻/CD8⁺ or TCRγδ + cell counts, and anti-Mycoplasma antibody levels were chosen to perform a blood transcriptome analysis with a porcine generic array enriched with immunity-related genes. Differentially expressed (DE) genes for in vitro production of IL2 and IL10, phagocytosis activity and CD4⁻/CD8⁺ cell counts were identified. Gene set enrichment analysis revealed a significant over-representation of immune response functions. To validate the microarray-based results, a subset of DE genes was confirmed by RT-qPCR. An independent set of 74 animals was used to validate the covariation between gene expression levels and ITs. Five potential gene biomarkers were found for prediction of IL2 (RALGDS), phagocytosis (ALOX12) or CD4⁻/CD8⁺ cell count (GNLY, KLRG1 and CX3CR1). On average, these biomarkers performed with a sensitivity of 79% and a specificity of 86%. CONCLUSIONS: Our results confirmed that gene expression profiling in blood represents a relevant molecular phenotype to refine ITs in pig and to identify potential biomarkers that can provide new insights into immune response analysis.

Analysis of porcine adipose tissue transcriptome reveals differences in de novo fatty acid synthesis in pigs with divergent muscle fatty acid composition.

BACKGROUND: In pigs, adipose tissue is one of the principal organs involved in the regulation of lipid metabolism. It is particularly involved in the overall fatty acid synthesis with consequences in other lipid-target organs such as muscles and the liver. With this in mind, we have used massive, parallel high-throughput sequencing technologies to characterize the porcine adipose tissue transcriptome architecture in six Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition (three per group). RESULTS: High-throughput RNA sequencing was used to generate a whole characterization of adipose tissue (backfat) transcriptome. A total of 4,130 putative unannotated protein-coding sequences were identified in the 20% of reads which mapped in intergenic regions. Furthermore, 36% of the unmapped reads were represented by interspersed repeats, SINEs being the most abundant elements. Differential expression analyses identified 396 candidate genes among divergent animals for intramuscular fatty acid composition. Sixty-two percent of these genes (247/396) presented higher expression in the group of pigs with higher content of intramuscular SFA and MUFA, while the remaining 149 showed higher expression in the group with higher content of PUFA. Pathway analysis related these genes to biological functions and canonical pathways controlling lipid and fatty acid metabolisms. In concordance with the phenotypic classification of animals, the major metabolic pathway differentially modulated between groups was de novo lipogenesis, the group with more PUFA being the one that showed lower expression of lipogenic genes. CONCLUSIONS: These results will help in the identification of genetic variants at loci that affect fatty acid composition traits. The implications of these results range from the improvement of porcine meat quality traits to the application of the pig as an animal model of human metabolic diseases.

The genome sequencing of an albino Western lowland gorilla reveals inbreeding in the wild.

BACKGROUND: The only known albino gorilla, named Snowflake, was a male wild born individual from Equatorial Guinea who lived at the Barcelona Zoo for almost 40 years. He was diagnosed with non-syndromic oculocutaneous albinism, i.e. white hair, light eyes, pink skin, photophobia and reduced visual acuity. Despite previous efforts to explain the genetic cause, this is still unknown. Here, we study the genetic cause of his albinism and making use of whole genome sequencing data we find a higher inbreeding coefficient compared to other gorillas. RESULTS: We successfully identified the causal genetic variant for Snowflake's albinism, a non-synonymous single nucleotide variant located in a transmembrane region of SLC45A2. This transporter is known to be involved in oculocutaneous albinism type 4 (OCA4) in humans. We provide experimental evidence that shows that this amino acid replacement alters the membrane spanning capability of this transmembrane region. Finally, we provide a comprehensive study of genome-wide patterns of autozygogosity revealing that Snowflake's parents were related, being this the first report of inbreeding in a wild born Western lowland gorilla. CONCLUSIONS: In this study we demonstrate how the use of whole genome sequencing can be extended to link genotype and phenotype in non-model organisms and it can be a powerful tool in conservation genetics (e.g., inbreeding and genetic diversity) with the expected decrease in sequencing cost.

Gastrointestinal Microbial Ecology of Weaned Piglets Fed Diets with Different Levels of Glyphosate.

Glyphosate possesses antimicrobial properties, and the present study investigated potential effects of feed glyphosate on piglet gastrointestinal microbial ecology. Weaned piglets were allocated to four diets (glyphosate contents [mg/kg feed]: 0 mg/kg control [CON; i.e., basal diet with no glyphosate added], 20 mg/kg as Glyphomax commercial herbicide [GM20], and 20 mg/kg [IPA20] and 200 mg/kg [IPA200] as glyphosate isopropylamine [IPA] salt). Piglets were sacrificed after 9 and 35 days of treatment, and stomach, small intestine, cecum, and colon digesta were analyzed for glyphosate, aminomethylphosphonic acid (AMPA), organic acids, pH, dry matter content, and microbiota composition. Digesta glyphosate contents reflected dietary levels (on day 35, 0.17, 16.2, 20.5, and 207.5 mg/kg colon digesta, respectively). Overall, we observed no significant glyphosate-associated effects on digesta pH, dry matter content, and-with few exceptions-organic acid levels. On day 9, only minor gut microbiota changes were observed. On day 35, we observed a significant glyphosate-associated decrease in species richness (CON, 462; IPA200, 417) and in the relative abundance of certain Bacteroidetes genera: CF231 (CON, 3.71%; IPA20, 2.33%; IPA200, 2.07%) and g_0.24 (CON, 3.69%; IPA20, 2.07%; IPA200, 1.75%) in cecum. No significant changes were observed at the phylum level. In the colon, we observed a significant glyphosate-associated increase in the relative abundance of Firmicutes (CON, 57.7%; IPA20, 69.4%; IPA200, 66.1%) and a decrease in Bacteroidetes (CON, 32.6%; IPA20, 23.5%). Significant changes were only observed for few genera, e.g., g_0.24 (CON, 7.12%; IPA20, 4.59%; IPA200, 4.00%). In conclusion, exposing weaned piglets to glyphosate-amended feed did not affect gastrointestinal microbial ecology to a degree that was considered actual dysbiosis, e.g., no potential pathogen bloom was observed. IMPORTANCE Glyphosate residues can be found in feed made from genetically modified glyphosate-resistant crops treated with glyphosate or from conventional crops, desiccated with glyphosate before harvest. If these residues affect the gut microbiota to an extent that is unfavorable to livestock health and productivity, the widespread use of glyphosate on feed crops may need to be reconsidered. Few in vivo studies have been conducted to investigate potential impact of glyphosate on the gut microbial ecology and derived health issues of animals, in particular livestock, when exposed to dietary glyphosate residues. The aim of the present study was therefore to investigate potential effects on the gastrointestinal microbial ecology of newly weaned piglets fed glyphosate-amended diets. Piglets did not develop actual gut dysbiosis when fed diets, containing a commercial herbicide formulation or a glyphosate salt at the maximum residue level, defined by the European Union for common feed crops, or at a 10-fold-higher level.

Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition.

BACKGROUND: New advances in high-throughput technologies have allowed for the massive analysis of genomic data, providing new opportunities for the characterization of the transcriptome architectures. Recent studies in pigs have employed RNA-Seq to explore the transcriptome of different tissues in a reduced number of animals. The main goal of this study was the identification of differentially-expressed genes in the liver of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq. RESULTS: The liver transcriptomes of two female groups (H and L) with phenotypically extreme intramuscular fatty acid composition were sequenced using RNA-Seq. A total of 146 and 180 unannotated protein-coding genes were identified in intergenic regions for the L and H groups, respectively. In addition, a range of 5.8 to 7.3% of repetitive elements was found, with SINEs being the most abundant elements. The expression in liver of 186 (L) and 270 (H) lncRNAs was also detected. The higher reproducibility of the RNA-Seq data was validated by RT-qPCR and porcine expression microarrays, therefore showing a strong correlation between RT-qPCR and RNA-Seq data (ranking from 0.79 to 0.96), as well as between microarrays and RNA-Seq (r=0.72). A differential expression analysis between H and L animals identified 55 genes differentially-expressed between groups. Pathways analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and fatty acid metabolism. In concordance with the phenotypic classification, the pathways analysis inferred that linolenic and arachidonic acids metabolism was altered between extreme individuals. In addition, a connection was observed among the top three networks, hence suggesting that these genes are interconnected and play an important role in lipid and fatty acid metabolism. CONCLUSIONS: In the present study RNA-Seq was used as a tool to explore the liver transcriptome of pigs with extreme phenotypes for intramuscular fatty acid composition. The differential gene expression analysis showed potential gene networks which affect lipid and fatty acid metabolism. These results may help in the design of selection strategies to improve the sensorial and nutritional quality of pork meat.

Caecal microbiota composition of experimental inbred MHC-B lines infected with IBV differs according to genetics and vaccination.

Interactions between the gut microbiota and the immune system may be involved in vaccine and infection responses. In the present study, we studied the interactions between caecal microbiota composition and parameters describing the immune response in six experimental inbred chicken lines harboring different MHC haplotypes. Animals were challenge-infected with the infectious bronchitis virus (IBV), and half of them were previously vaccinated against this pathogen. We explored to what extent the gut microbiota composition and the genetic line could be related to the immune response, evaluated through flow cytometry. To do so, we characterized the caecal bacterial communities with a 16S rRNA gene amplicon sequencing approach performed one week after the IBV infectious challenge. We observed significant effects of both the vaccination and the genetic line on the microbiota after the challenge infection with IBV, with a lower bacterial richness in vaccinated chickens. We also observed dissimilar caecal community profiles among the different lines, and between the vaccinated and non-vaccinated animals. The effect of vaccination was similar in all the lines, with a reduced abundance of OTU from the Ruminococcacea UCG-014 and Faecalibacterium genera, and an increased abundance of OTU from the Eisenbergiella genus. The main association between the caecal microbiota and the immune phenotypes involved TCRϒδ expression on TCRϒδ+ T cells. This phenotype was negatively associated with OTU from the Escherichia-Shigella genus that were also less abundant in the lines with the highest responses to the vaccine. We proved that the caecal microbiota composition is associated with the IBV vaccine response level in inbred chicken lines, and that the TCRϒδ+ T cells (judged by TCRϒδ expression) may be an important component involved in this interaction, especially with bacteria from the Escherichia-Shigella genus. We hypothesized that bacteria from the Escherichia-Shigella genus increased the systemic level of bacterial lipid antigens, which subsequently mitigated poultry γδ T cells.

Comparative analysis of the caecal tonsil transcriptome in two chicken lines experimentally infected with Salmonella Enteritidis.

Managing Salmonella enterica Enteritidis (SE) carriage in chicken is necessary to ensure human food safety and enhance the economic, social and environmental sustainability of chicken breeding. Salmonella can contaminate poultry products, causing human foodborne disease and economic losses for farmers. Both genetic selection for a decreased carriage and gut microbiota modulation strategies could reduce Salmonella propagation in farms. Two-hundred and twenty animals from the White Leghorn inbred lines N and 61 were raised together on floor, infected by SE at 7 days of age, transferred into isolators to prevent oro-fecal recontamination and euthanized at 12 days post-infection. Caecal content DNA was used to measure individual Salmonella counts (ISC) by droplet digital PCR. A RNA sequencing approach was used to measure gene expression levels in caecal tonsils after infection of 48 chicks with low or high ISC. The analysis between lines identified 7516 differentially expressed genes (DEGs) corresponding to 62 enriched Gene Ontology (GO) Biological Processes (BP) terms. A comparison between low and high carriers allowed us to identify 97 DEGs and 23 enriched GO BP terms within line 61, and 1034 DEGs and 288 enriched GO BP terms within line N. Among these genes, we identified several candidate genes based on their putative functions, including FUT2 or MUC4, which could be involved in the control of SE infection, maybe through interactions with commensal bacteria. Altogether, we were able to identify several genes and pathways associated with differences in SE carriage level. These results are discussed in relation to individual caecal microbiota compositions, obtained for the same animals in a previous study, which may interact with host gene expression levels for the control of the caecal SE load.

Genetic and functional evaluation of MITF as a candidate gene for cutaneous melanoma predisposition in pigs.

Cutaneous melanoma arises from transformed melanocytes and is caused mainly by environmental effects such as ultraviolet radiation and to a lesser extent by predisposing genetic variants. Only a few susceptibility genes for cutaneous melanoma have been identified so far in human; therefore, animal models represent a valuable alternative for genetic studies of this disease. In a previous quantitative trait locus (QTL) study, several susceptibility regions were identified in a swine biomedical model, the MeLiM (Melanoblastoma-bearing Libechov minipig) pigs. This article details the fine-mapping of a QTL located on SSC13 (Sus scrofa chromosome 13) through an increase in marker density. New microsatellites were used to confirm the results of the first analysis, and MITF (microphthalmia-associated transcription factor) was selected as a candidate gene for melanoma development. A single-marker association analysis was performed with single-nucleotide polymorphisms (SNPs) spread over the locus, but it did not reveal a significant association with diverse melanoma-related traits. In parallel, MITF alternative transcripts were characterized and their expression was investigated in different porcine tissues. The obtained results showed a complex transcriptional regulation concordant with the one present in other mammals. Notably, the ratio between MITF+ and MITF- isoforms in melanoma samples followed the same pattern as in human tumors, which highlights the adequacy of the MeLiM pig as a model for human melanoma. In conclusion, although MITF does not seem to be the causal gene of the QTL initially observed, we do not exclude a prominent role of its transcription and function in the outbreak and evolution of the tumors observed in pigs.