Skip to main content

Temporal transcriptome profiling of floating apical out chicken enteroids suggest stability and reproducibility


Enteroids are miniature self-organising three-dimensional (3D) tissue cultures which replicate much of the complexity of the intestinal epithelium. We recently developed an apical-out leukocyte-containing chicken enteroid model providing a novel physiologically relevant in vitro tool to explore host–pathogen interactions in the avian gut. However, the replicate consistency and culture stability have not yet been fully explored at the transcript level. In addition, causes for the inability to passage apical-out enteroids were not determined. Here we report the transcriptional profiling of chicken embryonic intestinal villi and chicken enteroid cultures using bulk RNA-seq. Comparison of the transcriptomes of biological and technical replicate enteroid cultures confirmed their high level of reproducibility. Detailed analysis of cell subpopulation and function markers revealed that the mature enteroids differentiate from late embryonic intestinal villi to recapitulate many digestive, immune and gut-barrier functions present in the avian intestine. These transcriptomic results demonstrate that the chicken enteroid cultures are highly reproducible, and within the first week of culture they morphologically mature to appear similar to the in vivo intestine, therefore representing a physiologically-relevant in vitro model of the chicken intestine.


Organoids bridge the gap between traditional cell-lines and in vivo studies and facilitate research into multicellular biological mechanisms. The transcriptional analysis of enteroids has provided a greater understanding of the maturation stage, functions and responses of these in vitro culture models, as well as providing opportunities to study rare cell types [1]. Several RNA-seq datasets describing mammalian gel-embedded enteroids have already been published [2,3,4], but detailed analysis of the transcriptomes of avian enteroid cultures is lacking. Most RNA-seq studies of the chicken intestine have focused on functionality-related genes involved in digestion and absorption of energy and nutrients [5,6,7] and host–pathogen responses [8, 9] in mature birds. However, the late-embryonic to post-hatch gut undergoes immense changes in growth, metabolism and development. Analysis of the transcriptional changes in these critical life-stage processes could help identify novel methods to improve nutrition and gut-barrier function in the post-hatch chick [10,11,12].

Typically, gel-embedded enteroid crypts feed into a central functional lumen and are lined by a single layer of highly polarised epithelial cells. The apical brush borders of the epithelial cells face the lumen, with the basolateral surfaces in contact with the extracellular matrix scaffold [3]. This “basal-out” orientation means the luminal surface, which interacts with nutrients, toxins and microorganisms, is not directly accessible. Co et al. [13] demonstrated that by removing the extracellular matrix they could reverse the epithelial polarity of enteroids, allowing easy access to the apical surface. Their “apical-out” method involved an initial embedding step where the human and murine enteroids were established within Matrigel (a solubilized extracellular matrix secreted by Engelbreth-Holm-Swarm mouse sarcoma cells) for 7–20 days before moving them into a suspension. The first step appears necessary for mammalian enteroids to successfully retain their structural integrity in suspension. The mammalian apical-out enteroids were not described as surviving past 5 days in suspension culture and the ability to passage enteroids in this orientation does not appear possible. This technique has also been repeated in porcine enteroids where, similar to the human and mouse cultures, ~20% of enteroids did not reverse their polarisation [14]. This raises concerns with reproducibility of these suspension cultures and, alongside the loss of enteroid viability after a few days, restricts longer term use of this technique (reviewed in [15]).

We recently published a protocol to culture floating apical-out avian enteroids directly from embryonic intestinal villi tips and showed they reflect the cellular diversity and barrier function of the chicken intestinal epithelium [16]. The epithelium of the tips seals over in culture, enclosing an inner core of lamina propria, which likely supplies many of the stem cell niche factors. This inner core also contributes to the unique presence of functional intra-epithelial and lamina propria leukocytes. These avian enteroids develop more rapidly than their mammalian counterparts and appear to remain viable for up to 2 weeks in culture. We have also developed enteroids from mature birds however, as they have reduced longevity and budding compared to those developed from embryonic tissue, we have not explored these cultures in detail [16].

To demonstrate that apical-out chicken enteroids are a robust and reproducible cell culture model it is necessary to provide evidence of transcriptional homogeneity between biological and experimental replicates. Use of starting material (embryos) from the same breed and age and consistent culture conditions helps to reduce biological variation amongst enteroid cultures. However, it is uncertain whether inherent (epi) genetic and environmental variation among individual embryos renders individual biological replicates highly heterogeneous. Since these enteroid cultures could represent an important in vitro system for a diverse array of nutritional, pharmaceutical and host–pathogen studies, it is also necessary to have a detailed understanding of their cellular and transcriptional similarity to the avian intestine in vivo. Just like apical-out mammalian enteroid cultures, and despite many passage technique and supplementary growth factor trials, the floating chicken enteroids were unable to be indefinitely propagated. A detailed chronological study of the enteroid transcriptome could also identify the molecular mechanisms that limit their long-term culture.

In the current study, using floating apical-out avian enteroids derived from embryonic day 18 (ED18) chicken intestinal villi, we provide an in-depth analysis of their transcriptional responses over a 7 day culture period. The initial aim was to assess for homogeneity between replicates in order to demonstrate the highly reproducible consistency of this culture method. These data were also used to determine the period of stable gene expression in the cultures and identify the optimal time-window for experimental studies. We also aimed to describe the functional maturation of the enteroids compared to the in vivo intestine post-hatch at the transcript level. Our final aim was to identify putative transcriptional changes in floating apical-out enteroids that may underlie their inability to be indefinitely propagated. The results suggest that the chicken enteroid cultures are transcriptionally highly reproducible, and within the first week of culture the enteroids morphologically mature similar to the chick’s intestine. Therefore the apical-out chicken enteroids represent a physiologically-relevant in vitro model of the chicken intestine.

Materials and methods


Enteroid experiments were performed using ED18 Hy-Line Brown chickens (Gallus gallus; 3 days before hatching) obtained from the National Avian Research Facility, Edinburgh, UK. In addition, 8 week old Hy-Line Brown chickens were used for comparative studies. This breed and age of adult chicken were used in line with the 3Rs as they were concurrently in use by other groups. Ethical approvals were obtained from The Roslin Institute’s and University of Edinburgh’s Animal Welfare Ethics Review Board. All animals were housed in premises licensed under a UK Home Office Establishment License in full compliance with the requirements of the Animals (Scientific Procedures) Act 1986 and with approval from The Roslin Institute Animal Welfare Ethics Review Board.

Villus isolation and enteroid culturing

Isolation of avian ED18 intestinal villi containing stem cells and generation of floating chicken enteroids was performed as described previously [16]. In brief, the small intestine (from duodenum to ileocaecal junction) was removed post-mortem, cut into 5 mm sections and washed in Ca2+- and Mg2+-free Phosphate-Buffered Saline (PBS). The tissue was digested in Dulbecco’s Modified Eagle’s Medium (DMEM) (Thermo Fisher Scientific (TFS)) with 0.2 mg/mL Collagenase from Clostridium histolyticum Type IA (Merck) at 37 °C for 50 min. The tube was shaken vigorously, supernatant collected and these steps were repeated to generate 4 fractions. Fractions were centrifuged at 100 g for 4 min and crypts/villi were plated out in Advanced DMEM/F12 (TFS) supplemented with 10 mM HEPES (TFS), 2 mM L-Glutamine (TFS), 50 U/mL Penicillin/Streptomycin (Merck) and 2% B27 supplement (50X; TFS). Culture of avian enteroids occurred at 37 °C, 5% CO2.

RNA-seq experimental design

A total of 3 biological replicates (each culture composed of 3 × ED18 Hy-Line chickens) were used for the RNA-seq analysis. This was needed to allow enough ED18 enteroids to be available for sampling for each technical replicate at each time point. The isolated villi from each biological replicate were seeded (at 100 villi/well) in two identical 24-well tissue culture plates which were cultured in identical conditions. Each culture plate was treated as a technical replicate, therefore there were 2 technical replicates for each biological replicate. For sample collection, 8 wells per plate were collected at each of the defined time points (freshly isolated villi defined as 0 h, 1 day, 3 day, 4 day and 7 days of culture) and pooled to form one sample for transcriptomic analysis.

RNA sequencing and analysis

Enteroid total RNA was extracted and the quality and concentration assessed. Technical validation of the samples by tape-station confirmed they had good RNA integrity (RIN 6.1–9.6) and concentration (104–427 ng/µL) (Additional file 1). The libraries were sequenced, reads trimmed and aligned as previously described [16]. In brief, obtained reads were trimmed for quality and to remove adaptor sequences using Cutadapt. Reads after trimming were required to have a minimum length of 50 bases. Paired-end reads from Illumina sequencing were aligned to the Gallus gallus genome (Gallus_gallus-5.0) using STAR. The annotation used for counting was the standard GTF-format annotation for that reference (annotation v91). Raw counts for each annotated gene were obtained using the featureCounts software (v1.5.2).

The Principal Component Analysis (PCA) was constructed using normalised and filtered expression data from the unmerged technical replicates. P-values were obtained from an analysis of variance (ANOVA) test which assessed the associations between continuous value ranges in principal components and categorical experimental variables (group, pool and time). The differential gene expression analysis was performed as previously described [16], in brief differential gene expression analysis was performed within the Bioconductor edgeR package (v3.16.565). Comparison of the embryonic enteroid transcriptome at 0–7 days post cultivation revealed that there were no differentially expressed genes between the technical replicates (FDR < 0.05), demonstrating the consistency and reproducibility of the enteroid system. The sumTechReps function in EdgeR was used to merge technical replicates. The raw counts table was filtered to remove genes consisting predominantly of near-zero counts, filtering on counts per million (CPM) to avoid artefacts due to library depth. Statistical assessment of differential expression was carried out with the likelihood-ratio test. Differentially expressed genes were defined as those with FDR < 0.05 and logFC > 2. Heatmaps were constructed in R using the pheatmap package (v1.0.1066).

The annotated freshly isolated villi (0 h) and enteroid RNA seq normalised data sets were imported into the bioinformatics tool Graphia (Kajeka, Edinburgh, UK) [17,18,19] and a pairwise transcript-to-transcript Pearson’s correlation matrix was calculated as previously described [2]. The technical replicates were merged (average rather than sum due to the absence of the COB_2_168 hr sample) for this part of the analysis. A Pearson correlation coefficient cut-off threshold of r ≥ 0.95 was selected. The resulting graph was clustered using the Markov Clustering algorithm at an inflation value (which determines cluster granularity) of 2.2. The first 50 clusters (by number of genes) were manually grouped into similar gene expression profiles in relation to time points of culture, and each cluster location was identified on the network graph.

To enable functional annotation and interpretation of the clusters, the top 20 significantly overrepresented gene ontology (GO) within each expression profile group were identified using the Molecular Signatures Database [20]. To be accepted as significant, an over-represented GO terms needed a hypergeometric P < 0.05, and FDR < 0.05. Where the same genes within clusters represented several GO terms, only the most significant were selected for further consideration. Individual genes associated within the RNA-seq normalised dataset were also assessed for cellular functions and activities using a combination of literature review, Kyoto Encyclopaedia of Genes and Genomes (KEGG) Pathway database and GeneCards database. Heatmaps were constructed in R using the pheatmap package (v. 1.0.10) [21].


Enteroid total RNA was extracted and the quality and concentration assessed. First-strand cDNA synthesis was performed using a SuperScript IV reverse transcription kit (Invitrogen) according to the manufacturer’s instructions. Primer design for qPCR and qPCR was carried as described in Borowska et al. [22]. Primers are described in Additional file 2.

Transmission electron microscopy

Enteroids and villi were fixed in 3% glutaraldehyde in 0.1 M sodium cacodylate buffer, pH 7.3, for 2 h and processed as described previously [16]. Ultrathin sections(60 nm) were stained in uranyl acetate and lead citrate and imaged using a JEOL JEM-1400 Plus TEM and analysed using ImageJ (v1.52e, Fiji).

Whole mount and immunohistochemical staining

Details of the sources, clone numbers and concentrations of the primary and secondary antibodies used for immunohistochemistry are provided in Additional file 3. Enteroids were fixed, blocked, stained and mounted as previously described [16]. Samples were visualised using a LSM710 or LSM880 Confocal Microscope (Zeiss) and processed using ImageJ v1.52e.

Proliferation was assessed in a 5-ethynyl-2′-deoxyuridine (EdU) incorporation assay as previously described [16].

Statistics and reproducibility

The figure legends provide most details of sample sizes, replicates, number of repeats and statistics. Unless otherwise stated, experiment data is representative of at least 3 independent cultures each containing 2–3 embryos. Four pooled wells with ~100 enteroids/well were used for each histology sample. All measurements were recorded from distinct samples.

In the RNA-seq data each RNA sample was comprised of 8 pooled wells containing ~100 enteroids/well. An ANOVA test was performed to compare the associations between continuous value ranges in principal components and categorical experimental factors. Statistical assessment of the RNA-seq differential expression data was carried out using a quasi-likelihood F-test which accounts for uncertainty in estimating the dispersion for each gene. A Pearson’s correlation coefficient was used to measure the strength and direction of the linear relationship between genes/transcripts in the RNA-seq network analysis.

Data records

The mRNA expression datasets for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB37491 and PRJEB51227.


Transcriptomic analysis suggests enteroid cultures are highly reproducible

Enteroid cultures were prepared from isolated villi of an individual embryo as described in [16], with three independent cultures pooled to make one biological replicate, and each culture split into two technical replicates after villus isolation (Figure 1). Samples for transcriptome analysis were collected at 0 h (freshly isolated ED18 intestinal villi) then at 1, 3, 4 and 7 days of culture. The 1 day time point was selected as it represents a time of initial rapid and major morphological change in the enteroids, days 3 and 4 were chosen as the cultures appear more morphologically stable at this time, and 7 day enteroids were of interest as this time point is typically used for passage of basal-out enteroids.

Figure 1
figure 1

Overview of experimental design and RNA-seq data analysis pipeline. RNA was isolated from enteroids at 1, 3, 4 and 7 days and the freshly isolated ED18 intestinal villi (0 h). The validated RNA was sequenced on a Illumina NovaSeq 6000 system. All RNA-seq reads were preprocessed for a quality assessment. The filtered transcriptome reads were aligned to the genome and the expression profile was analysed.

To assess the reproducibility of the cultures, the degree of heterogeneity in the biological and technical replicate global transcriptional profiles was compared using a PCA. This showed that the biological and technical replicates tightly clustered by time point (Figure 2), confirming that the global transcriptomes of the enteroid culture replicates were highly consistent. The PC1 and 2 plot, where technical replicate per time point and time were the main drivers for the differences seen among the samples, showed clustering of the biological and technical replicates for each time point to independent locations on the plot, except for the 3 and 4 day replicates which occupied the same area, suggesting they shared a similar global transcriptional profile. However, defined clustering of ED18 intestinal villi (0 h) on the left side and the enteroid time points (1–7 day) on the right side of the PC1 axis is apparent which may represent the initial change from in vivo growth to in vitro culture.

Figure 2
figure 2

Principle Component Analysis (PCA) of RNA-seq expression data. Plot of the first and second principal components (PCs) from PCA using all samples. In PC1 and PC2, technical replicate per time point and time were the significant drivers for differences. Sample “COA_1_0hr” describes “Chicken Organoid, biological replicate A, technical replicate 1, 0 h time point” etc. The % values after the PCs indicate their contribution to the variability in the data.

Early enteroid cultures continue the rapid development as seen in the chick

A gene-to-gene correlation network graph was created next using the bioinformatics tool Graphia in order to explore the functional maturation of the enteroids during culture. A gene-to-gene correlation matrix was calculated for each transcript across each of the mRNA-seq data sets using a Pearson correlation threshold of r ≥ 0.95. The resulting network graph is presented in Figure 3 and contains groups of genes (clusters) with similar co-expression profiles across all data sets at r ≥ 0.95. The genes within the top 50 clusters (containing ≥ 17 genes) are presented in Additional file 4. Clusters with similar mean expression profiles typically occupied similar regions of the network graph. The clusters predominantly expressed at 0 h (Figure 3A) are associated with cell cycle and organ growth, whereas clusters predominantly expressed at higher levels in 1–7 day enteroids (Figure 3B) are associated with intestinal epithelial-cell specific functions. The clusters associated with cell energy production, cell metabolism are predominantly expressed at 1–4 days in culture, and clusters predominantly expressed at 3–7 days in culture cell are associated with secretion and lipid metabolism were found from 1 to 7 days in culture (Figures 3C–E). Clusters predominantly expressed in 7 day enteroids (Figure 3F) are mainly associated with cell differentiation.

Figure 3
figure 3

Network analysis of ED18 intestinal villi-enteroid RNA-seq expression data. Main component of the network graph derived from all five data sets samples. Here, the nodes represent transcripts (genes) and the edges represent correlations between individual expression profiles above r ≥ 0.95. Clusters that contain genes which are predominantly expressed at similar levels in all enteroids occupy different regions in the network graph. Coloured broken lines delineate these regions containing clusters of genes with similar mean expression profiles (AF). The graph x axis shows the samples ordered from time of cultures. The y axis shows the mean expression intensity (transcripts/million reads, TPM) for the cluster.

The mean gene expression profiles and representative human GO term enrichment annotations for the largest 50 clusters from the network graph are provided in Additional files 5 and 6. A summary of GO term enrichment annotations from select clusters was created for ease of reference which highlight processes and functions for cell cycle, tissue development and cell functions within the enteroid culture time points (Tables 1 and 2).

Table 1 GO term enrichment annotations for network clusters grouped into similar time-point expression profiles.
Table 2 GO term enrichment annotations for ED18 enteroid network clusters grouped into similar time-point expression profiles.

The 0 h time point (freshly isolated villi) denotes a stage of rapid development in the late embryonic intestine and accordingly GO terms associated with tissue development were predominant in the villus-related clusters including regulation of cell differentiation, organ growth and cell cycle (Table 1). Analysing cell cycle related cluster 1 in more detail, multiple individual genes were also represented in the KEGG Gallus gallus cell cycle pathway map e.g. cyclin dependent kinase 1 (CDK1), WEE1, cell division cycle 25A (CDC25A), eukaryotic translation termination factor 1 (ERF1) (Figure 4, Additional file 4). GO terms associated with cell cycle were also found in 1–4 day enteroid clusters (Table 2). In addition, many of the genes enriched in 0 h and 1 day enteroids were also represented in GO terms associated with transcription and translational regulation including translational initiation, peptide biosynthetic process and RNA processing (Table 1).

Figure 4
figure 4

Cell cycle related genes highlighted in the ED18 intestinal villi. KEGG pathway: Cell cycle Gallus gallus (chicken). Purple and red circled genes were contained in embryonic tissue analysis Cluster 1 whose mean gene expression profile was high in ED18 intestinal villi compared to enteroids. Associated GO terms in Cluster 1 were predominantly expressed in the cell cycle. Red circled genes have their mean expression profiles demonstrated in adjacent graphs. The graph x axis shows the samples ordered from time of cultures with 0 h representing isolated ED18 intestinal villi. The y axis shows the mean expression intensity for the cluster (transcripts/million reads, TPM).

Cluster 3, which was highly expressed in 1 day enteroids, contained genes with GO terms associated with cell energy production such as generation of precursor metabolites and energy and adenosine triphosphate (ATP) metabolic process (Table 2). This transcriptional analysis, confirms as anticipated, that cell proliferation is especially prevalent in the first couple of days of enteroid cultivation, a timeframe also reflective of the rapid development of the post-hatch chick intestine in vivo. Consistent with the transcriptional analysis, EdU + staining of enteroid cultures showed the percentage of proliferating cells dropped from almost 10% at 14 h to 5% after 7 days of culture (Additional file 7, [16]).

Enteroid cultures suggest prolonged expression of intestinal stem cell homeostasis pathways

Intestinal stem cell (ISC) maintenance is essential for enteroid cultures and is largely dependent on the Wingless-related integration site (WNT) signalling pathway. Immediately after isolation of the villi, the expression of some ISC-related genes was higher compared to 1 day enteroids most likely reflecting the relative abundance of ISC in ED18 villi (Figure 5A, Additional file 8). In addition, many ISC and WNT signalling pathway genes were steadily expressed in the enteroids including WNT4 and WNT5a, both of which are expressed by mesenchymal cells in mice (Figure 5B, Additional file 8, [23]). There was also a 7–17 fold downregulation of negative WNT regulators e.g. Zinc and ring finger 3 H1 (ZNRF3H1), ZNRF3H2, Ring finger protein 43 (RNF43) in the enteroids compared to 0 h (Figures 5B). However, several ISC markers such as OLFM4 and homeostatic factors such as WNT6, typically produced by Paneth cells in murine crypts [24], were significantly downregulated in the enteroids compared to 0 h (Figures 5A and B; Additional file 9). This reduced expression most likely reflects a dilution of the villi stem cell population as the enteroids are maturing in culture but could also suggest a reduction in activity. However, the expression of most of these ISC genes and related homeostatic factors subsequently remained stable throughout the later culture time points.

Figure 5
figure 5

Single gene analyses related to stem cell homeostasis and epithelial cell differentiation. Expression of A stem cell markers, B Wnt signalling pathway genes, and C intestinal epithelial cell differentiation related genes at 0 h and enteroids at 1, 3, 4 and 7 days of culture were analysed by RNA sequencing. Heat maps show the mean relative expression levels (log2 counts per million reads) of a range of mammalian gene-sets. RNA sequencing data is representative of 3 independent experiments, each comprising of 2 technical replicates, each containing 3 embryos.

Enteroid cultures reproduce the cellular diversity of the intestinal epithelium

Expression of gene sets characteristically associated with mammalian Paneth cells, goblet cells, enterocytes, enteroendocrine, transit amplifying, intestinal stem, and tuft cells were identified in the enteroids throughout the culture period. Furthermore, the transcriptomes for most epithelial cell-types remained relatively stable during this period. However, there was a reduction in gene expression associated with differentiation to enteroendocrine cells including neurogenin 3 (NEUROG3; 15-fold) and forkhead box A2 (FOXA2; two-fold) (Figure 5C). In addition, the expression of the epithelial differentiation marker alkaline phosphatase (ALPP) decreased over time. In summary, the transcriptome analysis suggests that the cellular diversity in the enteroid cultures is similar to the embryonic intestine.

Brush border develops in the enteroid cultures

In order to confirm if this in vitro intestinal model represented the typical array of in vivo functions, the transcriptome was also investigated for evidence of absorption and digestion at the gene level. The gene expression clusters that were predominantly expressed in the enteroids data sets compared to 0 h (freshly isolated villi) were significantly enriched in genes with GO terms related to cellular lipid metabolic process, chylomicron assembly, cellular protein catabolic process and glycerolipid metabolic process (Table 2). Looking at small intestinal digestion related genes in the enteroids in further detail, there was a relatively steady expression of many genes over the sample time points including brush border enzymes e.g. lactase (LCT) and solute carrier family SLC13A1, and those involved in fat metabolism e.g. apolipoprotein APOA1 and APOA4 (lipid transport), fatty acid binding protein FABP2 (fatty acid transport), and sterol carrier protein SCP2 (cholesterol uptake and transport) (Figure 6A) [25, 26].

Figure 6
figure 6

Single gene analyses of development and function of the enteroid cultures. Expression of A digestion related genes in ED18 intestinal villi (0 h) and enteroids at 1, 3, 4 and 7 days of culture were analysed by RNA sequencing. Heat maps show the mean relative expression levels (log2 counts per million reads) of a range of mammalian gene-sets. RNA sequencing data is representative of 3 independent experiments, each comprising of 2 technical replicates, each containing 3 embryos. B Quantification of microvilli length in ED18 intestine, 3 and 7 day enteroids and 8 week chicken intestine measured from TEM images; averaged for 3 cultures containing 2–3 embryos. *P < 0.05, unpaired Student’s t-test compared to ED18. C Representative TEM images demonstrating microvilli length in ED18 intestine, 2 and 7 day chicken enteroids, and 8 week chicken intestine. Scale bar: 4 µm.

The expression of genes related to the GO term “cell projection organisation” in the 7 day enteroids may reflect microvilli development, essential structures which increase the intestinal epithelial surface area for both secretion and absorption. To confirm whether the microvilli of apical-out enteroids continue to develop in vitro, their lengths across the culture period were measured in transmission electron microscopy images. This analysis revealed that the microvilli on the enterocyte surfaces doubled in the enteroids after 7 days of culture when compared to the ED18 intestine (Figures 6B and C). This morphological analysis alongside the transcriptome data suggests that the enteroids develop post-hatch digestive functions during cultivation.

Development of a comprehensive lamina propria component

ED18 villi isolation and floating culture conditions resulted in an apical-out enteroid phenotype with a lamina propria core [16]. In vivo the lamina propria is composed of multiple elements including immunologically competent cells, components of the enteric nervous system (glial cells and neuron axons) and mesenchymal cells. Consistent with this, GO terms related to immune system development and function were represented during the entirety of the culture period (Tables 1 and 2). Genes representative of several populations of lymphocytes, including T cells, B cells, macrophages, dendritic cells and NK cells were present at all time points (Additional file 10), consistent with the detection of CD45 + lymphocytes in the enteroids (Additional file 10). Genes associated with innate immune function were also relatively stable expressed over time, although some genes increased and some genes decreased over time in culture (Figure 7A; Additional file 9), e.g. pattern recognition receptors (Toll-like receptor 3 (TLR3) increased, whereas cytokine receptor gene IL20RB decreased, IL1B expression varied over time and antimicrobial enzyme avidin (AVD) increased.

Figure 7
figure 7

Cell types identified in the chicken enteroid lamina propria. Expression of A innate immune function related genes, B glial cell, E mesenchymal cell, J smooth muscle genes and K intestinal development gene sets in freshly isolated villi (0 h), 1, 3, 4 and 7 day chicken enteroids was compared by RNA sequencing analysis. Heat maps show the mean relative expression levels (log2 counts per million reads) of a range of mammalian gene-sets. RNA sequencing data is representative of 3 independent experiments, each comprising of 2 technical replicates, each containing 3 embryos. CD, FI Confocal images of chicken enteroids stained for lamina propria cell population markers (green) and counterstained with DAPI (blue) and Phalloidin (red). Enteroids at 2 days of culture stained for (C, D) HuC/HuD (green, glial cells), (F, G) vimentin (green, mesenchymal cells), (H, I) βIII-tubulin (green, mesenchymal cells). Scale bar: 20 µm. Images are representative of data from at least 3 independent cultures each containing 2–3 embryos.

Development of both the enteric nervous system and the muscularis mucosa was indicated in the 0 h transcriptomes with the expression of many genes with GO terms related to neuron projection, synapse and muscle cell proliferation, as would be expected for this late stage of embryonic maturation (Table 1). Further analysis of the genes within villi-associated cluster 20 confirmed the presence of many genes encoding proteins associated with synaptic signalling, neurotransmission, neuron growth as well as several genes encoding ligand- or voltage-gated ion channels (Additional File 1, cluster 20). Glial cell marker genes displayed steady expression in the enteroids (Figure 7B) and the presence of the enteric nervous system was further confirmed by staining for HuC/HuD, a pan-neuronal marker which labels neurons and glial cells. Multiple HuC/HuD + cells were identified scattered individually (Figure 7C) or in clusters (Figure 7D) throughout the central core of the 2 day enteroids. Since the presence of neuronal cell bodies would not be expected in the lamina propria, these are likely glial populations, reflecting characteristic cellular and tissue spatial distribution patterns [27].

Vimentin is the major cytoskeletal component of mesenchymal cells and, alongside other mesenchymal markers and highly expressed in the enteroid cultures (Figure 7E). Vimentin positive stained cells were accordingly restricted to the enteroid lamina propria core (Figures 7F and G). βIII-tubulin is generally used as a neuronal differentiation marker but it has also been shown to be expressed by other cell types including mesenchymal stem cells and perivascular cells [28, 29]. Since gut neuronal bodies are located in the submucosal ganglia [30], the scattered distribution of βIII-tubulin + cells in the enteroid core are most likely representative of mesenchymal stem cells with typical perinuclear staining (Figures 7H and I).

Smooth muscle fibres, blood and lymph vessels are typically present within a villus core (Figure 7J). GO terms related to blood vessel morphogenesis and epithelial or endothelial tube development were highly enriched on day 7 of culture (Table 2). A major contributor to angiogenesis, vasculogenesis and endothelial cell growth is vascular endothelial growth factor A (VEGFA), which is produced by many cell types including macrophages, and its transcript levels were steadily expressed over the enteroid cultures (Figure 7K). Analysis of genes specifically related to intestinal development showed many have steady transcript levels through the enteroid cultures (Figure 7K) [26]. This work provides a thorough transcriptomics validation of the enteroid core and provides a foundation for future studies to determine the full diversity of enteroid cell populations on a transcriptomic, proteomic and functional level.

No indication of epithelial-mesenchymal transition in 7 days enteroids

Although the floating chicken enteroid cultures develop profuse villus-like structures without external supplementation, they are not viable long-term. We previously showed these enteroid cultures display no general increase in cell stress-related genes indicating that the culture conditions are not causing excessive apoptosis or necrosis [16]. Therefore a possible cause of the reduced enteroid longevity could be due to a change in cell type and function. During epithelial-mesenchymal transition (EMT), epithelial cells acquire mesenchymal fibroblast-like properties, losing apical-basal polarity, weakening cell–cell junctions, rearranging the cytoskeleton and acquiring increased motility (reviewed in [31]). EMT will occur in the basal-out enteroids of some mammalian species unless there is specific inhibition of transforming growth factor β (TGFβ) and p38 mitogen activated protein kinase (MAPK) pathways to block this activity. We therefore compared the differentially expressed genes (FDR < 0.05, minimum logFC > 1) in 7 day versus 4 day enteroids to determine whether there was significant upregulation of genes associated with EMT and epithelial migration. This analysis showed a ~six-fold increase in fibroblast activation protein (FAP), hepatocyte growth factor (HGF) and integrin subunit A4 (ITGA4) (Additional file 9). Analysis of genes generally upregulated during EMT showed a steady expression across the enteroid time points and although some were significant differentially expressed the fold changes between time points were low (Figure 8A and Additional file 9). A similar steady expression was found of EMT-associated transcription factors (Figure 8A, Additional file 9), and evidence of ongoing epithelial development, including microvilli lengthening, and digestive function was still apparent at 7 days of culture (Figures 5C and 7A–C). This suggests that it is unlikely that the lack of long term viability and ability to passage the enteroids is mediated by EMT.

Figure 8
figure 8

EMT and growth factor single gene analysis of villi and enteroid cultures. Expression of A EMT and B growth factor associated genes in ED18 intestinal villi (0 h) and enteroids analysed by RNA sequencing. Heat maps show the mean relative expression levels (log2 counts per million reads) of a range of mammalian described gene-sets. RNA sequencing data is representative of 3 independent experiments, each comprising of 2 technical replicates, each containing 3 embryos.

To ascertain if deficiencies in supplemented growth factors and/or inhibited molecules, typically targeted in mammalian cultures, could be the cause of reduced culture longevity, the expression of genes encoding key factors was compared between 0 h (freshly isolated villi) and chicken enteroids (Figure 8B). Most genes had consistent to increasing expression across the cultures except for epidermal growth factor (EGF), WNT6 and R-spondin 2 (RSPO2) which were downregulated in the enteroids compared to 0 h. Of note, MAPK14 and MAPK11 which are targeted by p38MAPK-inhibitors, and TGFB1 and Rho associated coiled-coil containing protein kinase 1 (ROCK1) which are targeted by TGFβ-inhibitors were also steadily expressed. In summary, both EMT and growth factor deficiencies are unlikely contributors to the inability to passage apical-out enteroids, and further studies are required to ascertain both the cause and potential solutions to alleviate this longevity problem.


Transcriptome profiling analyses were performed on data from floating avian enteroid and ED18 intestinal villi (0 h) samples to demonstrate the reproducibility of gene expression within and between cultures, show the maturation of embryonic enteroids to post-hatch functions in culture, and explore the propagation issues facing apical-out enteroid cultures.

To support reproducibility of cell based experiments, a desirable trait for any tissue culture model is culture consistency. This is of particular importance in multicellular complex organoid models. The enteroid RNA-seq analysis demonstrated that the main cause for differences between sample transcriptomes was time of culture, with a high degree of similarity in global transcriptional activity between technical and biological replicates. As would be expected of biologically distinct samples with natural variation, there were slight differences in biological replicate transcriptomes identified within the PCA [32]. Individual variability in donor genotype, microbiome and disease status has been shown to affect the epithelial phenotype of enteroid cultures [33,34,35]. Compared to wild chicken breeds, commercial layer chicken lines have a fairly low genotypic diversity therefore the effect of variability of genotype on the Hy-Line chicken enteroid cultures should be small [36]. The use of late-embryonic villi to develop the enteroid cultures also reduces the chance that infection of the donor will influence the culture phenotype, although some avian viruses can be vertically transmitted [37]. The influence of microbiome on cultures will be limited since the avian intestine is mostly colonised post-hatch, enabling the use of enteroids for microbiome-host interaction studies. However, ensuring adequate numbers of biological replicates for enteroid experiments will improve the efficiency of statistical testing.

The divergence of ED18 intestinal villi and enteroid samples to opposite sides of the PC1 axis highlights a marked difference in their global gene expression profiles. This disparity may in part be due to adaptions to the in vitro culture environment, which do not reproduce all in vivo factors that affect intestinal development and maturation, including a developing microbiome. In addition, it may in part be due to the dramatic and rapid development of budding enteroids from embryonic villi. In the chicken during the first 24 h post-hatch there are striking changes in the mucosal morphology of the small intestine, including intensive cryptogenesis and a dramatic increase in villus length [10, 38]. Similar morphological changes are demonstrated by the growing chicken enteroids where there is rapid closure of the villi fragments within hours of plating out and development of multiple buds by 1 day which continue to lengthen over the following days [16]. In conjunction with the rapid growth, the 1 day enteroid transcriptome denotes upregulation of many genes involved in development, proliferation and differentiation. The reproducibility of in vitro models is important for a robust analysis, and the enteroids display steady transcriptional profiles for cell differentiation and most ISC homeostasis gene sets in the first 7 days of culture, with 3–4 days highlighted as a particularly stable timeframe [39]. However, although steady expression during the culture period was observed, there is a substantial decrease in gene expression of ISC genes when freshly isolated villi are compared with cultured enteroids. Whether this is due to the isolation step, dilution of the gene expression during proliferation or decrease of stem cells needs further investigation. Zhang and Wong [40] also demonstrated an increase in OLFM4 expression in post-hatch chicks over 0–7 days of age indicate this ISC marker may play a role in intestinal maturation.

As a precocial species, in parallel with the morphological changes, the ability of the chick gut to digest and absorb nutrients has been shown to increase quickly enabling chickens to have an almost fully mature intestine shortly after hatch [11, 38, 41]. For example, in vivo brush border enzyme activities have been shown to increase rapidly in the post-hatch chicken and the enteroids mirror this with an increase in expression of selected brush border enzyme genes [38]. Nutrient digestion and absorption in the enteroids was significantly centered in the lipids but not the other macro-nutrients. This mirrors the post-hatch chick intestine where there is a selective hierarchy at this age for fatty acid uptake over glucose or amino acids due to the yolk being an important source of nutrients for up to 72 h post-hatch [42, 43]. Although transcript differences do not necessarily equate with functional significance, that the chicken enteroids develop strong expression of digestion-related genes and associated pathways strongly indicates that the in vitro conditions contain the important cues for post-hatch maturation and correspond to post hatch transcriptomic analysis of the small intestine [44, 45].

Despite the capability existing to consecutively passage gel-embedded chicken enterospheres (organoid nomenclature for the poorly-differentiated spherical structures observed using this technique), currently there is no evidence for the ability to passage floating enteroid cultures [13, 14, 16, 46, 47]. Various trials have been performed in the floating chicken enteroids, but continued bud growth after passage has not been recognised [16]. Blockade of p38MAPK and TGF-β pathways in human, bovine and porcine gel-embedded enteroid cultures is required, in part, to help to protect the cells within them against EMT [4, 48,49,50]. However, trials of p38MAPK and TGF-β molecular inhibitors in floating chicken enteroid media showed no improvement in post-passage growth or budding [16]. Murine enteroid cultures require active supplementation with TGF-β in order to induce an EMT model indicating species-specific differences in these pathways which may also be present in the chicken [51].

Using transcriptome-wide comparisons of avian apical-out enteroid and ED18 intestinal villi samples we have constructed a detailed picture of gene expression patterns in developing chicken enteroids and the late-embryonic chicken intestine. This has demonstrated the consistency of this culture model as well as identifying hallmarks of multiple key intestinal functional attributes such as to development of digestive function and organ development. The enteroids have been demonstrated to have a transcriptionally stable timeframe where they are morphologically mature, indicating useful culture stages to conduct physiologically-relevant experiments. Although the enteroids do not have the ability to be propagated, these results indicate that this unique model is an adequate in vitro tool for post-hatch chicken intestinal studies into e.g. host–pathogen interactions, growth of previously unculturable micro-organisms, breed phenotype and nutrition studies, and feed additive and heat stress trials (reviewed in [52]). Introduction of the apical-out chicken enteroid and ED18 chicken intestinal villi transcriptome to the wider scientific community will provide a valuable resource for future avian intestinal research.

Availability of data and materials

The mRNA expression datasets for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB37491 and PRJEB51227. The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.


  1. Grun D, Lyubimona A, Kester L, Wiebrands K, Basak O, Sasaki N, Clevers H, Van Oudenaarden A (2015) Single-cell messenger RNA sequencing reveals rare intestinal cell types. Nature 525:251–255

    Article  PubMed  Google Scholar 

  2. Hamilton CA, Young R, Jayaraman S, Sehgal A, Paxton E, Thomson S, Katzer F, Hope J, Innes E, Morrison LJ, Mabbott NA (2018) Development of in vitro enteroids derived from bovine small intestinal crypts. Vet Res 49:54

    Article  PubMed  PubMed Central  Google Scholar 

  3. Sato T, Stange DE, Ferrante M, Vries RGJ, Van Es JH, Van Den Brink S, Van Houdt WJ, Pronk A, Van Gorp J, Siersema PD, Clevers H (2011) Long-term expansion of epithelial organoids from human colon, adenoma, adenocarcinoma, and Barrett’s epithelium. Gastroenterology 141:1762–1772

    Article  CAS  PubMed  Google Scholar 

  4. Sato T, de Vries RG, Snippert HJ, van de Wetering M, Barker N, Stange DE, Van Es JH, Abo A, Kujala P, Peters P, Clevers H (2009) Single Lgr5 stem cells build crypt–villus structures in vitro without a mesenchymal niche. Nature 459:262–265

    Article  CAS  PubMed  Google Scholar 

  5. Liu R, Liu J, Zhao G, Li W, Zheng M, Wang J, Li Q, Cui H, Wen J (2019) Relevance of the intestinal health-related pathways to broiler residual feed intake revealed by duodenal transcriptome profiling. Poult Sci 98:1102–1110

    Article  CAS  PubMed  Google Scholar 

  6. Sabino M, Cappelli K, Capomaccio S, Pascucci L, Biasato I, Verini-Supplizi A, Valiani A, Trabalza-Marinucci M (2018) Dietary supplementation with olive mill wastewaters induces modifications on chicken jejunum epithelial cell transcriptome and modulates jejunum morphology. BMC Genomics 19:576

    Article  PubMed  PubMed Central  Google Scholar 

  7. Xiao C, Deng J, Zeng L, Sun T, Yang Z, Yang X (2021) Transcriptome analysis identifies candidate genes and signaling pathways associated with feed efficiency in Xiayan chicken. Front Genet 12:607719

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Wang H, Li S, Qi K, Xue T, Tu J, Zhou X, Liu H (2016) Application of RNA-Seq technology for screening differentially expressed genes in chicken small intestine cells affected by avian pathogenic Escherichia coli. Chinese J Prev Vet Med 38:111–115

    Google Scholar 

  9. Truong AD, Hong YH, Lillehoj HS (2015) High-throughput sequencing reveals differing immune responses in the intestinal mucosa of two inbred lines afflicted with necrotic enteritis. Vet Immunol Immunopathol 166:116–124

    Article  CAS  PubMed  Google Scholar 

  10. Geyra A, Uni Z, Sklan D (2001) Enterocyte dynamics and mucosal development in the posthatch chick. Poult Sci 80:776–782

    Article  CAS  PubMed  Google Scholar 

  11. Uni Z, Tako E, Gal-Garber O, Sklan D (2003) Morphological, molecular, and functional changes in the chicken small intestine of the late-term embryo. Poult Sci 82:1747–1754

    Article  CAS  PubMed  Google Scholar 

  12. Uni Z, Geyra A, Ben-Hur H, Sklan D (2000) Small intestinal development in the young chick: crypt formation and enterocyte proliferation and migration. Br Poult Sci 41:544–551

    Article  CAS  PubMed  Google Scholar 

  13. Co JY, Margalef-Catala M, Li X, Mah AT, Kuo CJ, Monack DM, Amieva MR (2019) Controlling epithelial polarity: a human enteroid model for host-pathogen interactions. Cell Rep 26:2509–2520

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Li Y, Yang N, Chen J, Huang X, Zhang N, Yang S, Liu G, Liu G (2020) Next-generation porcine intestinal organoids: an apical-out organoid model for swine enteric virus infection and immune response investigations. J Virol 94:e01006-e1020

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Beaumont M, Blanc F, Cherbuy C, Egidy G, Giuffra E, Lacroix-Lamande S, Wiedemann A (2021) Intestinal organoids in farm animals. Vet Res 52:33

    Article  PubMed  PubMed Central  Google Scholar 

  16. Nash TJ, Morris KM, Mabbott NA, Vervelde L (2021) Inside-out chicken enteroids with leukocyte component as a model to study host-pathogen interactions. Comm Biol 4:377

    Article  CAS  Google Scholar 

  17. Freeman TC, Goldovsky L, Brosch M, Van Dongen S, Maziere P, Grocock RJ, Freilich S, Thornton J, Enright AJ (2007) Construction, visualisation, and clustering of transcription networks from microarray expression data (network analysis of transcriptomics data). PLoS Comp Biol 3:2032–2042

    Article  CAS  Google Scholar 

  18. Athanasios T, van Stijn D, Anton JE, Tom CF (2009) Network visualization and analysis of gene expression data using BioLayout Express3D. Nat Protocols 4:1535

    Article  Google Scholar 

  19. Mabbott NA, Baillie JK, Brown H, Freeman TC, Hume DA (2013) An expression atlas of human primary cells: inference of gene function from coexpression networks. BMC Genomics 14:632

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Aravind S, Pablo T, Vamsi KM, Sayan M, Benjamin LE, Michael AG, Amanada P, Scott LP, Todd RG, Eric SL, Jill PM (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 102:15545

    Article  Google Scholar 

  21. Kolde R. (2019) pheatmap: Pretty Heatmaps

  22. Borowska D, Kuo R, Bailey RA, Watson KA, Kaiser P, Vervelde L, Stevens MP (2019) Highly multiplexed quantitative PCR-based platform for evaluation of chicken immune responses. PLoS One 14:e0225658

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lickert H, Kispert A, Kutsch S, Kemler R (2001) Expression patterns of Wnt genes in mouse gut development. Mech Dev 105:181–184

    Article  CAS  PubMed  Google Scholar 

  24. Flanagan DJ, Austin CR, Vincan E, Phesse TJ (2018) Wnt signalling in gastrointestinal epithelial stem cells. Genes 9:178

    Article  PubMed  PubMed Central  Google Scholar 

  25. Workman MJ, Mahe MM, Trisno S, Poling HM, Watson CL, Sundaram N, Chang C-F, Schiesser J, Aubert P, Stanley EG, Elefanty AG, Miyaoka Y, Mandegar MA, Conklin BR, Neunlist M, Brugmann SA, Helmrath MA, Wells JM (2017) Engineered human pluripotent-stem-cell-derived intestinal tissues with a functional enteric nervous system. Nat Med 23:49–59

    Article  CAS  PubMed  Google Scholar 

  26. Finkbeiner SR, Hill DR, Altheim CH, Dedhia PH, Taylor MJ, Tsai Y-H, Chin AM, Mahe MM, Watson CL, Freeman JJ, Nattiv R, Thomson M, Klein OD, Shroyer NF, Helmrath MA, Teitelbaum DH, Dempsey PJ, Spence JR (2015) Transcriptome-wide analysis reveals hallmarks of human intestine development and maturation in vitro and in vivo. Stem Cell Rep 4:1140–1155

    Article  CAS  Google Scholar 

  27. Boesmans W, Hao MM, Fung C, Li Z, Van Den Haute C, Tack J, Pachnis V, Vanden Berghe P (2019) Structurally defined signaling in neuro-glia units in the enteric nervous system. Glia 67:1167–1178

    Article  PubMed  PubMed Central  Google Scholar 

  28. Lee MK, Tuttle JB, Rebhun LI, Cleveland DW, Frankfurter A (1990) The expression and posttranslational modification of a neuron-specific beta-tubulin isotype during chick embryogenesis. Cell Motil Cytoskeleton 17:118–132

    Article  CAS  PubMed  Google Scholar 

  29. Powell DW, Pinchuk IV, Saada JI, Chen X, Mifflin RC (2011) Mesenchymal cells of the intestinal lamina propria. Annu Rev Physiol 73:213–237

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Foudah D, Monfrini M, Donzelli E, Niada S, Brini AT, Orciani M, Tredici G, Miloso M (2014) Expression of neural markers by undifferentiated mesenchymal-like stem cells from different sources. J Immunol Res 2014:987678

    Article  PubMed  PubMed Central  Google Scholar 

  31. Kalluri R, Weinberg RA (2009) The basics of epithelial-mesenchymal transition. J Clin Invest 119:1420–1428

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Blainey P, Krzywinski M, Altman N (2014) Points of significance: replication. Nat Methods 11:879–880

    Article  CAS  PubMed  Google Scholar 

  33. Noben M, Verstockt B, De Bruyn M, Hendriks N, Van Assche G, Vermeire S, Verfaillie C, Ferrante M (2017) Epithelial organoid cultures from patients with ulcerative colitis and Crohn’s disease: a truly long-term model to study the molecular basis for inflammatory bowel disease? Gut 66:2193–2195

    Article  CAS  PubMed  Google Scholar 

  34. Abo H, Chassaing B, Harusato A, Quiros M, Brazil JC, Ngo VL, Viennois E, Merlin D, Gewirtz A, Nusrat A, Denning TL (2020) Erythroid differentiation regulator-1 induced by microbiota in early life drives intestinal stem cell proliferation and regeneration. Nat Commun 11:513

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Van der Hee B, Madsen O, Vervoort J, Smidt H, Wells JM (2020) Congruence of transcription programs in adult stem cell-derived jejunum organoids and original tissue during long-term culture. Front Cell Dev Biol 8:375

    Article  PubMed  PubMed Central  Google Scholar 

  36. Malomane DK, Simianer H, Weigend A, Reimer C, Schmitt AO, Weigend S (2019) The SYNBREED chicken diversity panel: a global resource to assess chicken diversity at high genomic resolution. BMC Genomics 20:345

    Article  PubMed  PubMed Central  Google Scholar 

  37. Guy JS (1998) Virus infections of the gastrointestinal tract of poultry. Poult Sci 77:1166–1175

    Article  CAS  PubMed  Google Scholar 

  38. Uni Z, Noy Y, Sklan D (1999) Posthatch development of small intestinal function in the poult. Poult Sci 78:215–222

    Article  CAS  PubMed  Google Scholar 

  39. Van de Flier LG, Haegebarth A, Stang DE, Van de Wetering M, Clevers H (2009) OLFM4 is a robust marker for stem cells in human intestine and marks a subset of colorectal cancer cells. Gastroenterology 137:15–17

    Article  PubMed  Google Scholar 

  40. Zhang H, Wong E (2018) Identification of cells expressing OLFM4 and LGR5 mRNA by in situ hybridization in the yolk sac and small intestine of embryonic and early post-hatch chicks. Poult Sci 97:628–633

    Article  CAS  PubMed  Google Scholar 

  41. Sklan D (2001) Development of the digestive tract of poultry. World’s Poultry Sci J 57:415–428

    Article  Google Scholar 

  42. Noy Y, Uni Z, Sklan D (1996) Routes of yolk utilisation in the newly-hatched chick. Br Poult Sci 37:987–996

    Article  CAS  PubMed  Google Scholar 

  43. Noble R, Moore J (1964) Studies on the lipid metabolism of the chick embryo. Can J Biochem 42:1729–1741

    Article  CAS  PubMed  Google Scholar 

  44. Gong H, Wang T, Chu Q, Wu M, Lang W, Lan H, Zhu L, Zhou Y, Wen Q, Zheng X (2022) Transcriptome profiling reveals morphogenesis-related candidate genes and pathways in the chick embryonic small intestine. Br Poult Sci 63:194–201

    Article  CAS  PubMed  Google Scholar 

  45. Schokker DJ, Hoekman AJW, Smits MA, Rebel JMJ (2009) Gene expression patterns associated with chicken jejunal development. Dev Comp Immunol 33:1156–1164

    Article  CAS  PubMed  Google Scholar 

  46. Oost MJ, Ijaz A, van Haarlem DA, Summeren KV, Velkers FC, Kraneveld AD, Venema K, Jansen CA, Pieters RHH, Ten Klooster JP (2022) Chicken-derived RSPO1 and WNT3 contribute to maintaining longevity of chicken intestinal organoid cultures. Sci Rep 12:10563

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Powell RH, Behnke MS (2017) WRN conditioned media is sufficient for in vitro propagation of intestinal organoids from large farm and small companion animals. Biol Open 6:698–705

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Xie BY, Wu AW (2016) Organoid culture of isolated cells from patient-derived tissues with colorectal cancer. Chin Med J 129:2469–2475

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Li H, Liu X, Chen F, Zuo K, Wu C, Yan Y, Chen W, Lin W, Xie Q (2018) Avian influenza virus subtype H9N2 affects intestinal microbiota, barrier structure injury, and inflammatory intestinal disease in the chicken ileum. Viruses 10:270

    Article  PubMed  PubMed Central  Google Scholar 

  50. Moustakas A, Heldin CH (2007) Signaling networks guiding epithelial-mesenchymal transitions during embryogenesis and cancer progression. Cancer Sci 98:1512–1520

    Article  CAS  PubMed  Google Scholar 

  51. Hahn S, Nam M-O, Noh JH, Lee DH, Han HW, Kim DH, Hahm KB, Hong SP, Yoo J-H, Yoo J (2017) Organoid-based epithelial to mesenchymal transition (OEMT) model: from an intestinal fibrosis perspective. Sci Rep 7:2435

    Article  PubMed  PubMed Central  Google Scholar 

  52. Nash T, Vervelde L (2022) Advances, challenges and future applications of avian intestinal in vitro models. Avian Pathol 51:317–329

    Article  CAS  PubMed  Google Scholar 

Download references


We thank Dr A. de Groof and Dr P. Vermeij (MSD Animal Health, The Netherlands) for their support and fruitful discussions. The RNA sequencing, quality filtering, PCA, ANOVA test and differential analysis were performed by Edinburgh Genomics (Edinburgh, UK).


This research was supported by an iCase studentship from the Biotechnology and Biological Sciences Research Council (BBSRC) in collaboration with MSD Animal Health (BB/MO14819), and Institute Strategic Program Grant funding from the BBSRC (BBS/E/D/10002071). The funding bodies had no role in the design of the study and collection, analysis of the data, and in writing the manuscript.

Author information

Authors and Affiliations



LV devised the project. TJN, NAM & LV contributed to the design and analysis of the study. TJN performed the experiments, collected the data and prepared the figures. KMM and TJN analysed the RNA-seq data with support from NAM, TJN and LV wrote the paper with input from NAM and KMM. LV and NAM obtained funding. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lonneke Vervelde.

Ethics declarations

Ethics approval and consent to participate

Ethical approvals were obtained from The Roslin Institute’s and University of Edinburgh’s Animal Welfare Ethics Review Board. All animals were housed in premises licensed under a UK Home Office Establishment License in full compliance with the requirements of the Animals (Scientific Procedures) Act 1986 and with approval from The Roslin Institute Animal Welfare Ethics Review Board.

Consent for publication

For the purpose of open access, the author has applied a CC-BY public copyright licence to any author accepted manuscript version arising from this submission.

Competing interests

PhD research support and decision to publish by MSD Animal Health. MSD-AH did not play a role in the study design, data collection and analysis. Pending patent; University of Edinburgh, Prof L. Vervelde and Dr T. Nash, PCT/GB2020/051607 (drawing priority from 1909655.1), application submitted July 2020 and titled “Avian Enteroids”.

Additional information

Handling editor: Vincent Béringue

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. RNA samples collected for the transcriptome sequencing.

Description of the ED18 intestinal villi and chicken enteroid RNA samples submitted for RNA sequencing.


Additional file 2. Primers and probes used in RT-qPCR.


Additional file 3. Primary antibodies used for immunohistochemistry.

Additional file 4. Genes found in the largest 50 co-expression cluster derived from the RNA-seq analysis data

. Table providing list of genes found in each co-expression cluster in the network graph derived from the RNA-seq analysis data.

Additional file 5. Mean expression profiles of the genes in each of the largest 50 co-expression clusters.

Individual mean expression profiles of the genes in each of the largest 50 co-expression clusters derived from the network graph. The x axis shows the samples ordered by time of cultures. The y axis shows the mean expression intensity (transcripts/million reads, TPM) for the cluster.

Additional file 6. GO term enrichment annotations for the genes in the largest 50 co-expression clusters.

Table listing the representative GO term enrichment annotations for the genes in the largest 50 co-expression clusters derived from the network graph. Note, the Entrez Gene IDs are for human orthologs.

Additional file 7. Proliferation in the enteroids.

Graph demonstrating the percentage of proliferating (EdU+) cells within enteroids over 7 days of culture, alongside z-stack image of proliferating (EdU+) cells within a 14 h enteroid.

Additional file 8. Cell-type gene expression in chicken enteroid cultures.

Expression of stem cell, Paneth cell, enteroendocrine cell, goblet cell and enterocyte genes in 1 day chicken enteroid cultures. Confocal images of chicken enteroids at 2 days of culture stained for.

Additional file 9. Differentially expressed gene summary.

Table listing the DEGs (logFC, P-value, FDR) from adjacent culture time points.

Additional file 10. Leucocytes identified in the chicken enteroid lamina propria.

Expression of immune cell population gene sets analysed by RNA sequencing analysis. Confocal images of chicken enteroids at 2 days of culture stained for CD45.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nash, T.J., Morris, K.M., Mabbott, N.A. et al. Temporal transcriptome profiling of floating apical out chicken enteroids suggest stability and reproducibility. Vet Res 54, 12 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Chicken
  • 3D organoid
  • apical-out
  • transcriptome
  • intestine
  • stem cell