Sequencing and genomic annotation of the chicken (Gallus gallus) Hox clusters, and mapping of evolutionarily conserved regions

Unraveling genomic associations with feed efficiency and body weight traits in chickens through an integrative approach

Background

Feed efficiency and growth rate have been targets for selection to improve chicken production. The incorporation of genomic tools may help to accelerate selection. We genotyped 529 individuals using a high-density SNP chip (600 K, Affymetrix®) to estimate genomic heritability of performance traits and to identify genomic regions and their positional candidate genes associated with performance traits in a Brazilian F₂ Chicken Resource population. Regions exhibiting selection signatures and a SNP dataset from resequencing were integrated with the genomic regions identified using the chip to refine the list of positional candidate genes and identify potential causative mutations.

Results

Feed intake (FI), feed conversion ratio (FC), feed efficiency (FE) and weight gain (WG) exhibited low genomic heritability values (i.e. from 0.0002 to 0.13), while body weight at hatch (BW1), 35 days-of-age (BW35), and 41 days-of-age (BW41) exhibited high genomic heritability values (i.e. from 0.60 to 0.73) in this F₂ population. Twenty unique 1-Mb genomic windows were associated with BW1, BW35 or BW41, located on GGA1–4, 6–7, 10, 14, 24, 27 and 28. Thirty-eight positional candidate genes were identified within these windows, and three of them overlapped with selection signature regions. Thirteen predicted deleterious and three high impact sequence SNPs in these QTL regions were annotated in 11 positional candidate genes related to osteogenesis, skeletal muscle development, growth, energy metabolism and lipid metabolism, which may be associated with body weight in chickens.

Conclusions

The use of a high-density SNP array to identify QTL which were integrated with whole genome sequence signatures of selection allowed the identification of candidate genes and candidate causal variants. One novel QTL was detected providing additional information to understand the genetic architecture of body weight traits. We identified QTL for body weight traits, which were also associated with fatness in the same population. Our findings form a basis for further functional studies to elucidate the role of specific genes in regulating body weight and fat deposition in chickens, generating useful information for poultry breeding programs.

Genome-wide profiling of chicken dendritic cell response to infectious bursal disease

BACKGROUND

Avian infectious bursal disease virus (IBDV) is a highly contagious, immunosuppressive disease of young chickens, which causes high mortality rates and large economic losses in the poultry industry. Dendritic cells (DCs), which are antigen-presenting cells, have the unique ability to induce both innate and acquired immune responses and may significantly influence virus pathogenicity. To understand the interaction between IBDV and DCs, a microarray was used to analyse the response of DCs infected by IBDV.

RESULTS

IBDV infection induced 479 upregulated and 466 downregulated mRNAs in chicken DCs. Analysis of Gene Ontology suggested that transcription from the RNA polymerase II promoter and the RNA biosynthetic process were enriched, and pathway analyses suggested that oxidative phosphorylation, as well as the T cell receptor and Interleukin-17 (IL-17) signalling pathways might be activated by IBDV infection. Moreover, microRNA (miRNA) and long non-coding RNA (lncRNA) alterations in IBDV-infected chicken DCs were observed. A total of 18 significantly upregulated or downregulated miRNAs and 441 significantly upregulated or downregulated lncRNAs were identified in IBDV-stimulated DCs. We constructed 42 transcription factor (TF)-miRNA-mRNA interactions involving 1 TF, 3 miRNAs, and 42 mRNAs in IBDV-stimulated DCs. Finally, we predicted the target genes of differentially expressed lncRNAs, and constructed lncRNA-mRNA regulatory networks.

CONCLUSIONS

The results of this study suggest a mechanism to explain how IBDV infection triggers an effective immune response in chicken DCs.

Comparative spatiotemporal analysis of Hox gene expression in early stages of intermediate mesoderm formation

BACKGROUND

Hox genes are key players in AP patterning of the vertebrate body plan and are necessary for organogenesis. Several studies provide evidence for the role Hox genes play during kidney development and especially regarding metanephros initiation and formation. However, the role Hox genes play during early stages of kidney development is largely unknown. A recent study in our lab revealed the role Hoxb4 plays in conferring the competence of intermediate mesodermal cells to respond to kidney inductive signals and express early kidney regulators.

RESULTS

As a first step in understanding the role Hox genes play in setting the formation of the pronephros morphogenetic field and the expression of early regulators of kidney development, we studied in detail the expression pattern of 10 Hox genes in relation to the 6th somite axial level, the anterior sharp border of the kidney field. Despite the idea of spatial co-linearity as exemplified in the Hox gene expression pattern in late developmental stages, a very dynamic spatio-temporal expression of these genes was found in early stages.

CONCLUSIONS

Since mesodermal patterning occurs at gastrula stages, the relevance of a "Hox code" at early stages is questioned in this study.

The crest phenotype in chicken is associated with ectopic expression of HOXC8 in cranial skin

The Crest phenotype is characterised by a tuft of elongated feathers atop the head. A similar phenotype is also seen in several wild bird species. Crest shows an autosomal incompletely dominant mode of inheritance and is associated with cerebral hernia. Here we show, using linkage analysis and genome-wide association, that Crest is located on the E22C19W28 linkage group and that it shows complete association to the HOXC-cluster on this chromosome. Expression analysis of tissues from Crested and non-crested chickens, representing 26 different breeds, revealed that HOXC8, but not HOXC12 or HOXC13, showed ectopic expression in cranial skin during embryonic development. We propose that Crest is caused by a cis-acting regulatory mutation underlying the ectopic expression of HOXC8. However, the identification of the causative mutation(s) has to await until a method becomes available for assembling this chromosomal region. Crest is unfortunately located in a genomic region that has so far defied all attempts to establish a contiguous sequence.

Hox clusters of the bichir (Actinopterygii, Polypterus senegalus) highlight unique patterns of sequence evolution in gnathostome phylogeny

Teleost fishes have extra Hox gene clusters owing to shared or lineage-specific genome duplication events in rayfinned fish (actinopterygian) phylogeny. Hence, extrapolating between genome function of teleosts and human or even between different fish species is difficult. We have sequenced and analyzed Hox gene clusters of the Senegal bichir (Polypterus senegalus), an extant representative of the most basal actinopterygian lineage. Bichir possesses four Hox gene clusters (A, B, C, D); phylogenetic analysis supports their orthology to the four Hox gene clusters of the gnathostome ancestor. We have generated a comprehensive database of conserved Hox noncoding sequences that include cartilaginous, lobe-finned, and ray-finned fishes (bichir and teleosts). Our analysis identified putative and known Hox cis-regulatory sequences with differing depths of conservation in Gnathostoma. We found that although bichir possesses four Hox gene clusters, its pattern of conservation of noncoding sequences is mosaic between outgroups, such as human, coelacanth, and shark, with four Hox gene clusters and teleosts, such as zebrafish and pufferfish, with seven or eight Hox gene clusters. Notably, bichir Hox gene clusters have been invaded by DNA transposons and this trend is further exemplified in teleosts, suggesting an as yet unrecognized mechanism of genome evolution that may explain Hox cluster plasticity in actinopterygians. Taken together, our results suggest that actinopterygian Hox gene clusters experienced a reduction in selective constraints that surprisingly predates the teleost-specific genome duplication.

A general scenario of Hox gene inventory variation among major sarcopterygian lineages

BACKGROUND

Hox genes are known to play a key role in shaping the body plan of metazoans. Evolutionary dynamics of these genes is therefore essential in explaining patterns of evolutionary diversity. Among extant sarcopterygians comprising both lobe-finned fishes and tetrapods, our knowledge of the Hox genes and clusters has largely been restricted in several model organisms such as frogs, birds and mammals. Some evolutionary gaps still exist, especially for those groups with derived body morphology or occupying key positions on the tree of life, hindering our understanding of how Hox gene inventory varied along the sarcopterygian lineage.

RESULTS

We determined the Hox gene inventory for six sarcopterygian groups: lungfishes, caecilians, salamanders, snakes, turtles and crocodiles by comprehensive PCR survey and genome walking. Variable Hox genes in each of the six sarcopterygian group representatives, compared to the human Hox gene inventory, were further validated for their presence/absence by PCR survey in a number of related species representing a broad evolutionary coverage of the group. Turtles, crocodiles, birds and placental mammals possess the same 39 Hox genes. HoxD12 is absent in snakes, amphibians and probably lungfishes. HoxB13 is lost in frogs and caecilians. Lobe-finned fishes, amphibians and squamate reptiles possess HoxC3. HoxC1 is only present in caecilians and lobe-finned fishes. Similar to coelacanths, lungfishes also possess HoxA14, which is only found in lobe-finned fishes to date. Our Hox gene variation data favor the lungfish-tetrapod, turtle-archosaur and frog-salamander relationships and imply that the loss of HoxD12 is not directly related to digit reduction.

CONCLUSIONS

Our newly determined Hox inventory data provide a more complete scenario for evolutionary dynamics of Hox genes along the sarcopterygian lineage. Limbless, worm-like caecilians and snakes possess similar Hox gene inventories to animals with less derived body morphology, suggesting changes to their body morphology are likely due to other modifications rather than changes to Hox gene numbers. Furthermore, our results provide basis for future sequencing of the entire Hox clusters of these animals.

The use and abuse of -omes

The diverse fields of Omics research share a common logical structure combining a cataloging effort for a particular class of molecules or interactions, the underlying -ome, and a quantitative aspect attempting to record spatiotemporal patterns of concentration, expression, or variation. Consequently, these fields also share a common set of difficulties and limitations. In spite of the great success stories of Omics projects over the last decade, much remains to be understood not only at the technological, but also at the conceptual level. Here, we focus on the dark corners of Omics research, where the problems, limitations, conceptual difficulties, and lack of knowledge are hidden.

Re-annotation is an essential step in systems biology modeling of functional genomics data

One motivation of systems biology research is to understand gene functions and interactions from functional genomics data such as that derived from microarrays. Up-to-date structural and functional annotations of genes are an essential foundation of systems biology modeling. We propose that the first essential step in any systems biology modeling of functional genomics data, especially for species with recently sequenced genomes, is gene structural and functional re-annotation. To demonstrate the impact of such re-annotation, we structurally and functionally re-annotated a microarray developed, and previously used, as a tool for disease research. We quantified the impact of this re-annotation on the array based on the total numbers of structural- and functional-annotations, the Gene Annotation Quality (GAQ) score, and canonical pathway coverage. We next quantified the impact of re-annotation on systems biology modeling using a previously published experiment that used this microarray. We show that re-annotation improves the quantity and quality of structural- and functional-annotations, allows a more comprehensive Gene Ontology based modeling, and improves pathway coverage for both the whole array and a differentially expressed mRNA subset. Our results also demonstrate that re-annotation can result in a different knowledge outcome derived from previous published research findings. We propose that, because of this, re-annotation should be considered to be an essential first step for deriving value from functional genomics data.

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome-deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)-(Nk4/tinman)-(Nk3/bagpipe)-(Lbx/ladybird)-(Tlx/c15)-(Nk7)-(Nk6/hgtx)-(Nk1/slouch)-(Nk5/Hmx). Of these genes, only NKX2.6-NKX3.1, LBX1-TLX1 and LBX2-TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have "missed" certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1-Tlx1, Lbx2-Tlx2 and Nkx3.1-Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2-Hmx3 genes are also tightly linked. Finally, we show that Lbx1-Tlx1 and Hmx2-Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2-Tlx2 and Nkx3.1-Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.

Atypical relaxation of structural constraints in Hox gene clusters of the green anole lizard

Hox genes control many aspects of embryonic development in metazoans. Previous analyses of this gene family revealed a surprising diversity in terms of gene number and organization between various animal species. In vertebrates, Hox genes are grouped into tightly organized clusters, claimed to be devoid of repetitive sequences. Here, we report the genomic organization of the four Hox loci present in the green anole lizard and show that they have massively accumulated retrotransposons, leading to gene clusters larger in size when compared to other vertebrates. In addition, similar repeats are present in many other development-related gene-containing regions, also thought to be refractory to such repetitive elements. Transposable elements are major sources of genetic variations, including alterations of gene expression, and hence this situation, so far unique among vertebrates, may have been associated with the evolution of the spectacular realm of morphological variations in the body plans of Squamata. Finally, sequence alignments highlight some divergent evolution in highly conserved DNA regions between vertebrate Hox clusters, which may coincide with the emergence of mammalian-specific features.

Establishment of Hox vertebral identities in the embryonic spine precursors

The vertebrate spine exhibits two striking characteristics. The first one is the periodic arrangement of its elements-the vertebrae-along the anteroposterior axis. This segmented organization is the result of somitogenesis, which takes place during organogenesis. The segmentation machinery involves a molecular oscillator-the segmentation clock-which delivers a periodic signal controlling somite production. During embryonic axis elongation, this signal is displaced posteriorly by a system of traveling signaling gradients-the wavefront-which depends on the Wnt, FGF, and retinoic acid pathways. The other characteristic feature of the spine is the subdivision of groups of vertebrae into anatomical domains, such as the cervical, thoracic, lumbar, sacral, and caudal regions. This axial regionalization is controlled by a set of transcription factors called Hox genes. Hox genes exhibit nested expression domains in the somites which reflect their linear arrangement along the chromosomes-a property termed colinearity. The colinear disposition of Hox genes expression domains provides a blueprint for the regionalization of the future vertebral territories of the spine. In amniotes, Hox genes are activated in the somite precursors of the epiblast in a temporal colinear sequence and they were proposed to control their progressive ingression into the nascent paraxial mesoderm. Consequently, the positioning of the expression domains of Hox genes along the anteroposterior axis is largely controlled by the timing of Hox activation during gastrulation. Positioning of the somitic Hox domains is subsequently refined through a crosstalk with the segmentation machinery in the presomitic mesoderm. In this review, we focus on our current understanding of the embryonic mechanisms that establish vertebral identities during vertebrate development.

Unexpectedly large number of conserved noncoding regions within the ancestral chordate Hox cluster

The single amphioxus Hox cluster contains 15 genes and may well resemble the ancestral chordate Hox cluster. We have sequenced the Hox genomic complement of the European amphioxus Branchiostoma lanceolatum and compared it to the American species, Branchiostoma floridae, by phylogenetic footprinting to gain insights into the evolution of Hox gene regulation in chordates. We found that Hox intergenic regions are largely conserved between the two amphioxus species, especially in the case of genes located at the 3' of the cluster, a trend previously observed in vertebrates. We further compared the amphioxus Hox cluster with the human HoxA, HoxB, HoxC, and HoxD clusters, finding several conserved noncoding regions, both in intergenic and intronic regions. This suggests that the regulation of Hox genes is highly conserved across chordates, consistent with the similar Hox expression patterns in vertebrates and amphioxus.