Search

BackgroundThe accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.ResultsA total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.ConclusionsThese validated results show significant improvement over current bovine genome annotations.

BackgroundStructural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence.ResultsWe present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies.ConclusionsThis SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

IntroductionInvestigating the genetic footprints of historical temperature selection can get insights to the local adaptation and feasible influences of climate change on long-term population dynamics.ObjectChicken is a significative species to study genetic adaptation on account of its similar domestication track related to human activity with the most diversified varieties. Yet, few studies have demonstrated the genetic signatures of its adaptation to naturally tropical and frigid environments.MethodHere, we generated whole genome resequencing of 119 domesticated chickens in China including the following breeds which are in order of breeding environmental temperature from more tropical to more frigid: Wenchang chicken (WCC), green-shell chicken (GSC), Tibetan chicken (TBC), and Lindian chicken (LDC).ResultsOur results showed WCC branched off earlier than LDC with an evident genetic admixture between WCC and LDC, suggesting their closer genetic relationship. Further comparative genomic analyses solute carrier family 33 member 1 (SLC33A1) and thyroid stimulating hormone receptor (TSHR) genes exhibited stronger signatures for positive selection in the genome of the more tropical WCC. Furthermore, genotype data from about 3,000 African local ecotypes confirmed that allele frequencies of single nucleotide polymorphisms (SNPs) in these 2 genes appeared strongly associated with tropical environment adaptation. In addition, the NADH:ubiquinone oxidoreductase subunit S4 (NDUFS4) gene exhibited a strong signature for positive selection in the LDC genome, and SNPs with marked allele frequency differences indicated a significant relationship with frigid environment adaptation.ConclusionOur findings partially clarify how selection footprints from environmental temperature stress can lead to advantageous genomic adaptions to tropical and frigid environments in poultry and provide a valuable resource for selective breeding of chickens.

The objective of this study was to investigate the effects of dietary lysine restriction on the global gene expression profile of skeletal muscle in growing pigs. Twelve crossbred (Yorkshire × Landrace) barrows (initial BW 22.6 ± 2.04 kg) were randomly assigned to two dietary treatments (LDD: a lysine-deficient diet; LAD: a lysine-adequate diet) according to a completely randomized experiment design (n = 6). After feeding for 8 weeks, skeletal muscle was sampled from the longissimus dorsi of individual pigs. The muscle total RNA was isolated and cDNA libraries were prepared for RNA sequencing (RNA-Seq) analysis. The RNA-Seq data obtained was then analyzed using the CLC Genomics Workbench to identify differentially expressed genes (DEGs). A total of 80 genes (padj ≤ 0.05) were differentially expressed in the longissimus dorsi muscle of the pigs fed LDD vs. LAD, of which 46 genes were downregulated and 34 genes were upregulated. Gene Ontology (GO) analysis of the DEGs (padj ≤ 0.05) for functional annotation identified those GO terms that are mostly associated with the molecular functions of structural molecules and metabolic enzymes (e.g., oxidoreductase and endopeptidase), biological process of acute-phase response, and amino acid metabolism including synthesis and degradation in the extracellular matrix region. Collectively, the results of this study have provided some novel insight regarding the molecular mechanisms of muscle growth that are associated with dietary lysine supply.

Unlike PIWI-interacting RNA (piRNA) in other species that mostly target transposable elements (TEs), >80% of piRNAs in adult mammalian testes lack obvious targets. However, mammalian piRNA sequences and piRNA-producing loci evolve more rapidly than the rest of the genome for unknown reasons. Here, through comparative studies of chickens, ducks, mice, and humans, as well as long-read nanopore sequencing on diverse chicken breeds, we find that piRNA loci across amniotes experience: (1) a high local mutation rate of structural variations (SVs, mutations ≥ 50 bp in size); (2) positive selection to suppress young and actively mobilizing TEs commencing at the pachytene stage of meiosis during germ cell development; and (3) negative selection to purge deleterious SV hotspots. Our results indicate that genetic instability at pachytene piRNA loci, while producing certain pathogenic SVs, also protects genome integrity against TE mobilization by driving the formation of rapid-evolving piRNA sequences.

Growth differentiation factor-10 (GDF10) with its methylation trait has recently been found to play a crucial regulatory and communication role in cancers. This investigation aims to identify GDF10 methylation site-associated genes that are closely associated with endometrial cancer (EC) patients' survival based on normal and UCEC samples from the UCSC Xena database. Our study revealed for the first time that EC exhibited significantly higher levels of GDF10 promoter methylation in comparison with normal tissues. Multiple differentiated methylation sites, which have prognostic value due to their apparent survival differences, were found in the GDF10 promoter region. We performed weighted gene coexpression network analysis (WGCNA) on EC tissues and paraneoplastic tissues while using these differentially methylated sites as phenotypes for selecting the most correlated key modules and their internal genes. To obtain a gene set, the key module genes and differentially expressed genes (DEGs) of EC were intersected. The least absolute shrinkage and selection operator (LASSO) regression along with multivariate Cox regression were performed from the gene set and we screened out the key genes B4GALNT3, DNAJC22, and GREB1. Finally, a prognostic model was validated for effectiveness based on these genes. Additionally, Kaplan-Meier analysis and time-dependent receiver operating characteristics (ROC) were applied to assess and verify the model, and they showed good prognosis prediction. Moreover, the differences in risk scores were statistically significant with age, tumor stage, and grade. They may be related to the immune infiltration of tumors as well. In conclusion, based on the methylation-related genes associated with GDF10, we developed a prognosis model for EC patients. It might provide a fresh view for further research and treatment of EC.

Genomics encompasses the entire tree of life, both extinct and extant, and the evolutionary processes that shape this diversity. To date, genomic research has focused on humans, a small number of agricultural species, and established laboratory models. Fewer than 18,000 of ∼2,000,000 eukaryotic species (

BackgroundCross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes.ResultsHere, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues.ConclusionsIn summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.

To identify and annotate transcript isoforms in the chicken genome, we generated Nanopore long-read sequencing data from 68 samples that encompassed 19 diverse tissues collected from experimental adult male and female White Leghorn chickens. More than 23.8 million reads with mean read length of 790 bases and average quality of 18.2 were generated. The annotation and subsequent filtering resulted in the identification of 55,382 transcripts at 40,547 loci with mean length of 1,700 bases. We predicted 30,967 coding transcripts at 19,461 loci, and 16,495 lncRNA transcripts at 15,512 loci. Compared to existing reference annotations, we found ∼52% of annotated transcripts could be partially or fully matched while ∼47% were novel. Seventy percent of novel transcripts were potentially transcribed from lncRNA loci. Based on our annotation, we quantified transcript expression across tissues and found two brain tissues (i.e., cerebellum and cortex) expressed the highest number of transcripts and loci. Furthermore, ∼22% of the transcripts displayed tissue specificity with the reproductive tissues (i.e., testis and ovary) exhibiting the most tissue-specific transcripts. Despite our wide sampling, ∼20% of Ensembl reference loci were not detected. This suggests that deeper sequencing and additional samples that include different breeds, cell types, developmental stages, and physiological conditions, are needed to fully annotate the chicken genome. The application of Nanopore sequencing in this study demonstrates the usefulness of long-read data in discovering additional novel loci (e.g., lncRNA loci) and resolving complex transcripts (e.g., the longest transcript for the TTN locus).

BackgroundBeef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus.ResultsComparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity.ConclusionsWe present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating 47 human genome-wide association studies, we demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases.

Infectious bronchitis virus (IBV) induces respiratory and urogenital disease in chickens. Although IBV replicates in the gastrointestinal tract, enteric lesions are uncommon. We have reported a case of runting-stunting syndrome in commercial broilers from which an IBV variant was isolated from the intestines. The isolate, CalEnt, demonstrated an enteric tissue tropism in chicken embryos and SPF chickens experimentally. Here, we determined the full genome of CalEnt and compared it to other IBV strains, in addition to comparing the pathobiology of CalEnt and M41 in commercial broilers. Despite the high whole-genome identity to other IBV strains, CalEnt is rather unique in its nucleotide composition. The S gene phylogenetic analyses showed great similarity between CalEnt and Cal 99. Clinically, vent staining was slightly more frequent in CalEnt-infected birds than those challenged with M41. Furthermore, IBV IHC detection was more evident and the viral shedding in feces was overall higher with the CalEnt challenge compared with M41. Despite underlying intestinal lesions caused by coccidiosis and salmonellosis vaccination, microscopic lesions in CalEnt-infected chickens were more severe than in M41-infected chickens or controls, supporting the enteric tropism of CalEnt. Further studies in SPF chickens are needed to determine the pathogenesis of the virus, its molecular mechanisms for the enteric tropism, and its influence in intestinal health.

Characterizing transcription start sites is essential for understanding the regulatory mechanisms that control gene expression. Recently, a new bovine genome assembly (ARS-UCD1.2) with high continuity, accuracy, and completeness was released; however, the functional annotation of the bovine genome lacks precise transcription start sites and contains a low number of transcripts in comparison to human and mouse. By using the RAMPAGE approach, this study identified transcription start sites at high resolution in a large collection of bovine tissues. We found several known and novel transcription start sites attributed to promoters of protein-coding and lncRNA genes that were validated through experimental and in silico evidence. With these findings, the annotation of transcription start sites in cattle reached a level comparable to the mouse and human genome annotations. In addition, we identified and characterized transcription start sites for antisense transcripts derived from bidirectional promoters, potential lncRNAs, mRNAs, and pre-miRNAs. We also analyzed the quantitative aspects of RAMPAGE to produce a promoter activity atlas, reaching highly reproducible results comparable to traditional RNA-seq. Coexpression networks revealed considerable use of tissue-specific promoters, especially between brain and testicle, which expressed several genes in common from alternate loci. Furthermore, regions surrounding coexpressed modules were enriched in binding factor motifs representative of each tissue. The comprehensive annotation of promoters in such a large collection of tissues will substantially contribute to our understanding of gene expression in cattle and other mammalian species, shortening the gap between genotypes and phenotypes.

Highly pathogenic avian influenza viruses (HPAIVs) in gallinaceous poultry are associated with viral infection of the endothelium, the induction of a 'cytokine storm, and severe disease. In contrast, in Pekin ducks, HPAIVs are rarely endothelial tropic, and a cytokine storm is not observed. To date, understanding these species-dependent differences in pathogenesis has been hampered by the absence of a pure culture of duck and chicken endothelial cells. Here, we use our recently established in vitro cultures of duck and chicken aortic endothelial cells to investigate species-dependent differences in the response of endothelial cells to HPAIV H5N1 infection. We demonstrate that chicken and duck endothelial cells display a different transcriptional response to HPAI H5N1 infection in vitro-with chickens displaying a more pro-inflammatory response to infection. As similar observations were recorded following in vitro stimulation with the viral mimetic polyI:C, these findings were not specific to an HPAIV H5N1 infection. However, similar species-dependent differences in the transcriptional response to polyI:C were not observed in avian fibroblasts. Taken together, these data demonstrate that chicken and duck endothelial cells display a different response to HPAIV H5N1 infection, and this may help account for the species-dependent differences observed in inflammation in vivo.

A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome - or any mammalian transcriptome - appears significantly more feasible.

Newcastle disease (ND) has a great impact on poultry health and welfare with its most virulent (velogenic) strain. In addition, issues exacerbated by the increase in global temperatures necessitates a greater understanding of the host immune response when facing a combination of biotic and abiotic stress factors in poultry production. Previous investigations have revealed that the host immune response is tissue-specific. The goal of this study was to identify genes and/or signaling pathways associated with immune response to NDV (Newcastle disease virus) in the trachea, an essential organ where NDV replicate after the infection, by profiling the tissue specific transcriptome response in two genetically distinct inbred chicken lines when exposed to both abiotic and biotic stressors. Fayoumis appear to be able to respond more effectively (lower viral titer, higher antibody levels, immune gene up-regulation) and earlier than Leghorns. Our results suggest NDV infection in Fayoumis appears to elicit proinflammatory processes, and pathways such as the inhibition of cell viability, cell proliferation of lymphocytes, and transactivation of RNA, more rapidly than in Leghorns. These differences in immune response converge at later timepoints which may indicate that Leghorns eventually regulate its immune response to infection. The profiling of the gene expression response in the trachea adds to our understanding of the chicken host response to NDV infection and heat stress on a whole genome level and provides potential candidate genes and signaling pathways for further investigation into the characterization of the time-specific and pathway specific responses in Fayoumis and Leghorns.

In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources.

How have complex brains evolved from simple circuits? Here we investigated brain region evolution at cell-type resolution in the cerebellar nuclei, the output structures of the cerebellum. Using single-nucleus RNA sequencing in mice, chickens, and humans, as well as STARmap spatial transcriptomic analysis and whole-central nervous system projection tracing, we identified a conserved cell-type set containing two region-specific excitatory neuron classes and three region-invariant inhibitory neuron classes. This set constitutes an archetypal cerebellar nucleus that was repeatedly duplicated to form new regions. The excitatory cell class that preferentially funnels information to lateral frontal cortices in mice becomes predominant in the massively expanded human lateral nucleus. Our data suggest a model of brain region evolution by duplication and divergence of entire cell-type sets.

BackgroundAlthough considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues.ResultsOverall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions.ConclusionsThe lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

Newcastle disease virus (NDV) replication depends on the translation machinery of the host cell; therefore, the eukaryotic translation initiation factor 2 (eIF2) gene family is a likely candidate for control of viral replication. We hypothesized that differential expression of host genes related to translation and innate immune response could contribute to differential resistance to NDV in inbred Fayoumi and Leghorn lines. The expression of twenty-one genes related to the interferon signaling pathway and the eIF2 family was evaluated at two- and six-days post infection (dpi) in the spleen from both lines, either challenged by NDV or nonchallenged. Higher expression of OASL in NDV challenged versus nonchallenged spleen was observed in Leghorns at 2 dpi. Lower expression of EIF2B5 was found in NDV challenged than nonchallenged Fayoumis and Leghorns at 2 dpi. At 2 dpi, NDV challenged Fayoumis had lower expression of EIF2B5 and EIF2S3 than NDV challenged Leghorns. At 6 dpi, NDV challenged Fayoumis had lower expression of EIF2S3 and EIF2B4 than NDV challenged Leghorns. The genetic line differences in expression of eIF2-related genes may contribute to their differential resistance to NDV and also to understanding the interaction between protein synthesis shut-off and virus control in chickens.

Interferon regulatory factor 7 (IRF7) is known as the master transcription factor of the type I interferon response in mammalian species along with IRF3. Yet birds only have IRF7, while they are missing IRF3, with a smaller repertoire of immune-related genes, which leads to a distinctive immune response in chickens compared to in mammals. In order to understand the functional role of IRF7 in the regulation of the antiviral response against avian influenza virus in chickens, we generated IRF7-/- chicken embryonic fibroblast (DF-1) cell lines and respective controls (IRF7wt) by utilizing the CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9) system. IRF7 knockout resulted in increased viral titers of low pathogenic avian influenza viruses. Further RNA-sequencing performed on H6N2-infected IRF7-/- and IRF7wt cell lines revealed that the deletion of IRF7 resulted in the significant down-regulation of antiviral effectors and the differential expression of genes in the MAPK (mitogen-activated protein kinase) and mTOR (mechanistic target of rapamycin) signaling pathways. Dynamic gene expression profiling of the host response between the wildtype and IRF7 knockout revealed potential signaling pathways involving AP1 (activator protein 1), NF-κB (nuclear factor kappa B) and inflammatory cytokines that may complement chicken IRF7. Our findings in this study provide novel insights that have not been reported previously, and lay a solid foundation for enhancing our understanding of the host antiviral response against the avian influenza virus in chickens.

Two environmental factors, Newcastle disease and heat stress, are concurrently negatively impacting poultry worldwide and warrant greater attention into developing genetic resistance within chickens. Using two genetically distinct and highly inbred layer lines, Fayoumi and Leghorn, we explored how different genetic backgrounds affect the bursal response to a treatment of simultaneous Newcastle disease virus (NDV) infection at 6 days postinfection (dpi) while under chronic heat stress. The bursa is a primary lymphoid organ within birds and is crucial for the development of B cells. We performed RNA-seq and ChIP-seq targeting histone modifications on bursa tissue. Differential gene expression revealed that Leghorn, compared to Fayoumi, had significant down-regulation in genes involved in cell proliferation, cell cycle, and cell division. Interestingly, we also found greater differences in histone modification levels in response to treatment in Leghorns than Fayoumis, and biological processes enriched in associated target genes of H3K27ac and H3K4me1 were similarly associated with cell cycle and receptor signaling of lymphocytes. Lastly, we found candidate variants between the two genetic lines within exons of differentially expressed genes and regulatory elements with differential histone modification enrichment between the lines, which provides a strong foundation for understanding the effects of genetic variation on NDV resistance under heat stress. This study provides further understanding of the cellular mechanisms affected by NDV infection under heat stress in chicken bursa and identified potential genes and regulatory regions that may be targets for developing genetic resistance within chickens.

Heat stress results in reduced productivity, anorexia, and mortality in chickens. The objective of the study was to identify genes and signal pathways associated with heat stress and Newcastle disease virus (NDV) infection in the liver of chickens through RNA-seq analysis, using two highly inbred chicken lines (Leghorn and Fayoumi). All birds were held in the same environment until 14 days of age. On day 14, half the birds were exposed to 38 °C with 50% relative humidity for 4 h, then 35 °C until the end of the experiment. The remaining birds were kept at 25 °C throughout the experiment. The heat-treated birds were inoculated at 21 days of age with 107 EID50 (One EID50 unit is the amount of virus that will infect 50 percent of inoculated embryos) NDV La Sota strain to investigate the effects of both heat stress and NDV infection. Physiological parameters were recorded as blood phenotypes at three stages: acute heat (AH), chronic heat (CH1), and chronic heat combined with NDV infection (CH&NDV), at 4 h, 7 days, and 10 days post-initiation of heat treatment, respectively. Our previous work revealed that the heat-resilient Fayoumi line maintained a more stable acid-base balance in their blood compared to the Leghorn line. Liver samples were harvested on both AH and CH&NDV to characterize the transcriptome profiles of these two inbred lines. Both genetic lines and treatments had large impact on the liver transcriptome. Fayoumi birds had more differentially expressed genes (DEGs) than Leghorn birds for both treatments. Metabolic and immune-related genes were on the DEG list, with Fayoumi having more immune-related DEGs than Leghorns, which was confirmed by gene functional enrichment analysis. Weighted correlation network analysis (WGCNA) indicated that the driver genes such as Solute Carrier Family genes could be very important for stabilizing the acid-base balance in Fayoumi birds during heat stress. Therefore, candidate genes such solute carrier family genes could be potential genetic targets that are regulated by Fayoumis to maintain physical hemostasis under heat stress. Differential gene expression showed that Leghorns mainly performed metabolic regulation in response to heat stress and NDV infection, while Fayoumis regulated both immune and metabolic functions. This study provides novel insights and enhances our understandings of liver response to heat stress of heat resilient and susceptible inbred chicken lines.

Newcastle disease (ND) is one of the most challenging infectious diseases affecting poultry production in Africa, causing major economic losses. To date, Newcastle disease virus isolates from several African countries have been grouped into class II NDV genotypes I, IV, V, VI, VII, XI, XIII, XIV, XVII, XVIII and XXI. Although ND is endemic in many African countries, information on circulating genotypes is still scarce. In Tanzania, outbreaks with genotypes V and XIII have been reported. In West and Central Africa, genotypes XIV, XVII, and XVIII are the most predominant. To investigate other genotypes circulating in Tanzania and Ghana, we performed molecular genotyping on isolates from Tanzania and Ghana using the MinION, a third-generation portable sequencing device from Oxford Nanopore Technologies. Using the MinION, we successfully sequenced the NDV F gene hypervariable region of 24 isolates from Tanzania and four samples from Ghana. In Tanzania, genotypes V, VII and XIII were detected. All isolates from Ghana belonged to genotype XVIII. The data obtained in this study reflect the genetic diversity of NDV in Africa and highlight the importance of surveillance for monitoring the distribution of NDV genotypes and viral evolution.

Newcastle disease (ND) is a global threat to domestic poultry, especially in rural areas of Africa and Asia, where the loss of entire backyard local chicken flocks often threatens household food security and income. To investigate the genetics of Ghanaian local chicken ecotypes to Newcastle disease virus (NDV), in this study, three popular Ghanaian chicken ecotypes (regional populations) were challenged with a lentogenic NDV strain at 28 days of age. This study was conducted in parallel with a similar study that used three popular Tanzanian local chicken ecotypes and after two companion studies in the United States, using Hy-line Brown commercial laying birds. In addition to growth rate, NDV response traits were measured following infection, including anti-NDV antibody levels [pre-infection and 10 days post-infection (dpi)], and viral load (2 and 6 dpi). Genetic parameters were estimated, and two genome-wide association study analysis methods were used on data from 1,440 Ghanaian chickens that were genotyped on a chicken 600K Single Nucleotide Polymorphism (SNP) chip. Both Ghana and Tanzania NDV challenge studies revealed moderate to high (0.18 - 0.55) estimates of heritability for all traits, except viral clearance where the heritability estimate was not different from zero for the Tanzanian ecotypes. For the Ghana study, 12 quantitative trait loci (QTL) for growth and/or response to NDV from single-SNP analyses and 20 genomic regions that explained more than 1% of genetic variance using the Bayes B method were identified. Seven of these windows were also identified as having at least one significant SNP in the single SNP analyses for growth rate, anti-NDV antibody levels, and viral load at 2 and 6 dpi. An important gene for growth during stress, CHORDC1 associated with post-infection growth rate was identified as a positional candidate gene, as well as other immune related genes, including VAV2, IL12B, DUSP1, and IL17B. The QTL identified in the Ghana study did not overlap with those identified in the Tanzania study. However, both studies revealed QTL with genes vital for growth and immune response during NDV challenge. The Tanzania parallel study revealed an overlapping QTL on chromosome 24 for viral load at 6 dpi with the US NDV study in which birds were challenged with NDV under heat stress. This QTL region includes genes related to immune response, including TIRAP, ETS1, and KIRREL3. The moderate to high estimates of heritability and the identified QTL suggest that host response to NDV of local African chicken ecotypes can be improved through selective breeding to enhance increased NDV resistance and vaccine efficacy.

Disease has large negative impacts on poultry production. A more comprehensive understanding of host-pathogen interaction can lead to new and improved strategies to maintain health. In particular, host genetic factors can lead to a more effective response to pathogens, hereafter termed resistance. Fayoumi and Leghorn chicken lines have demonstrated relative resistance and susceptibility, respectively, to the Newcastle disease virus (NDV) vaccine strain and many other pathogens. This biological model was used to better understand the host response to a vaccine strain of NDV across three tissues and time points, using RNA-seq. Analyzing the Harderian gland, trachea, and lung tissues together using weighted gene co-expression network analysis (WGCNA) identified important genes that were co-expressed and associated with parameters including: genetic line, days post-infection (dpi), challenge status, sex, and tissue. Pathways and driver genes, such as EIF2AK2, MPEG1, and TNFSF13B, associated with challenge status, dpi, and genetic line were of particular interest as candidates for disease resistance. Overall, by jointly analyzing the three tissues, this study identified genes and gene networks that led to a more comprehensive understanding of the whole animal response to lentogenic NDV than that obtained by analyzing the tissues individually.

As a major infectious disease in chickens, Newcastle disease virus (NDV) causes considerable economic losses in the poultry industry, especially in developing countries where there is limited access to effective vaccination. Therefore, enhancing resistance to the virus in commercial chickens through breeding is a promising way to promote poultry production. In this study, we investigated gene expression changes at 2 and 6 days post inoculation (dpi) at day 21 with a lentogenic NDV in a commercial egg-laying chicken hybrid using RNA sequencing analysis. By comparing NDV-challenged and non-challenged groups, 526 differentially expressed genes (DEGs) (false discovery rate (FDR) < 0.05) were identified at 2 dpi, and only 36 at 6 dpi. For the DEGs at 2 dpi, Ingenuity Pathway Analysis predicted inhibition of multiple signaling pathways in response to NDV that regulate immune cell development and activity, neurogenesis, and angiogenesis. Up-regulation of interferon induced protein with tetratricopeptide repeats 5 (IFIT5) in response to NDV was consistent between the current and most previous studies. Sprouty RTK signaling antagonist 1 (SPRY1), a DEG in the current study, is in a significant quantitative trait locus associated with virus load at 6 dpi in the same population. These identified pathways and DEGs provide potential targets to further study breeding strategy to enhance NDV resistance in chickens.

One of the primary aims of the Functional Annotation of ANimal Genomes (FAANG) initiative is to characterize tissue-specific regulation within animal genomes. To this end, we used chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map four histone modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in eight prioritized tissues collected as part of the FAANG equine biobank from two thoroughbred mares. Data were generated according to optimized experimental parameters developed during quality control testing. To ensure that we obtained sufficient ChIP and successful peak-calling, data and peak-calls were assessed using six quality metrics, replicate comparisons, and site-specific evaluations. Tissue specificity was explored by identifying binding motifs within unique active regions, and motifs were further characterized by gene ontology (GO) and protein-protein interaction analyses. The histone marks identified in this study represent some of the first resources for tissue-specific regulation within the equine genome. As such, these publicly available annotation data can be used to advance equine studies investigating health, performance, reproduction, and other traits of economic interest in the horse.

One of the primary aims of the Functional Annotation of ANimal Genomes (FAANG) initiative is to characterize tissue-specific regulation within animal genomes. To this end, we used chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map four histone modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in eight prioritized tissues collected as part of the FAANG equine biobank from two thoroughbred mares. Data were generated according to optimized experimental parameters developed during quality control testing. To ensure that we obtained sufficient ChIP and successful peak-calling, data and peak-calls were assessed using six quality metrics, replicate comparisons, and site-specific evaluations. Tissue specificity was explored by identifying binding motifs within unique active regions, and motifs were further characterized by gene ontology (GO) and protein-protein interaction analyses. The histone marks identified in this study represent some of the first resources for tissue-specific regulation within the equine genome. As such, these publicly available annotation data can be used to advance equine studies investigating health, performance, reproduction, and other traits of economic interest in the horse.

Newcastle disease virus (NDV) is a highly contagious avian pathogen that poses a tremendous threat to poultry producers in endemic zones due to its epidemic potential. To investigate host genetic resistance to NDV while under the effects of heat stress, a genome-wide association study (GWAS) was performed on Hy-Line Brown layer chickens that were challenged with NDV while under high ambient temperature to identify regions associated with host viral titer, circulating anti-NDV antibody titer, and body weight change. A single nucleotide polymorphism (SNP) on chromosome 1 was associated with viral titer at two days post-infection (dpi), while 30 SNPs spanning a quantitative trait loci (QTL) on chromosome 24 were associated with viral titer at 6 dpi. Immune related genes, such as CAMK1d and CCDC3 on chromosome 1, associated with viral titer at 2 dpi, and TIRAP, ETS1, and KIRREL3, associated with viral titer at 6 dpi, were located in two QTL regions for viral titer that were identified in this study. This study identified genomic regions and candidate genes that are associated with response to NDV during heat stress in Hy-Line Brown layer chickens. Regions identified for viral titer on chromosome 1 and 24, at 2 and 6 dpi, respectively, included several genes that have key roles in regulating the immune response.

Avian influenza has caused significant damage to the poultry industry globally. Consequently, efforts have been made to elucidate the disease mechanisms as well as the mechanisms of disease resistance. Here, by investigating two chicken breeds with distinct responses to avian influenza virus (AIV), Leghorn GB2 and Fayoumi M43, we compared their genome, methylation, and transcriptome differences. MX1, HSP90AB1, and HSP90B1 exhibited high degrees of genetic differentiation (FST) between the two species. Except for the MX1-involved direct anti-virus mechanism, we found that at the methylation and transcriptome levels, the more AIV-resistant breed, Fayoumi, exhibited less variation compared with Leghorn after AIV inoculation, which included change trends in differentially expressed regions, top-fold change genes with FDR-corrected p < 0.05, immune response related genes, and housekeeping genes. Fayoumi also showed better consistency regarding changes in methylation and changes at the transcriptome level. Our results suggest a homeostasis hypothesis for avian influenza resistance, with Fayoumi maintaining superior homeostasis at both the epigenetic and gene expression levels. Three candidate genes-MX1, HSP90AB1, and HSP90B1-showed genetic differentiation and altered gene expression, methylation, and protein expression, which merit attention in further functional studies.

Newcastle Disease (ND) is a continuing global threat to domestic poultry, especially in developing countries, where severe outbreaks of velogenic ND virus (NDV) often cause major economic losses to households. Local chickens are of great importance to rural family livelihoods through provision of high-quality protein. To investigate the genetic basis of host response to NDV, three popular Tanzanian chicken ecotypes (regional populations) were challenged with a lentogenic (vaccine) strain of NDV at 28 days of age. Various host response phenotypes, including anti-NDV antibody levels (pre-infection and 10 days post-infection, dpi), and viral load (2 and 6 dpi) were measured, in addition to growth rate. We estimated genetic parameters and conducted genome-wide association study analyses by genotyping 1399 chickens using the Affymetrix 600K chicken SNP chip. Estimates of heritability of the evaluated traits were moderate (0.18-0.35). Five quantitative trait loci (QTL) associated with growth and/or response to NDV were identified by single-SNP analyses, with some regions explaining ≥1% of genetic variance based on the Bayes-B method. Immune related genes, such as ETS1, TIRAP, and KIRREL3, were located in regions associated with viral load at 6 dpi. The moderate estimates of heritability and identified QTL indicate that NDV response traits may be improved through selective breeding of chickens to enhance increased NDV resistance and vaccine efficacy in Tanzanian local ecotypes.

BackgroundNewcastle disease virus, in its most pathogenic form, threatens the livelihood of rural poultry farmers where there is a limited infrastructure and service for vaccinations to prevent outbreaks of the virus. Previously reported studies on the host response to Newcastle disease in chickens have not examined the disease under abiotic stressors, such as heat, which commonly experienced by chickens in regions such as Africa. The objective of this study was to elucidate the underlying biological mechanisms that contribute to disease resistance in chickens to the Newcastle disease virus while under the effects of heat stress.ResultsDifferential gene expression analysis identified genes differentially expressed between treated and non-treated birds across three time points (2, 6, and 10 days post-infection) in Fayoumi and Leghorn birds. Across the three time points, Fayoumi had very few genes differentially expressed between treated and non-treated groups at 2 and 6 days post-infection. However, 202 genes were differentially expressed at 10 days post-infection. Alternatively, Leghorn had very few genes differentially expressed at 2 and 10 days post-infection but had 167 differentially expressed genes at 6 days post-infection. Very few differentially expressed genes were shared between the two genetic lines, and pathway analysis found unique signaling pathways specific to each genetic line. Fayoumi had significantly lower viral load, higher viral clearance, higher anti-NDV antibody levels, and fewer viral transcripts detected compared to Leghorns. Fayoumis activated immune related pathways including SAPK/JNK and p38 MAPK signaling pathways at earlier time points, while Leghorn would activate these same pathways at a later time. Further analysis revealed activation of the GP6 signaling pathway that may be responsible for the susceptible Leghorn response.ConclusionsThe findings in this study confirmed our hypothesis that the Fayoumi line was more resistant to Newcastle disease virus infection compared to the Leghorn line. Within line and interaction analysis demonstrated substantial differences in response patterns between the two genetic lines that was not observed from the within line contrasts. This study has provided novel insights into the transcriptome response of the Harderian gland tissue during Newcastle disease virus infection while under heat stress utilizing a unique resistant and susceptible model.

BackgroundNumerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs.ResultsComprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen.ConclusionsLncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

: Newcastle disease virus (NDV) is a devastating worldwide poultry pathogen with major implications for global food security. In this study, two highly inbred and genetically distinct chicken lines, Fayoumis and Leghorns, were exposed to a lentogenic strain of NDV, while under the effects of heat stress, in order to understand the genetic mechanisms of resistance during high ambient temperatures. Fayoumis, which are relatively more resistant to pathogens than Leghorns, had larger numbers of differentially expressed genes (DEGs) during the early stages of infection when compared to Leghorns and subsequently down-regulated their immune response at the latter stages to return to homeostasis. Leghorns had very few DEGs across all observed time points, with the majority of DEGs involved with metabolic and glucose-related functions. Proteomic analysis corroborates findings made within Leghorns, while also identifying interesting candidate genes missed by expression profiling. Poor correlation between changes observed in the proteomic and transcriptomic datasets highlights the potential importance of integrative approaches to understand the mechanisms of disease response. Overall, this study provides novel insights into global protein and expression profiles of these two genetic lines, and provides potential genetic targets involved with NDV resistance during heat stress in poultry.

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.

During mammalian viral infections, interferon regulatory factor 7 (IRF7) partners with IRF3 to regulate the type I interferon response. In chickens, however, it is still unclear how IRF7 functions in the host innate immune response, especially given that IRF3 is absent. To further elucidate the functional role of chicken IRF7 during avian influenza virus (AIV) infection, we generated inducible IRF7 overexpression DF-1 cell lines and performed in vitro infection using low pathogenic AIVs (LPAIVs). Overexpression of IRF7 resulted in higher viral replication of H6N2 and H10N7 LPAIVs compared to empty vector control cells regardless of IRF7 expression level. In addition, a high rate of induced cell death was observed due to elevated level of IRF7 upon viral infection. RNA-seq and subsequent transcriptome analysis of IRF7 overexpression and control cells discovered candidate genes possibly controlled by chicken IRF7. Functional annotation revealed potential pathways modulated by IRF7 such as TGF-beta signaling pathway, FoxO signaling pathway and cell structural integrity related pathways. Next, we analyzed the host response alteration due to the IRF7 overexpression and additionally discovered the possible connection of chicken IRF7 and JAK-STAT signaling pathway. These findings suggest that chicken IRF7 could modulate a wide range of cellular processes in the host innate immune response thus meticulous control of IRF7 expression is crucial to the host in response to AIV infection.

Enhancing genetic resistance of chickens to Newcastle Disease Virus (NDV) provides a promising way to improve poultry health, and to alleviate poverty and food insecurity in developing countries. In this study, two inbred chicken lines with different responses to NDV, Fayoumi and Leghorn, were challenged with LaSota NDV strain at 21 days of age. Through transcriptome analysis, gene expression in spleen at 2 and 6 days post-inoculation was compared between NDV-infected and control groups, as well as between chicken lines. At a false discovery rate

In low income countries, chickens play a vital role in daily life. They provide a critical source of protein through egg production and meat. Newcastle disease, caused by avian paramyxovirus type 1, has been ranked as the most devastating disease for scavenging chickens in Africa and Asia. High mortality among flocks infected with velogenic strains leads to a devastating loss of dietary protein and buying power for rural households. Improving the genetic resistance of chickens to Newcastle Disease virus (NDV), in addition to vaccination, is a practical target for improvement of poultry production in low income countries. Because response to NDV has a component of genetic control, it can be influenced through selective breeding. Adding genomic information to a breeding program can increase the amount of genetic progress per generation. In this study, we challenged a commercial egg-laying line with a lentogenic strain of NDV, measured phenotypic responses, collected genotypes, and associated genotypes with phenotypes. Collected phenotypes included viral load at 2 and 6 days post-infection (dpi), antibody levels pre-challenge and 10 dpi, and growth rates pre- and post-challenge. Six suggestive QTL associated with response to NDV and/or growth were identified, including novel and known QTL confirming previously reported associations with related traits. Additionally, previous RNA-seq analysis provided support for several of the genes located in or near the identified QTL. Considering the trend of negative genetic correlation between antibody and Newcastle Disease tolerance (growth under disease) and estimates of moderate to high heritability, we provide evidence that these NDV response traits can be influenced through selective breeding. Producing chickens that perform favorably in challenging environments will ultimately increase the supply of quality protein for human consumption.

Behind each eye of the chicken resides a unique lymph tissue, the Harderian gland, for which RNA sequencing (RNA-seq) analysis is novel. We characterized the response of this tissue to Newcastle disease virus (NDV) in two inbred lines with different susceptibility to NDV across three time points. Three-week-old relatively resistant (Fayoumi) and relatively susceptible (Leghorn) birds were inoculated with a high-titered (107EID50) La Sota strain of NDV via an oculonasal route. At 2, 6, and 10 days post infection (dpi) Harderian glands were collected and analyzed via RNA-seq. The Fayoumi had significantly more detectable viral transcripts in the Harderian gland at 2 dpi than the Leghorn, but cleared the virus by 6 dpi. At all three time points, few genes were declared differentially expressed (DE) between the challenged and nonchallenged birds, except for the Leghorns at 6 dpi, and these DE genes were predicted to activate an adaptive immune response. Relative to the Leghorn, the Fayoumi was predicted to activate more immune pathways in both challenged and nonchallenged birds suggesting a more elevated immune system in the Fayoumis under homeostatic conditions. Overall, this study helped characterize the function of this important tissue and its response to NDV.

BackgroundNewcastle disease virus (NDV) is a threat to poultry production worldwide. A better understanding of mechanisms of resistance and susceptibility to this virus will improve measures for NDV prevention and control. Males and females from resistant Fayoumi and susceptible Leghorn lines were either challenged with a lentogenic strain of the virus or given a mock infection at 3 weeks of age. The lung transcriptomes generated by RNA-seq were studied using contrasts across the challenged and nonchallenged birds, the two lines, and three time points post-infection, and by using Weighted Gene Co-expression Network Analysis (WGNCA).ResultsGenetic line and sex had a large impact on the lung transcriptome. When contrasting the challenged and nonchallenged birds, few differentially expressed genes (DEG) were identified within each line at 2, 6, and 10 days post infection (dpi), except for the more resistant Fayoumi line at 10 dpi, for which several pathways were activated and inhibited at this time. The interaction of challenge and line at 10 dpi significantly impacted 131 genes (False Discovery Rate (FDR)

To select metabolic biomarkers and differentially expressed genes (DEGs) associated with resistant-ascites syndrome (resistant-AS), we used innovative techniques such as metabolomics and transcriptomics to comparatively examine resistant-AS chickens and AS controls. Metabolomic evaluation of chicken serum using ultra-performance liquid chromatography-quadruple time-of-flight high-sensitivity mass spectrometry (UPLC-QTOF/HSMS) showed significantly altered lysoPC(18:1), PE(18:3/16:0), PC(20:1/18:3), DG(24:1/22:6/0:0), PS(18:2/18:0), PI(16:0/16:0), PS(18:0/18:1), PS(14:1/14:0), dihydroxyacetone, ursodeoxycholic acid, tryptophan, L-valine, cycloserine, hypoxanthine, and 4-O-Methylmelleolide concentrations on day 21 and LysoPC(18:0), LysoPE(20:1/0:0), LysoPC(16:0), LysoPE(16:0/0:0), hypoxanthine, dihydroxyacetone, 4-O-Methylmelleolide, LysoPC(18:2), and PC(14:1/22:1) concentrations on day 35, between the susceptible and resistant groups. Compared to the susceptible group, transcriptomic analysis of liver samples using RNA-seq revealed 413 DEGs on day 21 and 214 DEGs on day 35 in the resistant group. Additional evaluations using gene ontology (GO) indicate that significant enrichment occurred in the oxygen transportation, defensive reactions, and protein modifications of the decreased DEGs as well as in the cell morphological formation, neural development, and transforming growth factor (TGF)-beta signalling of the increased DEGs on day 21. Oxygen transportation was also significantly enriched for downregulated DEGs on day 35. The combinatory evaluation of the metabolome and the transcriptome suggests the possible involvement of glycerophospholipid metabolism in the development of resistant-AS in broilers.

The Functional Annotation of Animal Genomes (FAANG) Consortium recently held a Gathering On FAANG (GO-FAANG) Workshop in Washington, DC on October 7-8, 2015. This consortium is a grass-roots organization formed to advance the annotation of newly assembled genomes of domesticated and non-model organisms (www.faang.org). The workshop gathered together from around the world a group of 100+ genome scientists, administrators, representatives of funding agencies and commodity groups to discuss the latest advancements of the consortium, new perspectives, next steps and implementation plans. The workshop was streamed live and recorded, and all talks, along with speaker slide presentations, are available at www.faang.org. In this report, we describe the major activities and outcomes of this meeting. We also provide updates on ongoing efforts to implement discussions and decisions taken at GO-FAANG to guide future FAANG activities. In summary, reference datasets are being established under pilot projects; plans for tissue sets, morphological classification and methods of sample collection for different tissues were organized; and core assays and data and meta-data analysis standards were established.

We describe the organization of a nascent international effort, the Functional Annotation of Animal Genomes (FAANG) project, whose aim is to produce comprehensive maps of functional elements in the genomes of domesticated animal species.

BackgroundDNA cytosine methylation is an important epigenetic modification that has significant effects on a variety of biological processes in animals. Avian species hold a crucial position in evolutionary history. In this study, we used whole-genome bisulfite sequencing (MethylC-seq) to generate single base methylation profiles of lungs in two genetically distinct and highly inbred chicken lines (Fayoumi and Leghorn) that differ in genetic resistance to multiple pathogens, and we explored the potential regulatory role of DNA methylation associated with immune response differences between the two chicken lines.MethodsThe MethylC-seq was used to generate single base DNA methylation profiles of Fayoumi and Leghorn birds. In addition, transcriptome profiling using RNA-seq from the same chickens and tissues were obtained to interrogate how DNA methylation regulates gene transcription on a genome-wide scale.ResultsThe general DNA methylation pattern across different regions of genes was conserved compared to other species except for hyper-methylation of repeat elements, which was not observed in chicken. The methylation level of miRNA and pseudogene promoters was high, which indicates that silencing of these genes may be partially due to promoter hyper-methylation. Interestingly, the promoter regions of more recently evolved genes tended to be more highly methylated, whereas the gene body regions of evolutionarily conserved genes were more highly methylated than those of more recently evolved genes. Immune-related GO (Gene Ontology) terms were significantly enriched from genes within the differentially methylated regions (DMR) between Fayoumi and Leghorn, which implicates DNA methylation as one of the regulatory mechanisms modulating immune response differences between these lines.ConclusionsThis study establishes a single-base resolution DNA methylation profile of chicken lung and suggests a regulatory role of DNA methylation in controlling gene expression and maintaining genome transcription stability. Furthermore, profiling the DNA methylomes of two genetic lines that differ in disease resistance provides a unique opportunity to investigate the potential role of DNA methylation in host disease resistance. Our study provides a foundation for future studies on epigenetic modulation of host immune response to pathogens in chickens.

We describe the organization of a nascent international effort, the Functional Annotation of Animal Genomes (FAANG) project, whose aim is to produce comprehensive maps of functional elements in the genomes of domesticated animal species.

In order to develop novel strategies to protect against increasingly virulent bird-linked pathogens, a better understanding of the avian antiviral response mechanism is essential. Type I interferons (IFNs) are recognized as the first line of defense in a host's antiviral response; and it has been suggested that IRF7, a member of the IFN regulatory factor (IRF) family, plays an important role in modulating the immune response to avian influenza virus infection in chickens. The objective of this study was to identify candidate genes and pathways associated with IRF7 regulation at the transcriptome level as a first step towards elucidating the underlying cellular mechanisms of IRF7 modulation in the chicken antiviral response. IRF7 overexpression and knockdown DF-1 cell lines were established and stimulated by various pathogen-associated molecular patterns. Significant IRF7 and type I IFN expression changes were observed in both the IRF7 overexpression cell line and the IRF7 knockdown cell line upon exposure to the double stranded RNA (dsRNA) analog poly(I:C). Using RNA-seq based transcriptome analysis, we identified potential novel genes that IRF7 may help regulate as part of the host immune response to dsRNA; potential biomarkers and therapeutic targets revealed as a result of this study warrant further investigation. Based on our results, we suggest that IRF7 may have conserved functional activity in the avian antiviral response, and plays a crucial role in type I IFN regulation.

Avian gastrointestinal (GI) tracts are highly populated with a diverse array of microorganisms that share a symbiotic relationship with their hosts and contribute to the overall health and disease state of the intestinal tract. The microbiome of the young chick is easily prone to alteration in its composition by both exogenous and endogenous factors, especially during the early posthatch period. The genetic background of the host and exposure to pathogens can impact the diversity of the microbial profile that consequently contributes to the disease progression in the host. The objective of this study was to profile the composition and structure of the gut microbiota in young chickens from two genetically distinct highly inbred lines. Furthermore, the effect of the Salmonella Enteritidis infection on altering the composition makeup of the chicken microbiome was evaluated through the 16S rRNA gene sequencing analysis. One-day-old layer chicks were challenged with S. Enteritidis and the host cecal microbiota profile as well as the degree of susceptibility to Salmonella infection was examined at 2 and 7 days post infection. Our result indicated that host genotype had a limited effect on resistance to S. Enteritidis infection. Alpha diversity, beta diversity, and overall microbiota composition were analyzed for four factors: host genotype, age, treatment, and postinfection time points. S. Enteritidis infection in young chicks was found to significantly reduce the overall diversity of the microbiota population with expansion of Enterobacteriaceae family. These changes indicated that Salmonella colonization in the GI tract of the chickens has a direct effect on altering the natural development of the GI microbiota. The impact of S. Enteritidis infection on microbial communities was also more substantial in the late stage of infection. Significant inverse correlation between Enterobacteriaceae and Lachnospiraceae family in both non-infected and infected groups, suggested possible antagonistic interaction between members of these two taxa, which could potentially influences the overall microbial population in the gut. Our results also revealed that genetic difference between two lines had minimal effect on the establishment of microbiota population. Overall, this study provided preliminary insights into the contributing role of S. Enteritidis in influencing the overall makeup of chicken's gut microbiota.