Dynamics of DNA methylation during osteogenic differentiation of porcine synovial membrane mesenchymal stem cells from two metabolically distinct breeds.

Mesenchymal stem cells (MSCs), with the ability to differentiate into osteoblasts, adipocytes, or chondrocytes, show evidence that the donor cell's metabolic type influences the osteogenic process. Limited knowledge exists on DNA methylation changes during osteogenic differentiation and the impact of diverse donor genetic backgrounds on MSC differentiation. In this study, synovial membrane mesenchymal stem cells (SMSCs) from two pig breeds (Angeln Saddleback, AS; German Landrace, DL) with distinct metabolic phenotypes were isolated, and the methylation pattern of SMSCs during osteogenic induction was investigated. Results showed that most differentially methylated regions (DMRs) were hypomethylated in osteogenic-induced SMSC group. These DMRs were enriched with genes of different osteogenic signalling pathways at different time points including Wnt, ECM, TGFB and BMP signalling pathways. AS pigs consistently exhibited a higher number of hypermethylated DMRs than DL pigs, particularly during the peak of osteogenesis (day 21). Predicting transcription factor motifs in regions of DMRs linked to osteogenic processes and donor breeds revealed influential motifs, including KLF1, NFATC3, ZNF148, ASCL1, FOXI1, and KLF5. These findings contribute to understanding the pattern of methylation changes promoting osteogenic differentiation, emphasizing the substantial role of donor the metabolic type and epigenetic memory of different donors on SMSC differentiation.

Transcriptomic Response of Differentiating Porcine Myotubes to Thermal Stress and Donor Piglet Age.

Climate change is a current concern that directly and indirectly affects agriculture, especially the livestock sector. Neonatal piglets have a limited thermoregulatory capacity and are particularly stressed by ambient temperatures outside their optimal physiological range, which has a major impact on their survival rate. In this study, we focused on the effects of thermal stress (35 °C, 39 °C, and 41 °C compared to 37 °C) on differentiating myotubes derived from the satellite cells of Musculus rhomboideus, isolated from two different developmental stages of thermolabile 5-day-old (p5) and thermostable 20-day-old piglets (p20). Analysis revealed statistically significant differential expression genes (DEGs) between the different cultivation temperatures, with a higher number of genes responding to cold treatment. These DEGs were involved in the macromolecule degradation and actin kinase cytoskeleton categories and were observed at lower temperatures (35 °C), whereas at higher temperatures (39 °C and 41 °C), the protein transport system, endoplasmic reticulum system, and ATP activity were more pronounced. Gene expression profiling of HSP and RBM gene families, which are commonly associated with cold and heat responses, exhibited a pattern dependent on temperature variability. Moreover, thermal stress exhibited an inhibitory effect on cell cycle, with a more pronounced downregulation during cold stress driven by ADGR genes. Additionally, our analysis revealed DEGs from donors with an undeveloped thermoregulation capacity (p5) and those with a fully developed thermoregulation capacity (p20) under various cultivation temperature. The highest number of DEGs and significant GO terms was observed under temperatures of 35 °C and 37 °C. In particular, under 35 °C, the DEGs were enriched in insulin, thyroid hormone, and calcium signaling pathways. This result suggests that the different thermoregulatory capacities of the donor piglets determined the ability of the primary muscle cell culture to differentiate into myotubes at different temperatures. This work sheds new light on the underlying molecular mechanisms that govern piglet differentiating myotube response to thermal stress and can be leveraged to develop effective thermal management strategies to enhance skeletal muscle growth.

RNA-Seq-based discovery of genetic variants and allele-specific expression of two layer lines and broiler chicken.

Recent advances in the selective breeding of broilers and layers have made poultry production one of the fastest-growing industries. In this study, a transcriptome variant calling approach from RNA-seq data was used to determine population diversity between broilers and layers. In total, 200 individuals were analyzed from three different chicken populations (Lohmann Brown (LB), n = 90), Lohmann Selected Leghorn (LSL, n = 89), and Broiler (BR, n = 21). The raw RNA-sequencing reads were pre-processed, quality control checked, mapped to the reference genome, and made compatible with Genome Analysis ToolKit for variant detection. Subsequently, pairwise fixation index (F ST) analysis was performed between broilers and layers. Numerous candidate genes were identified, that were associated with growth, development, metabolism, immunity, and other economically significant traits. Finally, allele-specific expression (ASE) analysis was performed in the gut mucosa of LB and LSL strains at 10, 16, 24, 30, and 60 weeks of age. At different ages, the two-layer strains showed significantly different allele-specific expressions in the gut mucosa, and changes in allelic imbalance were observed across the entire lifespan. Most ASE genes are involved in energy metabolism, including sirtuin signaling pathways, oxidative phosphorylation, and mitochondrial dysfunction. A high number of ASE genes were found during the peak of laying, which were particularly enriched in cholesterol biosynthesis. These findings indicate that genetic architecture as well as biological processes driving particular demands relate to metabolic and nutritional requirements during the laying period shape allelic heterogeneity. These processes are considerably affected by breeding and management, whereby elucidating allele-specific gene regulation is an essential step towards deciphering the genotype to phenotype map or functional diversity between the chicken populations. Additionally, we observed that several genes showing significant allelic imbalance also colocalized with the top 1% of genes identified by the FST approach, suggesting a fixation of genes in cis-regulatory elements.

Exploration of genotype-by-environment interactions affecting gene expression responses in porcine immune cells.

As one of the keys to healthy performance, robustness of farm animals is gaining importance, and with this comes increasing interest in genetic dissection of genotype-by-environment interactions (G×E). Changes in gene expression are among the most sensitive responses conveying adaptation to environmental stimuli. Environmentally responsive regulatory variation thus likely plays a central role in G×E. In the present study, we set out to detect action of environmentally responsive cis-regulatory variation by the analysis of condition-dependent allele specific expression (cd-ASE) in porcine immune cells. For this, we harnessed mRNA-sequencing data of peripheral blood mononuclear cells (PBMCs) stimulated in vitro with lipopolysaccharide, dexamethasone, or their combination. These treatments mimic common challenges such as bacterial infection or stress, and induce vast transcriptome changes. About two thirds of the examined loci showed significant ASE in at least one treatment, and out of those about ten percent exhibited cd-ASE. Most of the ASE variants were not yet reported in the PigGTEx Atlas. Genes showing cd-ASE were enriched in cytokine signaling in immune system and include several key candidates for animal health. In contrast, genes showing no ASE featured cell-cycle related functions. We confirmed LPS-dependent ASE for one of the top candidates, SOD2, which ranks among the major response genes in LPS-stimulated monocytes. The results of the present study demonstrate the potential of in vitro cell models coupled with cd-ASE analysis for the investigation of G×E in farm animals. The identified loci may benefit efforts to unravel the genetic basis of robustness and improvement of health and welfare in pigs.

The dynamics of molecular, immune and physiological features of the host and the gut microbiome, and their interactions before and after onset of laying in two hen strains.

Aggregation of data, including deep sequencing of mRNA and miRNA data in jejunum mucosa, abundance of immune cells, metabolites, or hormones in blood, composition of microbiota in digesta and duodenal mucosa, and production traits collected along the lifespan, provides a comprehensive picture of lifelong adaptation processes. Here, respective data from two laying hen strains (Lohmann Brown-Classic (LB) and Lohmann LSL-Classic (LSL) collected at 10, 16, 24, 30, and 60 wk of age were analyzed. Data integration revealed strain- and stage-specific biosignatures, including elements indicative of molecular pathways discriminating the strains. Although the strains performed the same, they differed in the activity of immunological and metabolic functions and pathways and showed specific gut-microbiota-interactions in different production periods. The study shows that both strains employ different strategies to acquire and maintain their capabilities under high performance conditions, especially during the transition phase. Furthermore, the study demonstrates the capacity of such integrative analyses to elucidate molecular pathways that reflect functional biodiversity. The bioinformatic reduction of the multidimensional data provides good guidance for further manual review of the data.

Glucocorticoid receptor hypersensitivity enhances inflammatory signaling and inhibits cell cycle progression in porcine PBMCs.

The consequences of glucocorticoid receptor (GR) hypersensitivity during infection have so far received little attention. We previously discovered that a natural gain-of-function Ala610Val substitution in the porcine GR aggravates response of pigs to lipopolysaccharide (LPS)-induced endotoxemia, which can be alleviated by dexamethasone (DEX) pretreatment. In this work, we investigated the relevant molecular basis of these phenotypes by transcriptomic profiling of porcine peripheral blood mononuclear cells (PBMCs) carrying different GR genotypes, in unstimulated conditions or in response to DEX and/or LPS in vitro. The Val allele differentially regulated abunda+nt genes in an additive-genetic manner. A subset of more than 200 genes was consistently affected by the substitution across treatments. This was associated with upregulation of genes related i.a. to endo-lysosomal system, lipid and protein catabolism, and immune terms including platelet activation, and antigen presentation, while downregulated genes were mainly involved in cell cycle regulation. Most importantly, the set of genes constitutively upregulated by Val includes members of the TLR4/LPS signaling pathway, such as LY96. Consequently, when exposing PBMCs to LPS treatment, the Val variant upregulated a panel of additional genes related to TLR4 and several other pattern recognition receptors, as well as cell death and lymphocyte signaling, ultimately amplifying the inflammatory responses. In contrast, when stimulated by DEX treatment, the Val allele orchestrated several genes involved in anti-inflammatory responses during infection. This study provides novel insights into the impact of GR hypersensitivity on the fate and function of immune cells, which may be useful for endotoxemia therapy.

DNA methylation landscapes from pig's limbic structures underline regulatory mechanisms relevant for brain plasticity.

Epigenetic dynamics are essential for reconciling stress-induced responses in neuro-endocrine routes between the limbic brain and adrenal gland. CpG methylation associates with the initiation and end of regulatory mechanisms underlying responses critical for survival, and learning. Using Reduced Representation Bisulfite Sequencing, we identified methylation changes of functional relevance for mediating tissue-specific responses in the hippocampus, amygdala, hypothalamus, and adrenal gland in pigs. We identified 4186 differentially methylated CpGs across all tissues, remarkably, enriched for promoters of transcription factors (TFs) of the homeo domain and zinc finger classes. We also detected 5190 differentially methylated regions (DMRs, 748 Mb), with about half unique to a single pairwise. Two structures, the hypothalamus and the hippocampus, displayed 860 unique brain-DMRs, with many linked to regulation of chromatin, nervous development, neurogenesis, and cell-to-cell communication. TF binding motifs for TFAP2A and TFAP2C are enriched amount DMRs on promoters of other TFs, suggesting their role as master regulators, especially for pathways essential in long-term brain plasticity, memory, and stress responses. Our results reveal sets of TF that, together with CpG methylation, may serve as regulatory switches to modulate limbic brain plasticity and brain-specific molecular genetics in pigs.

Prenatal transcript levels and metabolomics analyses reveal metabolic changes associated with intrauterine growth restriction and sex.

The metabolic changes associated with intrauterine growth restriction (IUGR) particularly affect the liver, which is a central metabolic organ and contributes significantly to the provision of energy and specific nutrients and metabolites. Therefore, our aim was to decipher and elucidate the molecular pathways of developmental processes mediated by miRNAs and mRNAs, as well as the metabolome in fetal liver tissue in IUGR compared to appropriate for gestational age groups (AGA). Discordant siblings representing the extremes in fetal weight at day 63 post conception (dpc) were selected from F2 fetuses of a cross of German Landrace and Pietrain. Most of the changes in the liver of IUGR at midgestation involved various lipid metabolic pathways, both on transcript and metabolite levels, especially in the category of sphingolipids and phospholipids. Differentially expressed miRNAs, such as miR-34a, and their differentially expressed mRNA targets were identified. Sex-specific phenomena were observed at both the transcript and metabolite levels, particularly in male. This suggests that sex-specific adaptations in the metabolic system occur in the liver during midgestation (63 dpc). Our multi-omics network analysis reveals interactions and changes in the metabolic system associated with IUGR and identified an important biosignature that differs between IUGR and AGA piglets.

Inhibition of KRAS, MEK and PI3K Demonstrate Synergistic Anti-Tumor Effects in Pancreatic Ductal Adenocarcinoma Cell Lines.

Kirsten rat sarcoma virus (KRAS) mutations are widespread in pancreatic ductal adenocarcinoma (PDAC) and contribute significantly to tumor initiation, progression, tumor relapse/resistance, and prognosis of patients. Although inhibitors against KRAS mutations have been developed, this therapeutic approach is not routinely used in PDAC patients. We investigated the anti-tumor efficacy of two KRAS inhibitors BI-3406 (KRAS::SOS1 inhibitor) and sotorasib (KRAS G12C inhibitor) alone or in combination with MEK1/2 inhibitor trametinib and/or PI3K inhibitor buparlisib in seven PDAC cell lines. Whole transcriptomic analysis of combined inhibition and control groups were comparatively analyzed to explore the corresponding mechanisms of inhibitor combination. Both KRAS inhibitors and corresponding combinations exhibited cytotoxicity against specific PDAC cell lines. BI-3406 enhance the efficacy of trametinib and buparlisib in BXPC-3, ASPC-1 and MIA PACA-2, but not in CAPAN-1, while sotorasib enhances the efficacy of trametinib and buparlisib only in MIA PACA-2. The whole transcriptomic analysis demonstrates that the two triple-inhibitor combinations exert antitumor effects by affecting related cell functions, such as affecting the immune system, cell adhesion, cell migration, and cytokine binding. As well as directly involved in RAF/MEK/ERK pathway and PI3K/AKT pathway affect cell survival. Our current study confirmed inhibition of KRAS and its downstream pathways as a potential novel therapy for PDAC and provides fundamental data for in vivo evaluations.

Allele-biased expression of the bovine APOB gene associated with the cholesterol deficiency defect suggests cis-regulatory enhancer effects of the LTR retrotransposon insertion.

The insertion of an endogenous retroviral long terminal repeat (LTR) sequence into the bovine apolipoprotein B (APOB) gene is causal to the inherited genetic defect cholesterol deficiency (CD) observed in neonatal and young calves. Affected calves suffer from developmental abnormalities, symptoms of incurable diarrhoea and often die within weeks to a few months after birth. Neither the detailed effects of the LTR insertion on APOB expression profile nor the specific mode of inheritance nor detailed phenotypic consequences of the mutation are undisputed. In our study, we analysed German Holstein dairy heifers at the peak of hepatic metabolic load and exposed to an additional pathogen challenge for clinical, metabolic and hepatic transcriptome differences between wild type (CDF) and heterozygote carriers of the mutation (CDC). Our data revealed that a divergent allele-biased expression pattern of the APOB gene in heterozygous CDC animals leads to a tenfold higher expression of exons upstream and a decreased expression of exons downstream of the LTR insertion compared to expression levels of CDF animals. This expression pattern could be a result of enhancer activity induced by the LTR insertion, in addition to a previously reported artificial polyadenylation signal. Thus, our data support a regulatory potential of mobile element insertions. With regard to the phenotype generated by the LTR insertion, heterozygote CDC carriers display significantly differential hepatic expression of genes involved in cholesterol biosynthesis and lipid metabolism. Phenotypically, CDC carriers show a significantly affected lipomobilization compared to wild type animals. These results reject a completely recessive mode of inheritance for the CD defect, which should be considered for selection decisions in the affected population. Exemplarily, our results illustrate the regulatory impact of mobile element insertions not only on specific host target gene expression but also on global transcriptome profiles with subsequent biological, functional and phenotypic consequences in a natural in-vivo model of a non-model mammalian organism.

Genetic regulation and variation of expression of miRNA and mRNA transcripts in fetal muscle tissue in the context of sex, dam and variable fetal weight.

BACKGROUND: Impaired skeletal muscle growth in utero can result in reduced birth weight and pathogenesis of intrauterine growth restriction. Fetal and placental growth is influenced by many factors including genetic, epigenetic and environmental factors. In fact, the sex and genotype of the fetus itself, as well as the mother providing it with a suitable environment, influence the growth of the fetus. Hence, our goal was to decipher and elucidate the molecular pathways of developmental processes mediated by miRNAs and mRNAs in fetal muscle tissue in the context of sex, dam, and fetal weight. Therefore, we analyse the variation of miRNA and mRNA expression in relation to these factors. In addition, the coincidence of genetic regulation of these mRNAs and miRNAs, as revealed by expression quantitative trait loci (eQTL) analyses, with sex-, mother- and weight-associated expression was investigated. METHODS: A three-generation pig F2 population (n = 118) based on reciprocal crossing of German Landrace (DL) and Pietrain (Pi) was used. Genotype information and transcriptomic data (mRNA and miRNA) from longissimus dorsi muscle (LDM) of pig fetuses sampled at 63 days post-conception (dpc) were used for eQTL analyses. RESULTS: The transcript abundances of 13, 853, and 275 probe-sets were influenced by sex, dam and fetal weight at 63 dpc, respectively (FDR < 5%). Most of significant transcripts affected by sex were located on the sex chromosomes including KDM6A and ANOS1 or autosomes including ANKS1B, LOC100155138 and miR-153. The fetal muscle transcripts associated with fetal weight indicated clearer metabolic directions than maternally influenced fetal muscle transcripts. Moreover, coincidence of genetic regulation (eQTL) and variation in transcript abundance due to sex, dam and fetal weight were identified. CONCLUSIONS: Integrating information on eQTL, sex-, dam- and weight-associated differential expression and QTL for fetal weight allowed us to identify molecular pathways and shed light on the basic biological processes associated with differential muscle development in males and females, with implications for adaptive fetal programming.

Multi-Omics Reveals Different Strategies in the Immune and Metabolic Systems of High-Yielding Strains of Laying Hens.

Lohmann Brown (LB) and Lohmann Selected Leghorn (LSL) are two commercially important laying hen strains due to their high egg production and excellent commercial suitability. The present study integrated multiple data sets along the genotype-phenotype map to better understand how the genetic background of the two strains influences their molecular pathways. In total, 71 individuals were analyzed (LB, n = 36; LSL, n = 35). Data sets include gut miRNA and mRNA transcriptome data, microbiota composition, immune cells, inositol phosphate metabolites, minerals, and hormones from different organs of the two hen strains. All complex data sets were pre-processed, normalized, and compatible with the mixOmics platform. The most discriminant features between two laying strains included 20 miRNAs, 20 mRNAs, 16 immune cells, 10 microbes, 11 phenotypic traits, and 16 metabolites. The expression of specific miRNAs and the abundance of immune cell types were related to the enrichment of immune pathways in the LSL strain. In contrast, more microbial taxa specific to the LB strain were identified, and the abundance of certain microbes strongly correlated with host gut transcripts enriched in immunological and metabolic pathways. Our findings indicate that both strains employ distinct inherent strategies to acquire and maintain their immune and metabolic systems under high-performance conditions. In addition, the study provides a new perspective on a view of the functional biodiversity that emerges during strain selection and contributes to the understanding of the role of host-gut interaction, including immune phenotype, microbiota, gut transcriptome, and metabolome.

Brain Transcriptome Responses to Dexamethasone Depending on Dose and Sex Reveal Factors Contributing to Sex-Specific Vulnerability to Stress-Induced Disorders.

BACKGROUND: Glucocorticoid (GC) receptor (GR) signaling in the hypothalamus (Hyp) and in the superordinate limbic structures, such as the hippocampus (Hip), conveys feedback regulation of the neuroendocrine stress response and acts upon other neurobiological functions that ultimately influence mental health. These responses are strongly influenced by sex, but the molecular causes are still largely unexplored. METHODS: To investigate GR targets and their GC sensitivity in the Hyp and Hip, we treated juvenile male and female piglets with 10 (D10) or 60 (D60) µg/kg dexamethasone (DEX), a selective GR agonist, and analyzed transcriptome responses compared to a saline control group using RNA sequencing. RESULTS: Both doses influenced similar biological functions, including cellular response to lipid and immune cell-related functions, but the transcriptional response to D10 was considerably weaker, particularly in the Hip. Weighted Gene Co-expression Network Analysis revealed a network of genes coordinately regulated by DEX in both structures, among which the alpha-arrestin ARRDC2 takes a central position. Distinct functional groups of genes were differentially regulated by DEX between sexes depending on the dose; at D10, these included particularly mitochondrial genes, whereas at D60 interferon signaling and lipid homeostasis genes were enriched. The general and sex-specific transcriptional responses to DEX highlight microglia as the prominent target. Several key marker genes of disease-associated microglia were regulated by DEX depending on sex, such as TREM2 and LPL. CONCLUSION: The discovered expression signatures suggest that DEX induced a dysfunctional state of microglia in males, while in females microglia were primed, which could entail predisposition for different mental disorders.

Erratum: Ali et al. In Utero Fetal Weight in Pigs Is Regulated by microRNAs and Their Target Genes. Genes 2021, 12, 1264.

The authors wish to make the following correction to their paper published in Genes [...].

Multi-Transcript Level Profiling Revealed Distinct mRNA, miRNA, and tRNA-Derived Fragment Bio-Signatures for Coping Behavior Linked Haplotypes in HPA Axis and Limbic System.

The molecular basis of porcine coping behavior (CB) relies on a sophisticated interplay of genetic and epigenetic features. Deep sequencing technologies allowed the identification of a plethora of new regulatory small non-coding RNA (sncRNA). We characterized mRNA and sncRNA profiles of central parts of the physiological stress response system including amygdala, hippocampus, hypothalamus and adrenal gland using systems biology for integration. Therefore, ten each of high- (HR) and low- (LR) reactive pigs (n = 20) carrying a CB associated haplotype in a prominent QTL-region on SSC12 were selected for mRNA and sncRNA expression profiling. The molecular markers related to the LR group included ATP1B2, MPDU1, miR-19b-5p, let-7g-5p, and 5'-tiRNA Leu in the adrenal gland, miR-194a-5p, miR-125a-5p, miR-7-1-5p, and miR-107-5p in the hippocampus and CBL and PVRL1 in the hypothalamus. Interestingly, amygdalae of the LR group showed 5'-tiRNA and 5'-tRF (5'-tRF Lys , 5'-tiRNA Lys , 5'-tiRNA Cys , and 5'-tiRNA Gln ) enrichment. Contrarily, molecular markers associated with the HR group encompassed miR-26b-5p, tRNA Arg , tRNA GlyiF in the adrenal gland, IGF1 and APOD in the amygdala and PBX1, TOB1, and C18orf1 in the hippocampus and miR-24 in the hypothalamus. In addition, hypothalami of the HR group were characterized by 3'-tiRNA enrichment (3'-tiRNAGln, 3'-tiRNA Asn , 3'-tiRNA Val , 3'-tRF Pro , 3'-tiRNA Cys , and 3'-tiRNA Ala ) and 3'-tRFs enrichment (3'-tRF Asn , 3'-tRF Glu , and 3'-tRF Val ). These evidence suggest that tRNA-derived fragments and their cleavage activity are a specific marker for coping behavior. Data integration revealed new bio-signatures of important molecular interactions on a multi-transcript level in HPA axis and limbic system of pigs carrying a CB-associated haplotype.

In Utero Fetal Weight in Pigs Is Regulated by microRNAs and Their Target Genes.

Impaired skeletal muscle growth in utero can result in reduced birth weight and poor carcass quality in pigs. Recently, we showed the role of microRNAs (miRNAs) and their target genes in prenatal skeletal muscle development and pathogenesis of intrauterine growth restriction (IUGR). In this study, we performed an integrative miRNA-mRNA transcriptomic analysis in longissimus dorsi muscle (LDM) of pig fetuses at 63 days post conception (dpc) to identify miRNAs and genes correlated to fetal weight. We found 13 miRNAs in LDM significantly correlated to fetal weight, including miR-140, miR-186, miR-101, miR-15, miR-24, miR-29, miR-449, miR-27, miR-142, miR-99, miR-181, miR-199, and miR-210. The expression of these miRNAs decreased with an increase in fetal weight. We also identified 1315 genes significantly correlated to fetal weight at 63 dpc, of which 135 genes were negatively correlated as well as identified as potential targets of the above-listed 13 miRNAs. These miRNAs and their target genes enriched pathways and biological processes important for fetal growth, development, and metabolism. These results indicate that the transcriptomic profile of skeletal muscle can be used to predict fetal weight, and miRNAs correlated to fetal weight can serve as potential biomarkers of prenatal fetal health and growth.

Single-cell RNA sequencing of freshly isolated bovine milk cells and cultured primary mammary epithelial cells.

Bovine mammary function at molecular level is often studied using mammary tissue or primary bovine mammary epithelial cells (pbMECs). However, bulk tissue and primary cells are heterogeneous with respect to cell populations, adding further transcriptional variation in addition to genetic background. Thus, understanding of the variation in gene expression profiles of cell populations and their effect on function are limited. To investigate the mononuclear cell composition in bovine milk, we analyzed a single-cell suspension from a milk sample. Additionally, we harvested cultured pbMECs to characterize gene expression in a homogeneous cell population. Using the Drop-seq technology, we generated single-cell RNA datasets of somatic milk cells and pbMECs. The final datasets after quality control filtering contained 7,119 and 10,549 cells, respectively. The pbMECs formed 14 indefinite clusters displaying intrapopulation heterogeneity, whereas the milk cells formed 14 more distinct clusters. Our datasets constitute a molecular cell atlas that provides a basis for future studies of milk cell composition and gene expression, and could serve as reference datasets for milk cell analysis.

Transcriptome analysis of porcine PBMCs reveals lipopolysaccharide-induced immunomodulatory responses and crosstalk of immune and glucocorticoid receptor signaling.

The current level of knowledge on transcriptome responses triggered by endotoxins and glucocorticoids in immune cells in pigs is limited. Therefore, in the present study, we treated porcine peripheral blood mononuclear cells (PBMCs) with lipopolysaccharide (LPS) and dexamethasone (DEX) separately or combined for 2 hours. The resultant transcriptional responses were examined by mRNA sequencing. We found that the LPS treatment triggered pronounced inflammatory responses as evidenced by upregulation of pro-inflammatory cytokines, chemokines, and related signaling pathways like NF-κB. Concurrently, a series of downregulated pro-inflammatory and upregulated anti-inflammatory molecules were identified. These are involved in the inhibition of TLR, NF-κB, and MAPK cascades and activation of signaling mediated by Tregs and STAT3, respectively. These findings suggested that LPS initiated also an anti-inflammatory process to prevent an overwhelming inflammatory response. The transcriptome responses further revealed substantial crosstalk of immune responses and glucocorticoid receptor (GR) signaling. This was apparent in four aspects: constitutive inhibition of T cell signaling by DEX through a subset of genes showing no response to LPS; inhibition of LPS-induced inflammatory genes by DEX; attenuation of DEX action by LPS paralleled by the regulation of genes implicated in cytokine and calcium signaling; and DEX-induced changes in genes associated with the activation of pro-inflammatory TLR, NF-κB, iNOS, and IL-1 signaling. Consequently, our study provides novel insights into inflammatory and GR signaling in pigs, as well as an understanding of the application of glucocorticoid drugs for the treatment of inflammatory disorders.

Dietary phosphorus and calcium in feed affects miRNA profiles and their mRNA targets in jejunum of two strains of laying hens.

Phosphorus (P) and calcium (Ca) are critical for egg production in laying hens. Most of P in plant-based poultry diet is bound as phytic acid and needs to be hydrolysed before absorption. To increase P bioavailability, exogenous phytases or bioavailable rock phosphate is added in feed. There is growing evidence of the importance of miRNAs as the epicentre of intestinal homeostasis and functional properties. Therefore, we demonstrated the expression of miRNA profiles and the corresponding target genes due to the different levels of P (recommended vs. 20% reduction) and/or Ca (recommended vs. 15% reduction) in feed. Jejunal miRNA profiles of Lohmann Selected Leghorn (LSL) and Lohmann Brown (LB) laying hens strains were used (n = 80). A total of 34 and 76 miRNAs were differentially expressed (DE) in the different diet groups within LSL and LB strains respectively. In LSL, the DE miRNAs and their targets were involved in calcium signaling pathway, inositol phosphate metabolism, and mitochondrial dysfunction. Similarly, in LB miRNAs targets were enriched in metabolic pathways such as glutathione metabolism, phosphonate metabolism and vitamin B6 metabolism. Our results suggest that both strains employ different intrinsic strategies to cope with modulated P and Ca supply and maintain mineral homeostasis.