GOing from functional genomics to biological significance

The chicken genome is sequenced and this, together with microarray and other functional genomics technologies, makes post-genomic research possible in the chicken. At this time, however, such research is hindered by a lack of genomic structural and functional annotations. Bio-ontologies have been developed for different annotation requirements, as well as to facilitate data sharing and computational analysis, but these are not yet optimally utilized in the chicken. Here we discuss genomic annotation and bio-ontologies. We focus specifically on the Gene Ontology (GO), chicken GO annotations and how these can facilitate functional genomics in the chicken. The GO is the most developed and widely used bio-ontology. It is the de facto standard for functional annotation. Despite its critical importance in analyzing microarray and other functional genomics data, relatively few chicken gene products have any GO annotation. When these are available, the average quality of chicken gene products annotations (defined using evidence code weight and annotation depth) is much less than in mouse. Moreover, tools allowing chicken researchers to easily and rapidly use the GO are either lacking or hard to use. To address all of these problems we developed ChickGO and AgBase. Chicken GO annotations are provided by complementary work at MSU-AgBase and EBI-GOA. The GO tools pipeline at AgBase uses GO to derive functional and biological significance from microarray and other functional genomics data. Not only will improved genomic annotation and tools to use these annotations benefit the chicken research community but they will also facilitate research in other avian species and comparative genomics.     

Cross-species overgo hybridization and comparative physical mapping within avian genomes

The chicken genome sequence facilitates comparative genomics within other avian species. We performed cross-species hybridizations using overgo probes designed from chicken genomic and zebra finch expressed sequence tags (ESTs) to turkey and zebra finch BAC libraries. As a result, 3772 turkey BACs were assigned to 336 markers or genes, and 1662 zebra finch BACs were assigned to 164 genes. As expected, cross-hybridization was more successful with overgos within coding sequences than within untranslated region, intron or flanking sequences and between chicken and turkey, when compared with chicken-zebra finch or zebra finch-turkey cross-hybridization. These data contribute to the comparative alignment of avian genome maps using a 'one sequence, multiple genomes' strategy.     

Application of chicken microarrays for gene expression analysis in other avian species

Background

With the threat of emerging infectious diseases such as avian influenza, whose natural hosts are thought to be a variety of wild water birds including duck, we are armed with very few genomic resources to investigate large scale immunological gene expression studies in avian species. Multiple options exist for conducting large gene expression studies in chickens and in this study we explore the feasibility of using one of these tools to investigate gene expression in other avian species.

Results

In this study we utilised a whole genome long oligonucleotide chicken microarray to assess the utility of cross species hybridisation (CSH). We successfully hybridised a number of different avian species to this array, obtaining reliable signals. We were able to distinguish ducks that were infected with avian influenza from uninfected ducks using this microarray platform. In addition, we were able to detect known chicken immunological genes in all of the hybridised avian species.

Conclusion

Cross species hybridisation using long oligonucleotide microarrays is a powerful tool to study the immune response in avian species with little available genomic information. The present study validated the use of the whole genome long oligonucleotide chicken microarray to investigate gene expression in a range of avian species.

     

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.     

The grapevine gene nomenclature system

Background

Grapevine (Vitis vinifera L.) is one of the most important fruit crops in the world and serves as a valuable model for fruit development in woody species. A major breakthrough in grapevine genomics was achieved in 2007 with the sequencing of the Vitis vinifera cv. PN40024 genome. Subsequently, data on structural and functional characterization of grape genes accumulated exponentially. To better exploit the results obtained by the international community, we think that a coordinated nomenclature for gene naming in species with sequenced genomes is essential. It will pave the way for the accumulation of functional data that will enable effective scientific discussion and discovery. The exploitation of data that were generated independently of the genome release is hampered by their heterogeneous nature and by often incompatible and decentralized storage. Classically, large amounts of data describing gene functions are only available in printed articles and therefore remain hardly accessible for automatic text mining. On the other hand, high throughput “Omics” data are typically stored in public repositories, but should be arranged in compendia to better contribute to the annotation and functional characterization of the genes.

Results

With the objective of providing a high quality and highly accessible annotation of grapevine genes, the International Grapevine Genome Project (IGGP) commissioned an international Super-Nomenclature Committee for Grape Gene Annotation (sNCGGa) to coordinate the effort of experts to annotate the grapevine genes. The goal of the committee is to provide a standard nomenclature for locus identifiers and to define conventions for a gene naming system in this paper.

Conclusions

Learning from similar initiatives in other plant species such as Arabidopsis, rice and tomato, a versatile nomenclature system has been developed in anticipation of future genomic developments and annotation issues. The sNCGGa’s first outreach to the grape community has been focused on implementing recommended guidelines for the expert annotators by: (i) providing a common annotation platform that enables community-based gene curation, (ii) developing a gene nomenclature scheme reflecting the biological features of gene products that is consistent with that used in other organisms in order to facilitate comparative analyses.     

Evolutionary constraint in flanking regions of avian genes

An important comprehension from comparative genomic analysis is that sequence conservation beyond neutral expectations is frequently found outside protein-coding regions, indicating important functional roles of noncoding DNA. Understanding the causes of constraint on noncoding sequence evolution forms an important area of research, not least in light of the importance for understanding the evolution of gene expression. We aligned all orthologous genes of chicken and zebra finch together with 5 kb of their upstream and downstream noncoding sequences, to study the evolution of gene flanking sequences in the avian genome. Using ancestral repeats as a neutral reference, we detected significant evolutionary constraint in the 3' flanking region, highest directly after termination (60%) and then gradually decreasing to about 20% 5 kb downstream. Constraint was higher in annotated 3' untranslated regions (UTRs) than in non-UTRs at the same distance from the stop codon and higher in sequences annotated as microRNA (miRNA)-binding sites than in non-miRNA-binding sites within 3' UTRs. Constraint was also higher when estimated for a smaller data set of genes from more closely related songbird species, indicating turnover of functional elements during avian evolution. On the 5' flanking side constraint was readily seen within the first 125 bp immediately upstream of the start codon (34%) and was about 10% for remaining sequence within 5 kb upstream. Analysis of chicken polymorphism data gave further support for the highest constraint directly before and after the translated region. Finally, we found that genes evolving under the highest constraint measured by d(N)/d(S) also had the highest level of constraint in the 3' flanking region. This study broadens the insights into gene flanking sequence evolution by adding new findings from a vertebrate lineage other than mammals.     

Using the chicken genome sequence in the development and mapping of genetic markers in the turkey (Meleagris gallopavo)

The efficacy of employing the chicken genome sequence in developing genetic markers and in mapping the turkey genome was studied. Eighty previously uncharacterized microsatellite markers were identified for the turkey using BLAST alignment to the chicken genome. The chicken sequence was then used to develop primers for polymerase chain reaction where the turkey sequence was either unavailable or insufficient. A total of 78 primer sets were tested for amplification and polymorphism in the turkey, and informative markers were genetically mapped. Sixty-five (83%) amplified turkey genomic DNA, and 33 (42%) were polymorphic in the University of Minnesota/Nicholas Turkey Breeding Farms mapping families. All but one marker genetically mapped to the position predicted from the chicken genome sequence. These results demonstrate the usefulness of the chicken sequence for the development of genomic resources in other avian species.     

稻属植物的基因组进化

水稻所在的稻属(Oryza)共有24个左右的物种。由于野生稻含有大量的优良农艺性状基因, 在水稻遗传学研究中日益受到重视。随着国际稻属基因组计划的开展, 越来越多的稻属基因组序列被测定, 稻属成为进行比较、功能和进化基因组学研究的模式系统。近期开展的一系列研究对稻属不同基因组区段以及全基因组序列的比较分析, 揭示了稻属在基因组大小、基因移动、多倍体进化、常染色质到异染色质的转化以及着丝粒区域的进化等方面的分子机制。转座子的活性以及转座子因非均等重组或非法重组而造成的删除, 对稻属基因组的扩增和收缩具有重要作用。DNA双链断裂修复介导的基因移动, 特别是非同源末端连接, 是稻属基因组非共线性基因形成的主要来源。稻属基因组从常染色质到异染色质的转换过程, 伴随着转座子的大量扩增、基因片段的区段性和串联重复以及从基因组其他位置不断捕获异染色质基因。对稻属不同物种间基因拷贝数、特异基因和重要农艺性状基因的进化等研究, 可揭示稻属不同物种间表型和适应性差异的分子基础, 将加速水稻的育种和改良。     

Multiple complex families of endogenous retroviruses are highly conserved in the genus Gallus.

We have analyzed the genome of the domestic chicken for the presence of genetic sequences related to the envelope protein-encoding genes of avian sarcoma/leukosis retroviruses to determine the organization, structure, potential functionality, and distribution of such sequences. We have previously identified in the genus Gallus an extensive group of endogenous avian retroviruses termed EAV-0. Southern blot and sequence analysis presented here of EAV-0 elements revealed that the majority of the EAV-0 elements in the domestic chicken genome have large deletions in their env genes. Screening of a line 0 chicken genomic DNA library for potential full-length env gene-containing endogenous elements yielded three provirus clones of a previously unrecognized group of endogenous retroviruses. These three clones, E13, E33, and E51, are more closely related to each other (80% or more sequence identity) than to other avian retroviruses (70% or less sequence identity). The E13 element has a large deletion in env, but the E51 element has full-length and highly divergent SU- and TM-coding domains. Complete sequence analysis of the E51 env gene region revealed a defective SU-coding domain and an intact TM-coding domain. Sequence analysis of the E51, E33, and E13 3' termini revealed highly distinctive long terminal repeats of approximately 360 bp which appear to be the products, in part, of long terminal repeat domain shuffling. Hybridization analysis with E51 and E33 env gene probes indicated that they are members of an extensive group of elements present in all Gallus species, and at least one element, E51, could be shown by polymerase chain reaction amplification and direct sequencing to have integrated prior to Gallus speciation.     

The family of Deg/HtrA proteases in plants

ABSTRACT: BACKGROUND: The Deg/HtrA family of ATP-independent serine endopeptidases is present in nearly all organisms from bacteria to human and vascular plants. In recent years, multiple deg/htrA protease genes were identified in various plant genomes. During genome annotations most proteases were named according to the order of discovery, hence the same names were sometimes given to different types of Deg/HtrA enzymes in different plant species. This can easily lead to false inference of individual protease functions based solely on a shared name. Therefore, the existing names and classification of these proteolytic enzymes does not meet our current needs and a phylogeny-based standardized nomenclature is required. RESULTS: Using phylogenetic and domain arrangement analysis, we improved the nomenclature of the Deg/HtrA protease family, standardized protease names based on their well-established nomenclature in Arabidopsis thaliana, and clarified the evolutionary relationship between orthologous enzymes from various photosynthetic organisms across several divergent systematic groups, including dicots, a monocot, a moss and a green alga. Furthermore, we identified a "core set" of eight proteases shared by all organisms examined here that might provide all the proteolytic potential of Deg/HtrA proteases necessary for a hypothetical plant cell. CONCLUSIONS: In our proposed nomenclature, the evolutionarily closest orthologs have the same protease name, simplifying scientific communication when comparing different plant species and allowing for more reliable inference of protease functions. Further, we proposed that the high number of Deg/HtrA proteases in plants is mainly due to gene duplications unique to the respective organism.     

RS-SNP: a random-set method for genome-wide association studies

Background

Common carp is one of the most important aquaculture teleost fish in the world. Common carp and other closely related Cyprinidae species provide over 30% aquaculture production in the world. However, common carp genomic resources are still relatively underdeveloped. BAC end sequences (BES) are important resources for genome research on BAC-anchored genetic marker development, linkage map and physical map integration, and whole genome sequence assembling and scaffolding.

Result

To develop such valuable resources in common carp (Cyprinus carpio), a total of 40,224 BAC clones were sequenced on both ends, generating 65,720 clean BES with an average read length of 647 bp after sequence processing, representing 42,522,168 bp or 2.5% of common carp genome. The first survey of common carp genome was conducted with various bioinformatics tools. The common carp genome contains over 17.3% of repetitive elements with GC content of 36.8% and 518 transposon ORFs. To identify and develop BAC-anchored microsatellite markers, a total of 13,581 microsatellites were detected from 10,355 BES. The coding region of 7,127 genes were recognized from 9,443 BES on 7,453 BACs, with 1,990 BACs have genes on both ends. To evaluate the similarity to the genome of closely related zebrafish, BES of common carp were aligned against zebrafish genome. A total of 39,335 BES of common carp have conserved homologs on zebrafish genome which demonstrated the high similarity between zebrafish and common carp genomes, indicating the feasibility of comparative mapping between zebrafish and common carp once we have physical map of common carp.

Conclusion

BAC end sequences are great resources for the first genome wide survey of common carp. The repetitive DNA was estimated to be approximate 28% of common carp genome, indicating the higher complexity of the genome. Comparative analysis had mapped around 40,000 BES to zebrafish genome and established over 3,100 microsyntenies, covering over 50% of the zebrafish genome. BES of common carp are tremendous tools for comparative mapping between the two closely related species, zebrafish and common carp, which should facilitate both structural and functional genome analysis in common carp.     

A review of catfish genomics: progress and perspectives

Catfish is one of the lower teleosts whose genome research is important for evolutionary genomics. As the major aquaculture species in the USA, its genome research also has practical and economical implications. Much progress has been made in recent years, including the development of large numbers of molecular markers, the construction of framework genetic linkage maps, the identification of putative markers involved in performance traits, and the development of genomic resources. Repetitive elements have been identified and characterized in the catfish genome that should facilitate physical analysis of the catfish genome. A large number of genes or full-length cDNAs have been analysed using genomic approaches, providing information on gene structure, gene evolution and gene expression in relation to functions. Catfish genome research has come to a stage when physical mapping through BAC contig construction is greatly demanded, in order to develop regional markers for QTL analysis and for large-scale comparative mapping. The current effort in large-scale EST analysis and type I marker mapping should further enhance research efficiency through comparative mapping. Candidate gene identification is being accelerated through the use of cDNA microarrays.     

假基因的组成、分布及其分子进化

假基因(pseudogene)是指基因组中与正常基因序列相似, 但是缺乏功能的DNA 序列。通过序列同源性搜索, 可以收集基因组中假基因的群体特性、染色体分布和同源家族等特性。假基因很好地保留了数百万年前基因组中祖先基因的分子记录, 被视为“基因化石”, 因此假基因在进化和比较基因组学中是重要的资源。应用假基因和基因比较体系, 可以探究生物基因的进化史和基因组稳定性。如: 用Ka/Ks比值确定假基因的自然选择压、物种亲缘关系和进化距离, 分析假基因自身的进化趋势, 探讨DNA 突变的成因等。     

Construction and preliminary characterization of a river buffalo bacterial artificial chromosome library

River buffalo genome analyses have advanced significantly in the last decade, and the genome sequence of Bubalus bubalis will be available shortly. Nonetheless, large-insert DNA library resources such as bacterial artificial chromosomes (BAC) are still required for validation and accurate assembly of the genome sequence. We constructed a river buffalo BAC library containing 52,224 clones with an average insert size of 97 kb, representing 1.7 × coverage of the genome. This genomic resource for river buffalo will facilitate further studies in this economically important species allowing for instance, whole genome physical mapping and isolation of genes and gene clusters, contributing to the elucidation of gene organization and identification of regulatory elements.     

假基因的组成、分布及其分子进化

假基因（pseudogene）是指基因组中与正常基因序列相似,但是缺乏功能的DNA序列.通过序列同源性搜索,可以收集基因组中假基因的群体特性、染色体分布和同源家族等特性.假基因很好地保留了数百万年前基因组中祖先基因的分子记录,被视为＂基因化石＂,因此假基因在进化和比较基因组学中是重要的资源.应用假基因和基因比较体系,可以探究生物基因的进化史和基因组稳定性.如：用Ka/Ks比值确定假基因的自然选择压、物种亲缘关系和进化距离,分析假基因自身的进化趋势,探讨DNA突变的成因等.