A novel method for genome-wide profiling of dynamic host-pathogen interactions using 3' end enriched RNA-seq

Profiling of alternative polyadenylation and gene expression in PEDV-infected IPEC-J2 cells

Since 2010, porcine epidemic diarrhea virus (PEDV) has received global attention with the emergence of variant strains characterized with high pathogenicity. The pathogen-host interaction after PEDV infection is still unclear. To investigate this issue, high-throughput-based sequencing technology is one of the optimal choices. In this study, we used in vitro transcription sequencing alternative polyadenylation sites (IVT-SAPAS) method, which allowed accurate profiling of gene expression and alternative polyadenylation (APA) sites to profile APA switching genes and differentially expressed genes (DEGs) in IPEC-J2 cells during PEDV variant strain infection. We found 804 APA switching genes, including switching in tandem 3' UTRs and switching between coding region and 3' UTR, and 1,677 DEGs in host after PEDV challenge. These genes participated in variety of biological processes such as cellular process, metabolism and immunity reactions. Moreover, 413 genes, most of which are the "focus" genes in interaction networks, were found to be involved in both APA switching genes and DEGs, suggesting these genes were synchronously regulated by different mechanisms. In summary, our results gave a relatively comprehensive insight into dynamic host-pathogen interactions in the regulation of host gene transcripts during PEDV infection.

Association of Candidate Genes with Response to Heat and Newcastle Disease Virus

Newcastle disease is considered the number one disease constraint to poultry production in low and middle-income countries, however poultry that is raised in resource-poor areas often experience multiple environmental challenges. Heat stress has a negative impact on production, and immune response to pathogens can be negatively modulated by heat stress. Candidate genes and regions chosen for this study were based on previously reported associations with response to immune stimulants, pathogens, or heat, including: TLR3, TLR7, MX, MHC-B (major histocompatibility complex, gene complex), IFI27L2, SLC5A1, HSPB1, HSPA2, HSPA8, IFRD1, IL18R1, IL1R1, AP2A2, and TOLLIP. Chickens of a commercial egg-laying line were infected with a lentogenic strain of NDV (Newcastle disease virus); half the birds were maintained at thermoneutral temperature and the other half were exposed to high ambient temperature before the NDV challenge and throughout the remainder of the study. Phenotypic responses to heat, to NDV, or to heat + NDV were measured. Selected SNPs (single nucleotide polymorphisms) within 14 target genes or regions were genotyped; and genotype effects on phenotypic responses to NDV or heat + NDV were tested in each individual treatment group and the combined groups. Seventeen significant haplotype effects, among seven genes and seven phenotypes, were detected for response to NDV or heat or NDV + heat. These findings identify specific genetic variants that are associated with response to heat and/or NDV which may be useful in the genetic improvement of chickens to perform favorably when faced with pathogens and heat stress.

Alternative polyadenylation analysis in animals and plants: newly developed strategies for profiling, processing and validation

Alternative polyadenylation is an essential RNA processing event that contributes significantly to regulation of transcriptome diversity and functional dynamics in both animals and plants. Here we review newly developed next generation sequencing methods for genome-wide profiling of alternative polyadenylation (APA) sites, bioinformatics pipelines for data processing and both wet and dry laboratory approaches for APA validation. The library construction methods LITE-Seq (Low-Input 3'-Terminal sequencing) and PAC-seq (PolyA Click sequencing) tag polyA+ cDNA, while BAT-seq (BArcoded, three-prime specific sequencing) and PAPERCLIP (Poly(A) binding Protein-mediated mRNA 3'End Retrieval by CrossLinking ImmunoPrecipitation) enrich polyA+ RNA. Interestingly, only WTTS-seq (Whole Transcriptome Termini Site sequencing) targets both polyA+ RNA and polyA+ cDNA. Varieties of bioinformatics pipelines are well established to pursue read quality control, mapping, clustering, characterization and pathway analysis. The RHAPA (RNase H alternative polyadenylation assay) and 3'RACE-seq (3' rapid amplification of cDNA end sequencing) methods directly validate APA sites, while WTSS-seq (whole transcriptome start site sequencing), RNA-seq (RNA sequencing) and public APA databases can serve as indirect validation methods. We hope that these tools, pipelines and resources trigger huge waves of interest in the research community to investigate APA events underlying physiological, pathological and psychological changes and thus understand the information transfer events from genome to phenome relevant to economically important traits in both animals and plants.