- Research article
- Open Access
Deciphering alternative splicing and nonsense-mediated decay modulate expression in primary lymphoid tissues of birds infected with avian pathogenic E. coli (APEC)
BMC Genetics volume 18, Article number: 21 (2017)
Avian pathogenic E. coli (APEC) can lead to a loss in millions of dollars in poultry annually because of mortality and produce contamination. Studies have verified that many immune-related genes undergo changes in alternative splicing (AS), along with nonsense mediated decay (NMD), to regulate the immune system under different conditions. Therefore, the splicing profiles of primary lymphoid tissues with systemic APEC infection need to be comprehensively examined.
Gene expression in RNAseq data were obtained for three different immune tissues (bone marrow, thymus, and bursa) from three phenotype birds (non-challenged, resistant, and susceptible birds) at two time points. Alternative 5′ splice sites and exon skipping/inclusion were identified as the major alternative splicing events in avian primary immune organs under systemic APEC infection. In this study, we detected hundreds of differentially-expressed-transcript-containing genes (DETs) between different phenotype birds at 5 days post-infection (dpi). DETs, PSAP and STT3A, with NMD have important functions under systemic APEC infection. DETs, CDC45, CDK1, RAG2, POLR1B, PSAP, and DNASE1L3, from the same transcription start sites (TSS) indicate that cell death, cell cycle, cellular function, and maintenance were predominant in host under systemic APEC.
With the use of RNAseq technology and bioinformatics tools, this study provides a portrait of the AS event and NMD in primary lymphoid tissues, which play critical roles in host homeostasis under systemic APEC infection. According to this study, AS plays a pivotal regulatory role in the immune response in chicken under systemic APEC infection via either NMD or alternative TSSs. This study elucidates the regulatory role of AS for the immune complex under systemic APEC infection.
Avian pathogenic E. coli (APEC) can cause numerous diseases, such as colibacillosis, septicemia, pericarditis, and airsacculitis, leading to significant economic loss in the poultry industry because of mortality and produce contamination [1, 2]. Researchers have also indicated that APEC and human extraintestinal pathogenic E. coli (ExPEC) shared similar structures, indicating the zoonotic risk attributed to APEC strains [3–5]. Therefore, retail contaminated chicken can be considered an important reservoir for APEC to cause ExPEC infections in humans.
Studies have demonstrated that numerous immune-related genes are undergoing changes in alternative splicing (AS) to regulate the immune system under different conditions [6, 7], such as CD3, CTLA4, CD44, FYN, and VAV1 [6–10]. In addition, non-sense mediated decay (NMD), another post-transcriptional mechanism for gene expression regulation, mainly occurs under pathogen-induced stress . Kalyna et al.  reported on the combination of AS and NMD to regulate gene expression, suggesting that AS-NMD can control a number of splicing proteins. Lewis et al.  also revealed that AS-NMD is a widely used post-transcriptional regulatory strategy. Thus, studying the two major sources of transcriptome diversity - differential splicing and NMD - is critical to understanding the mechanisms of AS and its regulation among different conditions, uncovering structural and functional diversity.
Research on AS in different chicken tissues under various conditions is rarely conducted. Chacko and Ranganathan  indicated that about 23% of chicken genes undergo AS. Thus, the underlying post-transcriptional stage molecular mechanisms of hosts need to be elucidated to improve APEC infection response, as well as the control and prevention of the disease. In the present study, we determined the landscape of AS and NMD in the primary lymphoid tissues (bone marrow, thymus, and bursa) of birds with systemic APEC infection to offer insights into the regulatory elements that contribute to the pathobiology of APEC infection.
All libraries were sequenced by using Illumina® HiSeq 2000 with 100 bp single-end reads. RNAseq obtained approximately 21–28 million raw sequence read for each treatment for each of the three immune tissues except treatment non-challenged birds at 1 day post-infection (dpi) in bursa (Additional file 1: Figure S4). Trimmed reads that passed the quality filter were about 21–27 million (Additional file 1: Figure S4). On average, 80.52% (bone marrow), 80.93% (bursa), and 78.62% (thymus) were uniquely mapped to the chicken reference transcriptome (Additional file 1: Figure S4). Moreover, the spliced reads comprised an average of 26,16, 20.18, and 20.54% of the unique mapped reads for bone marrow, bursa, and thymus, respectively (Additional file 1: Figure S4). These results showed that more spliced reads were detected in the bone marrow than in the bursa and the thymus. Moreover, both resistant and susceptible birds had more spliced reads compared with non-challenged birds in the bone marrow and the thymus at both time post-infection (Additional file 1: Figure S4). Cufflinks was used to estimate transcript abundance with an expression-level threshold of FPKM ≥ 0.5. The sensitivity and specificity were about 90 and 45%, respectively, at the transcript level for each sample by using the Cuffcompare software.
Overall, the percent of alternative splicing (AS) events was similar in each of the primary lymphoid tissues of different phenotype birds. The most common AS events were alternative 5′ splice sites (A5SS), exon skipping/inclusion (ESI), intron retention (IR), and alternative 3′ splice sites (A3SS), with 20.24, 18.89, 14.03, and 10.59%, respectively, in the bone marrow (Fig. 1a). Similar phenomenon was also detected in the bursa and the thymus (Fig. 1b and c). Collectively, the A5SS and ESI were the major AS events for chicken primary immune organs under systemic APEC infection.
Nine 2-way contrasts were generated among the six experimental groups in each of the three tissues (Additional file 1: Figure S1-3). However, only the 5 dpi contrasts included large numbers of differentially expressed (DE) isoforms in genes in all the three tissues (Additional file 1: Figure S1-3): susceptible vs. non-challenged birds at 5 dpi and susceptible vs. resistant birds at 5 dpi.
The contrasting susceptible vs. non-challenged birds at 5 dpi had 174, 184, and 67 significant DE isoforms belonging to 145, 176, and 62 significant DE genes, generating 70, 56, and 22 novel isoforms as well as 9, 4, and 2 non-sense mediated decays (NMDs) in the bone marrow, bursa, and thymus, respectively (Additional file 1: Figure S1-3, Additional file 2: Table S1-3). For this contrast, 10, 2, and 12 DE genes from the same transcript start site (TSS) were found in the bone marrow, bursa, and thymus, respectively (Fig. 2a and Additional file 2: Table S4-6); 7, 6, and 2 DE genes from different TSSs were found in the bone marrow, bursa, and thymus (Fig. 3a and Additional file 2: Table S7-9).
When susceptible birds were contrasted with resistant birds at 5 dpi, 127, 24, and 610 significant DE isoforms belonging to 112, 23, and 577 significant DE genes, producing 46, 5, and 157 novel isoforms, as well as 2, 0, and 21 NMDs, were found in the bone marrow, bursa, and thymus, respectively (Additional file 1: Figure S3-5, Additional file 2: Table S10-12). In this contrast, 4 and 19 DE genes from the same TSS were found in the bone marrow and thymus, respectively (Fig. 2b and Additional file 2: Table S13-14). Meanwhile, 6, 1, and 13 DE genes came from different TSSs in the bone marrow, bursa, and thymus in the contrasting susceptible vs. resistant birds at 5 dpi (Fig. 3b and Additional file 2: Table S15-17).
The phagosome and lysosome pathways were the common pathways identified in susceptible birds that responded to systemic APEC infection in the bone marrow and thymus. Cell adhesion molecules (CAMs) were greatly induced in the bursa in contrasting susceptible vs. non-challenged birds at 5 dpi and in the thymus in the contrasting susceptible vs. resistant birds at 5 dpi (Table 1, Additional file 2: Table S18). The thymus in susceptible birds also expressed specific pathways to respond to systemic APEC infection: cytokine-cytokine receptor interaction, cell cycle, p53 signaling pathway, and apoptosis (Table 1, Additional file 2: Table S18). All of these detected significant pathways were consistent with the results obtained in previous single-tissue studies [15–17].
Two novel isoforms PASP (CUFF.22759.5) and PASP (CUFF.22759.2) involved in the lysosome pathway were detected in the present study. The novel isoform PASP (CUFF.22759.5) was also exhibited NMD. Validation of PSAP transcript CUFF.22759.5 was displayed in Fig. 2. Higher expression level of transcript CUFF.22759.5 was identified in challenged birds compared with that in non-challenged birds (Fig. 2b). In the cell cycle pathway, the isoform STT3A (ENSGALT00000001434) was found to exhibit NMD (Additional file 2: Table S18). Genes CDK1 and CDC45 participated in cell cycle and p53 signaling pathway were detected in several isoforms from the same TSS. The important isoforms (especially related to immune and growth response) from the same and different TSSs are presented in Figs. 3 and 4, such as RAG2, POLR1B, DNASE1L3, TXND5, and Ig lambda chain V-1 region.
This study provides a comprehensive description of the whole transcriptomic alternative splicing (AS) changes and non-sense mediated decay (NMD) in the primary lymphoid tissues of resistant and non-challenged birds compared with susceptible birds by using RNAseq data analysis. Thus, this study represents a powerful resource that enables further investigation into the molecular mechanisms of different host phenotypes in response to systemic APEC infection. AS events were generally ubiquitous in the transcriptome analysis of primary lymphoid tissues in chickens of different phenotypes. Exon skipping/inclusion (ESI) and alternative 5′ splice sites (A5SS) are predominant specific AS events in primary lymphoid tissues.
NMD is one of the conserved RNA-level post-transcriptional mechanisms for regulating gene expression pathway involved in development and stress response to exhibit RNA degradation and translational repression with premature termination [18, 19]. Both Rayson et al.  and Shi et al.  groups demonstrated the closer association of NMD with pathogen- or wounding-induced stresses. In the present study, two important genes (PSAP and STT3A) were involved in the significant pathways with NMD in the thymus of contrasting susceptible vs. resistant birds at 5 dpi. These two genes were significant DE genes and DE isoforms, referred to as differentially-expressed-transcript-containing genes (DETs).
In the current study on AS results, the PSAP gene gives rise to three different isoforms that have various structures involved in the lysosome pathway (Additional file 2: Table S18). PSAP is the precursor of lysosomal activator molecules as a cofactor to degradation of glycosphingolipids by lysosomal hydrolases [21–23]. Lysosomal PSAP breaks down proteins acting as biological catalysts and removes cells from internal and external waste. Mutation of the first half of the PASP C-terminus causes a problem in its transportation to lysosomes in the immune response . PSAP had a crucial role in preventing cell death or apoptosis and promote cell survival [25, 26].
PASP also showed various responses in biological processes , such as the differentiation of male reproductive organs , spermatogenesis , and fertilization . PSAP deficiency in human and mice can lead to death [30, 31]. Both Liu et al.  and Kramer et al.  reported the single nucleotide polymorphism (SNP) of PSAP was associated with spleen bacterial load under the Salmonella enteritidis challenge, indicating the PSAP is a positional candidate gene in response to pathogen-induced stress. In present study, two of the three PSAP isoforms were novel: CUFF.22759.2 (length = 4204 & FC = 7.56) and CUFF.22759.5 (length = 4143 & FC = 19.31). The isoform CUFF.22759.5 was found with NMD in the thymus of contrasting susceptible vs. resistant birds at 5 dpi (Additional file 2: Table S18).
Another DETs with NMD is STT3A involved in the significant cell cycle pathway in the thymus of contrasting susceptible vs. resistant birds at 5 dpi (Additional file 2: Table S18). This transcript ENSGALT00000001434 (STT3A) was significantly down-regulated (length = 1798 & FC = 3.36) in the thymus (Additional file 2: Table S18). STT3A plays an important role in DNA damage and replication checkpoints against spontaneous chromosome missegregation , which can protect the host against tunmorigenesis . Modified STT3A was normally used for antimitotic chemotherapeutic drugs to kill tumor cells [35, 36]. Koiwa et al.  demonstrated that the STT3A subunit isoform has an important function control adaptive response to stress via reduction cell division in plant. On the basis of current results, DET STT3A was more highly expressed in the thymus of resistant birds than that in susceptible birds. Resistant birds could use the important gene STT3A to resist systemic APEC infection via NMD modification.
In addition to NMD, TSS is a crucial step to control gene expression activity under different conditions . TSS was useful in understanding the transcription architecture and regulatory elements. Pal et al.  reported that the alternative TSS selection has a more important function in mRNA isoform diversity, compared with alternative splicing in some mammalian tissues because of TSS can turn on/off switches. Normally, the core promoters have a wide range of alternative TSS to start transcribing different isoforms [40–42].
Many interesting DETs with the same TSS were found in the thymus of contrasting susceptible vs. resistant birds at 5 dpi: CDC45, CDK1, RAG2, POLR1B, PSAP, and DNASE1L3 (Fig. 2b and Additional file 2: Table S14). CDC45 and CDK1 were involved in the significant pathways cell cycle and p53 signaling pathway in this study (Additional file 2: Table S18). Recombination activating gene 2 (RAG2) has an important function in chicken B cells undergoing immunoglobulin gene conversion . Polymerase (RNA) I polypeptide B (POLR1B) is a necessary enzyme for transcription. Deoxyribonuclease I-like 3 (DNASE1L3) is involved in immunoglobulin V gene diversification . All aformentioned DETs were examined for down-regulation in the susceptible birds compared with the resistant birds at 5 dpi (Fig. 2b), indicating that DETs from cell death, cell cycle, cellular function, and maintenance were predominant under systemic APEC infection. Moreover, it is might be a key mechanism for resistant birds to resist APEC disease.
In the contrast of susceptible vs. resistant birds at 5 dpi, DETs, TXNDC5 and the Ig lambda chain V-1 region (ENSGALG00000021139), were significantly up-regulated with different TSSs in the bone marrow (Fig. 3b and Additional file 2: Table S17). TXNDC5 is involved in cell proliferation and migration  and is a potential pharmacological target . In addition, the Ig lambda chain V-1 region contributes to the innate immune response . In different tissues, alternative use of TSS is regulated in a specific pattern [48–50]. Under different conditions, altering the TSS choice can regulate protein synthesis to exhibit various biological effects . Under systemic APEC infection, susceptible birds seem to have actively responded to the disease via properly expression of protein isoforms with different TSSs. In the present study, comprehensive detection of TSSs in different tissues can help elucidate the gene regulatory mechanisms and improve the prediction of gene regulation under systemic APEC infection.
RNAseq and bioinformatics tools were used in this study to provide a portrait of the alternative splicing (AS) event and non-sense mediated decay (NMD) in primary lymphoid tissues, which can significantly influence host homeostasis under different conditions. This study indicates that AS plays a pivotal regulatory role in the immune response of chicken under systemic APEC infection via either NMD or alternative transcription start sites (TSSs). Two differentially-expressed-transcript-containing genes (DETs), PSAP and STT3A, were particularly important in this study because of their immune and growth response function. Many interesting DETs from the same and different TSS were detected, such as CDC45, CDK1, RAG2, POLR1B, DNASE1L3, TXNDC5 and the Ig lambda chain V-1 region. Alternative TSSs seem to have an important function in pathogen-responsive expression patterns by regulating protein synthesis to exhibit various biological effects. The results of this study can be used to explore immune complexity and sheds light on the regulatory role of AS and NMD in systemic APEC infection, improving our understanding of the transcription architecture, genomic structure, and regulatory elements.
Experiment and RNAseq data
Briefly, one group of broilers was challenged with APEC O1 and the other group was non-challenged as control. Birds were euthanized with carbon dioxide gas followed by cervical dislocation at 1 day post-infection (dpi) and 5 dpi and bone marrow, thymus, and bursa were collected. At necropsy, a veterinary pathologist classified challenged birds as susceptible (total lesion score = 0–1) and resistant birds (total lesion score = 6–7) by using the standard pathology scoring system : pericardium (0–2), liver (0–2), and air sacs (0–3). In total, six treatment groups were generated: susceptible birds at 1 and 5 dpi, resistant birds at 1 and 5 dpi, and non-challenged birds 1 and 5 dpi.
RNAseq mapping reads and transcript assemble
Firstly, adapter for each read was removed via Fastx toolkit software v0.0.13. FastQC software v0.10.1 was used to evaluate reads quality for all samples with cutoff Phred score 32. NGS QC toolkit v2.2 was conducted to improve reads quality if the Phred score was lower than 32. After trimming adaptor and quality control, available reads of each sample were aligned to Ensemble gallus gallus 4.0 reference genome by using TopHat2 v2.0.9 with default parameters. TopHat2 is a popular and efficient tool to detect splice junction between exons [53, 54]. A local software of the Integrative Genomics Viewer v 2.2.17  was performed to check the mapping results. The output of TopHat2 for each sample was in bam format, which was processed directly in the follow-up assembler program cufflinks v2.1.1 with default parameters to generate a transcriptome assembly. Cufflinks can use a known reference annotation to detect novel transcripts by assembling the reads into transcripts . Fragments Per Kilobase of exon per Million fragments mapped (FPKM) was used to measured the transcript abundances through Cufflinks. After Cufflinks assemble, the assemblies of all samples for each tissue were merged together by using Cuffmerge v2.1.1 to provide a uniform basis for calculating gene and transcript expression. Moreover, the cuffcompare program was performed to test transcripts that are identical to the reference chicken genome (Ensembl gallus gallus 4.0 reference genome).
Differential analysis and further biological pathway analysis
The Cuffmerge file along with alignment bam files produced by TopHat2 was used to conduct differential analysis between two treatment groups via Cuffdiff2 v2.1.1 with default parameters. Many simple tabular output files were generated by Cuffdiff2: FPKM tracking files (transcript FPKMs file, gene FPKMs file, primary transcript FPKMs file, and coding sequence FPKMs file), differential expression tests files (primary transcript differential FPKM file, transcript differential FPKM file, gene differential FPKM file, and coding sequence differential FPKM file), differential splicing tests file, differential coding output file, and differential promoter use file. CummeRbund v2.10.0 , a powerful R based visualization package, was used to create volcano, scatter, and box plots for the Cuffdiff2 output. In R package environment, CummeRbund can reorganize the output files of Cuffdiff2 to create a SQLite database about the relationships between genes, transcripts, transcription start sites (TSS), and coding sequences (CDS) regions. The cutoff for significant DE genes was q value < 0.05 and fold change (FC) was more than 1.5. Significant differentially expressed (DE) genes in Cuffdiff2 were used to identify the significantly changed pathways (adjusted p value < 0.05) by using GOseq.
Alternative splicing detection
ASTALAVISTA v3.0  was used to identify and classify the percent of the types of alternative splicing (AS) for each tissue using the merged assembled transcriptome gtf file. spliceR v1.10.0, a powerful R package freely available from the Bioconductor repository, was used to classify the alternative splicing and predict the coding potential from RNAseq data . This package can detect single and multiple exon skipping/inclusion (ESI and MESI), intron retention/inclusion (ISI), alternative 5′/3′ splice sites (A5SS/A3SS), alternative transcription start site (ATSS), transcription termination site (ATTS), mutually exclusive exons (MEE), as well as the nonsense mediated decay (NMD) sensitivity of transcripts based on stop codon position. Also, this software can detect the different isoforms within a gene whether were from the same TSS or different TSS. Visualization and the generation of GTF files for genome browsers can be also obtained from spliceR. Moreover, the significant DE isoform was q value < 0.05 and fold change > 1.5.
Validation of alternative splicing by RT-qPCR
To validate the alternative splicing (AS) results, independent samples (challenged and non-challenged birds) were used to detect the relative expression level of transcript CUFF.22759.5 (PSAP) in thymus. After isolating the total RNA, RT-qPCR was used to amplify the splicing isoform by using one primer in the alternative exon and an opposing primer in a constitutive exon (Fig. 2a). Moreover, we used a boundary-spanning primer for the sequence encompassing the exon–exon junction with the opposing primer in a constitutive exon in this validated study.
Total RNA was extracted from different phenotype birds using TRIZOL Reagent according to manufacturer’s instructions (Life Technologies, CA, US). Isolated RNA was treated with RNase-free Dnase I to remove DNA contamination. A nanodrop spectrophotometer was used to measure the quantity and quality of purified RNA. The primers sequence of PSAP (CUFF.22759.5) was: F, 5- GGCTGTTGACCCTGCTT-3; R, 5-GATATGGCTGTGTTCAATAATGCA-3. QuantiTect SYBR Green RT-PCR kit (Qiagen Inc., Valencia, CA) was used to test transcript expression. 28S rRNA, a housekeeping gene, was served as internal control. Gene assay was run in triplicate for the same individual samples in RT-qPCR. The reaction products were resolved on 2% TAE/agarose gels. All of the candidate fragments were cloned for sequencing to check transcript CUFF.22759.5. The adjusted cycle threshold (Ct) values were calculated using the equation: 40 – [Ct sample gene mean + (Ct 28S median – Ct 28S mean) (slope of sample gene/slope of 28S)]. The qPCR data were analyzed using JMP 8.0.2 statistical software (SAS Institute Inc., Cary, NC).
Alternative 3′ splice sites
Alternative 5′ splice sites
Avian pathogenic Escherichia coli
Alternative transcription start site
Transcription termination site
Cell adhesion molecules
Extraintestinal pathogenic Escherichia coli
Fragments Per Kilobase of exon per Million fragments mapped
Gene expression omnibus
Mutually exclusive exons
Transcription start site
JanBen T, Schwarz C, Preikschat P, Voss M, Philipp HC, Wieler LH. Virulence-associated genes in avian pathogenic Escherichia coli (APEC) isolated from internal organs of poultry having died from colibacillosis. Int J Med Microbiol. 2001;291:371–8.
Stordeur P, Bree A, Mainil J, Moulin-Schouleur M. Pathogenicity of pap-negative avian Escherichia coli isolated from septicaemic lesions. Microbes Infect. 2004;6:637–45.
Johnson TJ, Kariyawasam S, Wannamuehler Y, Mangiamele P, Johnson SJ, Doetkott C, Skyberg JA, Lynne AM, Johnson JR, Nolan LK. The genome sequence of avian pathogenic Escherichia coli strain O1:K1:H7 shares strong similarities with human extraintestinal pathogenic E. coli genomes. J Bacteriol. 2007;189:3228–36.
Johnson TJ, Wannamuehler Y, Doetkott C, Johnson SJ, Rosenberger SC, Nolan LK. Identification of minimal predictors of avian pathogenic Escherichia coli virulence for use as a rapid diagnostic tool. J Clin Microbiol. 2008;46:3987–96.
Johnson TJ, Logue CM, Wannamuehler Y, Kariyawasam S, Doetkott C, DebRoy C, White DG, Nolan LK. Examination of the source and extended virulence genotypes of Escherichia coli contaminating retail poultry meat. Foodborne Pathog Dis. 2009;6:657–67.
Moulton VR, Tsokos GC. Alternative splicing factor/splicing factor 2 regulates the expression of the zeta subunit of the human T cell receptor-associated CD3 complex. J Biol Chem. 2010;285:12490–6.
Ip JY, Tong A, Pan Q, Topp JD, Blencowe BJ, Lynch KW. Global analysis of alternative splicing during T-cell activation. RNA. 2007;13:563–72.
Magistrelli G, Jeannin P, Herbault N, Benoit De Coignac A, Gauchat JF, Bonnefoy JY, Delneste Y. A soluble form of CTLA-4 generated by alternative splicing is expressed by nonstimulated human T cells. Eur J Immunol. 1999;29:3596–602.
Arch R, Wirth K, Hofmann M, Ponta H, Matzku S, Herrlich P, Zöller M. Participation in normal immune responses of a metastasis-inducing splice variant of CD44. Science. 1992;257:682–5.
Davidson D, Viallet J, Veillette A. Unique catalytic properties dictate the enhanced function of p59fynT, the hemopoietic cell-specific isoform of the Fyn tyrosine protein kinase, in T cells. Mol Cell Biol. 1994;14:4554–64.
Shi C, Baldwin IT, Wu J. Arabidopsis plants having defects in nonsense-mediated mRNA decay factors UPF1, UPF2, and UPF3 show photoperiod-dependent phenotypes in development and stress responses. J Integr Plant Biol. 2012;54:99e114.
Kalyna M, Simpson C, Syed NH, Lewandowska D, Marquez Y, Kusenda B, Marshall J, Dinh HQ, Barta A, Brown JW. Alternative splicing and non-mediated decay modulate expression of important regulatory genes in Arabidopsis. Nucleic Acids Res. 2012;40:2454–69. doi:10.1093/nar/gkr932.
Lewis BP, Green RE, Brenner SE. Evidence for the widespread coupling of alternative splicing and nonsense-mediated mRNA decay in humans. Proc Natl Acad Sci U S A. 2003;100:189–92.
Chacko E, Ranganathan S. Comprehensive splicing graph analysis of alternative splicing patterns in chicken, compared to human and mouse. BMC Genomics. 2009;10 Suppl 1:S5. doi:10.1186/1471-2164-10-S1-S5.
Sun H, Liu P, Nolan LK, Lamont SJ. Avian pathogenic Escherichia coli (APEC) infection alters bone marrow transcriptome in chickens. BMC Genomics. 2015;16:690. doi:10.1186/s12864-015-1850-4. PMID: 26369556.
Sun H, Liu P, Nolan LK, Lamont SJ. Novel pathways revealed in bursa of Fabricius transcriptome in response to extraintestinal pathogenic Escherichia coli (ExPEC) infection. PLoS ONE. 2015;10:e0142570. doi:10.1371/journal.pone.0142570.
Sun H, Liu P, Nolan LK, Lamont SJ. Thymus transcriptome reveals novel pathways in response to avian pathogenic Escherichia coli infection. Poult Sci. 2016;95(12):2803–14. doi:10.3382/ps/pew202.
Zhang Z, Xin D, Wang P, Zhou L, Hu L, Kong X, Hurst LD. Noisy splicing, more than expression regulation, explains why some exons are subject to nonsense-mediated mRNA decay. BMC Biol. 2009;7:23. doi:10.1186/1741-7007-7-23.
Kurihara Y, Matsui A, Hanada K, Kawashima M, Ishida J, Morosawa T, Tanaka M, Kaminuma E, Mochizuki Y, Matsushima A, Toyoda T, Shinozaki K, Seki M. Genome-wide suppression of aberrant mRNA-like noncoding RNAs by NMD in Arabidopsis. Proc Natl Acad Sci U S A. 2009;106:2453–8.
Rayson S, Arciga-Reyes L, Wootton L, De Torres Zabala M, Truman W, Graham N, Grant M, Davies B. A role for nonsense-mediated mRNA decay in plants: pathogen responses are induced in Arabidopsis thaliana NMD mutants. PLoS ONE. 2012;7:e31917.
Nakano T, Sandhoff K, Stqmper J, Christomanou H, Suzuki K. Structure of full-length cDNA coding for sulfatide activator, a Co-beta-glucosidase and two other homologous proteins: two alternate forms of the sulfatide activator. J Biochem. 1989;105:152–4.
O’Brien JS, Kretz KA, Dewji N, Wenger DA, Esch F, Fluharty AL. Coding of two sphingolipid activator proteins (SAP-1 and SAP-2) by same genetic locus. Science. 1988;241:1098–101.
Rorman EG, Grabowski GA. Molecular cloning of a human co-betaglucosidase cDNA: evidence that four sphingolipid hydrolase activator proteins are encoded by single genes in humans and rats. Genomics. 1989;5:486–92.
Carvelli L, Libin Y, Morales CR. Prosaposin: a protein with differential sorting and multiple functions. Histol Histopathol. 2015;30:647–60.
Misasi R, Sorice M, Di Marzio L, Campana WM, Molinari S, Cifone MG, Pavan A, Pontieri GM, O’Brien JS. Prosaposin treatment induces PC12 entry in the S phase of the cell cycle and prevents apoptosis: activation of ERKs and sphingosine kinase. FASEB J. 2001;15:467–74.
Misasi R, Garofalo T, Di Marzio L, Mattei V, Gizzi C, Hiraiwa M, Pavan A, Grazia Cifone M, Sorice M. Prosaposin: a new player in cell death prevention of U937 monocytic cells. Exp Cell Res. 2004;298:38–47.
Morales CR, Zhao Q, Lefrancois S, Ham D. Role of prosaposin in the male reproductive system: effect of prosaposin inactivation on the testis, epididymis, prostate, and seminal vesicles. Arch Androl. 2000;44:173–86.
Sylvester SR, Morales C, Oko R, Griswold MD. Sulfated glycoprotein-1 (saposin precursor) in the reproductive tract of the male rat. Biol Reprod. 1989;41:941–8.
Hammerstedt RH, Cramer PG, Barbato GF, Amann RP, O’Brien JS, Griswold MD. A fragment of prosaposin (SGP-1) from rooster sperm promotes sperm-egg binding and improves fertility in chickens. J Androl. 2001;22:361–75.
Elleder M, Jerabkova M, Befekadu A, Hreblcek M, Berna L, Ledvinova J, Hulkova H, Rosewich H, Schymik N, Paton BC, Harzer K. Prosaposin deficiency - a rarely diagnosed, rapidly progressing, neonatal neurovisceral lipid storage disease. Report of a further patient. Neuropediatrics. 2005;36:171–80.
Fujita N, Suzuki K, Vanier MT, Popko B, Maeda N, Klein A, Henseler M, Sandhoff K, Nakayasu H, Suzuki K. Targeted disruption of the mouse sphingolipid activator protein gene: A complex phenotype, including severe leukodystrophy and wide-spread storage of multiple sphingolipids. Hum Mol Genet. 1996;5:711–25.
Liu W, Lamont SJ. Candidate gene approach: potentional association of caspase-1, inhibitor of apoptosis protein-1, and prosaposin gene polymorphisms with response to Salmonella enteritidis challenge or vaccination in young chicks. Anim Biotechnol. 2003;14:61–76.
Kramer J, Malek M, Lamont SJ. Association of twelve candidate gene polymorphisms and response to challenge with Salmonella enteritidis in poultry. Anim Genet. 2003;34:339–48.
Zachos G, Rainey MD, Gillespie DA. Chk1-deficient tumor cells are viable but exhibit multiple checkpoint and survival defects. EMBO J. 2003;22:713–23.
Kops GJ, Weaver BA, Cleveland DW. On the road to cancer: aneuploidy and the mitotic checkpoint. Nat Rev Cancer. 2005;5:773–85.
Zachos G, Black EJ, Walker M, Scott MT, Vagnarelli P, Earnshaw WC, Gillespie DA. Chk1 is required for spindle checkpoint function. Dev Cell. 2007;12:247–60.
Koiwa H, Li F, McCully MG, Mendoza I, Koizumi N, Manabe Y, Nakagawa Y, Zhu J, Rus A, et al. The STT3A subunit isoform of the Arabidopsis oligosaccharyltransferase controls adaptive responses to salt/osmotic stress. Plant Cell. 2003;15:2273–84.
Hochheimer A, Tjian R. Diversified transcription initiation complexes expand promoter selectivity and tissue-specific gene expression. Genes Dev. 2003;17:1309–20.
Pal S, Gupta R, Kim H, Wickramasinghe P, Baubet V, Showe LC, Dahmane N, Davuluri RV. Alternative transcription exceeds alternative splicing in generating the transcriptome diversity of cerebellar development. Genome Res. 2011;21:1260–72.
Suzuki Y, Taira H, Tsunoda T, Mizushima-Sugano J, Sese J, Hata H, Ota T, Isogai T, Tanaka T, Morishita S, et al. Diverse transcriptional initiation revealed by fine, large-scale mapping of mRNA start sites. EMBO Rep. 2001;2:388–93.
Carninci P, Sandelin A, Lenhard B, Katayama S, Shimokawa K, Ponjavic J, Semple CA, Taylor MS, Engstrom PG, Frith MC, et al. Genome-wide analysis of mammalian promoter architecture and evolution. Nat Genet. 2006;38:626–35.
Smale ST, Kadonaga JT. The RNA polymerase II core promoter. Annu Rev Biochem. 2003;72:449–79.
Carlson LM, Oettinger MA, Schatz DG, Masteller EL, Hurley EA, McCormack WT, Baltimore D, Thompson CB. Selective expression of RAG-2 in chicken B cells undergoing immunoglobulin gene conversion. Cell. 1991;64:201–8.
Okamoto N, Okamoto M, Araki S, Arakawa H, Mizuta R, Kitamura D. Possible contribution of DNase gamma to immunoglobulin V gene diversification. Immunol Lett. 2009;125:22–30. doi:10.1016/j.imlet.2009.05.005.
Chang X, Xu B, Wang L, Wang Y, Wang Y, Yan S. Investigating a pathogenic role for TXNDC5 in tumors. Int J Oncol. 2013;43:1871–84. doi:10.3892/ijo.2013.2123.
Horna-Terrón E, Pradilla-Dieste A, Sánchez-de-Diego C, Osada J. TXNDC5, a newly discovered disulfide isomerase with a key role in cell physiology and pathology. Int J Mol Sci. 2014;15:23501–18.
Reynaud CA, Anquez V, Dahan A, Weill JC. A single rearrangement event generates most of the chicken immunoglobulin light chain diversity. Cell. 1985;40:283–91.
Kimura K, Wakamatsu A, Suzuki Y, Ota T, Nishikawa T, Yamashita R, Yamamoto J, Sekine M, Tsuritani K, Wakaguri H, et al. Diversification of transcriptional modulation: Large-scale identification and characterization of putative alternative promoters of human genes. Genome Res. 2006;16:55–65.
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456:470–6.
Yamashita R, Sathira NP, Kanai A, Tanimoto K, Arauchi T, Tanaka Y, Hashimoto S, Sugano S, Nakai K, Suzuki Y. Genome-wide characterization of transcriptional start sites in humans by integrative transcriptome analysis. Genome Res. 2011;21:775–89.
Rojas-Duran MF, Gilbert WV. Alternative transcription start site selection leads to large differences in translation activity in yeast. RNA. 2012;18:2299–305.
Peighambari SM, Julian RJ, Gyles CL. Experimental Escherichia coli respiratory infection in broilers. Avian Dis. 2000;44:759–69.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.
Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2012;14:178–92.
Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011;27:2325–9.
Goff L, Trapnell C. CummeRbund: analysis, exploration, manipulation, and visualization of Cufflinks high-throughput sequencing data. R package version 1.2.0. 2011.
Foissac S, Sammeth M. ASTALAVISTA: dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 2007;35 suppl 2:W297–9.
Vitting-Seerup K, Porse BT, Sandelin A, Waage J. spliceR: an R package for classification of alternative splicing and prediction of coding potential from RNAseq data. BMC Bioinf. 2014;15:81–8.
This work was supported by Yangzhou University Research Funding no. 137010990.
Availability of data and materials
RNAseq data of bone marrow, thymus, and bursa are available in Gene Expression Omnibus (GEO) with accession numbers GSE67302, GSE6901, and GSE70334, respectively.
HS designed the experiments, analyzed the data, and wrote the manuscript.
The author declares that there are no competing interests.
Consent for publication
For the RNAseq experiment, all animal care and experimental procedures were reviewed and approved by the Institutional Animal Care and Use Committee of Iowa State University (Log #11-07-6460-G). For independent validated experiment, all animal protocols were approved by the animal welfare committee of the Yangzhou University (permission number: SYXK(Su) IACUC 2012–0036).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. The significant differentially expressed (DE) isoforms and DE genes in the nine contrast of bone marrow. S5, susceptible birds at 5 days post-infection; S1, susceptible birds at 1 day post-infection; R5, resistant birds at 5 days post-infection; R1, resistant birds at 1 day post-infection; NC5, non-challenged birds at 5 days post-infection; NC1, non-challenged birds at 1 day post-infection. Pre-mature stop codon (PTC) is equal to non-sense mediated decay (NMD). Figure S2. The significant differentially expressed (DE) isoforms and DE genes in the nine contrast of bursa. S5, susceptible birds at 5 days post-infection; S1, susceptible birds at 1 day post-infection; R5, resistant birds at 5 days post-infection; R1, resistant birds at 1 day post-infection; NC5, non-challenged birds at 5 days post-infection; NC1, non-challenged birds at 1 day post-infection. Pre-mature stop codon (PTC) is equal to non-sense mediated decay (NMD). Figure S3. The significant differentially expressed (DE) isoforms and DE genes in the nine contrast of thymus. S5, susceptible birds at 5 days post-infection; S1, susceptible birds at 1 day post-infection; R5, resistant birds at 5 days post-infection; R1, resistant birds at 1 day post-infection; NC5, non-challenged birds at 5 days post-infection; NC1, non-challenged birds at 1 day post-infection. Pre-mature stop codon (PTC) is equal to non-sense mediated decay (NMD). Figure S4. Raw, clean, unique mapped, and splice reads distribution in each treatment in each of the three immune tissues. NC1, non-challenged birds at 1 day post-infection; NC5, non-challenged birds at 5 days post-infection; R1, resistant birds at 1 day post-infection; R5, resistant birds at 5 days post-infection; S1, susceptible birds at 1 day post-infection; S5, susceptible birds at 5 days post-infection (DOCX 2498 kb)
Table S1. Significant DEI and DEG in susceptible vs. non-challenged birds at 5 dpi in bone marrow. Table S2. Significant DEI and DEG in susceptible vs. non-challenged birds at 5 dpi in bursa. Table S3. Significant DEI and DEG in susceptible vs. non-challenged birds at 5 dpi in thymus. Table S4. DEG from the same TSS in susceptible vs. non-challenged birds at 5 dpi in bone marrow. Table S5. DEG from the same TSS in susceptible vs. non-challenged birds at 5 dpi in bursa. Table S6. DEG from the same TSS in susceptible vs. non-challenged birds at 5 dpi in thymus. Table S7. DEG from different TSS in susceptible vs. non-challenged birds at 5 dpi in bone marrow. Table S8. DEG from different TSS in susceptible vs. non-challenged birds at 5 dpi in bursa. Table S9. DEG from different TSS in susceptible vs. non-challenged birds at 5 dpi in thymus. Table S10. Significant DEI and DEG in susceptible vs. resistant birds at 5 dpi in bone marrow. Table S11. Significant DEI and DEG in susceptible vs. resistant birds at 5 dpi in bursa. Table S12. Significant DEI and DEG in susceptible vs. resistant birds at 5 dpi in thymus. Table S13. DEG from the same TSS in susceptible vs. resistant birds at 5 dpi in bone marrow. Table S14. DEG from the same TSS in susceptible vs. resistant birds at 5 dpi in thymus. Table S15. DEG from different TSS in susceptible vs. resistant birds at 5 dpi in bone marrow. Table S16. DEG from different TSS in susceptible vs. resistant birds at 5 dpi in bursa. Table S17. DEG from different TSS in susceptible vs. resistant birds at 5 dpi in thymus. Table S18. Significant DEI in the significantly changed pathways in different contrasts in different tissues. (XLSX 904 kb)
About this article
Cite this article
Sun, H. Deciphering alternative splicing and nonsense-mediated decay modulate expression in primary lymphoid tissues of birds infected with avian pathogenic E. coli (APEC). BMC Genet 18, 21 (2017) doi:10.1186/s12863-017-0488-4
- Alternative splicing
- Nonsense-mediated decay
- Bone marrow