Breeding and genetics of waterfowl

Comparative genomics of the waterfowl innate immune system

Animal species differ considerably in their ability to fight off infections. Finding the genetic basis of these differences is not easy, as the immune response is comprised of a complex network of proteins that interact with one another to defend the body against infection. Here, we used population- and comparative genomics to study the evolutionary forces acting on the innate immune system in natural hosts of the avian influenza virus (AIV). For this purpose, we used a combination of hybrid capture, next- generation sequencing and published genomes to examine genetic diversity, divergence, and signatures of selection in 127 innate immune genes at a micro- and macroevolutionary time scale in 26 species of waterfowl. We show across multiple immune pathways (AIV-, toll-like-, and RIG-I -like receptors signalling pathways) that genes involved genes in pathogen detection (i.e., toll-like receptors) and direct pathogen inhibition (i.e., antimicrobial peptides and interferon-stimulated genes), as well as host proteins targeted by viral antagonist proteins (i.e., mitochondrial antiviral-signaling protein, [MAVS]) are more likely to be polymorphic, genetically divergent, and under positive selection than other innate immune genes. Our results demonstrate that selective forces vary across innate immune signaling signalling pathways in waterfowl, and we present candidate genes that may contribute to differences in susceptibility and resistance to infectious diseases in wild birds, and that may be manipulated by viruses. Our findings improve our understanding of the interplay between host genetics and pathogens, and offer the opportunity for new insights into pathogenesis and potential drug targets.

Identification and profiling of microRNAs in the embryonic breast muscle of pekin duck

MicroRNAs (miRNAs) regulate gene expression by fully or partially binding to complementary sequences and play important roles in skeletal muscle development. However, the roles of miRNAs in embryonic breast muscle of duck are unclear. In this study, we analyzed the miRNAs profiling in embryonic breast muscle of Pekin duck at E13 (the 13(th) day of hatching), E19, and E27 by high-throughput sequencing. A total of 382 miRNAs including 359 preciously identified miRNAs 23 novel miRNA candidates were obtained. The nucleotide bias analysis of identified miRNAs showed that the miRNAs in Pekin duck was high conserved. The expression of identified miRNAs were significantly different between E13 and E19 as well as between E27 and E19. Fifteen identified miRNAs validated using stem-loop qRT-PCR can be divided into three groups: those with peak expression at E19, those with minimal expression at E19, and those with continuous increase from E11 to E27. Considering that E19 is the fastest growth stage of embryonic Pekin duck breast muscle, these three groups of miRNAs might be the potential promoters, the potential inhibitors, and the potential sustainer for breast muscle growth. Among the 23 novel miRNAs, novel-miRNA-8 and novel-miRNA-14 had maximal expression at some stages. The stem-loop qRT-PCR analysis of the two novel miRNAs and their two targets (MAP2K1 and PPARα) showed that the expression of novel-mir-8 and PPARα reached the lowest points at E19, while that of novel-mir-14 and MAP2K1 peaked at E19, suggesting novel-miRNA-8 and novel-miRNA-14 may be a potential inhibitor and a potential promoter for embryonic breast muscle development of duck. In summary, these results not only provided an overall insight into the miRNAs landscape in embryonic breast muscle of duck, but also a basis for the further investigation of the miRNAs roles in duck skeletal muscle development.

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.

Genetic diversity and population structure of 10 Chinese indigenous egg-type duck breeds assessed by microsatellite polymorphism

The genetic structure and diversity of 10 Chinese indigenous egg-type duck breeds were investigated using 29 microsatellite markers. The total number of animals examined were 569, on average 57 animals per breed were selected. The microsatellite marker set analysed provided 177 alleles (mean 6.1 alleles per locus, ranging from 3 to 10). All populations showed high levels of heterozygosity with the lowest estimate of 0.539 for the Jinding ducks, and the highest 0.609 observed for Jingjiang partridge ducks. The global heterozygote deficit across all populations (FIT) amounted to -0.363. About 10% of the total genetic variability originated from differences among breeds, with all loci contributing significantly. An unrooted consensus tree was constructed using the NeighborNet tree based on the Reynold's genetic distance. The structure software was used to assess genetic clustering of these egg-type duck breeds. Clustering analysis provided an accurate representation of the current genetic relations among the breeds. An integrated analysis was undertaken to obtain information on the population dynamics in Chinese indigenous egg-type duck breeds, and to better determine the conservation priorities.

Production benefits of the crossbreeding of indigenous and non-indigenous ducks--growing and laying period body weight and production performance

To evaluate different crosses and purebreds ducks in respect to various economic traits and to estimate different crossbreeding genetic parameters, a 3 x 3 complete diallel cross involving indigenous duck (DD), Khaki Campbell (KK) and White Pekin (WW) were used to produce three purebreds (DD, KK, WW) three crossbreds (DK, DW, KW) and three reciprocals (KD, WD, WK). A total of 609 ducklings produced were reared on deep litter and the females (316 in number) were evaluated for growing and laying period body weight along with the production performance traits. Different crossbreeding genetic parameters were estimated for different traits. All the traits in respect to body weight gain during growing and laying period and different production traits including laying house mortality rate showed significant (p <or= 0.05) difference between different genetic groups. In general, crossbreds perform better than the purebreds for most of the traits studied. General combining ability (GCA), specific combining ability (SCA) and reciprocal effect (RE) were significant (p <or= 0.01) for body weight and production traits. Egg weight showed significant (p <or= 0.01) difference in respect to GCA, SCA and RE for all the ages of measurement except RE for 30th week egg weight. Laying period mortality rate was only significant (p <or= 0.05) for SCA. Most of the crossbreds recorded heterosis rate in desirable direction for majority of the traits. Overall results revealed that the crossbreds perform well in respect to different traits than the purebreds and may be used to take advantage of heterosis. DW performs well in respect to majority of the traits measured and is of importance for commercial exploitation. Further, pure line selection with development of specialised sire and dam line followed by crossing may be of importance to enhance the performances in the crosses.

Eleven generations of selection for the duration of fertility in the intergeneric crossbreeding of ducks

A 12-generation selection experiment involving a selected line (S) and a control line (C) has been conducted since 1992 with the aim of increasing the number of fertile eggs laid by the Brown Tsaiya duck after a single artificial insemination (AI) with pooled Muscovy semen. On average, 28.9% of the females and 17.05% of the males were selected. The selection responses and the predicted responses showed similar trends. The average predicted genetic responses per generation in genetic standard deviation units were 0.40 for the number of fertile eggs, 0.45 for the maximum duration of fertility, and 0.32 for the number of hatched mule ducklings' traits. The fertility rates for days 2-8 after AI were 89.14% in the S line and 61.46% in the C line. Embryo viability was not impaired by this selection. The largest increase in fertility rate per day after a single AI was observed from d5 to d11. In G12, the fertility rate in the selected line was 91% at d2, 94% at d3, 92% at days 3 and 4 then decreased to 81% at d8, 75% at d9, 58% at d10 and 42% at d11. In contrast, the fertility rate in the control line showed an abrupt decrease from d4 (74%). The same tendencies were observed for the evolution of hatchability according to the egg set rates. It was concluded that selection for the number of fertile eggs after a single AI with pooled Muscovy semen could effectively increase the duration of the fertile period in ducks and that research should now be focused on ways to improve the viability of the hybrid mule duck embryo.

Duck (Anas platyrhynchos) linkage mapping by AFLP fingerprinting

Amplified fragment length polymorphism (AFLP) with multicolored fluorescent molecular markers was used to analyze duck (Anas platyrhynchos) genomic DNA and to construct the first AFLP genetic linkage map. These markers were developed and genotyped in 766 F2 individuals from six families from a cross between two different selected duck lines, brown Tsaiya and Pekin. Two hundred and ninety-six polymorphic bands (64% of all bands) were detected using 18 pairs of fluorescent TaqI/EcoRI primer combinations. Each primer set produced a range of 7 to 29 fragments in the reactions, and generated on average 16.4 polymorphic bands. The AFLP linkage map included 260 co-dominant markers distributed in 32 linkage groups. Twenty-one co-dominant markers were not linked with any other marker. Each linkage group contained three to 63 molecular markers and their size ranged between 19.0 cM and 171.9 cM. This AFLP linkage map provides important information for establishing a duck chromosome map, for mapping quantitative trait loci (QTL mapping) and for breeding applications.