Construction and characterization of a duck bacterial artificial chromosome library

A duck RH panel and its potential for assisting NGS genome assembly

Background

Owing to the low cost of the high throughput Next Generation Sequencing (NGS) technology, more and more species have been and will be sequenced. However, de novo assemblies of large eukaryotic genomes thus produced are composed of a large number of contigs and scaffolds of medium to small size, having no chromosomal assignment. Radiation hybrid (RH) mapping is a powerful tool for building whole genome maps and has been used for several animal species, to help assign sequence scaffolds to chromosomes and determining their order.

Results

We report here a duck whole genome RH panel obtained by fusing female duck embryonic fibroblasts irradiated at a dose of 6,000 rads, with HPRT-deficient Wg3hCl₂ hamster cells. The ninety best hybrids, having an average retention of 23.6% of the duck genome, were selected for the final panel. To allow the genotyping of large numbers of markers, as required for whole genome mapping, without having to cultivate the hybrid clones on a large scale, three different methods involving Whole Genome Amplification (WGA) and/or scaling down PCR volumes by using the Fluidigm BioMark^TM Integrated Fluidic Circuits (IFC) Dynamic Array^TM for genotyping were tested. RH maps of APL12 and APL22 were built, allowing the detection of intrachromosomal rearrangements when compared to chicken. Finally, the panel proved useful for checking the assembly of sequence scaffolds and for mapping EST located on one of the smallest microchromosomes.

Conclusion

The Fluidigm BioMark^TM Integrated Fluidic Circuits (IFC) Dynamic Array^TM genotyping by quantitative PCR provides a rapid and cost-effective method for building RH linkage groups. Although the vast majority of genotyped markers exhibited a picture coherent with their associated scaffolds, a few of them were discordant, pinpointing potential assembly errors. Comparative mapping with chicken chromosomes GGA21 and GGA11 allowed the detection of the first chromosome rearrangements on microchromosomes between duck and chicken. As in chicken, the smallest duck microchromosomes appear missing in the assembly and more EST data will be needed for mapping them. Altogether, this underlines the added value of RH mapping to improve genome assemblies.

Avian influenza rapidly induces antiviral genes in duck lung and intestine

Ducks are the natural reservoir of influenza A and survive infection by most strains. To characterize the duck immune response to influenza, we sought to identify innate immune genes expressed early in an infection. We used suppressive subtractive hybridization (SSH) to construct 3 libraries enriched in differentially expressed genes from lung RNA of a duck infected with highly pathogenic avian influenza virus A/Vietnam/1203/04 (H5N1), or lung and intestine RNA of a duck infected with low pathogenic avian influenza A/mallard/BC/500/05 (H5N2) compared to a mock-infected duck. Sequencing of 1687 clones identified a transcription profile enriched in genes involved in antiviral defense and other cellular processes. Major histocompatibility complex class I (MHC I), interferon induced protein with tricopeptide repeats 5 (IFIT5), and 2'-5' oligoadenylate synthetase-like gene (OASL) were increased more than 1000-fold in relative transcript abundance in duck lung at 1dpi with highly pathogenic VN1203. These genes were induced much less in lung or intestine following infection with low pathogenic BC500. The expression of these genes following infection suggests that ducks initiate an immediate and robust response to a potentially lethal influenza strain, and a minimal response to a low pathogenic strain.

Characterization of the MHC class II α-chain gene in ducks

In humans, classical MHC class II molecules include DQ, DR, and DP, which are similar in structure but consist of distinct α- and β-chains. The genes encoding these molecules are all located in the MHC class II gene region. In non-mammalian vertebrates such as chickens, only a single class II α-chain gene corresponding to the human DRA has been identified. Here, we report a characterization of the duck MHC class II α-chain (Anpl-DRA) encoding gene, which contains four exons encoding a typical signal peptide, a peptide-binding α1 domain, an immunoglobulin-like α2 domain, and Tm/Cyt, respectively. This gene is present in the duck genome as a single copy and is highly expressed in the spleen. Sequencing of cDNA and genomic DNA of the Anpl-DRA of different duck individuals/strains revealed low levels of genetic polymorphism, especially in the same strain, although most duck individuals have two different alleles. Otherwise, we found that the duck gene is located next to class II β genes, which is the same as in humans but different from the situation in chickens.

Molecular evolution of the vertebrate TLR1 gene family--a complex history of gene duplication, gene conversion, positive selection and co-evolution

Background

The Toll-like receptors represent a large superfamily of type I transmembrane glycoproteins, some common to a wide range of species and others are more restricted in their distribution. Most members of the Toll-like receptor superfamily have few paralogues; the exception is the TLR1 gene family with four closely related genes in mammals TLR1, TLR2, TLR6 and TLR10, and four in birds TLR1A, TLR1B, TLR2A and TLR2B. These genes were previously thought to have arisen by a series of independent gene duplications. To understand the evolutionary pattern of the TLR1 gene family in vertebrates further, we cloned the sequences of TLR1A, TLR1B, TLR2A and TLR2B in duck and turkey, constructed phylogenetic trees, predicted codons under positive selection and identified co-evolutionary amino acid pairs within the TLR1 gene family using sequences from 4 birds, 28 mammals, an amphibian and a fish.

Results

This detailed phylogenetic analysis not only clarifies the gene gains and losses within the TLR1 gene family of birds and mammals, but also defines orthologues between these vertebrates. In mammals, we predict amino acid sites under positive selection in TLR1, TLR2 and TLR6 but not TLR10. We detect co-evolution between amino acid residues in TLR2 and the other members of this gene family predicted to maintain their ability to form functional heterodimers. In birds, we predict positive selection in the TLR2A and TLR2B genes at functionally significant amino acid residues. We demonstrate that the TLR1 gene family has mostly been subject to purifying selection but has also responded to directional selection at a few sites, possibly in response to pathogen challenge.

Conclusions

Our phylogenetic and structural analyses of the vertebrate TLR1 family have clarified their evolutionary origins and predict amino acid residues likely to be important in the host's defense against invading pathogens.

Avian genomics in the 21st century

The chicken has long been an important model organism for developmental biology, as well as a major source of protein with billions of birds used in meat and egg production each year. Chicken genomics has been transformed in recent years, with the characterisation of large EST collections and most recently with the assembly of the chicken genome sequence. As the first livestock genome to be fully sequenced it leads the way for others to follow--with zebra finch later this year. The genome sequence and the availability of three million genetic polymorphisms are expected to aid the identification of genes that control traits of importance in poultry. As the first bird genome to be sequenced it is a model for the remaining 9,600 species thought to exist today. Many of the features of avian biology and organisation of the chicken genome make it an ideal model organism for phylogenetics and embryology, along with applications in agriculture and medicine. The availability of new tools such as whole-genome gene expression arrays and SNP panels, coupled with information resources on the genes and proteins are likely to enhance this position.