Census of orthologous genes and self-organizing maps of biologically relevant transcriptional patterns in chickens (Gallus gallus)

The launch of large-scale chicken expressed sequence tags (EST) projects has placed the chicken in the lead for the number of EST sequences in agriculturally important animals. More than 451,000 chicken ESTs derived from over 158 libraries have been deposited in the NCBI dbEST database as of December 2003. But how many genes these ESTs represent and how they are expressed in different chicken tissues/organs remain undetermined. In the present research, we developed a human gene-based strategy for census of chicken orthologous genes and identification of their expression patterns. Among 34,157 human coding genes used in the study, BLAST analysis revealed that 11,066 genes provisionally matched 248,628 chicken ESTs. Based on the average EST abundance of the orthologous genes, the current public repository of chicken ESTs could represent 20,000 provisional genes. Analysis of gene expression in 14 single tissues/organs showed that approximately 15% of genes were expressed exclusively in single tissue/organ whereas the remaining 85% of genes were co-expressed in two or more tissues/organs. A majority (91.15%) of genes expressed in chicken embryos were also expressed at post-hatch stages, indicating that most genes activated in chicken embryos could serve housekeeping functions. Self-organizing maps (SOM) analysis organized 8807 provisional genes in selected chicken tissues into 98 clusters with each cluster being indicative of common regulatory factors and pathways. A total of 969 provisional orthologous genes were identified as preferentially expressed genes (PEGs) in various chicken tissues/organs (LOD>3.0). No doubt, the present study on gene expression patterns will provide insight into dynamics of metabolic pathways and tissue/organ programming and reprogramming in chickens.     

Generation of EST and microarray resources for functional genomic studies on chicken intestinal health

Expressed sequenced tags (ESTs) and microarray resources have a great impact on the ability to study host response in mice and humans. Unfortunately, these resources are not yet available for domestic farm animals. The aim of this study was to provide genomic resources to study chicken intestinal health, in particular malabsorption syndrome (MAS), which affects mainly the intestine. Therefore a normalized and subtracted cDNA library containing more than 7000 clones was prepared. Randomly chosen clones were sequenced for control purposes. New ESTs were found and multiple ESTs not identified in the chicken intestine before were observed. The number of non-specific ESTs in this cDNA library was low. Based on this normalized and subtracted library a cDNA microarray was made. In a preliminary hybridization experiment with the microarray, genes were identified to be up- or downregulated in MAS infected chickens. This indicates that the generated resources are valuable tools to investigate chicken intestinal health by whole genome expression analysis approaches.     

Generation of EST and Microarray Resources for Functional Genomic Studies on Chicken Intestinal Health

Abstract

Expressed sequenced tags (ESTs) and microarray resources have a great impact on the ability to study host response in mice and humans. Unfortunately, these resources are not yet available for domestic farm animals. The aim of this study was to provide genomic resources to study chicken intestinal health, in particular malabsorption syndrome (MAS), which affects mainly the intestine. Therefore a normalized and subtracted cDNA library containing more than 7000 clones was prepared. Randomly chosen clones were sequenced for control purposes. New ESTs were found and multiple ESTs not identified in the chicken intestine before were observed. The number of non‐specific ESTs in this cDNA library was low. Based on this normalized and subtracted library a cDNA microarray was made. In a preliminary hybridization experiment with the microarray, genes were identified to be up‐ or downregulated in MAS infected chickens. This indicates that the generated resources are valuable tools to investigate chicken intestinal health by whole genome expression analysis approaches.     

Expressed sequence tag analysis of Eimeria-stimulated intestinal intraepithelial lymphocytes in chickens

Intraepithelial lymphocytes (IELs) play a critical role in protective immune response to intestinal pathogens such as Eimeria, the etiologic agent of avian coccidiosis. A list of genes expressed by intestinal IELs of Eimeria-infected chickens was compiled using the expressed sequence tag (EST) strategy. The 14,409 ESTs consisted of 1851 clusters and 7595 singletons, which revealed 9446 unique genes in the data set. Comparison of the sequence data with chicken DNA sequences in GenBank identified 125 novel clones. This EST library will provide a valuable resource for profiling global gene expression in normal and pathogen-infected chickens and identifying additional unique immune-related genes.     

Expressed sequence tag analysis of Eimeria-stimulated intestinal intraepithelial lymphocytes in chickens

Intraepithelial lymphocytes (IELs) play a critical role in protective immune response to intestinal pathogens such as Eimeria, the etiologic agent of avian coccidiosis. A list of genes expressed by intestinal IELs of Eimeria-infected chickens was compiled using the expressed sequence tag (EST) strategy. The 14,409 ESTs consisted of 1851 clusters and 7595 singletons, which revealed 9446 unique genes in the data set. Comparison of the sequence data with chicken DNA sequences in GenBank identified 125 novel clones. This EST library will provide a valuable resource for profiling global gene expression in normal and pathogen-infected chickens and identifying additional unique immune-related genes.     

Classification of sequences expressed during the primordial and basidiome stages of the cultivated mushroom Agaricus bisporus

Two complementary DNA (cDNA) libraries were constructed from tissues isolated from primordia and basidiomes of Agaricus bisporus to characterize genes involved in mushroom development. Using single-pass sequencing of 869 cDNA clones, we found 477 expressed sequence tags (ESTs), including 466 not previously described in the databases for A. bisporus. A BLASTX search revealed that 374 ESTs had similarities with protein sequences available from databases; 193 of these ESTs were categorized according to their putative function. Most ESTs were assigned to one of four roles: metabolism (23%), cell structure (15%), cell growth and division (12%), and protein destination and storage (10%). The remaining ESTs with putative homologues were classified in 10 additional categories. Many ESTs could not be functionally assigned. Based on redundancy levels, at least 4 ESTs were preferentially expressed in each tissue type. Sequence analysis also suggested the presence of paralog tyrosinase genes in the A. bisporus genome.     

Identification of genes related to Parkinson's disease using expressed sequence tags.

In a search for novel target genes related to Parkinson's disease (PD), two full-length cDNA libraries were constructed from a human normal substantia nigra (SN) and a PD patient's SN. An analysis of the gene expression profiles between them was done using the expressed sequence tags (ESTs) frequency. Data for the differently expressed genes were verified by quantitative real-time RT-PCR, immunohistochemical analysis and a cell death assay. Among the 76 genes identified with a significant difference (P > 0.9), 21 upregulated genes and 13 downregulated genes were confirmed to be differentially expressed in human PD tissues and/or in an MPTP-treated mice model by quantitative real-time RT-PCR. Among those genes, an immunohistochemical analysis using an MPTP mice model for alpha-tubulin including TUBA3 and TUBA6 showed that the protein levels are downregulated, as well as the RNA levels. In addition, MBP, PBP and GNAS were confirmed to accelerate cell death activity, whereas SPP1 and TUBA3 to retard this process. Using an analysis of ESTs frequency, it was possible to identify a large number of genes related to human PD. These new genes, MBP, PBP, GNAS, SPP1 and TUBA3 in particular, represent potential biomarkers for PD and could serve as useful targets for elucidating the molecular mechanisms associated with PD.     

Characterization of expressed sequence tags from a gallus gallus pineal gland cDNA library

The pineal gland is the circadian oscillator in the chicken, regulating diverse functions ranging from egg laying to feeding. Here, we describe the isolation and characterization of expressed sequence tags (ESTs) isolated from a chicken pineal gland cDNA library. A total of 192 unique sequences were analysed and submitted to GenBank; 6% of the ESTs matched neither GenBank cDNA sequences nor the newly assembled chicken genomic DNA sequence, three ESTs aligned with sequences designated to be on the Z_random, while one matched a W chromosome sequence and could be useful in cataloguing functionally important genes on this sex chromosome. Additionally, single nucleotide polymorphisms (SNPs) were identified and validated in 10 ESTs that showed 98% or higher sequence similarity to known chicken genes. Here, we have described resources that may be useful in comparative and functional genomic analysis of genes expressed in an important organ, the pineal gland, in a model and agriculturally important organism.     

Cloning and Tissue Expression of Chicken Heart Fatty Acid-Binding Protein and Intestine Fatty Acid-Binding Protein Genes

Fatty acid-binding proteins (FABPs) are members of a superfamily of lipid-binding proteins, occurring intracellularly in invertebrates and vertebrates. This study was designed to clone and characterize the genes of heart fatty acid-binding protein and intestine fatty acid-binding protein in the chicken. PCR primers were designed according to the chicken EST sequences to amplify cDNA of H-FABP and I-FABP genes from chicken heart and intestinal tissues. Analysis of sequence showed that the cDNA of the chicken H-FABP gene is 75 to 77% homologues to human, mouse, and pig H-FABP genes, and the chicken I-FABP gene is 71 to 72% homologues to human, mouse, and pig I-FABP genes. In addition, Northern blot analysis indicated that of the two genes, similar to the copartner of the mammal, H-FABP gene was expressed in a wide variety of tissues, and I-FABP gene was expressed only in intestinal tissues. The expression levels of the chicken H-FABP mRNA in heart and I-FABP mRNA in intestine had significant differences between the broilers from fat line and Bai'er layers at six weeks of age. The results of this study provided basic molecular information for studying the role of two FABPs in the regulation of fatty acid metabolism in avian species.     

Cloning and tissue expression of chicken heart fatty acid-binding protein and intestine fatty acid-binding protein genes

Fatty acid-binding proteins (FABPs) are members of a superfamily of lipid-binding proteins, occurring intracellularly in invertebrates and vertebrates. This study was designed to clone and characterize the genes of heart fatty acid-binding protein and intestine fatty acid-binding protein in the chicken. PCR primers were designed according to the chicken EST sequences to amplify cDNA of H-FABP and I-FABP genes from chicken heart and intestinal tissues. Analysis of sequence showed that the cDNA of the chicken H-FABP gene is 75 to 77% homologues to human, mouse, and pig H-FABP genes, and the chicken I-FABP gene is 71 to 72% homologues to human, mouse, and pig I-FABP genes. In addition, Northern blot analysis indicated that of the two genes, similar to the copartner of the mammal, H-FABP gene was expressed in a wide variety of tissues, and I-FABP gene was expressed only in intestinal tissues. The expression levels of the chicken H-FABP mRNA in heart and I-FABP mRNA in intestine had significant differences between the broilers from fat line and Bai'er layers at six weeks of age. The results of this study provided basic molecular information for studying the role of two FABPs in the regulation of fatty acid metabolism in avian species.     

Sequence of a Complete Chicken BG Haplotype Shows Dynamic Expansion and Contraction of Two Gene Lineages with Particular Expression Patterns

Many genes important in immunity are found as multigene families. The butyrophilin genes are members of the B7 family, playing diverse roles in co-regulation and perhaps in antigen presentation. In humans, a fixed number of butyrophilin genes are found in and around the major histocompatibility complex (MHC), and show striking association with particular autoimmune diseases. In chickens, BG genes encode homologues with somewhat different domain organisation. Only a few BG genes have been characterised, one involved in actin-myosin interaction in the intestinal brush border, and another implicated in resistance to viral diseases. We characterise all BG genes in B12 chickens, finding a multigene family organised as tandem repeats in the BG region outside the MHC, a single gene in the MHC (the BF-BL region), and another single gene on a different chromosome. There is a precise cell and tissue expression for each gene, but overall there are two kinds, those expressed by haemopoietic cells and those expressed in tissues (presumably non-haemopoietic cells), correlating with two different kinds of promoters and 5′ untranslated regions (5′UTR). However, the multigene family in the BG region contains many hybrid genes, suggesting recombination and/or deletion as major evolutionary forces. We identify BG genes in the chicken whole genome shotgun sequence, as well as by comparison to other haplotypes by fibre fluorescence in situ hybridisation, confirming dynamic expansion and contraction within the BG region. Thus, the BG genes in chickens are undergoing much more rapid evolution compared to their homologues in mammals, for reasons yet to be understood.     

Differences of Z chromosome and genomic expression between early- and late-feathering chickens

One duplicated segment on chicken Z chromosome is a causal mutation to the late-feathering phenotype. However, understanding biological process of the late-feathering formation is also of interest to chicken breeding and feather development theory. One hundred and thirty-seven valid single nucleotide polymorphisms (SNPs) from an SNP database were used to perform an association study of the Z chromosome in Xinghua chickens. Two SNPs, which were respectively on 9607480 bp and 10607757 bp, were significantly associated with feathering phenotypes. This result indicated the causal mutation of the late-feathering formation in Xinghua chickens was consistent with the previous report which showed the late-feathering locus ranged 9966364–10142688 bp on Z chromosome. Microarray expressions were implemented for six 1-day-old female Xinghua chicks. Compared to the early-feathering chicks, there were 249 and 83 upregulated and downregulated known genes in the late-feathering chicks. Forty-one genes were expressed in late-feathering chicks, but not in early-feathering ones. At least 14 significantly differentially expressed genes were directly related to keratin. In the region of the sex-linked feathering gene, only prolactin receptor (PRLR) gene was a significantly differentially expressed gene. Expression of PRLR in late-feathering chicks was 1.78-fold as that in early-feathering chicks. Late-feathering Wenchang chicks also had higher expression level of PRLR than early-feathering ones. This study suggested that increasing PRLR expression that resulted from the special variant on chicken Z chromosome caused the late-feathering phenotype.     

Skin Transcriptome Profiles Associated with Skin Color in Chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.