A Novel lncRNA Regulates the Toll-Like Receptor Signaling Pathway and Related Immune Function by Stabilizing FOS mRNA as a Competitive Endogenous RNA

Long non-coding RNAs (lncRNAs) have recently emerged as new regulatory molecules with diverse functions in regulating gene expression and significant roles in the immune response. However, the function of many unknown lncRNAs is still unclear. By studying the regulatory effect of daidzein (DA) on immunity, we identified a novel lncRNA with an immune regulatory function: lncRNA- XLOC_098131. In vivo, DA treatment upregulated the expression of lncRNA- XLOC_098131, FOS, and JUN in chickens and affected the expression of activator protein 1 (AP-1) to regulate MAPK signaling, Toll-like receptor signaling, and related mRNA expression. It also enhanced macrophage activity and increased the numbers of blood neutrophils and mononuclear cells, which can improve the body's ability to respond to stress and bacterial and viral infections. Furthermore, DA treatment also reduced B lymphocyte apoptosis and promoted the differentiation of B lymphocytes into plasma cells, which in turn resulted in the production of more immunoglobulins and the promotion of antigen presentation. In vitro, using HEK293FT cells, we demonstrated that mir-548s could bind to and decrease the expression of both FOS and lncRNA- XLOC_098131. LncRNA- XLOC_098131 served as a competitive endogenous RNA to stabilize FOS by competitively binding to miR-548s and thereby reducing its inhibitory effect of FOS expression. Therefore, we concluded that the novel lncRNA XLOC_098131 acts as a key regulatory molecule that can regulate the Toll-like receptor signaling pathway and related immune function by serving as a competitive endogenous RNA to stabilize FOS mRNA expression.


INTRODUCTION
As a template RNA molecule, mRNA plays an important role in DNA genetic information transmission and protein synthesis (1). However, with advancements in research technologies, researchers have identified that mRNA accounts for only ∼1.2% of transcripts, whereas non-coding RNA (ncRNA) in the transcriptome accounts for as much as 98%, which suggests that ncRNA might perform rich biological functions in organisms (2). Long non-coding RNAs (lncRNAs), which are non-coding RNA transcripts that are longer than 200 nucleotides and lack apparent open reading frames, play important roles in many cellular biological processes, including cell cycle progression, apoptosis, development, muscle differentiation, and immune regulation (3), and participate mainly in epigenetic, transcriptional, and posttranscriptional regulation (4). LncRNAs can regulate epigenetic processing by affecting chromatin reprogramming, DNA methylation and histone modification (5). In addition, lncRNAs mediate gene silencing by recruiting chromatin-modifying complexes toward target genes through sequence-specific binding (6). LncRNAs, such as Kcnq1ot1 and Air, which map to the Kcnq1 and Igf2r imprinted gene clusters, respectively, mediate the transcriptional silencing of multiple genes by interacting with chromatin and recruiting the chromatin-modifying machinery (7). Furthermore, lncRNAs can regulate gene expression at the transcriptional and posttranscriptional levels. The depletion of several ncRNAs leads to decreased expression of their neighboring protein-coding genes, including the master regulator of hematopoiesis, SCL (also called TAL1), Snai1 and Snai2 (8). After transcription, lncRNAs can regulate gene expression as competitive endogenous RNAs (ceRNAs) (5,9,10).
LncRNAs can participate in immune cell differentiation and the regulation of related processes. In CD8(+) T cells, some lncRNAs overlap with miRNAs and siRNAs, which suggests that lncRNAs act by being processed into smaller molecules that play a key role in the acquired immune response (11). The relationship between lncRNAs and dendritic cell differentiation was recently explained as follows: lnc-DC directly binds STAT3 in the cytoplasm and thus promotes STAT3 phosphorylation on tyrosine-705 by preventing STAT3 from binding to and being dephosphorylated by SHP1, and this phosphorylation activates the tyrosine kinase JAK/STAT signaling pathway and regulates dendritic cell differentiation (12). Recent studies revealed that in the immune system, in addition to directly regulating immune cells, lncRNAs can control inflammation by directly acting on inflammatory factors, and other studies found that lncRNAs can participate in the epigenetic modification of inflammation-related genes (13). Stimulation with lipopolysaccharide significantly increases the expression of lnc-IL7R in cells, and the expression of this lncRNA is associated with epigenetic regulation (14). In addition, lncRNA-Cox2 plays a key role in the regulation of IL-6 expression induced by Toll-like receptor stimulation (15). Multiple studies have shown that lncRNAs can be used as ceRNAs to indirectly regulate immunity-related mRNA expression through competitive binding to miRNA (16,17).
For example, let-7e regulates the inflammatory response of vascular endothelial cells by ceRNA crosstalk (18). The lncRNA XLOC_008466 regulates oncogenes in human lung cancer cells by targeting miR-874 (19), and lncRNAs play a positive role as competitive endogenous RNAs in gastric cancer (20). The involvement of lncRNAs in immune regulation is complicated, and many key immune regulatory lncRNAs have not yet been identified. Therefore, our study aimed to reveal some key novel lncRNAs in immune regulation.
Daidzein (4 ′ ,7-dihydroxyisoflavone) is a naturally occurring isoflavonic phytoestrogen belonging to the non-steroidal estrogen family and is mainly derived from leguminous plants (21). DA is also a major bioactive ingredient in the traditional Chinese medicine Gegen, which is frequently used for the treatment of fever, acute dysentery, diabetes, cardiac dysfunctions, and liver injury, among other conditions (22). The chemical structure of DA is similar to that of 17βestradiol (E2), and DA can selectively bind different estrogen receptors (ER) at different affinities to regulate the recruitment of co-repressors and co-activators and affect ER signaling (23). Therefore, DA exerts protective effects against some diseases that are linked to the regulation of estrogen, such as breast cancer, osteoporosis, diabetes, and cardiovascular diseases (24). DA also has several other biological activities that are independent of the ER, such as anti-inflammation and anticancer activities and protection of the skin and nerves. These beneficial effects are mainly due to the regulation of the immune response (25). DA has been shown to precisely orchestrate processes related to the regulation of the secretion of immune molecules, the proliferation and differentiation of immune cells, and immune signal pathways. Feeding mice DA after the onset of experimental allergic encephalomyelitis reduced the secretion of interferon-3 and interleukin-12, enhanced interleukin-10 production, suppressed lymphocyte proliferation, and decreased cytotoxicity (26). Studies have shown that under conditions of LPS-induced inflammation, DA could inhibit inflammatory hyperreactivity and eliminate inflammation by inhibiting the differentiation of B lymphocytes (27). In addition, mice treated with DA exhibited increased percentages of CD4(+) and CD28(+) T cells, and DA treatment regulated B lymphopoiesis and decreased the mRNA levels of RANKL in B220(+) cells (28). DA can affect the immune response by regulating the Toll-like and NF-κB pathways. We thus attempted to identify key lncRNAs involved in immune regulation by studying the daidzein-mediated changes in lncRNA expression associated with the regulation of the immune response (29,30).
Numerous studies have shown that plant extracts not only regulate the expression of coding genes but also exert regulatory effects on non-coding genes (31)(32)(33)(34). Although these studies highlight the relevance of lncRNAs in the immune regulation process, the precise molecular mechanisms remain largely unelucidated. Furthermore, the role of lncRNAs involved in the regulation of the immune response to DA remains unknown. In this study, we treated egg-laying hens with DA, screened for key mRNAs and lncRNAs related to immune regulation using transcriptome technology, and performed cell experiments to explore the specific molecular biological mechanisms underlying the regulation of immune responses by lncRNAs.

Materials
The DA used in this study was synthetically produced by the Kai Meng Co. (Xi An, Shanxi, China) Chemical Plant with a purity of 99.9%.

Feeding Experimental Design and Bird Management
The experimental animal procedures were approved by the China Agricultural University Animal Care and Use Committee (Beijing, China, permit number SYXK20130013). The experiment was performed with laying broiler breeder hens housed at a commercial farm (Zhuozhou, China) under standard conditions. After a 2-week acclimation period, a total of 480 57-week-old Ross 308 laying broiler breeder hens were allocated to two treatment groups: the DA-deficient group (DAD), which the diet is specially formulated lacking DA, and the DA-supplemented group (DS). Each treatment was replicated eight times, and each replicate included 30 broiler breeder hens. The hens were fed a nutritionally balanced corn-miscellaneous meal (CSCM) with DA added at 0 and 20 mg/kg for 8 weeks. The CSCM diets were formulated to meet the nutrient requirements of laying broiler breeders according to the NRC guidelines (1994) ( Table 1).

Sample Collection and Chemical Analysis
At the 4 and 8th weeks of the experiment, one chicken per replicate was selected and deprived of food for 8 h. One blood sample was collected from the wing vein of each replicate chick into vacuum blood collection tubes, and the serum was centrifuged at 3000× g for 15 min and stored at −20 • C until use for the detection of immunoglobulins (Igs). Another blood sample was collected from the wing vein of each replicate chick into vacuum blood collection tubes (with heparin sodium) for assessing lymphocyte proliferation and percentages. At the 8th week of the experiment, one chicken from each replicate was slaughtered. One liver sample was immediately collected from each replicate for measuring gene expression, and the samples were frozen in liquid nitrogen and stored in a freezer at −80 • C.

Serum Immunoglobulin Levels
The serum IgM and IgA levels were determined using a commercial ELISA kit (IDEXX laboratories Inc., Westbrook, Maine, USA) according to the manufacturer's recommended protocol.

Lymphocyte Classification and Proliferation
Peripheral blood mononuclear cells (PBMCs) were isolated by Ficoll density centrifugation. Briefly, heparinized blood was diluted with Hank's balanced salt solution at a ratio of 1:1 (no calcium and no magnesium, Life Technologies, Burlington, Vermont, USA) and carefully layered on top of Histopaque 1077

Blood Routine Examination
The routine blood examination was performed at Xiyuan Hospital, Beijing, China.

Next-Generation Sequencing (NGS)
Total RNA samples for sequencing were purified from 20 mg of tissue samples from eight chickens (four replicates in each treatment group) using the RNeasy Fibrous Tissue Mini messenger RNA (mRNA) extraction kit (Qiagen, Hilden, Germany) following the manufacturer's recommendations. The concentration and purity of total RNA were determined using a UV/Vis spectrophotometer (ACTGene, New Jersey, USA) at 260 nm, and the sample integrity was evaluated through a microfluidic assay using a Bioanalyzer system (Agilent Technologies, Inc., Santa Clara, California, USA  (36). Four replicates from each treatment were analyzed independently for library synthesis and sequencing, and the quality of the raw reads was assessed using FastQC (Version 0.10.1). Adapters, lowquality reads at the 3 ′ end, reads with fuzzy N bases, ribosomal RNA (rRNA), sequences shorter than 20 nt and low-quality reads (those with a Q < 20) were trimmed with the FASTX clipper (Version 0.0.13). All the double-end reads in eight samples from two treatment groups were separately aligned to the chicken reference genome (Gallus_gallus-5.0, version 81, Ensembl) using the spliced mapping algorithm in TopHat2 (version: 2.0.9) (37). Unless stated otherwise, all programs were run with the default parameters. The number of reads equivalent to mapped reads [reads per kilobase per million mapped reads (RPKM)] was used to normalize the expression of each gene. The quality of the obtained data was assessed based on the presence and abundance of contaminating sequences, the average read length, and the GC content. With the exception of the microarray design, the NGS experiment conformed to the MIAME guidelines (38).

Bioinformatics Analysis
All data processing steps for differential gene expression evaluation were performed using the Cuffdiff software package (

LncRNA Identification
To identify novel reliable lncRNA models, we considered only multiexonic transcripts and filtered them through the following highly stringent criteria: (1) size selection: only transcripts with ≥200 bp were kept; (2) read coverage threshold: transcripts with ≤3 reads were removed from our dataset; (3) open reading frame (ORF) filter: transcripts with a predicted ORF longer than 100 aa were removed; (4) known protein domain filter: transcripts were aligned to the Pfam and Swiss-Protein databases to eliminate transcripts with significant homology to known protein domains (41); and (5) protein-coding-score test: both the Coding-Non-Coding Index (CNCI) and the CPC were used to evaluate the coding potential of the candidate lncRNAs (42,43). Further analysis was conducted only with the lncRNAs that showed >1-log2 (fold change) or <-1-log2 (fold change) differential expression between the dietary groups, as demonstrated through t-tests with an adjusted P-value ≤ 0.05.

RNA Fluorescence in situ Hybridization (RNA FISH)
Cy3-labeled lncRNA-XLOC_098131 probes were obtained from RiboBio (Guangzhou, China). RNA FISH was performed using a fluorescent in situ hybridization kit (RiboBio) following the manufacturer's instructions.

′ -and 3 ′ -Rapid Amplification of cDNA Ends (RACE)
To determine the full-length sequence of lncRNA-XLOC_098131, RACE experiments were performed using the SMARTer RACE cDNA Amplification Kit (Clontech, Palo Alto, California, USA) following the manufacturer's recommended protocol. The gene-specific primers used for 5 ′ and 3 ′ RACE were GTTCTGCAGGAGAGCAGCAG and AAGTGGCCTGAGCTGGAGTC, respectively.

Quantitative Real-Time Reverse Transcription PCR
The expression levels of various genes were analyzed by qRT-PCR. The gene-specific primer sequences are shown in Table 2, and the analysis was performed on a 7500-fluorescence detection system (Applied Biosystems, Foster City, California, USA) using a commercial SYBR-Green PCR kit (Takara Bio Inc. Foster City, California, USA). According to the manufacturer's recommended protocol, the following PCR conditions were employed: 95 • C for 30 s and 40 cycles of 95 • C for 5 s and 60 • C for 34 s. In addition, melting curve analyses and subsequent agarose gel electrophoresis of the PCR products were conducted to confirm the amplification specificity. The relative gene expression data were analyzed using the 2 − Ct method.

Statistical Analysis
When a significant difference was observed between treatments, the individual treatment means were compared with Duncan's multiple comparison using SPSS Version 18.0. A P-value < 0.05 was considered to indicate statistical significance, and P < 0.1 was considered to indicate a trend toward statistical significance.  WBCP, white blood cell percentage (%); PLT, platelet count (10 9 /L); MPV, mean platelet volume (fL); PDW, platelet volume distribution width (%); PCT, platelet thrombocytocrit (%); NEUT (%), neutrophil ratio (%); lymph%, lymphocyte ratio (%); MONO%, monocyte ratio (%); EOS%, eosinophil ratio (%); LUC%, proportion of large unstained cells (%); BASO%, basophil ratio (%). The data are expressed as the means with SEMs and P-values. * Values with the same superscript within a column are not significantly different at P < 0.05, values with the different superscript (a,b) are significantly different at P < 0.05. mg/kg increased (P < 0.05) the ratio of B cells on the 4 and 8th weeks of the experiment (Figure 1C), but no significant effect (P > 0.05) on the percentage of CD3 ( Figure 1D) and the CD4/CD8 ratio was observed ( Figure 1E). We further examined the lipopolysaccharide (LPS)-and concanavalin A (CON A)-mediated stimulation of the lymphocyte proliferation rate and found that the DA treatment increased (P < 0.05) the proliferative responses of B cells on the 8th week of the experiment, as shown by changes in the LPS stimulation index (SI) values obtained in vitro using peripheral blood ( Figure 1F).

Serum Immunoglobulin Levels, Peripheral Blood Lymphocyte Classification, and Lymphocyte Proliferation in Hens
In contrast, no difference in lymphocyte proliferation was observed following CON A stimulation.

Routine Blood Examination
The routine blood examination results ( Table 3) showed that on the 4th week of the experiment, the platelet count (PLT) and pressure (PCT) in the DA-supplemented (DS) group were decreased (P < 0.05) compared with those found in the DAD group, and the DS group also exhibited a slight increase (P < 0.1) in the neutrophil ratio (NEUT%). Furthermore, DA significantly increased the percentage of monocytes (MONO%) and the proportion of large unstained cells (LUC%). On the 8th week of the experiment, the MONO and LUC% in the DS group were higher (P < 0.05) than those in the DAD group, and the NEUT% showed an increasing tendency (P < 0.1). In addition, the DA treatment did not affect the white blood cell percentage (WBCP), average platelet volume (MPV), platelet distribution width (PDW), lymphocyte ratio (lymph%), eosinophilia ratio (EOS%), or basophil ratio (BASO%).

RNA-Seq Statistics
In this study, we established eight cDNA libraries from the livers of the birds in the DAD and DS groups with four replicates per group. The RNA-Seq assay generated 79,846,624-93,885,866 raw reads per library, with an average of 88,374,633 and 83,729,690 paired-end reads in the DAD and DS groups, respectively. After filtering the low-quality reads, the average numbers of clean reads in the DAD and DS groups were 82,831,728 (93.73%) and 74,604,252 (89.10%), respectively. The clean reads were used in all further analyses. Approximately 83.0% of the reads in each library were uniquely mapped to the galGal4 assembly of the chicken genome, and the average mapping rates of the DAD and DS groups were 92.7 and 93%, respectively ( Table S1). As shown in Figure 2A, 146 differentially expressed genes (DEGs; 64 upregulated and 82 downregulated) were identified in the comparison of the DAD group with the DS group (FDR ≤ 0.05, fold change ≥1.5 or ≤0.7) using Cuffdiff software. The expression levels and fold changes of the DEGs are shown in Table S2, and the DEGs are also shown in a heat map ( Figure 2B) and volcano plot ( Figure 2C). Moreover, the transcriptome data of the DEGs were validated by RT-PCR (Figure 2E), and the results were consistent with the RNA-Seq results ( Figure 2D).

Functional Categorization of DEGs and Pathway Analysis
We performed a functional enrichment analysis using the 146 DEGs, and results showed revealed 1,853, 259, 407, and 21 terms associated with biological process regulation, cell components, molecular function, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively, and 498, 59, 105, and 3 significantly enriched terms in these categories, respectively ( Figure 2F). The DA treatment mainly affected the GO terms: oxidation-reduction process, fatty acid metabolic process, response to peptides, response to oxygen-containing compounds, antigen processing, presentation of peptide antigens via MHC class I and other biological processes; MHC class I protein complex, extracellular region portion, extracellular membrane-bound organelle, extracellular region, membrane-bound vesicle, and other cell components; and peptide antigen binding, oxidoreductase activity, antigen binding, receptor binding, and other molecular functions ( Figure 2G). A complete list of the different gene ontology (GO) terms identified in this analysis is provided in Table S3.
The 498 different biological process GO terms were further analyzed using OmicsBean software. According to the number  of enriched genes (Figures 3A,C), the DEGs identified from the comparison of the DS group with the DAD group were mainly enriched in multicellular organism development (28 genes), system development (27 genes), and oxidation-reduction process (22 genes). According to the different levels of the GO terms (Figure 3D), most differential genes were related to stimuli (level 2 GO), and continued subdivision of the grades showed that the differential genes were involved in the organism response to multiple stimuli (level 3-8 GO). A sorting of the different GO terms based on the P-value ( Figure 3E) revealed that these DEGs are most likely involved in biological processes related to antioxidation and the immune response. The PAS Z-score (Figure 3B) combines the P-value and fold change data and showed that compared with the DAD group, the DS group had significantly more immune response-related GO terms associated with enhanced immunity and significantly fewer terms related to lipid oxidation and other processes, which indicates that DS can be associated with an enhancement of the body's antioxidant capacity.
The analysis of the 59 different cell components according to the number of enriched genes (Figures 4A,C) revealed that the most abundant DEGs were enriched in intracellular organelles (60 genes) and cytoplasm (52 genes). In terms of the different GO classification levels (Figure 4B), the genes showing the highest differential expression are associated with organelles (level 2 GO), and further subdivisions showed that the DEGs are associated with membrane-bound organelles (level 3 GO). Sorting the different GO terms based on ascending P-values ( Figure 4D) revealed that the genes showing the greatest differential expression are likely enriched in the cellular component of the MHC protein complex.
The analysis of the 105 different molecular function-related GO terms showed that according to the number of enriched genes (Figures 5A,C), the most abundantly enriched DEGs are associated with cation binding (27 genes) and receptor binding (12 genes). According to the classification of the different GO levels (Figure 5B), most DEGs are related to continued subdivision identified receptor binding (level 4 GO) and peptide antigen binding (level 5 GO). By distinguishing the difference in the P-value of the GO terms (Figure 5D), we found that the molecular functions of these DEGs are peptide antigen binding and oxidoreductase activity.
Based on the above results, it can be hypothesized that the DA treatment had an impact on the immune response, and a statistical analysis of immune regulation-related GO terms (Table 4) showed that the differential genes are involved in the positive regulation of responses to stimulus, interleukin-8 endocytosis involved viral entry into host cell, chemical stimulus, drug stimulus, endogenous stimulus, abiotic stimulus, oxidative stress, reactive oxygen species, and other immune responses. The following immunity-related GO terms were also regulated: negative regulation of B cell apoptotic processes, antigen binding, interleukin-1-mediated signaling pathway, transport of virus, regulation of MAPK cascade, Through a KEGG pathway cluster analysis, we also found that the DEGs identified from the comparison of the DAD group with the DS group regulated the Toll-like receptor signaling pathway, salmonella infection pathway and herpes simplex infection pathway ( Figure 6A). A correlation analysis ( Figure 6B) showed that the DEGs are mainly involved in the Toll-like receptor signaling pathway, salmonella infection, herpes simplex infection, MAPK signaling pathway, pyruvate metabolism, cysteine and methionine metabolism, and cardiac muscle contraction, and the key genes regulating these processes are JUN, FOS, EGR1, IRS2, CTSK, and CYR61. Furthermore, an analysis of the most significant Toll-like receptor signaling pathway in KEGG (Figure 6C) revealed that the CTSK (upstream) gene in the pathway and the dimer AP-1, which is composed of proteins encoded by the JUN and FOS (downstream) genes in the pathway, play an important role in DA-mediated regulation of immunity. AP-1 can respond to a variety of stimuli, including cytokines, growth factors, stress, and bacterial and viral infections, by regulating the expression of downstream genes. This finding is consistent with the results of the GO analysis, which showed that DA addition mainly enhanced the body's ability to respond to multiple stimuli. Therefore, FOS and JUN are key genes that enable the DAmediated enhancement of the immune response and innate immunity in hens.

Differential lncRNA Selection
We performed a lncRNA screen of the DS and DAD groups and identified 1,095 lncRNAs (Table S4). An analysis of the transcript lengths and transcript exons of all lncRNAs (Figures 7A,B) revealed that the frequency of these lncRNA transcripts was mainly concentrated in 200-5,000 units and that the number of transcript exons mainly ranged from 2 to 4. LncRNAs that demonstrated differential expression >1-log2 (fold change) or <-1-log2 (fold change) between the two groups, as determined through t-tests with an adjusted P-value ≤ 0.05, were selected. In the DS group, nine and 22 lncRNAs were upregulated and downregulated, respectively compared with their levels in the DAD group ( Table 5).

Cis-Regulation Analysis of Differential lncRNAs
A screen of all mRNAs within 10 kb from the 31 differentially expressed lncRNAs yielded 392 mRNAs (Table S5). A comparing of these 392 mRNAs with transcriptome mRNA expression data did not result in a significant difference in the expression of these 392 mRNAs between the DS and DAD groups, suggesting that the differentially expressed lncRNAs identified from the comparison

Potential Binding Sites in the lncRNAs and the mRNAs
To explore the mechanism underlying the action of the differentially expressed lncRNAs and identify the targets, we analyzed the potential binding sites in the lncRNAs and the mRNAs of FOS and JUN, which are key genes that enable the DA-mediated enhancement in the immune response, using mirDB (http://www.mirdb.org/miRDB/). The results showed that lncRNA-XLOC_098131 has two sites for binding to miR-548s (miR-548e-5p), and the seed locations are 484 and 2,722 ( Figure 7C). In addition, lncRNA-XLOC_098131 has two binding sites for miR-548s (miR-548ah-5p), and the seed locations are 1,884 and 3,039 ( Figure 7D). Furthermore, FOS has a binding site for miR-548s (miR-548ay-3p) located at 1,544 (Figure 7E), and lncRNA-XLOC_098131 has two binding sites for miR-1180s, with seed locations of 961 and 2,501 ( Figure 7F). Moreover, JUN has a binding site for miR-1180s located at 1,506 ( Figure 7G). These results suggest that lncRNA-XLOC_098131 might regulate the expression of FOS and JUN by competitively binding to miR-548s and miR-1180s.

Characterization of the lncRNA-XLOC_098131 Sequence
The 5 ′ and 3 ′ rapid amplification of cDNA ends (RACE) analyses demonstrated that lncRNA-XLOC_098131 is a transcript consisting of 3,321 nucleotides, and similar to many lncRNAs, this transcript is polyadenylated. The full length of lncRNA-XLOC_098131 is shown in Table S6. The full length of lncRNA-XLOC_098131 was then used in a BLASTN search with the genomic sequences of chicken, human, mouse, pig, cow, and sheep using the online tool Ensemble (http:// asia.ensembl.org/) and considering an e-val<0.05. This search yielded more than 100 chicken alignment sequences, one human alignment sequence, one mouse alignment sequence, seven pig alignment sequences, one cow alignment sequence, and four sheep alignment sequences. The results of the first ten comparisons between chickens and other animals are shown in Table 6. For the chicken alignment, the sequence at genomic These results are consistent with the transcriptome results, which suggests that lncRNA-XLOC_098131 is located in 8:25819102-25826005. Furthermore, analyses of the high-scoring segment pair (HSP) distribution on the genome ( Figure 8A) and HSP distribution on the query sequence ( Figure 8B) revealed that lncRNA-XLOC_098131 is located on chromosome 8 in chicken. As shown in Figure 8C, lncRNA-XLOC_098131 is located on chromosome 8 at position 25819102-25821666 (illustrated by the red line), which is a non-coding area between the two protein-coding genes ENSGALG00000010818 and protein kinase AMP-activated catalytic subunit alpha 2 (PRKAA2) and near the ENSGALG00000044889 gene, which encodes complement C8 alpha chain (C8A), C8B, and DAB1. We performed a further BLASTN analysis, as shown in Figure 8D, and the BLASTN against Chicken Gallus_gallus-5.0 (Genomic sequence) showed no BLASTN pairing with lncRNA-XLOC_098131 and a forward strand orientation in the chicken genome sequence, but BLASTN revealed a high degree of matching with the reverse direction of the chicken genome sequence (lncRNA- XLOC_098131 is a reverse sequence). Furthermore, the BLASTN analysis of lncRNA-XLOC_098131 with the chicken brain Iso-Seq from the European Nucleotide Archive (https://www.ebi.ac. uk/ena) (44) resulted in a high degree of matching, which was slightly lower than that obtained with Chicken Gallus_gallus-5.0 (Genomic sequence); thus, the tissue specificity of lncRNA-XLOC_098131 is low in chickens. The BLASTN analysis of lncRNA-XLOC_098131 with chicken embryo Iso-Seq from the European Nucleotide Archive (https://www.ebi.ac.uk/ena) (44) resulted in a lower degree of matching, showing that the expression level of lncRNA-XLOC_098131 gradually increases with chicken development and highlighting the high temporal specificity of lncRNA-XLOC_098131. The analysis of lncRNA-XLOC_098131 sequence variants from all sources showed that the variability in the sequence of lncRNA-XLOC_098131 is low, indicating that lncRNA-XLOC_098131 has high sequence stability. Furthermore, the percentages of bases G and C in the sequence (%GC) are high, indicating that the sequence of lncRNA-XLOC_098131 has high stability. In addition, consistent with the fact that lncRNA-XLOC_098131 is a non-coding RNA, the Coding Potential Calculator (CPC) computational algorithm predicted that lncRNA-XLOC_098131 has a very low coding potential, similar to pouBW1, which is a well-documented Gallus lncRNA (45) (Figure 8E). These results indicate that lncRNA-XLOC_098131 is a real long non-coding reverse RNA sequence with good stability and low specificity between tissues but with a high temporal specificity and a certain degree of conservation among species.

Cellular Localization of lncRNA-XLOC_098131
To determine the cellular localization of the lncRNA-XLOC_098131 transcript, the nuclear, and cytosolic RNAs from chicken primary hepatocytes were isolated, and the expression of lncRNA-XLOC_098131 transcripts in both subcellular locations was measured. The qRT-PCR data showed that lncRNA-XLOC_098131 transcripts were more highly expressed in the cytosol than in the nucleus (Figure 8F). A FISH analysis of lncRNA-XLOC_098131 in hepatocytes also indicated that lncRNA-XLOC_098131 is mainly located in the cytoplasm (Figure 8G).

DA Treatment Can Upregulate the Expression of lncRNA-XLOC_098131 in vitro
Chicken primary hepatocytes were treated with 0, 1, 5, 10, 50, or 100 µm/L DA for 24 h, and the total RNA from the cells was then extracted for RT-PCR. The 10, 50, and 100 µm/L DA treatment groups exhibited significantly increased expression of lncRNA-XLOC_098131 compared with the control group ( Figure 8H). We then examined cultured cells treated with 10 µm/L DA for 0, 12, 18, 24, and 30 h and found that the cells treated with 10 µm/L DA for 12, 18, 24, or 30 h exhibited significantly increased expression of lncRNA-XLOC_098131 compared with the control group ( Figure 8I).

Binding Between lncRNA-XLOC_098131 and miR-548s
To assess the existence of direct binding between miR-548s and lncRNA-XLOC_098131 at endogenous levels, we constructed luciferase reporters containing lncRNA-XLOC_098131 that contained wild-type (WT) or mutated miR-548s-binding sites. We found that the overexpression of miR-548s reduced the luciferase activities of the WT reporter vector (Figure 9A) but not that of the empty vector ( Figure 9C) or that of the mutant reporter vector (Figure 9B). These data demonstrated that the predicted miR-548s could regulate luciferase expression with the 3 ′ UTR of lncRNA-XLOC_098131 (P < 0.05). Mutation of the binding site resulted in the disappearance of this regulation, confirming that miR-548s regulates luciferase through this binding site and that lncRNA-XLOC_098131 might be a target gene of miR-548s.

Binding Between FOS and miR-548s
To validate the hypothesized direct binding between miR-548s and FOS at endogenous levels, we constructed luciferase reporters containing FOS that contained WT or mutated miR-548s-binding sites. We found that miR-548s overexpression reduced the luciferase activities of the WT reporter vector ( Figure 9D) but not that of the empty vector ( Figure 9F) or that of mutant reporter vector ( Figure 9E). These data demonstrated that the predicted miR-548s could regulate luciferase expression with the 3 ′ UTR of FOS (P < 0.05). This regulation disappeared if the binding site was mutated, confirming that miR-548s regulates luciferase through this binding site and that FOS might be a target gene of miR-548s.

LncRNA-XLOC_098131 Upregulates FOS Expression by Serving as a Competitive Endogenous RNA
Because lncRNA-XLOC_098131 shares the regulatory miR-548s with FOS, we questioned whether lncRNA-XLOC_098131 could modulate FOS by regulating miR-548s. As shown in Figure 9G, the overexpression of lncRNA-XLOC_098131 WT, but not that of the mutant or lncRNA-XLOC_098131, increased the FOS transcript level. In addition, the co-transfection of the lncRNA-XLOC_098131 WT overexpression vector, FOS luciferase vector, and sufficient miR-548s resulted in no difference in the luciferase activity of the FOS luciferase vector compared with the control group, and the above treatments had no effect on the psiCHECK-2 luciferase activity.

Binding Between lncRNA-XLOC_098131 and miR-1180s
To validate the existence of direct binding between miR-1180s and lncRNA-XLOC_098131 at endogenous levels, we constructed luciferase reporters containing lncRNA-XLOC_098131 that contained WT or mutated miR-1180sbinding sites. However, compared with the control, the overexpression of miR-1180s could not reduce the luciferase activities of the WT reporter vector (Figure 9H), the mutant reporter vector (Figure 9I) or the empty vector ( Figure 9J). This result showed that lncRNA-XLOC_098131 cannot bind to miR-1180s through our predicted binding site. The regulation Genomic location, location of the hit on the genome; Overlapping gene(s), genes that overlap the genomic location hit; Orientation, strand of the genome the query sequence is found on; Query start and end, positions within the query sequence where the alignment begins and ends (1-based, inclusive); Length, full length of the alignment, including all gaps in either the query or the alignment; Score, BLAST score calculated from the alignment; E-val, probability that the alignment between the query sequence and subject sequence is due to chance; %ID, percentage of the aligned query sequence that is identical to the subject sequence.
Frontiers in Immunology | www.frontiersin.org  respond to stimuli are important manifestations of a good immune response (48,49). Studies investigating lymphocytes following DA treatment have shown that the ability of offspring broilers' B, but not T, lymphocytes to differentiate might be enhanced by adding DA to the breeder's diet at a dose of 20 mg/kg (35). In our study, 20 mg/kg DA increased the number of B lymphocytes in the peripheral blood in hens but had no significant effect on CD3(+) subset of T lymphocytes or on the CD4(+)/CD8(+) T cell ratio. Furthermore, the DA treatment resulted in an increase in the proliferative responses of B cells, as shown by the changes in the LPS SI values obtained in vitro using peripheral blood from hens on the 8th week of the experiment, but had no effect on the proliferative responses of T cells, which is consistent with the previous finding that 0.5 mg/kg DA has no significant effect on T lymphocyte proliferative responses in 5-to 6-week-old boars (50). These results suggest that the effect of DA on lymphocytes is mainly manifested in the regulation of B lymphocytes and that the effect on T lymphocytes is markedly lower. Enhanced B lymphocyte differentiation produces more plasma cells that can secrete more immunoglobulins, which is consistent with the finding that the DA treatment increased the serum IgA and IgM levels in hens. The transcriptome results showed that compared with the DAD group, the DS group exhibited significantly upregulated IRS2 gene expression, which is associated with the inhibition of B lymphocyte apoptosis. Additionally, DA treatment can promote IgA and IgM secretion by regulating the expression of antigen-presenting related genes, such as LOC100859408, ENSGALG00000004772, and LOC768350. Previous studies have found that soy isoflavones, including DA, can regulate the expression of the IRS2 gene, which might affect B lymphocytes (51). Therefore, treatment with 20 mg/kg DA can reduce B lymphocyte apoptosis, promote the differentiation of B lymphocytes into plasma cells that, in turn, secrete more immunoglobulins, promote the presentation of antigens and improve humoral immunity in hens. Properly elevated levels of neutrophils and mononuclear cells in the blood can increase the body's ability to respond to external stimuli and bacterial infections and enhance immunity (52,53). Based on our results, the DS group tended to exhibit an increase in the number of neutrophils in the blood of hens on the 4 and 8th weeks of the experiment (P < 0.1), and the number of monocytes in the DS group was significantly higher than that in the DAD group. The RNA-Seq results showed that the DEGs identified from the comparison of two groups such as EGRI1, CAV2, JUN, FOS, SOD3, and CDH2 are involved in the regulation of the interleukin-8 endocytosis viral entry into the host cell, the regulation of Toll-like receptor 21 signaling pathway and MAPK signaling pathway, which can modulate the neutrophils and macrophages. These results show that DA treatment can enhance the body's resistance to viruses and harmful bacteria, enhance the immune response to various stimuli, and improve the body's innate immunity. Previous studies have demonstrated that DA can regulate the content of monocytes and neutrophils and relieve toxic symptoms in mice (54) and that DA can play a positive role in regulating LPS-induced lung injury by regulating the neutrophil content (55).
The KEGG and correlation analyses performed in this study showed that the DS group had an enhanced ability to respond to multiple stimuli compared with the DAD group, and this enhanced ability was regulated by the MAPK signaling pathway, Toll-like pathway and salmonella-related pathway, which can regulate downstream genes to improve humoral and innate immunity, as discussed above. We found that the key genes regulating these processes are JUN, FOS, EGR1, IRS2, CTSK, and CYR61, and the dimer AP-1 composed of JUN and FOS plays a crucial regulatory role in the pathway. AP-1 can regulate the expression of multiple genes to respond to multiple stimuli, including cytokines, growth factors, stress, and bacterial and viral infections. Some studies have found that an appropriately high expression of AP-1 can enhance the body's immune and antibacterial ability and improve its ability to respond to multiple stimuli (56). The AP-1 family of transcription factors can activate Tolllike receptor agonists and positively regulate interleukin-4 (IL-4) to activate macrophages, thereby enhancing antibacterial activity (57). AP-1 plays a key role in regulating immunity as an intermediate regulatory factor; for example, CR3 and Dectin-1 synthesize macrophage cytokines through synergistic activation of the lipid and Syk-JNK-AP-1 pathways (58). Previous studies have shown that multiple plant estrogens can activate AP-1 expression (59), and the treatment of cells with DA significantly increases their expression of the FOS gene (60). DA can also increase FOS gene and protein expression in humans (61). Therefore, in our study, FOS and JUN were identified as the key genes controlling the immunity of hens following DA treatment. The addition of 20 mg/kg DA to the breeder hen diet can upregulate the expression of FOS, JUN, and AP-1 to regulate the Toll-like receptor signaling pathway, MAPK pathway, salmonella infection signal pathway, and related mRNA expression to enhance macrophage activity and increase blood neutrophils and mononuclear cells. This process improves the body's ability to respond to cytokines, growth factors, stress, and bacterial and viral infections and improves the innate immunity.
This study explored the role of lncRNAs in the regulation of key immunity-related genes following DA treatment. We sequenced lncRNAs that showed differential expression between the DA treatment and control groups, identified 1,096 lncRNAs, and found that the frequency of these lncRNA transcripts was mainly concentrated in the 200-5,000 unit range and that the number of transcript exons ranged from 2 to 4, which is consistent with the characteristics of lncRNAs (6). Furthermore, nine and 22 lncRNAs were significantly upregulated and downregulated after DA treatment. We further analyzed the relationship between the differentially expressed lncRNAs and key mRNAs. Previous studies have shown that lncRNAs can influence the expression of adjacent mRNAs that are located within 10 kb of the lncRNAs through cis-regulatory effects; for example, the lncRNA Malat1 plays a cis-regulatory role in the adult mouse (62), and the lncRNA GAPLINC regulates mRNA by cis-regulation in gastric cancer (63). In this study, we selected 31 differentially expressed lncRNAs for the analysis and found that 392 mRNAs are located near these lncRNAs. However, none of these mRNAs exhibited significantly differential expression between the DAD and DS groups, suggesting that the DA treatment could not regulate the expression of mRNAs by the cis-regulatory function of the lncRNAs that were found to be differentially expressed between the two groups. Many studies have found that lncRNAs function as competitive endogenous RNAs to regulate mRNA expression. For example, the lncRNA CCAT1 upregulates the proliferation and invasion of melanoma cells by inhibiting miR-33a (64), the lncRNA XIST regulates AKT expression by competitively binding miR-494 (65), and the lncRNA TUG1 binds to microRNA-9 to regulate the expression of Bcl2l11 (66). By predicting binding sites using an online tool, we found that lncRNA-XLOC_098131, which is the most highly upregulated lncRNA after DA treatment, might regulate FOS expression by competitively binding to miR-548s and regulate JUN expression by competitively binding to miR-1180s. Therefore, lncRNA-XLOC_098131 was selected for further experimentation. First, we also verified that DA treatment can upregulate lncRNA-XLOC_098131 expression in vitro. We then found that lncRNA-XLOC_098131 is a novel transcript located on chromosome 8 in chickens between 25819102 and 25826005, and a BLASTN analysis showed that this fragment is located in a non-coding sequence region, indicating that lncRNA-XLOC_098131 is a non-coding transcript. For further confirmation, we used the CPC computational algorithm and predicted that lncRNA-XLOC_098131 has a very low coding potential, similarly to pouBW1, a well-documented Gallus lncRNA. As several studies have revealed, lncRNAs have low tissue specificity in animals (67), and we found that a BLASTN analysis of lncRNA-XLOC_098131 with chicken brain Iso-Seq resulted in a high degree of matching, which was slightly lower than that obtained from the BLASTN analysis with the chicken genomic sequence. These findings demonstrate the low tissue specificity of lncRNA-XLOC_098131. However, further analysis revealed that lncRNA-XLOC_098131 has low expression in the chicken embryo and that its expression gradually increases with development, highlighting the high temporal specificity of lncRNA-XLOC_098131, and these results are consistent with previous research (68). LncRNAs are known to have high stability and a certain degree of conservation (69). We found that the variability of the sequence of lncRNA-XLOC_098131 is low and that the percentages of the bases G and C in the sequence (%GC) are high, indicating that the sequence of lncRNA-XLOC_098131 exhibits high stability. These results indicate that lncRNA-XLOC_098131 is a real long non-coding reverse RNA sequence with good stability and low specificity between tissues but with a high temporal specificity and a certain degree of conservation among species. LncRNAs can be distributed in the nucleus or cytoplasm, and lncRNAs in the cytoplasm can regulate the expression of mRNAs as ceRNAs (70). In this study, we conducted PCR and FISH experiments to show that lncRNA-XLOC_098131 is mainly located in the cytoplasm, where lncRNAs can bind to miRNAs through their binding sites and thereby eliminate the miRNA-mediated inhibition of mRNAs, which can also bind to miRNAs (71). We successfully demonstrated that both lncRNA-XLOC_098131 and FOS are target genes of miR-548s and that lncRNA-XLOC_098131 can bind to mir-548s through its binding sites to reduce the level of free miR-548s, which in turn reduces the binding of FOS and miR-548s, eliminating the miR-548smediated inhibition of FOS and upregulating the expression of FOS. Furthermore, when found at a sufficient high level, the binding of miR-548s to lncRNA-XLOC_098131 did not affect the intracellular free mir-548s content; thus, there was no change in the number of miR-548s that bound to FOS, and the expression of FOS was unchanged. These results confirmed that lncRNA-XLOC_098131 regulates FOS expression by competitively binding to miR-548s. However, subsequent experiments demonstrated that lncRNA-XLOC_098131 cannot regulate JUN gene expression by competitively binding to miR-1180s. LncRNAs can also affect mRNA expression through epigenetic or transcriptional regulation (6), and the molecular mechanism through which lncRNAs regulate JUN expression needs further investigation.
In summary, our study demonstrates that a novel lncRNA activated by DA can stabilize FOS mRNA by serving as a competitive endogenous RNA and can affect AP-1 with JUN to regulate the MAPK signaling, Toll-like receptor signaling, and salmonella signaling pathways in vivo. Dietary supplement DA can improve the numbers of blood neutrophils and mononuclear cells, and improve the innate immunity response. In addition, by increasing the content and differential ability of B lymphocytes and improving antigen presentation and immunoglobulin secretion, DA enhances the body's humoral immunity. In conclusion, our study shows that the novel lncRNA-XLOC_098131 plays a key role in the molecular mechanism underlying the immune function mediated by DA. This study provides new insight into the molecular mechanisms through which DA affects the Toll-like receptor signaling pathway-related immune response and the mechanisms through which lncRNAs regulate mRNA expression. As direct target genes, miR-548s plays important roles in mediating the regulation of FOS by the lncRNA-XLOC_098131 to affect immunity. These data constitute high-quality resources for future genome and functional research.

ETHICS STATEMENT
The experimental animal procedures were approved by the China Agricultural University Animal Care and Use Committee (Beijing, China).

AUTHOR CONTRIBUTIONS
HF designed and performed experiments, analyzed data, and wrote the paper. ZeL, LG, BZ, and ZhL performed some experiments, and analyzed some data. CN, JL, and MY analyzed some transcriptome data. BS, GL, and DT analyzed some data. JG, SY, and YW performed some experiments. YG initiated the study, designed animal experiments, analyzed data, and wrote the paper.

FUNDING
This study was financially supported by the National Natural Science Foundation of China: 31872375.

ACKNOWLEDGMENTS
We thank Professor JL for assisting with the transcriptome data analysis in this study. All authors read and approved the final version of the manuscript.