International Journal of Genomics
Volume 2019 (2019), Article ID 6905194, 11 pages
https://doi.org/10.1155/2019/6905194

Research Article

RNA-Seq-Based Gene Expression Pattern and Morphological Alterations in Chick Thymus during Postnatal Development

Zhouyiyuan Xue,¹ Abdur Rahman Ansari,² Xing Zhao,¹ Kun Zang,¹ Yu Liang,¹ Lei Cui,¹ Yafang Hu,¹ Ranran Cheng,¹ Xiaolong Zhang,¹ Juming Zhong,1,3 and Huazhen Liu¹

¹Department of Basic Veterinary Medicine, College of Animal Science and Veterinary Medicine, Huazhong Agricultural University, Wuhan, Hubei 430070, China
²Section of Anatomy and Histology, Department of Basic Sciences, College of Veterinary and Animal Sciences (CVAS) Jhang, University of Veterinary and Animal Sciences (UVAS), Lahore, Pakistan
³Department of Anatomy, Physiology and Pharmacology, College of Veterinary Medicine, Auburn University, Auburn, USA

Correspondence should be addressed to Huazhen Liu; lhz219@mail.hzau.edu.cn

Received 15 September 2018; Revised 20 December 2018; Accepted 25 February 2019; Published 21 April 2019

Academic Editor: Wilfred van IJcken

Copyright © 2019 Zhouyiyuan Xue et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The thymus is a lobulated unique lymphoid immune organ that plays a critical role in the selection, development, proliferation, and differentiation of T cells. The thymus of developing chickens undergoes continued morphological alterations; however, the biomolecular and transcriptional dynamics of the postnatal thymus in avian species is not clear yet. Therefore, the thymuses from chickens at different stages of development (at weeks 0, 1, 5, 9, 18, and 27) were used in the present study. The RNA-seq method was used to study the gene expression patterns. On average, 24120819 clean reads were mapped, differentially expressed genes (DEGs) were identified on the basis of log values (fold change), including 744 upregulated and 425 downregulated genes. The expression pattern revealed by RNA-seq was validated by quantitative real-time PCR (qPCR) analysis of four important genes, which are PCNA, CCNA2, CCNB2, and CDK1. Thus, the current study revealed that during postnatal development, the thymus undergoes severe atrophy. Thymus structure was damaged and gene expression changed dramatically, especially at the 27th week of age. Moreover, we found significant changes of several signaling pathways such as the cytokine-cytokine receptor interaction and cell cycle signaling pathways. Hence, it may be inferred that those signaling pathways might be closely related to the postnatal chicken thymus development.

1. Introduction

The thymus develops as a gland and is considered as a lobulated unique lymphoid organ [1, 2]. As a primary immune organ, it acts as the main player in the selection, development, proliferation, and differentiation of T cells [3, 4]. Although it has prime importance in immune functions [5, 6], it experiences a conserved evolutionary process of shrinkage with aging in almost all vertebrates that causes a decrease in thymic tissue mass as well as alters its architecture [7].

The development of T cells begins with the migration of multipotent progenitors from the bone marrow to the thymus [8]. The complicated microenvironment in the thymus provides the site for naïve T cells to mature under a specific cascade of transcriptional factors [9]. As a primary immune organ, it plays a crucial role in acquired immunity with the expression of several intrathymic transcripts [10]. The thymus morphology changes with age, and its cellular components may also play an important role in prenatal immunity [11]. Aging is described as the gradual loss of physiology and cellular morphology based on cell-to-cell transcriptional variability in many mammalian tissues [12]. Studies showed that with the increase in age, the thymus undergoes a severe involution, manifesting as increased cell death and fibrosis, concomitant with a loss of T cell production [13–15]. But the molecular mechanism and when this process initiates remain undetermined [16]. Since thymic morphology and cell types are evolutionarily preserved in several vertebrate species [17, 18], the thymus might have similar age-related changes in chickens [18].

The transcriptional signature of age-related thymic involution shows altered expression of certain genes and transcriptional factors [16]. The development of the thymus is regulated by many transcriptional factors such as the thymic epithelial cell- (TEC-) specific transcription factor Forkhead box N1 (Foxn1), transcription factor Fezf2, and thymic stromal lymphopoietin (TSLP)/TSLP Receptor (TSLPR) through several signaling pathways in mammalian species [19–21]. Although recent studies have provided much information about several key genes, signaling pathways functioning during thymus postnatal development need further studies. Furthermore, intrinsic mediators responsible for biological pressure and mechanisms of decreased thymic function as well as the molecular and transcriptional factors of thymus development and its atrophy are still unclear in avian species. Therefore, in the current study, we conducted RNA-seq to provide a macroscopic picture of the gene expression and explored the mechanism on how aging affects transcriptional dynamics using RNA sequencing of thymuses from developing chicks without immunological activation.

2. Materials and Methods

2.1. Ethics Statement

The current scientific investigation was conducted in accordance with the rules and regulations of the ethics committee for use of animals, Huazhong Agricultural University (HZAU), Wuhan, China.

2.2. Animals

One-day-old, healthy Cobb500 strain broiler chicks were purchased from the Wuhan Zhengda chicken breeding company. The newly hatched chicks were reared under conventional housing conditions. For transcriptome analysis of the chick thymus during aging, three chicks per time point (0, 1, 5, 9, 18, and 27 weeks) were sacrificed. Thymuses were harvested from each bird and weighed at certain defined time points. The thymus index was calculated using the following formula: thymus index of thymus (g)/body weight (kg) [22, 23]. For histomorphological analysis, parts of freshly harvested thymuses were fixed in 4% paraformaldehyde solution, embedded into paraffin wax, and stored at 4°C. For molecular studies and gene expression analysis, parts of freshly harvested thymuses were snap frozen in liquid nitrogen and stored at -80°C.

2.3. Hematoxylin and Eosin Staining

The paraffin wax thymic tissue blocks were cut into 4 μm thick tissue sections. After deparaffinization in xylene, the tissue sections were rehydrated in a graded series of alcohol to water. Staining of tissue sections was performed with Harris hematoxylin solution for 5 min and then differentiated in acid alcohol for 10 s. The tissue sections were rinsed under running tap water for 15 min and subsequently stained in eosin solution for 1 min. After dehydration in alcohol, the tissue sections were cleared in xylene and finally mounted with glass coverslips. Stained tissue sections were thoroughly examined by light microscopy (Olympus BX51, Tokyo, Japan) with an attached digital camera (DP72, Olympus).

2.4. RNA-Seq and Data Analysis

Tissue specimens were sent to BGI Genomics Co. Ltd. to conduct transcriptome sequencing. In this study, different time points were compared with the 0 week to identify DGEs. Three biological replicates were studied for each time point. Following construction of strand-specific libraries, sequencing was done on an Illumina HiSeq 4000 platform by BGI Genomics Co. Ltd. (Shenzhen, China). Differentially expressed genes (DEGs) were screened using the NOISeq method with absolute fold and divergent [24, 25]. The read length was 50 bp. Perl was used for handling of raw data, and the quality of the obtained data was carefully checked with FastQC 0.11.2 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Clean reads in this study were aligned to the specific chicken genome database (release: Gallus gallus 4.0) from NCBI using HISAT, by allowing one mismatch. The annotation of genome mapped data was done using the GFF3 file of the Gallus gallus 4.0 genome (ftp://ftp.ncbi.nlm.nih.gov/genomes/Gallus_gallus/). The expression of the transcript was represented in terms of FPKM (fragments per kilobase of exon per million fragments mapped). Gene expression profile analysis, GO analysis (http://geneontology.org/page/go-enrichment-analysis), and KEGG analysis were conducted depending on the data. In this study, the significantly enriched signaling pathways were sorted out when .

2.5. RNA-Seq Data Validation by Quantitative Real-Time PCR (qPCR)

For the validation of RNA-seq data, total RNA from thymic tissue was isolated and the genomic DNA was then removed by treating with RNase-free DNase I (Fermentas, St. Leon-Rot, Germany). Synthesis of first strand cDNA was accomplished using a RevertAid First Strand cDNA Synthesis Kit (Fermentas, St. Leon-Rot, Germany). The total volume (10 μL) of the qPCR reaction mixture contained 2 μL of each forward and reverse primer, 1 μL of template cDNA, and 5 μL of SYBR Select Master Mix for CFX (Applied Biosystems, Inc., USA). The qPCR reaction was accomplished by the following conditions: predenaturation at 95°C (5 min), followed by denaturation (40 cycles) at 95°C (30 s), annealing at 60°C (30 s), and finally elongation at 72°C (20 s) The reference gene for qPCR was β-actin. All the primer sequences used in this study are listed in Table 1. All the thymic tissue RNA samples were run in triplicate, and the expression levels of all desired genes were calculated by the ΔΔCt method [26].

Table 1: Primers used for real-time PCR.

Gene
Primer sequences (5 to 3)
Accession no.

β-Actin
f-TTGTTGACAATGGCTCCGGT
NM_205518.1

r-TCTGGGCTTCATCACCAACG

PCNA
f- TCTGAGGGCTTCGACACCTA
NM_204170.2

r- AACCTTTTCCTGATTTGGTGCTT

CDK1
f-TCTTCTGCCATTCAAGACGAGTTCTG
NM_205314.1

r-GATCCTAGCAGTACCTCTGGAGACC

CCNA2
f-GCCTTGCCTCATGGACCTTCAC
NM_205244.1

r-TCTGGTGCGTCAATAAGCGATACTG

CCNB2
f-TCCAGGTCCACTCAAGGTTCCAG
NM_001004369.1

r-CCACCAACTGAAGCCTCTTACGAG

2.6. Statistical Analysis

Statistical data were calculated as the deviation (SD). Statistical calculations and analyses as well as graphical representations were done with Prism software version 5.01 (GraphPad Software Inc., San Diego, USA). The statistical significance was accomplished with Bonferroni’s multiple comparison test after a one-way ANOVA test.

3. Results

3.1. The Changes in Thymus Weight, Index, and Morphology during Postnatal Development

The statistical analysis of chicken thymus weight and thymus index revealed that the thymus weight kept increasing with the increase of age and peaked at the 18th week and then decreased (Figure 1(a)), while the thymus index increased at the first week and then decreased (Figure 1(b)). HE staining showed that the chicken thymus was covered by a capsule which penetrated into the thymus and divided it into many niches. Each niche contained an outer cortex and an inner medulla. The thickness of the cortex decreased, and the clear corticomedullary junction became blurred with the advancement in age of birds from 0 to the 5th week and was extremely fuzzy at the 27th week. Furthermore, several round-shaped bubbles of dissolved adipose tissue appeared at the 5th week, and their number continuously increased with the postnatal growth of the birds (Figure 1(c)).

(a)

(b)

(c)

Figure 1: The changes of the thymus weight, index, and morphology during postnatal development. (a) The thymus weight kept increasing until the 18th week and then decreased. (b) The thymus index continued to decline. (c) Sections of thymuses were stained with hematoxylin and eosin (H&E) staining to analyze changes of the tissue structure. Light areas represent the medulla, and dark areas represent the cortex. The cortex was thick, and the boundary between the cortex and medulla was clear at the 0 week. With the increase in age, the cortex became thinner and the boundary blurred at the 5th week and was hardly distinguishable at the 27th week. In addition, many round-shaped bubbles were formed due to the dissolved adipose tissue, which appeared at the 5th week and kept increasing with the advancement of age. Scale bars = 200 μm. All the data are presented as . , .

3.2. Identification of DEGs in Chicken Thymus during Postnatal Development

The RNA-seq statistics pertaining to the average total reads, mapping reads, and unique matching for each of the 18 samples are shown in Figure 2(a). The differentially expressed genes (DEGs) were analyzed by using NOISeq software. Overall, there were 425 downregulated and 744 upregulated genes at all the five time points (0, 1, 5, 9, 18, and 27 weeks) (Supplementary File 1). Among these, we observed 114 DEGs (42 downregulated and 72 upregulated) at the 1st week, 151 DEGs (57 downregulated and 94 upregulated) at the 5th week, 183 DEGs (90 downregulated and 93 upregulated) at the 9th week, 328 DEGs (113 downregulated and 215 upregulated) at the 18th week, and 965 DEGs (345 downregulated and 620 upregulated) at the 27th week (Figure 2(b)). Maximum DEGs were observed at the 27th week; therefore, we selected the 27th week as a key time point to conduct further analysis. The scatter plots of all the expressed genes were drawn to represent the distribution of DEGs in screening threshold dimensions (Figure 2(c)). The top ten upregulated and downregulated DEGs at the 27th week are shown in Figure 2(d). Interestingly, the most of them are ncRNA (uncharacterized LOC Gallus gallus gene) and microRNA such as mir-7470, mir-1725, mir-29c, and mir-181a-1.

(a)

(b)

(c)

(d)

Figure 2: Analysis and validation of DEGs in the chicken thymus at different time points. (a) Alignment statistics of reads align to the reference genome. (b) The diagram shows total numbers of differentially expressed genes (DEGs) at 0 (1), 5, 9, 18, and 27 weeks. (c) Scatter plot diagram showing the log value of the gene expression of the 0-week thymus (X-axis) versus the log value of the gene expression of the 27th-week thymus (Y-axis). The blue color indicates downregulated genes, the orange color represents upregulated genes, and the brown color designates unchanged genes. Top legends on each figure show the statistics of screening threshold values. (d) Top 10 upregulated and downregulated DEGs at the 27th-week thymus. The positive numbers represented upregulated log values while negative numbers represent downregulated log values of DEGs.

3.3. Bioinformatics Analysis of DEGs

GO functional classification on DEGs for each pair was conducted following the previously described method [25]. All GO terms are grouped into three ontologies: biological process, cellular component, and molecular function. We found that most DEGs were enriched in GO terms at the 27th week, i.e., 1351 genes were annotated in biological processes, 845 genes in cellular components, and 313 genes in molecular function (Figure 3(a)). The top three enriched biological processes were the cellular process (190 genes), single-organism process (165 genes), and metabolic process (136 genes) (Figure 3(b)). The top three enriched cellular components were the cell (180 genes), cell part (180 genes), and organelle (133 genes). The top three enriched molecular function categories were the binding activity (174 genes), catalytic activity (64 genes), and transporter activity (15 genes) (Figure 3(b)). KEGG enrichment analysis of DEGs at the 27th week showed enrichment of 227 DEGs in cellular processes, 264 in the environmental information system, 79 in genetic information processing, 641 in human diseases, 224 in the metabolism, and 358 in organismal systems related to the signaling pathways (Figure 3(c)). Further analysis showed that DNA replication, cell cycle, cytokine-cytokine receptor, and mismatch repair pathways were significantly enriched (Figure 3(d)).

(a)

(b)

(c)

(d)

Figure 3: (a) Gene ontology (GO) analysis of DEGs at each time point showing enrichment of different commonly expressed genes in GO terms. (b) Diagram showing GO terms and their designated number of DEGs at the 27th-week thymus. The X-axis shows the number of DEGs while the Y-axis represents the GO terms. Three different colors designate three ontologies of GO terms: the blue color represents the biological process, the green color indicates cellular components, and the red color shows molecular functions. (c) Diagram showing the KEGG enrichment analysis of DEGs at the 27th-week thymus. (d) The X-axis shows the rich factor while the Y-axis represents different pathways. The colors represent the value. The size of the dot represents the DEG number in a certain pathway.

3.4. RNA-Seq Data Validation by Quantitative Real-Time PCR (qPCR)

The cell cycle signaling pathway is closely associated with cell proliferation and death; therefore, it may play an important role in postnatal development. PCNA, CCNA2, CCNB2, and CDK1 were chosen to conduct q-PCR. The expression patterns of these genes were similar to the results obtained from RNA-seq (Figure 4).

Figure 4: The expression patterns of four genes in chick thymuses were detected by qPCR. In order to test the reliability of the RNA-seq result, comparison analysis was performed for expression patterns of PCNA, CCNA2, CCNB2, and CDK1 genes.

4. Discussion

Gene expression can be measured by the application of the next generation sequencing (NGS) technologies such as RNA sequencing (RNA-seq) [27–29], since the changes in gene expression may play an important role in postnatal development of the thymus [30, 31]. For better understanding of the gene expression pattern, we conducted transcriptome analysis of the thymus at different developmental stages. A large body of scientific research has also described the basic histogenetical rule that the postnatal vertebrate thymus experiences distinct morphological and functional alterations [7, 32, 33]. Morphologically, the weight of the immune organ reflects its developmental speed [34], while the index represents its maturity [35]. Maximum DEGs were observed at the 27th week that dramatically coincided with alterations in the morphological structure of the thymus at the same time point as compared to the 0 week. In the current study, the increase in thymus weight before the 18th week might be the result of the increase in cell number and cell volume in the early postnatal developmental stage while the decline at the 27th week may be the result of the age-related atrophy. The shrinkage of the thymus with aging causes a decline in tissue mass and leads to certain structural changes [2, 7]. Interestingly, the most of the top ten upregulated and downregulated DEGs were ncRNA and microRNA at the 27th week. A lot of research experiments have showed that the developmental processes of organs are regulated by both protein-coding and noncoding regions of the genome [36, 37]. Many transcripts of both ncRNA and microRNA such as lncRNA205 and MiR-205 contribute to thymopoiesis and postnatal thymus development [2, 38]. Thymic involution is strongly mediated by microRNA along with adipose tissue and cytokine alterations with the advancement in age [39]. The study of normal thymic structure and its alteration offers a cornerstone for better understanding of its immune function [40]. Since the structure and the function are closely related, the microstructure of different aged thymuses was explored by HE staining that showed a downtrend, except for a small increase at the 3rd week, with significant alterations in ratio of the cortex to medulla. Similar changes in the aging thymus were also previously reported in rats [41–43], African ostrich chicks [44, 45], and cynomolgus macaques [46, 47]. In addition, we also found an accumulation of adipose tissue while aging. In a prior report, it was found that both adipocyte tissue and lipid-bearing cells increase within the thymic medulla with the advancement in age [48]. Thymic function during age-related atrophy is largely compromised of morphological alterations due to deposition of the adipose tissue within the immune organ together with the expansion of the perivascular spaces leading to greater risk of infections and autoimmune diseases [15, 49, 50]. Hence, all these studies proposed that the thymus of a developing chicken undergoes continued morphological alterations that must be driven by biomolecular and transcriptional dynamics in the immune organ microenvironment. Moreover, ncRNA and microRNA might be associated with age-associated thymic atrophy in the chicken at the 27th week of age at the transcriptional level.

Based upon the consistency in morphological and transcriptional findings, we performed GO and KEGG enrichment analyses using the most relevant DEGs at the 27th week of postnatal thymus development. In our study, DEGs annotated for GO analysis revealed the top three enriched molecular function categories as the binding activity, catalytic activity, and transporter activity. Several previous studies have reported that the catalytic activity plays an important role in thymic mass reduction [51–53] along with the transport and binding of fatty acids at adipose tissue deposition sites during age-related thymus atrophy [54, 55]. According to bioinformatics analysis, cell cycle, DNA replication, cytokine-cytokine receptor, and mismatch repair pathways were significantly enriched in the KEGG pathways in the chicken thymus at the 27th week of age. Almost similar transcripts participate in cell cycle regulation as implicated in DNA replication during transition of the G2-to-M or S phase [56, 57]. Cytokines mainly regulate immune response [58]. The balance between proinflammatory and anti-inflammatory cytokines is important for the body to maintain normal function [59, 60]. A previous study has shown that age-related thymic atrophy is associated with decreased expression of interleukin-7 (IL-7) [61]. Cyclin-A2 (CCNA2), cyclin-B2 (CCNB2), and cyclin-dependent kinase-1 (CDK1) genes showed similar expression patterns in RNA-seq and q-PCR results at the 27th week of age. CCNA2 is implicated in cell division and may exert antiproliferative effects and cell cycle arrest in the mammalian spleen during immunosuppression [62]. CDKs regulate the cell cycle functions and are involved in differentiation, transcription, and mRNA processing in eukaryotic cells [57, 63]. In agreement with our results, a similar RNA-seq study in murine found that several members of the cyclin family including CCNA2 and CCNB2 as well as CDKs are strongly associated with cell cycle regulation during the early stages of age-related thymic involution [64]. Thus, we concluded that the decrease in proliferation in the chicken thymus might be attributed to the arrest of the cell cycle, especially at the G1/S checkpoint by cyclin and cyclin-dependent kinases during age-associated thymic atrophy in chickens at the 27th week of age.

5. Conclusion

The current study revealed that during postnatal development, the thymus undergoes severe atrophy. The thymus structure was damaged and gene expression changed dramatically, especially at the 27th week of age. Moreover, we found significant activation of several signaling pathways such as the cytokine-cytokine receptor interaction and cell cycle signaling pathways. Hence, it may be inferred that those signaling pathways might be closely related to the postnatal chicken thymus development.

Data Availability

RNA-seq data supporting the current study has been deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under accession code PRJNA492316. The SRA data will be accessible with the following link after publication of the article: https://www.ncbi.nlm.nih.gov/sra/PRJNA492316, or before by contacting the corresponding author.

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the grants from the National Key Research and Development Program (2017YFE0113700), College Students’ Research Funds (2017048), Hubei Province Natural Science Foundation of China (2017CFB514), and National Natural Science Foundation of China (30800808).

Supplementary Materials

Venn diagrams showing downregulated (a) and upregulated (b) DGEs at different time points compared with the 0 week. H1, H2, and H3 denote three biological replicates. (Supplementary Materials)

References

I. Bódi, K. Minkó, D. Molnár, N. Nagy, and I. Oláh, “A novel aspect of the structure of the avian thymic medulla,” Cell and Tissue Research, vol. 359, no. 2, pp. 489–501, 2015.

A. R. Ansari and H. Liu, “Acute thymic involution and mechanisms for recovery,” Archivum Immunologiae et Therapiae Experimentalis, vol. 65, no. 5, pp. 401–420, 2017.

B. Bajoghli, P. Guo, N. Aghaallaei et al., “A thymus candidate in lampreys,” Nature, vol. 470, no. 7332, pp. 90–94, 2011.

J. Gameiro, P. Nagib, and L. Verinaud, “The thymus microenvironment in regulating thymocyte differentiation,” Cell Adhesion & Migration, vol. 4, no. 3, pp. 382–390, 2010.

J. F. A. P. Miller, “The discovery of thymus function and of thymus-derived lymphocytes,” Immunological Reviews, vol. 185, no. 1, pp. 7–14, 2002.

P. Samara, K. Ioannou, and O. E. Tsitsilonis, “Prothymosin alpha and immune responses: are we close to potential clinical applications?” Vitamins and Hormones, vol. 102, pp. 179–207, 2016.

D. P. Shanley, D. Aw, N. R. Manley, and D. B. Palmer, “An evolutionary perspective on the mechanisms of immunosenescence,” Trends in Immunology, vol. 30, no. 7, pp. 374–381, 2009.

V. W. C. Yu, B. Saez, C. Cook et al., “Specific bone cells produce DLL4 to generate thymus-seeding progenitors from bone marrow,” The Journal of Experimental Medicine, vol. 212, no. 5, pp. 759–774, 2015.

W. Seo and I. Taniuchi, “Transcriptional regulation of early T-cell development in the thymus,” European Journal of Immunology, vol. 46, no. 3, pp. 531–538, 2016.

T. Yamano, J. Nedjic, M. Hinterberger et al., “Thymic B cells are licensed to present self antigens for central T cell tolerance induction,” Immunity, vol. 42, no. 6, pp. 1048–1061, 2015.

C. O. Igbokwe and K. Ezenwaka, “Age-related morphological changes in the thymus of indigenous large white pig cross during foetal and postnatal development,” Anatomy, vol. 11, no. 1, pp. 12–20, 2017.

C. P. Martinez-Jimenez, N. Eling, H. C. Chen et al., “Aging increases cell-to-cell transcriptional variability upon immune stimulation,” Science, vol. 355, no. 6332, pp. 1433–1436, 2017.

V. D. Dixit, “Thymic fatness and approaches to enhance thymopoietic fitness in aging,” Current Opinion in Immunology, vol. 22, no. 4, pp. 521–528, 2010.

I. K. Chinn, C. C. Blackburn, N. R. Manley, and G. D. Sempowski, “Changes in primary lymphoid organs with aging,” Seminars in Immunology, vol. 24, no. 5, pp. 309–320, 2012.

D. B. Palmer, “The effect of age on thymic function,” Frontiers in Immunology, vol. 4, 2013.

S. Ki, D. Park, H. J. Selden et al., “Global transcriptional profiling reveals distinct functions of thymic stromal subsets and age-related changes during thymic involution,” Cell Reports, vol. 9, no. 1, pp. 402–415, 2014.

R. Rezzani, L. Nardo, G. Favero, M. Peroni, and L. F. Rodella, “Thymus and aging: morphological, radiological, and functional overview,” Age, vol. 36, no. 1, pp. 313–351, 2014.

A. B. Silva, D. Aw, and D. B. Palmer, “Evolutionary conservation of neuropeptide expression in the thymus of different species,” Immunology, vol. 118, no. 1, pp. 131–140, 2006.

L. Chen, S. Xiao, and N. R. Manley, “Foxn1 is required to maintain the postnatal thymic microenvironment in a dosage-sensitive manner,” Blood, vol. 113, no. 3, pp. 567–574, 2009.

H. Takaba, Y. Morishita, Y. Tomofuji et al., “Fezf2 orchestrates a thymic program of self-antigen expression for immune tolerance,” Cell, vol. 163, no. 4, pp. 975–987, 2015.

S.-F. Nie, L. F. Zha, Q. Fan et al., “Genetic regulation of the thymic stromal lymphopoietin (TSLP)/TSLP receptor (TSLPR) gene expression and influence of epistatic interactions between IL-33 and the TSLP/TSLPR axis on risk of coronary artery disease,” Frontiers in Immunology, vol. 9, 2018.

J. M. Sharma, J. E. Dohms, and A. L. Metz, “Comparative pathogenesis of serotype 1 and variant serotype 1 isolates of infectious bursal disease virus and their effect on humoral and cellular immune competence of specific-pathogen-free chickens,” Avian Diseases, vol. 33, no. 1, pp. 112–124, 1989.

H. Huang, A. Liu, H. Wu et al., “Transcriptome analysis indicated that Salmonella lipopolysaccharide-induced thymocyte death and thymic atrophy were related to TLR4-FOS/JUN pathway in chicks,” BMC Genomics, vol. 17, no. 1, p. 322, 2016.

S. Tarazona, F. García-Alcalde, J. Dopazo, A. Ferrer, and A. Conesa, “Differential expression in RNA-seq: a matter of depth,” Genome Research, vol. 21, no. 12, pp. 2213–2223, 2011.

A. R. Ansari, N.-Y. Li, Z.-J. Sun et al., “Lipopolysaccharide induces acute bursal atrophy in broiler chicks by activating TLR4-MAPK-NF-κB/AP-1 signaling,” Oncotarget, vol. 8, no. 65, pp. 108375–108391, 2017.

K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method,” Methods, vol. 25, no. 4, pp. 402–408, 2001.

Z. Wang, M. Gerstein, and M. Snyder, “RNA-seq: a revolutionary tool for transcriptomics,” Nature Reviews Genetics, vol. 10, no. 1, pp. 57–63, 2009.

M. L. Metzker, “Sequencing technologies—the next generation,” Nature Reviews Genetics, vol. 11, no. 1, pp. 31–46, 2010.

B. Li and C. N. Dewey, “RSEM: accurate transcript quantification from RNA-seq data with or without a reference genome,” BMC Bioinformatics, vol. 12, no. 1, p. 323, 2011.

J. Derbinski, A. Schulte, B. Kyewski, and L. Klein, “Promiscuous gene expression in medullary thymic epithelial cells mirrors the peripheral self,” Nature Immunology, vol. 2, no. 11, pp. 1032–1039, 2001.

M. Meredith, D. Zemmour, D. Mathis, and C. Benoist, “Aire controls gene expression in the thymic epithelium with ordered stochasticity,” Nature Immunology, vol. 16, no. 9, pp. 942–949, 2015.

B. Bodey, B. Bodey Jr., S. E. Siegel, and H. E. Kaiser, “Involution of the mammalian thymus, one of the leading regulators of aging,” In Vivo, vol. 11, no. 5, pp. 421–440, 1997.

N. R. Manley, E. R. Richie, C. C. Blackburn, B. G. Condie, and J. Sage, “Structure and function of the thymic microenvironment,” Frontiers in Bioscience, vol. 16, no. 1, pp. 2461–2477, 2011.

Z.-S. Bai, C.-Y. Che, Y.-X. Cheng, and Z.-H. Cai, “Analysis growth imbalance on immune organ of AA broiler,” Journal of Anhui Science and Technology University, vol. 2, 2017.

A. V. Khalyavkin and V. N. Krut’ko, “Early thymus involution – manifestation of an aging program or a program of development?” Biochemistry, vol. 80, no. 12, pp. 1622–1625, 2015.

T. Derrien, R. Johnson, G. Bussotti et al., “The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression,” Genome Research, vol. 22, no. 9, pp. 1775–1789, 2012.

A. Fatica and I. Bozzoni, “Long non-coding RNAs: new players in cell differentiation and development,” Nature Reviews Genetics, vol. 15, no. 1, pp. 7–21, 2014.

N. S. van Oers, A. R. Hoover, Q. Du et al., “A long noncoding RNA, lncRNA205, and an embedded microRNA, MiR-205 have overlapping and distinct contributions to thymopoiesis and development,” The Journal of Immunology, vol. 198, Supplement 1, 2017.

J. Dooley and A. Liston, “Molecular control over thymic involution: from cytokines and microRNA to aging and adipose tissue,” European Journal of Immunology, vol. 42, no. 5, pp. 1073–1079, 2012.

G. Pearse, “Normal structure, function and histology of the thymus,” Toxicologic Pathology, vol. 34, no. 5, pp. 504–514, 2006.

J. L. Tarry-Adkins, C. E. Aiken, T. J. Ashmore, D. S. Fernandez-Twinn, J.-H. Chen, and S. E. Ozanne, “A suboptimal maternal diet combined with accelerated postnatal growth results in an altered aging profile in the thymus of male rats,” The FASEB Journal, vol. 33, no. 1, pp. 239–253, 2019.

Y. Hamazaki, M. Sekai, and N. Minato, “Medullary thymic epithelial stem cells: role in thymic epithelial cell maintenance and thymic involution,” Immunological Reviews, vol. 271, no. 1, pp. 38–55, 2016.

S. Siddiqui, A. Lustig, A. Carter et al., “Genomic deletion of GIT2 induces a premature age-related thymic dysfunction and systemic immune system disruption,” Aging, vol. 9, no. 3, pp. 706–740, 2017.

J.-X. Wang, P. Li, X. T. Zhang, and L. X. Ye, “Distribution and morphology of ghrelin-immunopositive cells in the thymus of the African ostrich,” International Journal of Morphology, vol. 35, no. 2, pp. 541–546, 2017.

H.-B. Huang, K. Xiao, S. Lu et al., “Increased thymic cell turnover under boron stress may bypass TLR3/4 pathway in African ostrich,” PLoS One, vol. 10, no. 6, article e0129596, 2015.

P. W. Snyder, N. E. Everds, W. A. Craven, J. Werner, S. H. Tannehill-Gregg, and R. E. Guzman, “Maturity-related variability of the thymus in cynomolgus monkeys (Macaca fascicularis),” Toxicologic Pathology, vol. 44, no. 6, pp. 874–891, 2016.

M. S. Spoor, Z. A. Radi, and R. W. Dunstan, “Characterization of age- and gender-related changes in the spleen and thymus from control cynomolgus macaques used in toxicity studies,” Toxicologic Pathology, vol. 36, no. 5, pp. 695–704, 2008.

H. Yang, Y.-H. Youm, and V. D. Dixit, “Inhibition of thymic adipogenesis by caloric restriction is coupled with reduction in age-related thymic involution,” The Journal of Immunology, vol. 183, no. 5, pp. 3040–3052, 2009.

H. E. Lynch, G. L. Goldberg, A. Chidgey, M. R. M. van den Brink, R. Boyd, and G. D. Sempowski, “Thymic involution and immune reconstitution,” Trends in Immunology, vol. 30, no. 7, pp. 366–373, 2009.

E. Montecino-Rodriguez, B. Berent-Maoz, and K. Dorshkind, “Causes, consequences, and reversal of immune system aging,” The Journal of Clinical Investigation, vol. 123, no. 3, pp. 958–965, 2013.

K. Newton, D. L. Dugger, K. E. Wickliffe et al., “Activity of protein kinase RIPK3 determines whether cells die by necroptosis or apoptosis,” Science, vol. 343, no. 6177, pp. 1357–1360, 2014.

Y. Yin and W. Zhang, “The role of ghrelin in senescence: a mini-review,” Gerontology, vol. 62, no. 2, pp. 155–162, 2016.

A. V. Griffith, T. Venables, J. Shi et al., “Metabolic damage and premature thymus aging caused by stromal catalase deficiency,” Cell Reports, vol. 12, no. 7, pp. 1071–1079, 2015.

C. Pararasa, C. J. Bailey, and H. R. Griffiths, “Ageing, adipose tissue, fatty acids and inflammation,” Biogerontology, vol. 16, no. 2, pp. 235–248, 2015.

A. Bonen, M. L. Parolin, G. R. Steinberg et al., “Triacylglycerol accumulation in human obesity and type 2 diabetes is associated with increased rates of skeletal muscle fatty acid transport and increased sarcolemmal FAT/CD36,” The FASEB Journal, vol. 18, no. 10, pp. 1144–1146, 2004.

K. A. U. Gonzales, H. Liang, Y.-S. Lim et al., “Deterministic restriction on pluripotent state dissolution by cell-cycle pathways,” Cell, vol. 162, no. 3, pp. 564–579, 2015.

M. P. Swaffer, A. W. Jones, H. R. Flynn, A. P. Snijders, and P. Nurse, “CDK substrate phosphorylation and ordering the cell cycle,” Cell, vol. 167, no. 7, pp. 1750–1761.e16, 2016.

A. Iwasaki and R. Medzhitov, “Control of adaptive immunity by the innate immune system,” Nature Immunology, vol. 16, no. 4, pp. 343–353, 2015.

M. Gazdic, V. Volarevic, N. Arsenijevic, and M. Stojkovic, “Mesenchymal stem cells: a friend or foe in immune-mediated diseases,” Stem Cell Reviews and Reports, vol. 11, no. 2, pp. 280–287, 2015.

N. P. Singh, U. P. Singh, M. Rouse et al., “Dietary indoles suppress delayed-type hypersensitivity by inducing a switch from proinflammatory Th17 cells to anti-inflammatory regulatory T cells through regulation of microRNA,” The Journal of Immunology, vol. 196, no. 3, pp. 1108–1122, 2016.

D. Andrew and R. Aspinall, “Age-associated thymic atrophy is linked to a decline in IL-7 production,” Experimental Gerontology, vol. 37, no. 2-3, pp. 455–463, 2002.

K. A. Baken, J. L. A. Pennings, M. J. Jonker et al., “Overlapping gene expression profiles of model compounds provide opportunities for immunotoxicity screening,” Toxicology and Applied Pharmacology, vol. 226, no. 1, pp. 46–59, 2008.

A. R. Nebreda, “CDK activation by non-cyclin proteins,” Current Opinion in Cell Biology, vol. 18, no. 2, pp. 192–198, 2006.

H. Wu, X. Qin, H. Dai, and Y. Zhang, “Time-course transcriptome analysis of medullary thymic epithelial cells in the early phase of thymic involution,” Molecular Immunology, vol. 99, pp. 87–94, 2018.