<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Genet.</journal-id>
<journal-title>Frontiers in Genetics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Genet.</abbrev-journal-title>
<issn pub-type="epub">1664-8021</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">886875</article-id>
<article-id pub-id-type="doi">10.3389/fgene.2022.886875</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Genetics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Early-Life Exposure of Pigs to Topsoil Alters miRNA and mRNA Expression in Peripheral Blood Mononuclear Cells</article-title>
<alt-title alt-title-type="left-running-head">De Souza et al.</alt-title>
<alt-title alt-title-type="right-running-head">Topsoil Alters Piglet Gene Expression</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>De Souza</surname>
<given-names>M. M.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/630581/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Koltes</surname>
<given-names>D. A.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/473604/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Beiki</surname>
<given-names>H.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1314113/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Sales</surname>
<given-names>M. A.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/611179/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tsai</surname>
<given-names>T.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/597689/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Maxwell</surname>
<given-names>C. V.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/636769/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhao</surname>
<given-names>J.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/411174/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Koltes</surname>
<given-names>J. E.</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/227550/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Animal Science</institution>, <institution>Iowa State University</institution>, <addr-line>Ames</addr-line>, <addr-line>IA</addr-line>, <country>United States</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Animal Science</institution>, <institution>University of Arkansas-Division of Agriculture</institution>, <addr-line>Fayetteville</addr-line>, <addr-line>AR</addr-line>, <country>United States</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/214280/overview">Eveline M. Ibeagha-Awemu</ext-link>, Agriculture and Agri-Food Canada, Canada</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/637818/overview">Yaoyao Xia</ext-link>, South China Agricultural University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1189625/overview">Wang Bing</ext-link>, Tianjin Medical University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: J. E. Koltes, <email>jekoltes@iastate.edu</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Livestock Genomics, a section of the journal Frontiers in Genetics</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>08</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>13</volume>
<elocation-id>886875</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>03</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>06</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 De Souza, Koltes, Beiki, Sales, Tsai, Maxwell, Zhao and Koltes.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>De Souza, Koltes, Beiki, Sales, Tsai, Maxwell, Zhao and Koltes</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Exposure to less-hygienic conditions during early childhood has been associated with stimulation and development of the immune system. A recent study indicated that exposure of piglets to soil-borne microbes during lactation was related with modulation of gut microbiota and immune function. To identify the potential molecular mechanisms and pathways impacted by early-life topsoil exposure, we analyzed the messenger RNA (mRNA) and micro-RNA (miRNA) expression in peripheral blood mononuclear cells (PBMCs) from these piglets. Total RNA was extracted from the PBMCs of piglets exposed to topsoil only from d 4&#x2013;d 21 of life (mRNA <italic>n</italic> &#x3d; 6; miRNA <italic>n</italic> &#x3d; 5) or unexposed control pigs (mRNA <italic>n</italic> &#x3d; 6; miRNA <italic>n</italic> &#x3d; 8) at 11, 20, and 56&#xa0;days of age. Small RNA and mRNA were sequenced with 50-bp single-end reads using Illumina chemistry. Sequence data were quality checked with FASTQC software and aligned to the Sscrofa 11.1 genome with the STAR aligner for mRNA and mirDeep2 for miRNA. Differential expression (DE) analysis was performed using PROC Glimmix of SAS to evaluate changes in expression due to topsoil exposure over time with genes declared DE at a false discovery rate (FDR) of q &#x3c; 0.10. A total of 138 mRNA and 21 miRNAs were identified as DE for the treatment by age interaction. Ontology enrichment analysis of DE mRNA revealed Gene ontology (GO) terms directly involved in the connection between T-cell and antigen-presenting cells that are associated with T-cell activation. Key regulatory genes identified include <italic>PTPRJ, ITGB3, TRBV30, CD3D</italic>, mir-143, mir-29, and mir-148a. While these results require validation, this study provides data supporting the hypothesis that less-hygienic environments during early life may contribute to the development of the immune system.</p>
</abstract>
<kwd-group>
<kwd>immune system</kwd>
<kwd>peripheral blood mononuclear cells</kwd>
<kwd>livestock</kwd>
<kwd>microbiota</kwd>
<kwd>RNA sequencing</kwd>
<kwd>hygiene hypothesis</kwd>
</kwd-group>
<contract-sponsor id="cn001">Arkansas Biosciences Institute<named-content content-type="fundref-id">10.13039/100008231</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">National Institute of Food and Agriculture<named-content content-type="fundref-id">10.13039/100005825</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>Disease resistance of individuals exposed to less-hygienic conditions during early life has been called the hygiene hypothesis, first postulated by <xref ref-type="bibr" rid="B39">Strachan (1989)</xref>. It proposes that individuals exposed to less-hygienic conditions early in life may have altered immune system development which lowers the incidence of illness later in life (<xref ref-type="bibr" rid="B39">Strachan, 1989</xref>). While it is well understood that particular conditions, such as poor water sanitation, lead to the development and propagation of diseases, low environmental exposure has been associated with increased resistance to human chronic illnesses, such as asthma and allergies (<xref ref-type="bibr" rid="B21">Illi et al., 2001</xref>; <xref ref-type="bibr" rid="B30">MacNeill et al., 2013</xref>; <xref ref-type="bibr" rid="B19">Huang et al., 2021</xref>). It has been theorized that this low dose of exposure results in a mild stimulation of the immune system, particularly the adaptive immune system (<xref ref-type="bibr" rid="B40">Su et al., 2013</xref>).</p>
<p>A growing number of publications indicate that microbe exposure may be a significant contributor to the development of the immune system (reviewed by <xref ref-type="bibr" rid="B49">Zheng et al., 2020</xref>). <xref ref-type="bibr" rid="B14">Ennamorati et al. (2020)</xref> identified the gut microbiome as a major presenter of novel antigens to the immune system that influences the generation of regulatory T cells by the thymus in neonatal mice. The authors described how this process may facilitate homeostasis between microbiota and the host immune system in the host to prevent intestinal inflammation later in life. Studies in the chicken also indicate that intestinal microbiota may regulate T-cell production in the thymus in early life (<xref ref-type="bibr" rid="B6">Cheng et al., 2021</xref>). Moreover, postnatal microbial colonization of the body may act as a major stimulus of immune tissue development and modulation of miRNA expression (<xref ref-type="bibr" rid="B33">Masotti, 2012</xref>), which contributes to &#x201c;training&#x201d; of the immune system (<xref ref-type="bibr" rid="B18">Gury-BenAri et al., 2016</xref>). This early training of the immune system and development of microbial memory in early life may allow individuals to overcome an illness more quickly, thus, allowing them to spare energy to partition toward growth. In addition, development of a more robust immune system should reduce antibiotic use and subsequently antimicrobial resistance.</p>
<p>As livestock producers have strived to increase sustainability and growth efficiency of pigs, the environment and hygienic conditions have significantly changed over time. During early commercialized swine production, swine were reared outdoors and exposed to soil, weather, and predators. The move to modern, indoor housing has allowed for decreased exposure to parasites, fecal material, predators, and changes in weather conditions, which greatly improved the overall welfare of swine. It has been observed that early-life exposure of pigs to less-hygienic conditions (i.e., topsoil exposure) has increased growth (<xref ref-type="bibr" rid="B12">Dourmad et al., 2009</xref>; <xref ref-type="bibr" rid="B26">Lebret et al., 2011</xref>; <xref ref-type="bibr" rid="B46">Wenner et al., 2013</xref>; <xref ref-type="bibr" rid="B44">Vo et al., 2017</xref>) while decreasing feed intake (<xref ref-type="bibr" rid="B26">Lebret et al., 2011</xref>; <xref ref-type="bibr" rid="B46">Wenner et al., 2013</xref>). One explanation for these findings is that host exposure to microbes from the environment may influence host physiology.</p>
<p>Initial studies of neonatal pigs exposed to topsoil indicate microbial communities (<xref ref-type="bibr" rid="B44">Vo et al., 2017</xref>; <xref ref-type="bibr" rid="B45">Wen et al., 2021</xref>) and the immune cell population (<xref ref-type="bibr" rid="B41">Tsai et al., 2016</xref>; <xref ref-type="bibr" rid="B45">Wen et al., 2021</xref>) are altered later in life in response to early-life topsoil exposure. Both <xref ref-type="bibr" rid="B44">Vo et al. (2017)</xref> and <xref ref-type="bibr" rid="B45">Wen et al. (2021)</xref> identified higher levels of <italic>Prevotella</italic> at days 13 and 12, respectively, in the gut of piglets exposed to topsoil. In mice, perturbation of <italic>Prevotella</italic> in the gut microbiome resulted in exacerbated intestinal inflammation (<xref ref-type="bibr" rid="B20">Iljazovic et al., 2021</xref>). Moreover, <xref ref-type="bibr" rid="B44">Vo et al. (2017)</xref>, who used the same animals described in this study, identified a higher abundance of <italic>Prevotella</italic> at weaning (d 21) and at day 35 (mid-nursery). <xref ref-type="bibr" rid="B44">Vo et al. (2017)</xref> also identified higher levels of <italic>Bacteroidetes</italic> and Ruminococcaceae at day 13 and <italic>Clostridium</italic> and <italic>Clostridium cluster XI</italic>, <italic>Coprococcus</italic>, <italic>Campylobacter</italic>, <italic>Streptococcus</italic>, <italic>Mitsuokella</italic>, <italic>Ruminococcus</italic>, <italic>Dialister</italic>, and <italic>Sarcina</italic> during postweaning at the end of the nursery phase (d 56) in response to topsoil exposure. Interestingly, <italic>Clostridium</italic> species have been associated with colonic regulatory T-cell accumulation and improvement of the local and systemic immune responses (<xref ref-type="bibr" rid="B2">Atarashi et al., 2011</xref>). Although these studies demonstrate that the microbiome is altered in pigs exposed to topsoil, it is still unclear which specific microbes are directly impacting the animal&#x2019;s physiology. Little is known about the molecular mechanisms responsible for impacting animal growth or the immune system in this model. Understanding how the immune system is programmed in response to microbes and the environment is an important step in the development of technologies (e.g., probiotics) that may enhance efficiency in livestock production, and provide basic science knowledge about host&#x2013;microbiome interactions that are beneficial to improve human health.</p>
<p>The objective of this study was to identify changes in peripheral blood mononuclear cell (PBMC) gene expression and associated pathways in response to early-life topsoil exposure in piglets from a subset of piglets exposed to topsoil from previous studies (<xref ref-type="bibr" rid="B41">Tsai et al., 2016</xref>; <xref ref-type="bibr" rid="B44">Vo et al., 2017</xref>). Specifically, we sought to determine if there was evidence of early, enhanced, or altered immune system development based on mRNA and miRNA expression profiles using RNA-sequencing (RNA-seq) technology. Transcriptome profiling is a powerful tool to evaluate global changes in gene expression within a tissue. As miRNAs are well-established regulators of gene expression and may facilitate communication between the host and microbiome, we evaluated the relationship between miRNAs and known or predicted mRNA targets (<xref ref-type="bibr" rid="B33">Masotti, 2012</xref>; <xref ref-type="bibr" rid="B47">Yuan et al., 2018</xref>; <xref ref-type="bibr" rid="B17">Gebert and MacRae, 2019</xref>).</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>Materials and Methods</title>
<sec id="s2-1">
<title>Animal Population and Topsoil Exposure Protocol</title>
<p>A detailed description of the animal trial can be found in <xref ref-type="bibr" rid="B44">Vo et al. (2017)</xref>. Briefly, animal experiments were conducted at the University of Arkansas-Division of Agriculture Swine Research Unit in Fayetteville, Arkansas under protocols approved by the University of Arkansas Institute of Animal Care and Use Committee &#x23;13060. A total of 14 synthetic large white pigs were used for the experiment. Control piglets (control group; <italic>n</italic> &#x3d; 6 for mRNA and <italic>n</italic> &#x3d; 8 for miRNA analysis) were reared conventionally in farrowing crates, whereas treatment piglets (topsoil group; <italic>n</italic> &#x3d; 6 for mRNA and <italic>n</italic> &#x3d; 5 for miRNA analysis) were exposed daily to fresh topsoil from day 4 postpartum until the end of lactation (d 20). All piglets were weaned and transferred to a nursery facility without soil at d 20. Piglets were weighed and blood samples were collected at 11, 20, and 56&#xa0;days of age. Additional information about the individual pigs (e.g., sex, weight) is provided in <xref ref-type="sec" rid="s12">Supplementary Tables S1, S2</xref> for mRNA and miRNA experiments, respectively.</p>
</sec>
<sec id="s2-2">
<title>Blood Collection and Isolation of Peripheral Blood Mononuclear Cells</title>
<p>Whole blood samples (10&#xa0;ml) were collected through venipuncture via jugular vena cava into EDTA tubes (BD Vacutainer<sup>&#xae;</sup> K2EDTA-coated blood collection tube, Franklin Lakes, NJ) from piglets with median body weight (BW) in each litter at days 11, 20 (weaning), and 56 (the end of nursery) of age. PBMCs were isolated using Histopaque&#xae;-1,077 (MilliporeSigma, Burlington, MA) gradient centrifugation according to the manufacturer&#x2019;s protocol and stored at -80&#xb0;C until total RNA extraction.</p>
</sec>
<sec id="s2-3">
<title>RNA Isolation and RNA-Sequencing Analysis</title>
<p>Total RNA was extracted with TRIzol<sup>&#xae;</sup> using the Direct-Zol RNA Isolation Kit (Zymo Research Corp.), according to the manufacturer&#x2019;s protocol. RNA quality was determined using an Experion RNA analysis kit (Biorad, Hurcules, CA) to determine RNA degradation using the RNA quality index (RQI) and RNA quantity was measured by fluorometer using a Qubit (ThermoFisher, Waltham, MA). Total RNA samples that passed quantity (100 ng/uL) and quality criteria (RNA quality indicator, RQI &#x3d; 5.1&#x2013;9.8) were utilized for sequencing. The RNA-Seq analysis was conducted using Illumina HiSeq4000 chemistry at the Beijing Genomics Institute (BGI) Americas (BGI, Shenzhen, Guangdong, China). Small RNA and mRNA libraries were constructed from the total RNA samples for 50-bp single-end sequencing.</p>
</sec>
<sec id="s2-4">
<title>Bioinformatics Analyses</title>
<p>Sequence read quality was evaluated with FastQC version0.11.3 (<xref ref-type="bibr" rid="B1">Andrews and Babraham Bioinformatics, 2010</xref>) and adaptors and low-quality reads were trimmed using Trim Galore version 0.4.5 (<xref ref-type="bibr" rid="B24">Krueger, 2015</xref>). Sequence reads from mRNA libraries were aligned to the Sscrofa 11.1 reference genome using the STAR alignment software version 2.5.4a (<xref ref-type="bibr" rid="B11">Dobin et al., 2013</xref>). The number of read counts for each transcript was quantified with FeatureCounts version 2.0.3 (<xref ref-type="bibr" rid="B27">Liao et al., 2014</xref>). To determine miRNA expression levels, sequence reads were aligned to the Sscrofa11.1 reference genome using miRDeep2 (<xref ref-type="bibr" rid="B15">Friedl&#xe4;nder et al., 2008</xref>).</p>
</sec>
<sec id="s2-5">
<title>Differential Expression Analysis</title>
<p>Prior to analysis, data were filtered to remove unexpressed (i.e., genes with 0 read counts) and lowly expressed genes. Lowly expressed genes were defined as those not present in 24 of 36 mRNA or 17 of 39 miRNA samples, and with an average expression &#x3e;2 reads/sample. To account for differences in library size, all samples were normalized using the 75% Quantile method (<xref ref-type="bibr" rid="B5">Bullard et al., 2010</xref>). Differential expression analysis was performed with PROC Glimmix of SAS (SAS <xref ref-type="bibr" rid="B22">Institute Inc., 2013</xref>) using a negative binomial distribution to fit the model:<disp-formula id="equ1">
<mml:math id="m1">
<mml:mrow> <mml:msub>
<mml:mi>y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:msub>
<mml:mi>t</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>R</mml:mi>
<mml:mi>Q</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>g</mml:mi>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>l</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>z</mml:mi>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mi>x</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>a</mml:mi>
<mml:mi>g</mml:mi>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>l</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>z</mml:mi>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>where y &#x3d; normalized expression value for each sample k; treatment i &#x3d; topsoil or control; lane j &#x3d; sequencing lane 1, 2 or 3; RQI &#x3d; RNA quality score for each sample k from each individual l; parity m &#x3d; parity of the piglet&#x2019;s sow (1&#x2013;4); age &#x3d; age of piglet at sample k &#x3d; d 11, d 20, or d 56; sex n &#x3d; male or female, litter size &#x3d; number of piglets within piglet&#x2019;s litter (range: 5&#x2013;17), and e &#x3d; random residual error. In addition, pig was included as a random effect to account for the repeated measure across age. The same model was fit for miRNA, except that only one lane was used for sequencing, so no lane effect was needed. For both the mRNA and miRNA analyses, we focused on identifying genes that changed in expression overtime in response to early-life topsoil exposure (i.e., treatment <italic>vs</italic>. age interaction). To account for multiple testing, q-values were calculated in R using the q-value package (<xref ref-type="bibr" rid="B38">Storey et al., 2019</xref>) to determine the false discovery rate (FDR) for gene lists. Significant differentially expressed (DE) genes were declared at q &#x3c; 0.10.</p>
</sec>
<sec id="s2-6">
<title>Co-Expression Network Analysis</title>
<p>To identify pathways of genes that may be impacted by the early-life topsoil exposure treatment, co-expression network analysis was conducted. Co-expression networks were built separately for the control and topsoil mRNA and miRNA groups, respectively, using the Weighted Gene Co-expression Network Analysis (WGCNA) methods and associated R package (<xref ref-type="bibr" rid="B25">Langfelder and Horvath, 2008</xref>). The RNA-seq count data were normalized using the Trimmed Mean of M-values method (TMM) in the edgeR package (<xref ref-type="bibr" rid="B36">Robinson et al., 2009</xref>). Four network analyses were performed to compare the network modules (putative pathways) between control and topsoil treatments for both the mRNA and miRNA expression data. Spearman&#x2019;s correlation coefficients were used to obtain the pairwise similarity between gene expression profiles in a signed co-expression network. The soft threshold chosen was &#x3b2; &#x3d; 14 for all four analyses performed. The minimum module size chosen was 25 and 5 transcripts per module for mRNA and miRNA networks, respectively.</p>
</sec>
<sec id="s2-7">
<title>miRNA:mRNA Module Interactions</title>
<p>To determine the potential miRNA regulation of mRNA pathways, we identified miRNA modules statistically associated with mRNA modules. Spearman correlations were calculated between eigen values for mRNA and miRNA modules and correlation p-values were adjusted to control the FDR using the R package psych version 1.8.10 (<xref ref-type="bibr" rid="B35">Revelle, 2018</xref>). To discover individual mRNA or miRNA that may be regulating pathways, mRNAs and miRNA showing higher connectivity within a module were identified as possible hub genes based on module membership (MM) values. The MM values &#x3e;0.9 were selected as candidate hub (regulator) genes. We used TargetScan to retrieve the known target genes for each miRNA contained within the modules. Only the correlated target genes for the corresponding miRNA within in the correlated miRNA:mRNA module were kept for ontology enrichment analysis.</p>
</sec>
<sec id="s2-8">
<title>Functional Annotation (Ontology Enrichment) Analysis</title>
<p>Biological processes and molecular pathways that differ in response to early-life topsoil treatment were identified using gene ontology enrichment analysis. Functional analysis was carried out using the human orthologues genes since more complete annotation information is available for human genes. Biomart (<xref ref-type="bibr" rid="B13">Durinck et al., 2009</xref>) was used to retrieve <italic>Homo sapiens</italic> orthologues for all expressed genes, where orthologues were defined as DNA sequences with &#x3e;60% nucleotide similarity. Enriched pathways and biological processes were enriched for the DE gene list and miRNA target gene list using the Cytoscape plugins ClueGO v. 2.5.3 and CluePedia v. 1.5.3 (<xref ref-type="bibr" rid="B4">Bindea et al., 2009</xref>). The background (reference genes) used in the enrichment analysis were set as all expressed genes in PBMCs from our samples. Biological process gene ontology (GO) terms (GO tree level &#x2265;5) were considered significantly enriched within the gene list at a <italic>p</italic> &#x2264; 0.05 Benjamini and Hochberg corrected p-value. After identifying the gene-GO matrix, ClueGo clusters similar GO terms into functional groups based on the association strength between the ontology terms, using a kappa score (&#x2265;0.5). Group <italic>p</italic>-values were calculated using the number of unique genes found from the uploaded gene list that were statistically associated with the ontology terms included in the group and the total number of unique genes associated with those terms. <italic>p</italic>-value was adjusted using Benjamini and Hochberg (BH) (p-value &#x3c; 0.05). Enrichment of KEGG pathway terms from small RNA target genes was conducted using the DIANA-miRPath, v3.0 software (<xref ref-type="bibr" rid="B43">Vlachos et al., 2015</xref>). Predicted targets were identified using the DIANA-microT-CDS Targets of miRNAs analysis, with default setting (Micro T threshold &#x3d; 0.08, <italic>p</italic> &#x3c; 0.05) and enrichment analysis methods (Fisher&#x2019;s exact test).</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>Results</title>
<sec id="s3-1">
<title>Evaluation of Sequence Quality</title>
<p>A total of 16.5 million sequence reads were mapped, on average, per sample for the mRNA libraries to the Sscrofa 11.1 reference genome. A total of 13 million sequence reads were mapped, on average, per sample for the miRNA libraries. After applying zero and low expression filters, a total of 14,316 transcripts were available for differential expression analysis in SAS. All sequencing statistics are available in <xref ref-type="sec" rid="s12">Supplementary Tables S1, S2</xref>.</p>
</sec>
<sec id="s3-2">
<title>Differentially Expressed Genes (mRNA), Enriched Biological Process, and Pathway Annotations Identified in Response to Early-Life Topsoil Exposure</title>
<p>A total of 138 genes were identified as DE for the treatment by age interaction (q &#x3c; 0.10). The full list of DE genes and associated statistics are available in <xref ref-type="sec" rid="s12">Supplementary Table S3</xref>. A heat map displaying the expression levels of DE genes is provided in <xref ref-type="fig" rid="F1">Figure 1A</xref>. Trends in gene expression patterns observed over time when comparing the topsoil treatment to control for all DE genes are shown in <xref ref-type="sec" rid="s12">Supplementary Figures S1&#x2013;S6</xref> for all the DE genes. To better understand the function of the DE mRNAs, GO enrichment analysis was performed using the DE genes that had annotations for human orthologues (109 of 138 DE genes). A total of 59 GO terms were identified as statistically enriched for specific biological terms (p-adjusted &#x3c; 0.05; <xref ref-type="fig" rid="F1">Figure 1B</xref>, <xref ref-type="sec" rid="s12">Supplementary Table S4</xref>). Forty-three of these GO terms were clustered into 13 groups with similar GO term biological processes (BP) (<xref ref-type="fig" rid="F1">Figure 1B</xref>). The most statistically significant ontology term was platelet degranulation (adjusted <italic>p</italic> &#x3c; 0.01). The three largest clusters of enriched GO terms (Groups 1, 2, and 12; <xref ref-type="fig" rid="F1">Figure 1B</xref>) were involved in processes related to functionality or components of the immune system. Four additional gene ontology groups identified were also known to be involved in immune system, including the terms microtubule-based protein transport (group 4), integrin-mediated signaling pathways (group 8), positive regulation of peptidyl-tyrosine phosphorylation (group 9) and positive regulation of cell-substrate adhesion (group 13). The genes in these seven combined ontology term groups had higher levels of transcript abundance at day 20 in topsoil-treated animals (with the exception of <italic>IFT20</italic> and <italic>OPTN</italic>; <xref ref-type="fig" rid="F1">Figure 1C</xref>). A total of 15 out of 28 genes in these ontology groups exhibited higher levels of transcript abundance in the control piglets at days 11 and 56. Although the DE genes <italic>CD3D</italic> and <italic>TRBV30</italic> were not included in these ontology term groups, they also had higher expression in topsoil-treated piglets at day 20 in the same pattern observed in the seven ontology groups identified. Since <italic>CD3D</italic> and <italic>TRBV30</italic> play a major role during immune response as components of the T-cells receptor (TCR) complex, they are also highlighted in <xref ref-type="fig" rid="F1">Figure 1C</xref>. Other statistically significant ontology terms (adjusted <italic>p</italic> &#x3c; 0.05) included: negative regulation of insulin secretion, positive chemotaxis, BMP signaling, SMAD signaling, iron ion transport, and homeostasis and glucan metabolic process.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Gene expression and gene ontology (GO) enrichment analysis of the differentially expressed (DE) mRNA in response to topsoil treatment <bold>(A)</bold> A heatmap comparison of the expression of the DE mRNAs across time and treatment (control <italic>vs.</italic> topsoil) in peripheral blood mononuclear cells (PBMCs.). The three timepoints correspond to when pigs were still exposed to topsoil in the nursery (d 11) and at weaning (d 20) as well as after topsoil was removed at the postweaning (d56) timepoint. The heat map shows the clustered relative gene expression (least squares mean&#x2013;LSM) for control and topsoil treatments by age (day of life/age of piglet: 11, 20, and 56) <bold>(B)</bold> Visualization of the number and significance of grouped gene ontology (GO) terms identified using the DE genes in ClueGo software. The group p-values (provided on the <italic>x</italic>-axis) were calculated using the list of unique DE genes and adjusted using a Benjamini and Hochberg correction (BH). <italic>p</italic>-values were considered significant at BH adjusted p-value &#x3c; 0.05. The size of the bubbles represents the number of genes clustered in the ontology term group <bold>(C)</bold> A histogram showing the difference in transcript abundance between the control and topsoil treatment at days 11, 20, and 56 for selected genes with ontology terms related to functionality or components of the immune system. Differences in mRNA abundance were estimated using the LSM between control and topsoil treatments at each specific time point (piglet age). Positive differences (above zero) means higher abundance in the control group while negative differences (below zero) mean higher in topsoil group.</p>
</caption>
<graphic xlink:href="fgene-13-886875-g001.tif"/>
</fig>
</sec>
<sec id="s3-3">
<title>Differentially Expressed miRNA, Enriched Biological Process, and Pathway Annotations Identified in Response to Early-Life Topsoil Exposure</title>
<p>A total of 542 miRNAs were identified by miRDeep2 software. After filtering for annotations, a total of 254 miRNAs were used for the statistical analysis. Twenty-one unique miRNAs were identified as DE for the treatment by age interaction (q &#x3c; 0.10; miRNA DE results are listed in <xref ref-type="sec" rid="s12">Supplementary Table S5</xref>). A heat map displaying the expression levels of DE miRNA is provided in <xref ref-type="fig" rid="F2">Figure 2A</xref>. Trends in miRNA expression patterns observed over time when comparing treatments are shown in <xref ref-type="sec" rid="s12">Supplementary Table S7</xref>. The five most significantly DE miRNA included miR-1306, miR-660, miR-126, miR-1285, and miR-1468. Notably, six of the miRNAs identified as DE had gene annotations related to functionality or components of the immune system. To better understand the function of the DE miRNAs, enrichment analysis was performed using the KEGG ontology to identify enriched pathways that may be affected by these miRNAs (<xref ref-type="fig" rid="F2">Figure 2B</xref>, <xref ref-type="sec" rid="s12">Supplementary Table S6</xref>). Enrichment analysis identified 32 different pathways associated with these miRNAs. Of these 32 pathways, seven were directly associated with disease phenotypes and seven were associated with cellular signaling (i.e., neurotrophin signaling to Ras pathway signaling). Additional commonly enriched KEGG pathways were neurological signaling (4 pathways) and development (3 pathways). In total, these KEGG pathways included 11 of the DE miRNAs miR-126&#x2013;5p, miR-1306&#x2013;5p, miR-143&#x2013;3p, 148b-5p, miR199a-3p, miR199b-3p, miR-22&#x2013;3p, miR26b-5p, miR-296&#x2013;3p, miR545&#x2013;5p, miR-628&#x2013;5p. Most of these 11 miRNAs showed higher abundance at day 20 in piglets exposed to topsoil, except for miR-143&#x2013;3p, miR-296&#x2013;3p, and miR-22&#x2013;3p, as shown in <xref ref-type="fig" rid="F2">Figure 2C</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Gene expression and gene ontology (GO) enrichment analysis of the differentially expressed (DE) miRNA response to topsoil treatment <bold>(A)</bold> A heatmap comparison of the expression of the most highly DE miRNAs across time and treatment (control <italic>vs</italic>. topsoil) in peripheral blood mononuclear cells (PBMCs.) The three timepoints correspond to when pigs were still exposed to topsoil in the nursery (d 11) and at weaning (d 20) as well as after topsoil was removed at the postweaning (d56) timepoint. The heat map shows the clustered relative miRNA expression (least squares mean&#x2014;LSM) for control and topsoil treatments by age (day of life/age of piglet: 11, 20, and 56) <bold>(B)</bold> Visualization of the number and significance of KEGG pathway terms identified by DIANA pathway ontology enrichment analysis of DE miRNA. The size of the bubbles represents the number of target genes clustered in the KEGG ontology term group. Bubbles were colored by the number of miRNAs enriched within an ontology group <bold>(C)</bold> A histogram showing the difference in abundance of miRNA expression between control and topsoil treatment at each time points of miRNAs clustered in the KEGG ontology. Differences in miRNA abundance were estimated using the LSM between control and topsoil treatments at each specific time point (piglet age). Positive differences (above zero) mean higher abundance in the control group while negative differences (below zero) mean higher in topsoil group.</p>
</caption>
<graphic xlink:href="fgene-13-886875-g002.tif"/>
</fig>
</sec>
<sec id="s3-4">
<title>Identification of mRNA and miRNA Modules</title>
<p>To determine coordinated (i.e., correlated) changes in gene expression in response to early-life topsoil exposure, gene co-expression network modules were identified using 14,316 mRNA genes and 185 miRNAs, separately. A total of 99 control mRNA modules and 144 topsoil mRNA modules were identified (<xref ref-type="fig" rid="F3">Figures 3A,B</xref>). A total of 17 control miRNA modules and 22 topsoil miRNA modules were identified (<xref ref-type="fig" rid="F3">Figures 3C,D</xref>). We identified 22 DE genes as hub genes in control modules and 23 in topsoil modules (<xref ref-type="sec" rid="s12">Supplementary Table S3</xref>). For miRNA, four DE miRNA were hub genes in control modules while eight were hub in topsoil modules (<xref ref-type="sec" rid="s12">Supplementary Table S5</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Gene co-expression modules identified for the mRNA control group <bold>(A)</bold>, mRNA topsoil <bold>(B)</bold> miRNA control <bold>(C)</bold> and miRNA topsoil <bold>(D)</bold> groups. Co-expression modules are differentiated by color. Each group of co-expression modules is accompanied with the corresponding gene dendrogram.</p>
</caption>
<graphic xlink:href="fgene-13-886875-g003.tif"/>
</fig>
</sec>
<sec id="s3-5">
<title>Identification of miRNA:mRNA Module Interactions and Functional Analysis of Target Genes</title>
<p>To identify the potential regulation of mRNA by miRNA, miRNA:mRNA modules were identified. A total of 11 miRNA and mRNA module pairs were identified as significantly correlated (q-value &#x3c; 0.05; <xref ref-type="table" rid="T1">Table 1</xref>). In the control group, only one pair of negatively correlated modules was identified (miR-blue:mRNA-turquoise). Five negatively and five positively correlated modules were identified in the topsoil group with some modules identified as correlated across multiple modules. For example, the topsoil miR-blue module was correlated to three different mRNA modules. The black, brown and salmon miRNAs modules were each correlated with two different mRNA modules. Only one mRNA module (turquoise) was correlated to multiple miRNA modules in topsoil. <xref ref-type="table" rid="T2">Table 2</xref> summarizes the GO BP enrichment results for the four most highly correlated mRNAs modules with miRNA modules in the topsoil group. The complete list of enriched BPs is presented in <xref ref-type="sec" rid="s12">Supplementary Table S7</xref>. Only seven out of 11 miRNA:mRNA pairs enriched for GO terms. The miRblack:black interaction showed the highest number of enriched BP, totalizing 204&#xa0;BPs clustered in 36 groups, including pathways directly related to immune response, such as the groups lymphocyte homeostasis, natural killer cell-mediated cytotoxicity and regulation of T-cell receptor signaling pathways. The second highest enrichment was for the miRbrown:brown interaction with 36 BPs clustered in 10 groups including immunity-related pathways such as positive regulation of B-cell receptor signaling pathway and response to lipopolysaccharide. All the groups enriched for MiRblack:antiquewhite4 and miRblue:royalblue modules interaction were related to functionality or components of the immune system.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Identification of correlated (q &#x3c; 0.05) miRNA:mRNA modules. Modules are identified with colors names designated by WGCNA software within treatment group (topsoil or control). The number of miRNA (&#x23; miRNAs) and mRNA (&#x23; mRNA) within each module, correlation between module pairs (R), and significance of the correlation (q-value) are presented.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Treatment</th>
<th align="center">miRNA Module</th>
<th align="center">&#x23; miRNAs</th>
<th align="center">mRNA Module</th>
<th align="center">&#x23; mRNA</th>
<th align="center">q-value</th>
<th align="center">R</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Control</td>
<td align="left">miR-blue</td>
<td align="center">25</td>
<td align="left">turquoise</td>
<td align="center">1,480</td>
<td align="char" char=".">0.0018</td>
<td align="char" char=".">&#x2212;0.921</td>
</tr>
<tr>
<td rowspan="10" align="left">Topsoil</td>
<td align="left">miR-salmon</td>
<td align="center">8</td>
<td align="left">turquoise</td>
<td align="center">1,234</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">&#x2212;0.875</td>
</tr>
<tr>
<td align="left">miR-brown</td>
<td align="center">14</td>
<td align="left">brown</td>
<td align="center">739</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">&#x2212;0.839</td>
</tr>
<tr>
<td align="left">miR-blue</td>
<td align="center">15</td>
<td align="left">royalblue</td>
<td align="center">138</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">&#x2212;0.835</td>
</tr>
<tr>
<td align="left">miR-blue</td>
<td align="center">15</td>
<td align="left">sienna3</td>
<td align="center">91</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">&#x2212;0.835</td>
</tr>
<tr>
<td align="left">miR-black</td>
<td align="center">11</td>
<td align="left">black</td>
<td align="center">389</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">&#x2212;0.832</td>
</tr>
<tr>
<td align="left">miR-red</td>
<td align="center">12</td>
<td align="left">darkseagreen4</td>
<td align="center">60</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">0.825</td>
</tr>
<tr>
<td align="left">miR-salmon</td>
<td align="center">8</td>
<td align="left">thistle1</td>
<td align="center">68</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">0.825</td>
</tr>
<tr>
<td align="left">miR-blue</td>
<td align="center">15</td>
<td align="left">turquoise</td>
<td align="center">1,234</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">0.825</td>
</tr>
<tr>
<td align="left">miR-black</td>
<td align="center">11</td>
<td align="left">antiquewhite4</td>
<td align="center">58</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">0.853</td>
</tr>
<tr>
<td align="left">miR-brown</td>
<td align="center">14</td>
<td align="left">skyblue3</td>
<td align="center">87</td>
<td align="char" char=".">0.0485</td>
<td align="char" char=".">0.853</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Summary of enriched Gene Ontology (GO) biological processes (BP) for four miRNA:mRNA co-expression interaction networks in the topsoil group. miRNA targets contained in the correlated mRNA module were used for the enrichment.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="left">miRNA:mRNA module<xref ref-type="table-fn" rid="Tfn1">
<sup>a</sup>
</xref>,<xref ref-type="table-fn" rid="Tfn2">
<sup>b</sup>
</xref>
</th>
<th align="center">Target</th>
<th rowspan="2" align="center">Group GO-Term<xref ref-type="table-fn" rid="Tfn3">
<sup>c</sup>
</xref>
</th>
<th rowspan="2" align="center">Group PValue (BH)<xref ref-type="table-fn" rid="Tfn4">
<sup>d</sup>
</xref>
</th>
<th rowspan="2" align="center">N. BPs</th>
<th rowspan="2" align="center">N. Genes</th>
</tr>
<tr>
<th align="center">Genes<sup>$</sup>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="11" align="left">brown:brown</td>
<td rowspan="11" align="center">201</td>
<td align="left">Positive regulation of B-cell receptor signaling pathway</td>
<td align="center">0.0330</td>
<td align="center">3</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Response to lipopolysaccharide</td>
<td align="center">0.0079</td>
<td align="center">1</td>
<td align="center">12</td>
</tr>
<tr>
<td align="left">Response to progesterone</td>
<td align="center">0.0233</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Peptide catabolic process</td>
<td align="center">0.0198</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Negative regulation of DNA-dependent DNA replication</td>
<td align="center">0.0244</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Regulation of keratinocyte proliferation</td>
<td align="center">0.0233</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Protein glycosylation</td>
<td align="center">0.0062</td>
<td align="center">3</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Monocarboxylic acid catabolic process</td>
<td align="center">0.0180</td>
<td align="center">3</td>
<td align="center">7</td>
</tr>
<tr>
<td align="left">Negative regulation of intracellular protein transport</td>
<td align="center">0.0193</td>
<td align="center">4</td>
<td align="center">7</td>
</tr>
<tr>
<td align="left">Negative regulation of DNA biosynthetic process</td>
<td align="center">0.0055</td>
<td align="center">5</td>
<td align="center">9</td>
</tr>
<tr>
<td align="left">Response to peptide hormone</td>
<td align="center">0.0003</td>
<td align="center">13</td>
<td align="center">23</td>
</tr>
<tr>
<td align="left">black:antiquewhite4</td>
<td align="center">24</td>
<td align="left">Negative regulation of cytokine biosynthetic process</td>
<td align="center">2.1E-05</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td rowspan="3" align="left">blue:royalblue</td>
<td rowspan="3" align="center">64</td>
<td align="left">Interferon-gamma-mediated signaling pathway</td>
<td align="center">0.0006</td>
<td align="center">1</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">Defense response to virus</td>
<td align="center">7.2E-06</td>
<td align="center">3</td>
<td align="center">11</td>
</tr>
<tr>
<td align="left">Viral genome replication</td>
<td align="center">0.0004</td>
<td align="center">4</td>
<td align="center">4</td>
</tr>
<tr>
<td rowspan="37" align="left">black:black</td>
<td rowspan="37" align="center">176</td>
<td align="left">Lymphocyte homeostasis</td>
<td align="center">0.0213</td>
<td align="center">1</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">Natural killer cell-mediated cytotoxicity</td>
<td align="center">0.0241</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Regulation of T cell receptor signaling pathway</td>
<td align="center">0.0356</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Virus receptor activity</td>
<td align="center">0.0060</td>
<td align="center">7</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Virion assembly</td>
<td align="center">0.0006</td>
<td align="center">9</td>
<td align="center">14</td>
</tr>
<tr>
<td align="left">Protein serine/threonine phosphatase activity</td>
<td align="center">0.0349</td>
<td align="center">1</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">Blood coagulation, fibrin clot formation</td>
<td align="center">0.0121</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Regulation of glycoprotein metabolic process</td>
<td align="center">0.0382</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Vesicle targeting</td>
<td align="center">0.0201</td>
<td align="center">1</td>
<td align="center">5</td>
</tr>
<tr>
<td align="left">Activation of cysteine-type endopeptidase activity involved in apoptotic process</td>
<td align="center">0.0394</td>
<td align="center">1</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">Neuropeptide signaling pathway</td>
<td align="center">0.0099</td>
<td align="center">1</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">Regulation of G protein-coupled receptor signaling pathway</td>
<td align="center">0.0346</td>
<td align="center">1</td>
<td align="center">5</td>
</tr>
<tr>
<td align="left">Regulation of cell shape</td>
<td align="center">0.0112</td>
<td align="center">1</td>
<td align="center">7</td>
</tr>
<tr>
<td align="left">Antioxidant activity</td>
<td align="center">0.0357</td>
<td align="center">1</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">Embryonic digestive tract development</td>
<td align="center">0.0094</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Muscle filament sliding</td>
<td align="center">0.0127</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Platelet degranulation</td>
<td align="center">8.1E-06</td>
<td align="center">1</td>
<td align="center">12</td>
</tr>
<tr>
<td align="left">Cyclic nucleotide-dependent protein kinase activity</td>
<td align="center">0.0178</td>
<td align="center">1</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Artery development</td>
<td align="center">0.0267</td>
<td align="center">1</td>
<td align="center">4</td>
</tr>
<tr>
<td align="left">Actomyosin structure organization</td>
<td align="center">0.0129</td>
<td align="center">2</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Cardiac muscle cell apoptotic process</td>
<td align="center">0.0233</td>
<td align="center">2</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Renal water homeostasis</td>
<td align="center">0.0390</td>
<td align="center">2</td>
<td align="center">3</td>
</tr>
<tr>
<td align="left">Protein transport within lipid bilayer</td>
<td align="center">0.0020</td>
<td align="center">4</td>
<td align="center">7</td>
</tr>
<tr>
<td align="left">Myotube differentiation</td>
<td align="center">0.0202</td>
<td align="center">4</td>
<td align="center">6</td>
</tr>
<tr>
<td align="left">Calcium-mediated signaling using intracellular calcium source</td>
<td align="center">0.0197</td>
<td align="center">4</td>
<td align="center">9</td>
</tr>
<tr>
<td align="left">Negative regulation of potassium ion transport</td>
<td align="center">0.0116</td>
<td align="center">4</td>
<td align="center">9</td>
</tr>
<tr>
<td align="left">Platelet aggregation</td>
<td align="center">0.0007</td>
<td align="center">4</td>
<td align="center">13</td>
</tr>
<tr>
<td align="left">Regulation of focal adhesion assembly</td>
<td align="center">0.0201</td>
<td align="center">5</td>
<td align="center">5</td>
</tr>
<tr>
<td align="left">Regulation of muscle contraction</td>
<td align="center">0.0001</td>
<td align="center">6</td>
<td align="center">14</td>
</tr>
<tr>
<td align="left">Rab protein signal transduction</td>
<td align="center">0.0056</td>
<td align="center">6</td>
<td align="center">6</td>
</tr>
<tr>
<td align="left">Regulation of smooth muscle contraction</td>
<td align="center">0.0004</td>
<td align="center">7</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Cellular response to amino acid starvation</td>
<td align="center">0.0195</td>
<td align="center">8</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Chaperone-mediated protein complex assembly</td>
<td align="center">0.0116</td>
<td align="center">9</td>
<td align="center">9</td>
</tr>
<tr>
<td align="left">Regulation of transmembrane receptor protein serine/threonine kinase signaling</td>
<td align="center">0.0069</td>
<td align="center">12</td>
<td align="center">10</td>
</tr>
<tr>
<td align="left">Histamine secretion</td>
<td align="center">0.0006</td>
<td align="center">22</td>
<td align="center">20</td>
</tr>
<tr>
<td align="left">ADP metabolic process</td>
<td align="center">0.0056</td>
<td align="center">28</td>
<td align="center">15</td>
</tr>
<tr>
<td align="left">Positive regulation of hormone secretion</td>
<td align="center">4.4E-06</td>
<td align="center">42</td>
<td align="center">37</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="Tfn1">
<label>a</label>
<p>Individual miRNA, and mRNA, modules (i.e., highly correlated genes based on co-expression analysis) are identified with a unique color name.</p>
</fn>
<fn id="Tfn2">
<label>b</label>
<p>Only the four most highly correlated mRNAs, modules with miRNA, modules in topsoil group are showed in this summary table<sup>
<bold>$</bold>
</sup> For each mRNA, module, genes that were the target of the associated miRNA, modules were used for the enrichment analysis in ClueGO.</p>
</fn>
<fn id="Tfn3">
<label>c</label>
<p>Biological process ontology terms directly related to immune function are highlighted in red.</p>
</fn>
<fn id="Tfn4">
<label>d</label>
<p>The group p-value was calculated using the list of unique genes found within the uploaded gene list. <italic>p</italic>-values were adjusted using a Benjamini and Hochberg correction (BH) (p-value &#x3c; 0.05).</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>Discussion</title>
<p>The objective of this study was to identify potential mechanisms controlling how the immune system is programmed in response to microbes and the environment. To understand how early-life exposure to less-hygienic conditions alters immune training, piglets were exposed to topsoil only from d 4-d 21 of life, which resulted in the modulation of gut microbiota and growth performance (<xref ref-type="bibr" rid="B44">Vo et al., 2017</xref>). To determine changes in immune signaling, PBMC mRNA and miRNA gene expression were analyzed over the first 56&#xa0;days of life. The following passages provide potential explanations for the observed changes in PBMC gene expression in response to early-life topsoil exposure along with the potential limitations in the interpretation of this study.</p>
<p>Genes identified as DE indicate changes in the development of the immune system which under the parameters of this study may be in response to interactions with novel antigens or nutrient availability. This analysis identified DE genes associated with T-cell activation. The first two DE genes, <italic>PTPRJ</italic> and <italic>ITGB3,</italic> are known to be involved in the formation of the immune synapse or area of communication between a T cell and an antigen-presenting cell (APC). Both genes were upregulated at d 20 in the topsoil group (<xref ref-type="sec" rid="s12">Supplementary Figure S1</xref>, <xref ref-type="fig" rid="F1">Figure 1C</xref>) which may be in response to increased contact with microbes during the nursery phase from d 4 to d 21. Corroborating these results, <xref ref-type="bibr" rid="B28">Lin et al. (2004)</xref> found increased levels of <italic>PTPRJ</italic> in T cells upon activation in mice, indicating that this increase in the topsoil-exposed piglets may have elevated T-cell activation at d 20. Conversely, the lower levels of <italic>PTPRJ</italic> in na&#xef;ve T cells in mice (<xref ref-type="bibr" rid="B28">Lin et al., 2004</xref>) indicate that the lower expression of <italic>PTPRJ</italic> in the control piglets at d 20, when the piglets were in higher hygienic conditions, results from lower antigen contact and probably higher na&#xef;ve T cells.</p>
<p>Two of the identified DE genes makeup components of the T-cells receptor (TCR) complex, <italic>TRBV30</italic> and <italic>CD3D,</italic> are crucial to antigen recognition and initiation of cellular signaling. The diversity of paired alpha and beta chains of the TCR repertoire confers broad immune coverage against pathogens. Recent studies using next-generation sequencing have characterized the TCR repertoire of individuals. The TCR repertoire is composed of a diverse combination of alpha and beta TCR chains (<xref ref-type="bibr" rid="B3">Attaf et al., 2015</xref>), which may provide the pathogen exposure history and immunological memory of a given individual (<xref ref-type="bibr" rid="B10">DeWitt et al., 2018</xref>). <italic>TRBV30</italic> is one of the several genes encoding the V region of the variable domain of TCR beta chain which is responsible for recognizing specific peptides presented by APCs (<xref ref-type="bibr" rid="B9">Davis et al., 1995</xref>). The <italic>CD3D</italic> gene is associated with TCR chain activation after antigen recognition (<xref ref-type="bibr" rid="B8">Dave et al., 1997</xref>; <xref ref-type="bibr" rid="B16">Garcill&#xe1;n et al., 2021</xref>). Both <italic>TRBV30</italic> and <italic>CD3D</italic> exhibited similar trends of expression over time in this study. In piglets exposed to topsoil, expression levels increased from d 11 to d 20 and then decreased at d 56 (<xref ref-type="fig" rid="F1">Figure 1C</xref>). In control piglets, expression levels decreased from d 11 to d 20 and increased at d 56. The increased expression at d 20 followed by a decrease at d 56 in the topsoil group as well as the delayed peak in the control piglets may indicate earlier development of the immune system in topsoil-exposed piglets compared with the control. Early-life exposure to topsoil may have provided microbes or other substances that resulted in these differences in gene expression that allow for faster response to future exposure of these microbes. Once ingested by the piglets, microbes and other constituents of the topsoil would interact with Paneth cells and APCs that stimulate the immune system.</p>
<p>Other DE genes with known function in immune response (i.e., <italic>BMP8A, GPD2</italic>, <italic>PHKA1</italic> and <italic>PHKA2</italic>) (<xref ref-type="sec" rid="s12">Supplementary Figure S1</xref>) exhibited a similar pattern of expression, which may support the hypothesis that exposure to topsoil alters the timing or development of the immune system. Of particular interest, <italic>BMP8A</italic> is a transforming growth factor beta (TGF-&#x3b2;) superfamily member contained within a signaling pathway essential to activation and homeostasis of na&#xef;ve CD4<sup>&#x2b;</sup> T cells (<xref ref-type="bibr" rid="B32">Mart&#xed;nez et al., 2015</xref>).</p>
<p>The ontology enrichment analyses identified pathways involved in T-cell activation, specifically those involved in the interaction between the T cell and the APC, including the formation of the interactive space between the two also known as the immunological synapse. The DE genes <italic>STX3, PTPRJ</italic>, <italic>ITGB3, ITGA3</italic>, <italic>MYH9</italic>, and <italic>TLN1</italic> were clustered in several GO terms related to processes involved in many aspects of the immune synapse and downstream signaling. These genes function in vesicle-mediated transport to the plasma membrane (GO:0098876), protein transport within plasma membrane (GO:0099632), microtubule-based protein transport (GO:0099118), integrin-mediated signaling pathways (GO:0007229), positive regulation of peptidyl-tyrosine phosphorylation (GO:0050731), and positive regulation of cell-substrate adhesion (GO:0010811). After T-cell receptor stimulation, cytoskeleton reorganization (actin and microtubules) is needed to adapt the cell conformation and formation of the synapse architecture (<xref ref-type="bibr" rid="B42">Valitutti et al., 1995</xref>; <xref ref-type="bibr" rid="B31">Mart&#xed;n-C&#xf3;freces and S&#xe1;nchez-Madrid, 2018</xref>). Traffic through the membrane is important to regulate the activation of T cells by balancing the cellular localization of a wide type of components as receptors, signaling molecules or cytokines (reviewed by <xref ref-type="bibr" rid="B37">Soares et al., 2013</xref>). <italic>STX3</italic> has a putative role in the formation of the immune synapse of cytotoxic T cells as it is involved in protein trafficking involving vesicle fusion and exocytosis. STX3 co-localized with CD3 proteins which are components of TCR-CD3 complex on T-cell surfaces (<xref ref-type="bibr" rid="B34">Pattu et al., 2012</xref>). Two additional DE genes, <italic>PTPRJ</italic> and <italic>ITGB3,</italic> are known to be involved in the formation of the immune synapse and clustered in GO terms related to cell-substrate adhesion and regulation of tyrosine phosphorylation. The <italic>ITGB3</italic> gene plays a role in the immune synapse by encoding a subunit of specific integrins. Other DE genes involved in these processes include <italic>ITGA3</italic>, <italic>MYH9,</italic> and <italic>TLN1</italic> (<xref ref-type="fig" rid="F1">Figure 1C</xref>, <xref ref-type="sec" rid="s12">Supplementary Figure S1</xref>), all of which increase at d 20 in piglets exposed to topsoil suggesting that these genes support the formation of the immune synapse (<xref ref-type="fig" rid="F1">Figure 1C</xref>).</p>
<p>An important finding in relation to the development of the immune system was the identification of mir-143 as DE in response to topsoil exposure over time (<xref ref-type="fig" rid="F2">Figure 2C</xref>). MicroRNAs are important regulators of gene expression with important roles in mediating the function of the immune system (<xref ref-type="bibr" rid="B29">Ma et al., 2011</xref>; <xref ref-type="bibr" rid="B48">Zhang et al., 2018</xref>). Typically, each miRNA suppresses gene expression in tens to thousands of mRNA transcripts to coordinate the use of pathways within the transcriptome (<xref ref-type="bibr" rid="B17">Gebert and MacRae, 2019</xref>). Previous studies indicate that overexpression of mir-143 may enhance the conversion of cytotoxic T cells into memory cells by reprogramming cellular metabolism (<xref ref-type="bibr" rid="B48">Zhang et al., 2018</xref>). In addition, mir-143 was identified as a candidate regulator of a miRNA-target co-expression network in response to topsoil and a hub miRNA within the brown miRNA module. It was also negatively correlated with a group of mRNAs (brown mRNA module) only in the piglets exposed to topsoil during early life. Target mRNA genes within the miRNA:mRNA correlated modules containing mir-143 were enriched for biological processes and pathways such as response to lipopolysaccharide and positive regulation of B-cell receptor signaling pathway. Together, mir-143, the other hub miRNAs, and their co-expressed target genes in the brown module make up a miRNA-target network acting as a putative driver of the immune response in the animals exposed to topsoil in early life.</p>
<p>Analysis of mRNA-miRNA interactions also identified mir-29a-3p and mir-148&#x2013;5p as potential regulators of gene expression only in topsoil-exposed pigs. Upon bacterial infection, mir-29a-3p has been identified to be downregulated in innate immune cells, such as natural killer and adaptive immune cells such as helper T cells and cytotoxic T cells (<xref ref-type="bibr" rid="B29">Ma et al., 2011</xref>). Targets of mir-29a-3p were enriched for regulation of mTOR signaling which has a role in the innate and adaptive immune system by regulating metabolism. As reviewed by <xref ref-type="bibr" rid="B7">Chi. (2012)</xref>, mTOR signaling can act in different pathways affecting T-cell activation, differentiation, and immune homeostasis. The third interesting candidate miRNA-target network identified was driven by mir-148a-5p, which was negatively correlated to the royalblue module. The mRNA genes within the royalblue mRNA module were enriched for genes involved in response to type I interferon and interferon-gamma-mediated signaling pathway. Type I interferon is secreted by the innate immune cells after detection of the pathogen, but also can act in the adaptive immune response by stimulating antibody production of B cells or intensifying the T-cells function (Reviewed by <xref ref-type="bibr" rid="B23">Ivashkiv and Donlin, 2014</xref>).</p>
</sec>
<sec sec-type="conclusion" id="s5">
<title>Conclusion</title>
<p>Understanding the role hygienic environments have on host gene expression is important in our understanding of the development of the immune system. The findings of this study indicate that early-life exposure of piglets to topsoil impacts both mRNA and miRNA expression in PBMCs, particularly genes and miRNA involved in the development of T cells and their ability to interact with antigen-presenting cells. Key regulatory genes identified include: <italic>PTPRJ, ITGB3, TRBV30, CD3D</italic>, mir-143, mir-29, and mir-148a. While these results need to be interpreted with caution as they represent only a snapshot in time across a small window of time from a subset of peripheral immune cells, these results support the current hygiene hypothesis in that less-hygienic environments during early life contribute to the development of the immune system. However, additional studies are needed to characterize how topsoil exposure may impact gene expression within the gut and impact gut phenotypes critical to animal health.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data Availability Statement</title>
<p>The RNA-Seq data set generated and analyzed during the current study is available in the NCBI&#x2019;s Gene Expression Omnibus (GEO) under accession numbers GSE205321.</p>
</sec>
<sec id="s7">
<title>Ethics Statement</title>
<p>The animal study was reviewed and approved by the University of Arkansas Institute of Animal Care and Use Committee under protocol number 13060.</p>
</sec>
<sec id="s8">
<title>Author Contributions</title>
<p>CM, TT, and MS developed and conducted the experiment. MD performed bioinformatics analysis, interpreted results, and drafted the initial manuscript with JK and DK. HB conducted sample quality check, read trimming and mapping to genome, differential expression analysis, and gene co-expression network analysis. DK and MD conducted statistical analyses. JK oversaw the bioinformatics and statistical analysis of the data. MS and TT conducted animal and wet-lab experiments and coordinated RNA-seq analyses. JZ and CM provided funding for the experiment and analyses. All authors read, reviewed, and contributed to the writing of the manuscript.</p>
</sec>
<sec id="s9">
<title>Funding</title>
<p>This research was supported by a grant from the Arkansas Biosciences Institute (ABI) and the USDA-NIFA project &#x23;2018-67015-27479.</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ack>
<p>The authors thank Tim Johnson for his assistance in sample collection and Dr. Joan Lunney for helpful comments and feedback regarding early drafts of the manuscript.</p>
</ack>
<sec id="s12">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fgene.2022.886875/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fgene.2022.886875/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Presentation1.zip" id="SM1" mimetype="application/zip" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<sec id="s13">
<title>Abbreviations</title>
<p>APC, antigen-presenting cell; BP, biological process; DE, differentially expressed; FDR, false discovery rate; GO, gene ontology; MM, module membership; PBMC, peripheral blood mononuclear cells; RNA-Seq, RNA sequencing; RQI, RNA quality score; TCR, T-cells receptor; TMM, Trimmed Mean of M-values; WGCNA, the Weighted Gene Co-expression Network Analysis.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Andrews</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Bioinformatics</surname>
<given-names>Babraham</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>FastQC: A Quality Control Tool for High Throughput Sequence Data</article-title>. <comment>
<italic>Manual</italic>
</comment> http://www.bioinformatics.babraham.ac.uk/projects/. </citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Atarashi</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Tanoue</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Shima</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Imaoka</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kuwahara</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Momose</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Induction of Colonic Regulatory T Cells by Indigenous Clostridium Species</article-title>. <source>Science</source> <volume>331</volume>, <fpage>337</fpage>&#x2013;<lpage>341</lpage>. <pub-id pub-id-type="doi">10.1126/science.1198469</pub-id> </citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Attaf</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Huseby</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Sewell</surname>
<given-names>A. K.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>&#x391;&#x3b2; T Cell Receptors as Predictors of Health and Disease</article-title>. <source>Cell. Mol. Immunol.</source> <volume>12</volume>, <fpage>391</fpage>&#x2013;<lpage>399</lpage>. <pub-id pub-id-type="doi">10.1038/cmi.2014.134</pub-id> </citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bindea</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Mlecnik</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Hackl</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Charoentong</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Tosolini</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kirilovsky</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2009</year>). <article-title>ClueGO: A Cytoscape Plug-In to Decipher Functionally Grouped Gene Ontology and Pathway Annotation Networks</article-title>. <source>Bioinformatics</source> <volume>25</volume>, <fpage>1091</fpage>&#x2013;<lpage>1093</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btp101</pub-id> </citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bullard</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Purdom</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Hansen</surname>
<given-names>K. D.</given-names>
</name>
<name>
<surname>Dudoit</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Evaluation of Statistical Methods for Normalization and Differential Expression in mRNA-Seq Experiments</article-title>. <source>BMC Bioinforma.</source> <volume>11</volume>, <fpage>94</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2105-11-94</pub-id> </citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cheng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Thymic T-Cell Production Is Associated with Changes in the Gut Microbiota in Young Chicks</article-title>. <source>Front. Immunol.</source> <volume>12</volume>, <fpage>700603</fpage>. <pub-id pub-id-type="doi">10.3389/fimmu.2021.700603</pub-id> </citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chi</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Regulation and Function of mTOR Signalling in T Cell Fate Decisions</article-title>. <source>Nat. Rev. Immunol.</source> <volume>12</volume>, <fpage>325</fpage>&#x2013;<lpage>338</lpage>. <pub-id pub-id-type="doi">10.1038/nri3198</pub-id> </citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dave</surname>
<given-names>V. P.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Browne</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Alarcon</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Fernandez-Miguel</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Lafaille</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>1997</year>). <article-title>CD3delta Deficiency Arrests Development of the Alpha Beta but Not the Gamma Delta T Cell Lineage</article-title>. <source>EMBO J.</source> <volume>16</volume>, <fpage>1360</fpage>&#x2013;<lpage>1370</lpage>. <pub-id pub-id-type="doi">10.1093/emboj/16.6.1360</pub-id> </citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Davis</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>McHEYZER-Williams,a</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Chien</surname>
<given-names>Y.-H.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>T-cell Receptor V-Region Usage and Antigen Specificity</article-title>. <source>Ann. N. Y. Acad. Sci.</source> <volume>756</volume>, <fpage>1</fpage>&#x2013;<lpage>11</lpage>. <pub-id pub-id-type="doi">10.1111/j.1749-6632.1995.tb44477.x</pub-id> </citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>DeWitt</surname>
<given-names>W. S.</given-names>
</name>
<name>
<surname>Smith</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Schoch</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Hansen</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Matsen</surname>
<given-names>F. A.</given-names>
</name>
<name>
<surname>Bradley</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Human T Cell Receptor Occurrence Patterns Encode Immune History, Genetic Background, and Receptor Specificity</article-title>. <source>Elife</source> <volume>7</volume>, <fpage>e38358</fpage>. <pub-id pub-id-type="doi">10.7554/eLife.38358</pub-id> </citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dobin</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Davis</surname>
<given-names>C. A.</given-names>
</name>
<name>
<surname>Schlesinger</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Drenkow</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zaleski</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Jha</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>STAR: Ultrafast Universal RNA-Seq Aligner</article-title>. <source>Bioinformatics</source> <volume>29</volume>, <fpage>15</fpage>&#x2013;<lpage>21</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/bts635</pub-id> </citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dourmad</surname>
<given-names>J. Y.</given-names>
</name>
<name>
<surname>Hassouna</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Robin</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Guingand</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Meunier-Sala&#xfc;n</surname>
<given-names>M. C.</given-names>
</name>
<name>
<surname>Lebret</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Influence of Pig Rearing System on Animal Performance and Manure Composition</article-title>. <source>Animal</source> <volume>3</volume>, <fpage>606</fpage>&#x2013;<lpage>616</lpage>. <pub-id pub-id-type="doi">10.1017/S1751731108003601</pub-id> </citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Durinck</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Spellman</surname>
<given-names>P. T.</given-names>
</name>
<name>
<surname>Birney</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Huber</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Mapping Identifiers for the Integration of Genomic Datasets with the R/Bioconductor Package biomaRt</article-title>. <source>Nat. Protoc.</source> <volume>4</volume>, <fpage>1184</fpage>&#x2013;<lpage>1191</lpage>. <pub-id pub-id-type="doi">10.1038/nprot.2009.97</pub-id> </citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ennamorati</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Vasudevan</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Clerkin</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Halvorsen</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Verma</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ibrahim</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Intestinal Microbes Influence Development of Thymic Lymphocytes in Early Life</article-title>. <source>Proc. Natl. Acad. Sci. U.S.A.</source> <volume>117</volume>, <fpage>2570</fpage>&#x2013;<lpage>2578</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1915047117</pub-id> </citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Friedl&#xe4;nder</surname>
<given-names>M. R.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Adamidi</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Maaskola</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Einspanier</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Knespel</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2008</year>). <article-title>Discovering microRNAs from Deep Sequencing Data Using miRDeep</article-title>. <source>Nat. Biotechnol.</source> <volume>26</volume>, <fpage>407</fpage>&#x2013;<lpage>415</lpage>. <pub-id pub-id-type="doi">10.1038/nbt1394</pub-id> </citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garcill&#xe1;n</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Fuentes</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Marin</surname>
<given-names>A. V.</given-names>
</name>
<name>
<surname>Megino</surname>
<given-names>R. F.</given-names>
</name>
<name>
<surname>Chacon-Arguedas</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Mazariegos</surname>
<given-names>M. S.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>CD3G or CD3D Knockdown in Mature, but Not Immature, T Lymphocytes Similarly Cripples the Human TCR&#x3b1;&#x3b2; Complex</article-title>. <source>Front. Cell. Dev. Biol.</source> <volume>9</volume>, <fpage>608490</fpage>. <pub-id pub-id-type="doi">10.3389/fcell.2021.608490</pub-id> </citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gebert</surname>
<given-names>L. F. R.</given-names>
</name>
<name>
<surname>MacRae</surname>
<given-names>I. J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Regulation of microRNA Function in Animals</article-title>. <source>Nat. Rev. Mol. Cell. Biol.</source> <volume>20</volume>, <fpage>21</fpage>&#x2013;<lpage>37</lpage>. <pub-id pub-id-type="doi">10.1038/s41580-018-0045-7</pub-id> </citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gury-BenAri</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Thaiss</surname>
<given-names>C. A.</given-names>
</name>
<name>
<surname>Serafini</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Winter</surname>
<given-names>D. R.</given-names>
</name>
<name>
<surname>Giladi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Lara-Astiaso</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>The Spectrum and Regulatory Landscape of Intestinal Innate Lymphoid Cells Are Shaped by the Microbiome</article-title>. <source>Cell.</source> <volume>166</volume>, <fpage>1231</fpage>&#x2013;<lpage>1246</lpage>. <comment>e13</comment>. <pub-id pub-id-type="doi">10.1016/j.cell.2016.07.043</pub-id> </citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>C.-F.</given-names>
</name>
<name>
<surname>Chie</surname>
<given-names>W.-C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>I.-J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Effect of Environmental Exposures on Allergen Sensitization and the Development of Childhood Allergic Diseases: A Large-Scale Population-Based Study</article-title>. <source>World Allergy Organ. J.</source> <volume>14</volume>, <fpage>100495</fpage>. <pub-id pub-id-type="doi">10.1016/j.waojou.2020.100495</pub-id> </citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Iljazovic</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Roy</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>G&#xe1;lvez</surname>
<given-names>E. J. C.</given-names>
</name>
<name>
<surname>Lesker</surname>
<given-names>T. R.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Gronow</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Perturbation of the Gut Microbiome by Prevotella Spp. Enhances Host Susceptibility to Mucosal Inflammation</article-title>. <source>Mucosal Immunol.</source> <volume>14</volume>, <fpage>113</fpage>&#x2013;<lpage>124</lpage>. <pub-id pub-id-type="doi">10.1038/s41385-020-0296-4</pub-id> </citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Illi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>von Mutius</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Lau</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Bergmann</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Niggemann</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Sommerfeld</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2001</year>). <article-title>Early Childhood Infectious Diseases and the Development of Asthma up to School Age: A Birth Cohort Study</article-title>. <source>Br. Med. J.</source> <volume>322</volume>, <fpage>390</fpage>&#x2013;<lpage>395</lpage>. <pub-id pub-id-type="doi">10.1136/bmj.322.7283.390</pub-id> </citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Institute Inc</surname>
<given-names>S. A. S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Base SAS&#xae;9.4 Procedures Guide</article-title>. <source>Stat. Proced.</source> <pub-id pub-id-type="doi">10.1016/B978-0-444-59425-9.00014-7</pub-id> </citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ivashkiv</surname>
<given-names>L. B.</given-names>
</name>
<name>
<surname>Donlin</surname>
<given-names>L. T.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Regulation of Type I Interferon Responses</article-title>. <source>Nat. Rev. Immunol.</source> <volume>14</volume>, <fpage>36</fpage>&#x2013;<lpage>49</lpage>. <pub-id pub-id-type="doi">10.1038/nri3581</pub-id> </citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Krueger</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Trim Galore!: A Wrapper Tool Around Cutadapt and FastQC to Consistently Apply Quality and Adapter Trimming to FastQ Files</article-title>. <source>Babraham Inst.</source> <pub-id pub-id-type="doi">10.1002/maco.200603986</pub-id> </citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langfelder</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Horvath</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>WGCNA: An R Package for Weighted Correlation Network Analysis</article-title>. <source>BMC Bioinforma.</source> <volume>9</volume>, <fpage>559</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2105-9-559</pub-id> </citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lebret</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Prunier</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Bonhomme</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Foury</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Morm&#xe8;de</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Dourmad</surname>
<given-names>J. Y.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Physiological Traits and Meat Quality of Pigs as Affected by Genotype and Housing System</article-title>. <source>Meat Sci.</source> <volume>88</volume>, <fpage>14</fpage>&#x2013;<lpage>22</lpage>. <pub-id pub-id-type="doi">10.1016/j.meatsci.2010.11.025</pub-id> </citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Smyth</surname>
<given-names>G. K.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>FeatureCounts: An Efficient General Purpose Program for Assigning Sequence Reads to Genomic Features</article-title>. <source>Bioinformatics</source> <volume>30</volume>, <fpage>923</fpage>&#x2013;<lpage>930</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btt656</pub-id> </citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>J. W.</given-names>
</name>
<name>
<surname>Baker</surname>
<given-names>J. E.</given-names>
</name>
<name>
<surname>Weiss</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Regulated Expression of the Receptor-like Tyrosine Phosphatase CD148 on Hemopoietic Cells</article-title>. <source>J. Immunol.</source> <volume>173</volume>, <fpage>2324</fpage>&#x2013;<lpage>2330</lpage>. <pub-id pub-id-type="doi">10.4049/jimmunol.173.4.2324</pub-id> </citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>The microRNA miR-29 Controls Innate and Adaptive Immune Responses to Intracellular Bacterial Infection by Targeting Interferon-&#x3b3;</article-title>. <source>Nat. Immunol.</source> <volume>12</volume>, <fpage>861</fpage>&#x2013;<lpage>869</lpage>. <pub-id pub-id-type="doi">10.1038/ni.2073</pub-id> </citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>MacNeill</surname>
<given-names>S. J.</given-names>
</name>
<name>
<surname>Sozanska</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Danielewicz</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Debinska</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kosmeda</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Boznanski</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>Asthma and Allergies: Is the Farming Environment (Still) Protective in Poland? the GABRIEL Advanced Studies</article-title>. <source>Allergy</source> <volume>68</volume>, <fpage>771</fpage>&#x2013;<lpage>779</lpage>. <pub-id pub-id-type="doi">10.1111/all.12141</pub-id> </citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mart&#xed;n-C&#xf3;freces</surname>
<given-names>N. B.</given-names>
</name>
<name>
<surname>S&#xe1;nchez-Madrid</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Sailing to and Docking at the Immune Synapse: Role of Tubulin Dynamics and Molecular Motors</article-title>. <source>Front. Immunol.</source> <volume>9</volume>, <fpage>1174</fpage>. <pub-id pub-id-type="doi">10.3389/fimmu.2018.01174</pub-id> </citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mart&#xed;nez</surname>
<given-names>V. G.</given-names>
</name>
<name>
<surname>Saced&#xf3;n</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Hidalgo</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Valencia</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Fern&#xe1;ndez-Sevilla</surname>
<given-names>L. M.</given-names>
</name>
<name>
<surname>Hern&#xe1;ndez-L&#xf3;pez</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>The BMP Pathway Participates in Human Naive CD4&#x2b; T Cell Activation and Homeostasis</article-title>. <source>PLoS One</source> <volume>10</volume>, <fpage>e0131453</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0131453</pub-id> </citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Masotti</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Interplays between Gut Microbiota and Gene Expression Regulation by miRNAs</article-title>. <source>Front. Cell. Inf. Microbio.</source> <volume>2</volume>, <fpage>137</fpage>. <pub-id pub-id-type="doi">10.3389/fcimb.2012.00137</pub-id> </citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pattu</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Qu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Schwarz</surname>
<given-names>E. C.</given-names>
</name>
<name>
<surname>Strau&#xdf;</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Weins</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Bhat</surname>
<given-names>S. S.</given-names>
</name>
<etal/>
</person-group> (<year>2012</year>). <article-title>SNARE Protein Expression and Localization in Human Cytotoxic T Lymphocytes</article-title>. <source>Eur. J. Immunol.</source> <volume>42</volume>, <fpage>470</fpage>&#x2013;<lpage>475</lpage>. <pub-id pub-id-type="doi">10.1002/eji.201141915</pub-id> </citation>
</ref>
<ref id="B35">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Revelle</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2018</year>). <source>Psych: Procedures for Personality and Psychological Research</source>. <publisher-loc>Evanston, Illinois, USA</publisher-loc>: <publisher-name>Northwestern University</publisher-name>. <comment>Available at:</comment> https://CRAN.R-project.org/package&#x3d;psych. Version &#x3d; 1.8.10. </citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robinson</surname>
<given-names>M. D.</given-names>
</name>
<name>
<surname>McCarthy</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Smyth</surname>
<given-names>G. K.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>edgeR: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data</article-title>. <source>Bioinformatics</source> <volume>26</volume>, <fpage>139</fpage>&#x2013;<lpage>140</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btp616</pub-id> </citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Soares</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lasserre</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Alcover</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Orchestrating Cytoskeleton and Intracellular Vesicle Traffic to Build Functional Immunological Synapses</article-title>. <source>Immunol. Rev.</source> <volume>256</volume>, <fpage>118</fpage>&#x2013;<lpage>132</lpage>. <pub-id pub-id-type="doi">10.1111/imr.12110</pub-id> </citation>
</ref>
<ref id="B38">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Storey</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bass</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Dabney</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Robinson</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2019</year>). <source>Qvalue: Q-Value Estimation for False Discovery Rate Control</source>. <comment>R package version 2.18.0. Available at:</comment> http://github.com/jdstorey/qvalue. </citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Strachan</surname>
<given-names>D. P.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Hay Fever, Hygiene, and Household Size</article-title>. <source>Bmj</source> <volume>299</volume>, <fpage>1259</fpage>&#x2013;<lpage>1260</lpage>. <pub-id pub-id-type="doi">10.1136/bmj.299.6710.1259</pub-id> </citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Su</surname>
<given-names>L. F.</given-names>
</name>
<name>
<surname>Kidd</surname>
<given-names>B. A.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kotzin</surname>
<given-names>J. J.</given-names>
</name>
<name>
<surname>Davis</surname>
<given-names>M. M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Virus-Specific CD4&#x2b; Memory-Phenotype T Cells Are Abundant in Unexposed Adults</article-title>. <source>Immunity</source> <volume>38</volume>, <fpage>373</fpage>&#x2013;<lpage>383</lpage>. <pub-id pub-id-type="doi">10.1016/j.immuni.2012.10.021</pub-id> </citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tsai</surname>
<given-names>T. C.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Sales</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Erf</surname>
<given-names>G. F.</given-names>
</name>
<name>
<surname>Kegley</surname>
<given-names>E. B.</given-names>
</name>
<etal/>
</person-group> (<year>2016</year>). <article-title>259 Effect of Rearing Strategies during Lactation on Growth Performance and the Population of Innate and Adaptive Immune Cells in Pigs from Pre-weaning to Market</article-title>. <source>J. Anim. Sci.</source> <volume>94</volume>, <fpage>122</fpage>&#x2013;<lpage>123</lpage>. <pub-id pub-id-type="doi">10.2527/msasas2016-259</pub-id> </citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Valitutti</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Dessing</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Aktories</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Gallati</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lanzavecchia</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Sustained Signaling Leading to T Cell Activation Results from Prolonged T Cell Receptor Occupancy. Role of T Cell Actin Cytoskeleton</article-title>. <source>J. Exp. Med.</source> <volume>181</volume>, <fpage>577</fpage>&#x2013;<lpage>584</lpage>. <pub-id pub-id-type="doi">10.1084/jem.181.2.577</pub-id> </citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vlachos</surname>
<given-names>I. S.</given-names>
</name>
<name>
<surname>Zagganas</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Paraskevopoulou</surname>
<given-names>M. D.</given-names>
</name>
<name>
<surname>Georgakilas</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Karagkouni</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Vergoulis</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>DIANA-miRPath v3.0: Deciphering microRNA Function with Experimental Support</article-title>. <source>Nucleic Acids Res.</source> <volume>43</volume>, <fpage>W460</fpage>&#x2013;<lpage>W466</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkv403</pub-id> </citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vo</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Tsai</surname>
<given-names>T. C.</given-names>
</name>
<name>
<surname>Maxwell</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Carbonero</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Early Exposure to Agricultural Soil Accelerates the Maturation of the Early-Life Pig Gut Microbiota</article-title>. <source>Anaerobe</source> <volume>45</volume>, <fpage>31</fpage>&#x2013;<lpage>39</lpage>. <pub-id pub-id-type="doi">10.1016/j.anaerobe.2017.02.022</pub-id> </citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wen</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>van Dixhoorn</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Schokker</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Woelders</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Stockhofe-Zurwieden</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Rebel</surname>
<given-names>J. M. J.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Environmentally Enriched Housing Conditions Affect Pig Welfare, Immune System and Gut Microbiota in Early Life</article-title>. <source>Anim. Microbiome</source> <volume>3</volume>, <fpage>52</fpage>. <pub-id pub-id-type="doi">10.1186/s42523-021-00115-2</pub-id> </citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wenner</surname>
<given-names>B. A.</given-names>
</name>
<name>
<surname>Zerby</surname>
<given-names>H. N.</given-names>
</name>
<name>
<surname>Boler</surname>
<given-names>D. D.</given-names>
</name>
<name>
<surname>Gebreyes</surname>
<given-names>W. A.</given-names>
</name>
<name>
<surname>Moeller</surname>
<given-names>S. J.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Effect of Mannan Oligosaccharides (Bio-Mos) and Outdoor Access Housing on Pig Growth, Feed Efficiency and Carcass Composition</article-title>. <source>J. Animal Sci.</source> <volume>91</volume>, <fpage>4936</fpage>&#x2013;<lpage>4944</lpage>. <pub-id pub-id-type="doi">10.2527/jas.2013-6582</pub-id> </citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Burns</surname>
<given-names>M. B.</given-names>
</name>
<name>
<surname>Subramanian</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Blekhman</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Interaction between Host MicroRNAs and the Gut Microbiota in Colorectal Cancer</article-title>. <source>mSystems</source> <volume>3</volume>. <pub-id pub-id-type="doi">10.1128/msystems.00205-17</pub-id> </citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ping</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>miR-143 Regulates Memory T Cell Differentiation by Reprogramming T Cell Metabolism</article-title>. <source>J. I.</source> <volume>201</volume>, <fpage>2165</fpage>&#x2013;<lpage>2175</lpage>. <pub-id pub-id-type="doi">10.4049/jimmunol.1800230</pub-id> </citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Liwinski</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Elinav</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Interaction between Microbiota and Immunity in Health and Disease</article-title>. <source>Cell. Res.</source> <volume>30</volume>, <fpage>492</fpage>&#x2013;<lpage>506</lpage>. <pub-id pub-id-type="doi">10.1038/s41422-020-0332-7</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>