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.     

Differences in gene density on chicken macrochromosomes and microchromosomes

The chicken karyotype comprises six pairs of large macrochromosomes and 33 pairs of smaller microchromosomes^1. Cytogenetic evidence suggests that microchromosomes may be more gene-dense than macrochromosomes. In this paper, we compare the gene densities on macrochromosomes and microchromosomes based on sequence sampling of cloned genomic DNA, and from the distribution of genes mapped by genetic linkage and physical mapping. From these different approaches we estimate that microchromosomes are twice as gene-dense as macrochromosomes and show that sequence sampling is an effective means of gene discovery in the chicken. Using this method we have also detected a conserved linkage between the genes for serotonin 1D receptor ( HTR1D ) and the platelet-activating factor receptor protein gene ( PTAFR ) on chicken chromosome 5 and human chromosome 1p34 ·3. Taken together with its advantages as an experimental animal, and public access to genetic and physical mapping resources, the chicken is a useful model genome for studies on the structure, function and evolution of the vertebrate genome.     

Genome-Wide and Species-Wide In Silico Screening for Intragenic MicroRNAs in Human,Mouse and Chicken

MicroRNAs (miRNAs) are non-coding RNAs (ncRNAs) involved in regulation of gene expression. Intragenic miRNAs, especially those exhibiting a high degree of evolutionary conservation, have been shown to be coordinately regulated and/or expressed with their host genes, either with synergistic or antagonistic correlation patterns. However, the degree of cross-species conservation of miRNA/host gene co-location is not known and co-expression information is incomplete and fragmented among several studies. Using the genomic resources (miRBase and Ensembl) we performed a genome-wide in silico screening (GWISS) for miRNA/host gene pairs in three well-annotated vertebrate species: human, mouse, and chicken. Approximately half of currently annotated miRNA genes resided within host genes: 53.0% (849/1,600) in human, 48.8% (418/855) in mouse, and 42.0% (210/499) in chicken, which we present in a central publicly available Catalog of intragenic miRNAs (http://www.integratomics-time.com/miR-host/catalog). The miRNA genes resided within either protein-coding or ncRNA genes, which include long intergenic ncRNAs (lincRNAs) and small nucleolar RNAs (snoRNAs). Twenty-seven miRNA genes were found to be located within the same host genes in all three species and the data integration from literature and databases showed that most (26/27) have been found to be co-expressed. Particularly interesting are miRNA genes located within genes encoding for miRNA silencing machinery (DGCR8, DICER1, and SND1 in human and Cnot3, Gdcr8, Eif4e, Tnrc6b, and Xpo5 in mouse). We furthermore discuss a potential for phenotype misattribution of miRNA host gene polymorphism or gene modification studies due to possible collateral effects on miRNAs hosted within them. In conclusion, the catalog of intragenic miRNAs and identified 27 miRNA/host gene pairs with cross-species conserved co-location, co-expression, and potential co-regulation, provide excellent candidates for further functional annotation of intragenic miRNAs in health and disease.     

Genome-scale identification of membrane-associated human mRNAs

The subcellular localization of proteins is critical to their biological roles. Moreover, whether a protein is membrane-bound, secreted, or intracellular affects the usefulness of, and the strategies for, using a protein as a diagnostic marker or a target for therapy. We employed a rapid and efficient experimental approach to classify thousands of human gene products as either "membrane-associated/secreted" (MS) or "cytosolic/nuclear" (CN). Using subcellular fractionation methods, we separated mRNAs associated with membranes from those associated with the soluble cytosolic fraction and analyzed these two pools by comparative hybridization to DNA microarrays. Analysis of 11 different human cell lines, representing lymphoid, myeloid, breast, ovarian, hepatic, colon, and prostate tissues, identified more than 5,000 previously uncharacterized MS and more than 6,400 putative CN genes at high confidence levels. The experimentally determined localizations correlated well with in silico predictions of signal peptides and transmembrane domains, but also significantly increased the number of human genes that could be cataloged as encoding either MS or CN proteins. Using gene expression data from a variety of primary human malignancies and normal tissues, we rationally identified hundreds of MS gene products that are significantly overexpressed in tumors compared to normal tissues and thus represent candidates for serum diagnostic tests or monoclonal antibody-based therapies. Finally, we used the catalog of CN gene products to generate sets of candidate markers of organ-specific tissue injury. The large-scale annotation of subcellular localization reported here will serve as a reference database and will aid in the rational design of diagnostic tests and molecular therapies for diverse diseases.     

Fine mapping of the muscular dystrophy (AM) gene on chicken chromosome 2q

Our previous studies revealed that the genetic locus for chicken muscular dystrophy of abnormal muscle (AM) mapped to chromosome 2q, and that the region showed conserved synteny with human chromosome 8q11-24.3. In the current study, we mapped the chicken orthologues of genes from human chromosome 8q11-24 in order to identify the responsible gene. Polymorphisms in the chicken orthologues were identified in the parents of the resource family. Twenty-three genes and expressed sequence tags (ESTs) were mapped to chicken chromosome 2 by linkage analysis. The detailed comparative map shows a high conservation of synteny between chicken chromosome 2q and human chromosome 8q. The AM locus was mapped between [inositol(myo)-1(or4)-monophosphatase 1] (IMPA1) gene and [core-binding factor, runt domain, alpha-subunit 2; translocated to 1; cyclin D-related] (CBFA2T1) gene. The genes located between IMPA1 and CBFA2T1 are the most likely candidates for chicken muscular dystrophy.     

Screening of Rice Genes from the cDNA Catalog Using the Data Obtained by Protein Sequencing

The partial amino acid sequences of 121 rice proteins separated by two-dimensional gel electrophoresis (2D-PAGE), were determined for a protein sequence data file. In the Rice Genome Research Program (RGP), more than 20,000 cDNA clones randomly selected from rice cDNA libraries have been sequenced to construct a cDNA catalog. Complimentary DNAs encoding about 30% of proteins in the protein sequence data file could be identified in the catalog by computer search. It was deduced that 20,000–40,000 genes are present in the rice genome. Only half of about 20,000 cDNAs sequenced in the RGP, corresponding to 1/4–1/2 of genes present in the entire rice genome, should have unique sequences after considering gene redundancy. This is consistent with the fact that the cDNAs encoding about 30% of the sequenced proteins could be identified in the catalog. If the size of the cDNA catalog is enlarged further, cDNAs encoding all proteins separated by 2D-PAGE could be easily identified from the catalog by using the protein sequence data.     

用放射杂交板定位鸡的MC4R基因以及其在鸡和人染色体上同源区的比较分析

黑素皮质素受体(melanocortin-4 receptor,MC4R)基因的突变与猪、鼠和人等的食欲、肥胖和生长有关联性,然而对鸡的MC4R基因的功能却知之甚少。为了确定鸡的MC4R基因在染色体上的位置,使用鸡-仓鼠杂交板(ChickRH6)做了该基因的定位工作。通过扩增ChickRH6杂交板上的93个样品,然后经整合分析将mC4R基因定位在2号染色体上的标记MCW0062、BCL2和OVY附近,即2q12。这个连锁图上的5个标记基于两点分析与MC4R的LOD值都大于5。同时,以MC4R基因为标记做了鸡和人的染色体比较分析。结果显示鸡的2号染色体(GGA2)和人的18号染色体(HSA18)存在同源区,且基因BCL2和肥胖基因(obesity)位于MC4R基因附近。推测鸡的MC4R基因与人的MC4R基因可能具有相似的功能。该研究揭示了鸡和人MC4R基因的染色体分布,并用杂交放射板将鸡的MC4R基因定位在2号染色体的12区带。     

Composition of the lectin pathway of complement in Gallus gallus: absence of mannan-binding lectin-associated serine protease-1 in birds

The lectin pathway of complement is activated by multimolecular complexes that recognize and bind to microbial polysaccharides. These complexes comprise a multimeric carbohydrate recognition subunit (either mannan-binding lectin (MBL) or a ficolin), three MBL-associated serine proteases (MASP-1, -2, and -3), and MAp19 (a truncated product of the MASP-2 gene). In this study we report the cloning of chicken MASP-2, MASP-3, and MAp19 and the organization of their genes and those for chicken MBL and a novel ficolin. Mammals usually possess two MBL genes and two or three ficolin genes, but chickens have only one of each, both of which represent the undiversified ancestors of the mammalian genes. The primary structure of chicken MASP-2 is 54% identical with those of the human and mouse MASP-2, and the organization of its gene is the same as in mammals. MASP-3 is even more conserved; chicken MASP-3 shares approximately 75% of its residues with human and Xenopus MASP-3. It is more widely expressed than other lectin pathway components, suggesting a possible function of MASP-3 different from those of the other components. In mammals, MASP-1 and MASP-3 are alternatively spliced products of a single structural gene. We demonstrate the absence of MASP-1 in birds, possibly caused by the loss of MASP-1-specific exons during phylogeny. Despite the lack of MASP-1-like enzymatic activity in sera of chicken and other birds, avian lectin pathway complexes efficiently activate C4.     

Mapping and Expression Analysis of Chicken NDRG1 and NDRG3 Genes

N-myc downstream-regulated genes 1 and 3 (NDRG1 and NDRG3) are members of the alpha/beta hydrolase superfamily. Phylogenetic analysis of the family demonstrated that human NDRG1 and 3 belong to a subfamily. The mapping and gene expression patterns of these genes represent one step toward further investigation into their possible roles in the chicken (Gallus gallus). To map these genes in the chicken chromosome, a 6000 rads chicken-hamster radiation hybrid panel (ChickRH6) was used. Primers were designed according to the published human sequences for amplification of those two genes. We compared the corresponding human mRNA sequences with the predicted coding sequences of the chicken NDRG1 and 3 genes and found that the assembled contigs shared a high percentage of similarity with the human genes. PCR of samples from ChickRH6 revealed that the locations of NDRG1 and 3 are linked to the markers MYC (58 cRs away, LOD score 4.52) and SEQ0265 (10 cRs away, LOD score 17.81), respectively. This result adds two new markers to the chicken RH map, and it reinforces that the RH technique is indeed a powerful tool for mapping genes due to its rapidity, precision, convenience, and reproducibility. In addition, we detected the gene expression and distribution of chicken NDRG1 and 3 in seven tissues, including heart, liver, spleen, lung, muscle, brain, and thymus, by RT-PCR, and found that NDRG1 is relatively ubiquitously expressed in all the tested tissues and highly expressed in heart and liver, whereas NDRG3 is high in heart, muscle, and brain.     

A strategy to identify positional candidate genes conferring Marek''s disease resistance by integrating DNA microarrays and genetic mapping

Marker-assisted selection (MAS) to enhance genetic resistance to Marek's disease (MD), a herpesvirus-induced T cell cancer in chicken, is an attractive alternative to augment control with vaccines. Our earlier studies indicate that there are many quantitative trait loci (QTL) containing one or more genes that confer genetic resistance to MD. Unfortunately, it is difficult to sufficiently resolve these QTL to identify the causative gene and generate tightly linked markers. One possible solution is to identify positional candidate genes by virtue of gene expression differences between MD resistant and susceptible chicken using deoxyribonucleic acid (DNA) microarrays followed by genetic mapping of the differentially-expressed genes. In this preliminary study, we show that DNA microarrays containing approximately 1200 genes or expressed sequence tags (ESTs) are able to reproducibly detect differences in gene expression between the inbred ADOL lines 63 (MD resistant) and 72 (MD susceptible) of uninfected and Marek's disease virus (MDV)-infected peripheral blood lymphocytes. Microarray data were validated by quantitative polymerase chain reaction (PCR) and found to be consistent with previous literature on gene induction or immune response. Integration of the microarrays with genetic mapping data was achieved with a sample of 15 genes. Twelve of these genes had mapped human orthologues. Seven genes were located on the chicken linkage map as predicted by the human-chicken comparative map, while two other genes defined a new conserved syntenic group. More importantly, one of the genes with differential expression is known to confer genetic resistance to MD while another gene is a prime positional candidate for a QTL.     

Microarray analysis of genes differentially expressed in the liver of lean and fat chickens

Excessive accumulation of lipids in the adipose tissue is one of the main problems faced by the broiler industry nowadays. In chicken, lipogenesis occurs essentially in the liver, in which much of the triglycerides that accumulate in avian adipose tissue are synthesized. In order to better understand the gene expression and its regulation in chicken liver, the gene expression profiles of liver at developmental stages of chicken (1 week, 4 weeks and 7 weeks of age) were investigated and differentially expressed genes between lean and fat chicken lines divergently selected for abdominal fat content for eight generations were screened. Our data indicated that 4 weeks of age was a very important stage on chicken liver lipogenesis compared to 1 week and 7 weeks of age, and the glycometabolism in chicken liver could be related to lipid metabolism and the difference of glycometabolism could be another potential reason for the fat and lean phenotype occurrence besides the difference of lipogenesis in chicken liver. Our result have established groundwork for further study of the basic genetic control of chicken obesity and will benefit chicken research communities as well as researches that use chicken as a model organism for developmental biology and human therapeutics.     

Terrestrial vertebrates have two keratin gene clusters; striking differences in teleost fish

Keratins I and II form the largest subgroups of mammalian intermediate filament (IF) proteins and account as obligatory heteropolymers for the keratin filaments of epithelia. All human type I genes except for the K18 gene are clustered on chromosome 17q21, while all type II genes form a cluster on chromosome 12q13, that ends with the type I gene K18. Highly related keratin gene clusters are found in rat and mouse. Since fish seem to lack a keratin II cluster we screened the recently established draft genomes of a bird (chicken) and an amphibian (Xenopus). The results show that keratin I and II gene clusters are a feature of all terrestrial vertebrates. Because hair with its multiple hair keratins and inner root sheath keratins is a mammalian acquisition, the keratin gene clusters of chicken and Xenopus tropicalis have only about half the number of genes found in mammals. Within the type I clusters all genes have the same orientation. In type II clusters there is a rare gene of opposite orientation. Finally we show that the genes for keratins 8 and 18, which are the first expression pair in embryology, are not only adjacent in mammals, but also in Xenopus and three different fish. Thus neighboring K8 and K18 genes seem a feature shared by all vertebrates. In contrast to the two well defined keratin gene clusters of terrestrial vertebrates, three teleost fish show an excess of type I over type II genes, the lack of a keratin type II gene cluster and a striking dispersal of type I genes, that are probably the result of the teleost-specific whole genome duplication followed by a massive gene loss. This raises the question whether keratin gene clusters extend beyond the ancestral bony vertebrate to cartilage fish and lamprey. We also analyzed the complement of non-keratin IF genes of the chicken. Surprisingly, an additional nuclear lamin gene, previously overlooked by cDNA cloning, is documented on chromosome 10. The two splice variants closely resemble the lamin LIII a + b of amphibia and fish. This lamin gene is lost on the mammalian lineage.     

Fanconi anemia: genetic analysis of a human disease using chicken system

Fanconi anemia (FA) is a rare hereditary disorder characterized by skeletal abnormalities, bone marrow failure, and an increased incidence of cancer. The basic cellular abnormality in FA has been postulated to lie in the DNA repair mechanisms because cells from FA patients display chromosomal breakage, which is particularly remarkable following induction of DNA crosslinks. However, experimental evidence for this hypothesis has been lacking. To test whether DNA repair is really defective in FA cells, we disrupted several FA genes in chicken B cell line DT40. By measuring efficiency of gene conversion and hypermutation at the Immunoglobulin locus, we have shown that DT40 FA mutant cell lines exhibited defects in homologous DNA recombination, and possibly, translesion synthesis. However, levels of sister chromatid exchange, another important cellular event mediated by HR, were not reduced, possibly indicating the role of FA genes only in a subpathway of HR. Our results indicate that chicken DT40 cells could be highly useful in molecular dissection of basic biochemical processes, which are deficient in a human genetic disorder.