Normalized long read RNA sequencing in chicken reveals transcriptome complexity similar to human.

BACKGROUND: Despite the significance of chicken as a model organism, our understanding of the chicken transcriptome is limited compared to human. This issue is common to all non-human vertebrate annotations due to the difficulty in transcript identification from short read RNAseq data. While previous studies have used single molecule long read sequencing for transcript discovery, they did not perform RNA normalization and 5'-cap selection which may have resulted in lower transcriptome coverage and truncated transcript sequences. RESULTS: We sequenced normalised chicken brain and embryo RNA libraries with Pacific Bioscience Iso-Seq. 5' cap selection was performed on the embryo library to provide methodological comparison. From these Iso-Seq sequencing projects, we have identified 60 k transcripts and 29 k genes within the chicken transcriptome. Of these, more than 20 k are novel lncRNA transcripts with ~3 k classified as sense exonic overlapping lncRNA, which is a class that is underrepresented in many vertebrate annotations. The relative proportion of alternative transcription events revealed striking similarities between the chicken and human transcriptomes while also providing explanations for previously observed genomic differences. CONCLUSIONS: Our results indicate that the chicken transcriptome is similar in complexity compared to human, and provide insights into other vertebrate biology. Our methodology demonstrates the potential of Iso-Seq sequencing to rapidly expand our knowledge of transcriptomics.

A comparative analysis of host responses to avian influenza infection in ducks and chickens highlights a role for the interferon-induced transmembrane proteins in viral resistance.

BACKGROUND: Chickens are susceptible to infection with a limited number of Influenza A viruses and are a potential source of a human influenza pandemic. In particular, H5 and H7 haemagglutinin subtypes can evolve from low to highly pathogenic strains in gallinaceous poultry. Ducks on the other hand are a natural reservoir for these viruses and are able to withstand most avian influenza strains. RESULTS: Transcriptomic sequencing of lung and ileum tissue samples from birds infected with high (H5N1) and low (H5N2) pathogenic influenza viruses has allowed us to compare the early host response to these infections in both these species. Chickens (but not ducks) lack the intracellular receptor for viral ssRNA, RIG-I and the gene for an important RIG-I binding protein, RNF135. These differences in gene content partly explain the differences in host responses to low pathogenic and highly pathogenic avian influenza virus in chicken and ducks. We reveal very different patterns of expression of members of the interferon-induced transmembrane protein (IFITM) gene family in ducks and chickens. In ducks, IFITM1, 2 and 3 are strongly up regulated in response to highly pathogenic avian influenza, where little response is seen in chickens. Clustering of gene expression profiles suggests IFITM1 and 2 have an anti-viral response and IFITM3 may restrict avian influenza virus through cell membrane fusion. We also show, through molecular phylogenetic analyses, that avian IFITM1 and IFITM3 genes have been subject to both episodic and pervasive positive selection at specific codons. In particular, avian IFITM1 showed evidence of positive selection in the duck lineage at sites known to restrict influenza virus infection. CONCLUSIONS: Taken together these results support a model where the IFITM123 protein family and RIG-I all play a crucial role in the tolerance of ducks to highly pathogenic and low pathogenic strains of avian influenza viruses when compared to the chicken.

Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis.

A synergistic combination of two next-generation sequencing platforms with a detailed comparative BAC physical contig map provided a cost-effective assembly of the genome sequence of the domestic turkey (Meleagris gallopavo). Heterozygosity of the sequenced source genome allowed discovery of more than 600,000 high quality single nucleotide variants. Despite this heterozygosity, the current genome assembly (∼1.1 Gb) includes 917 Mb of sequence assigned to specific turkey chromosomes. Annotation identified nearly 16,000 genes, with 15,093 recognized as protein coding and 611 as non-coding RNA genes. Comparative analysis of the turkey, chicken, and zebra finch genomes, and comparing avian to mammalian species, supports the characteristic stability of avian genomes and identifies genes unique to the avian lineage. Clear differences are seen in number and variety of genes of the avian immune system where expansions and novel genes are less frequent than examples of gene loss. The turkey genome sequence provides resources to further understand the evolution of vertebrate genomes and genetic variation underlying economically important quantitative traits in poultry. This integrated approach may be a model for providing both gene and chromosome level assemblies of other species with agricultural, ecological, and evolutionary interest.

Systems analysis of immune responses in Marek's disease virus-infected chickens identifies a gene involved in susceptibility and highlights a possible novel pathogenicity mechanism.

Marek's disease virus (MDV) is a highly contagious oncogenic alphaherpesvirus that causes disease that is both a cancer model and a continuing threat to the world's poultry industry. This comprehensive gene expression study analyzes the host response to infection in both resistant and susceptible lines of chickens and inherent expression differences between the two lines following the infection of the host. A novel pathogenicity mechanism, involving the downregulation of genes containing HIC1 transcription factor binding sites as early as 4 days postinfection, was suggested from this analysis. HIC1 drives antitumor mechanisms, suggesting that MDV infection switches off genes involved in antitumor regulation several days before the expression of the MDV oncogene meq. The comparison of the gene expression data to previous QTL data identified several genes as candidates for involvement in resistance to MD. One of these genes, IRG1, was confirmed by single nucleotide polymorphism analysis to be involved in susceptibility. Its precise mechanism remains to be elucidated, although the analysis of gene expression data suggests it has a role in apoptosis. Understanding which genes are involved in susceptibility/resistance to MD and defining the pathological mechanisms of the disease gives us a much greater ability to try to reduce the incidence of this virus, which is costly to the poultry industry in terms of both animal welfare and economics.

Evolution of the chicken Toll-like receptor gene family: a story of gene gain and gene loss.

BACKGROUND: Toll-like receptors (TLRs) perform a vital role in disease resistance through their recognition of pathogen associated molecular patterns (PAMPs). Recent advances in genomics allow comparison of TLR genes within and between many species. This study takes advantage of the recently sequenced chicken genome to determine the complete chicken TLR repertoire and place it in context of vertebrate genomic evolution. RESULTS: The chicken TLR repertoire consists of ten genes. Phylogenetic analyses show that six of these genes have orthologs in mammals and fish, while one is only shared by fish and three appear to be unique to birds. Furthermore the phylogeny shows that TLR1-like genes arose independently in fish, birds and mammals from an ancestral gene also shared by TLR6 and TLR10. All other TLRs were already present prior to the divergence of major vertebrate lineages 550 Mya (million years ago) and have since been lost in certain lineages. Phylogenetic analysis shows the absence of TLRs 8 and 9 in chicken to be the result of gene loss. The notable exception to the tendency of gene loss in TLR evolution is found in chicken TLRs 1 and 2, each of which underwent gene duplication about 147 and 65 Mya, respectively. CONCLUSION: Comparative phylogenetic analysis of vertebrate TLR genes provides insight into their patterns and processes of gene evolution, with examples of both gene gain and gene loss. In addition, these comparisons clarify the nomenclature of TLR genes in vertebrates.

Analysis of the rdd locus in chicken: a model for human retinitis pigmentosa.

PURPOSE: To identify the locus responsible for the blind mutation rdd (retinal dysplasia and degeneration) in chickens and to further characterise the rdd phenotype. METHODS: The eyes of blind and sighted birds were subjected to ophthalmic, morphometric and histopathological examination to confirm and extend published observations. Electroretinography was used to determine age of onset. Birds were crossed to create pedigrees suitable for genetic mapping. DNA samples were obtained and subjected to a linkage search. RESULTS: Measurement of IOP, axial length, corneal diameter, and eye weight revealed no gross morphological changes in the rdd eye. However, on ophthalmic examination, rdd homozygotes have a sluggish pupillary response, atrophic pecten, and widespread pigmentary disturbance that becomes more pronounced with age. Older birds also have posterior subcapsular cataracts. At three weeks of age, homozygotes have a flat ERG indicating severe loss of visual function. Pathological examination shows thinning of the RPE, ONL, photoreceptors and INL, and attenuation of the ganglion cell layer. From 77 classified backcross progeny, 39 birds were blind and 38 sighted. The rdd mutation was shown to be sex-linked and not autosomal as previously described. Linkage analysis mapped the rdd locus to a small region of the chicken Z chromosome with homologies to human chromosomes 5q and 9p. CONCLUSIONS: Ophthalmic, histopathologic, and electrophysiological observations suggest rdd is similar to human recessive retinitis pigmentosa. Linkage mapping places rdd in a region homologous to human chromosomes 9p and 5q. Candidate disease genes or loci include PDE6A, WGN1, and USH2C. This is the first use of genetic mapping in a chicken model of human disease.

Genetic, ophthalmic, morphometric and histopathological analysis of the Retinopathy Globe Enlarged (rge) chicken.

PURPOSE: To identify the locus responsible for rge (retinopathy globe enlarged) in chickens and further characterise the rge phenotype. METHODS: A colony of chickens carrying the rge mutation was rederived from a single heterozygous animal of the original line. The eyes of blind, heterozygous and normal birds were subjected to ophthalmic, morphometric and histopathological examination to confirm and extend published observations. DNA samples were obtained and subjected to a whole genome linkage search. RESULTS: From 138 classified backcross progeny, 56 birds were blind and 82 sighted. Heterozygous birds were indistinguishable from wild type, but homozygotes had sluggish or unresponsive pupils, posterior sub-capsular lens opacities and an atrophic pecten. The fundus appeared normal with no significant pigmentary disturbance, but axial length and eye weight were increased. Pathology revealed focal retinal lesions. Linkage analysis placed the rge locus in a small region of chicken chromosome 1. CONCLUSIONS: rge is a severe recessive retinal dystrophy in chickens, with associated globe enlargement. Linkage mapping has highlighted chicken chromosome 1 in a region most probably homologous to human chromosomes 7q31-35, 21q21 or 22q12-21. Candidate disease loci include RP10 (IMPDH1) and uncharacterised Ushers (USH1E) and glaucoma (GLC1F) loci.

Pivotal Advance: Avian colony-stimulating factor 1 (CSF-1), interleukin-34 (IL-34), and CSF-1 receptor genes and gene products.

Macrophages are involved in many aspects of development, host defense, pathology, and homeostasis. Their normal differentiation, proliferation, and survival are controlled by CSF-1 via the activation of the CSF1R. A recently discovered cytokine, IL-34, was shown to bind the same receptor in humans. Chicken is a widely used model organism in developmental biology, but the factors that control avian myelopoiesis have not been identified previously. The CSF-1, IL-34, and CSF1R genes in chicken and zebra finch were identified from respective genomic/cDNA sequence resources. Comparative analysis of the avian CSF1R loci revealed likely orthologs of mammalian macrophage-specific promoters and enhancers, and the CSF1R gene is expressed in the developing chick embryo in a pattern consistent with macrophage-specific expression. Chicken CSF-1 and IL-34 were expressed in HEK293 cells and shown to elicit macrophage growth from chicken BM cells in culture. Comparative sequence and co-evolution analysis across all vertebrates suggests that the two ligands interact with distinct regions of the CSF1R. These studies demonstrate that there are two separate ligands for a functional CSF1R across all vertebrates.

Isolation and mapping the chicken zona pellucida genes: an insight into the evolution of orthologous genes in different species.

The avian oocyte is surrounded by a specialized extracellular glycoproteinaceous matrix, the perivitelline membrane, which is equivalent to the zona pellucida (ZP) in mammals and the chorion in teleosts. A number of related ZP genes encode the proteins that make up this matrix. These proteins play an important role in the sperm/egg interaction and may be involved in speciation. The human genome is known to contain ZP1, ZP2, ZP3, and ZPB genes, while a ZPAX gene has also been identified in Xenopus. The rapid evolution of these genes has confused the nomenclature and thus orthologous relationships across species. In order to clarify these homologies, we have identified ZP1, ZP2, ZPC, ZPB, and ZPAX genes in the chicken and mapped them to chromosomes 5, 14, 10, 6, and 3, respectively, establishing conserved synteny with human and mouse. The amino acid sequences of these genes were compared to the orthologous genes in human, mouse, and Xenopus, and have given us an insight into the evolution of these genes in a variety of different species. The presence of the ZPAX gene in the chicken has highlighted a pattern of probable gene loss by deletion in mouse and gene inactivation by deletion, and base substitution in human.

Comparative mapping of human Chromosome 19 with the chicken shows conserved synteny and gives an insight into chromosomal evolution.

Human Chromosome 19 (HSA19) is virtually completely sequenced. A complete physical contig map made up of BACs and cosmids is also available for this chromosome. It is, therefore, a rich source of information that we have used as the basis for a comparative mapping study with the chicken. Various orthologs of genes known to map to HSA19 have been mapped in the chicken. Five chicken microchromosomes (two of which were previously undefined) are seen to show conserved synteny with this chromosome, along with individual gene homologs on Chr 1 and another tiny microchromosome. Compared with the mouse, which has 12 chromosomal regions homologous to HSA19, the chicken genotype displays fewer evolutionary rearrangements. The ancestral nature of the chicken karyotype is demonstrated and may prove to be an excellent tool for studying genome evolution.