Yanqiang Wang

Genome-wide association study of body weight in chicken F2 resource population.

Chicken body weight is an economically important trait and great genetic progress has been accomplished in genetic selective for body weight. To identify genes and chromosome regions associated with body weight, we performed a genome-wide association study using the chicken 60 k SNP panel in a chicken F2 resource population derived from the cross between Silky Fowl and White Plymouth Rock. A total of 26 SNP effects involving 9 different SNP markers reached 5% Bonferroni genome-wide significance. A chicken chromosome 4 (GGA4) region approximately 8.6 Mb in length (71.6–80.2 Mb) had a large number of significant SNP effects for late growth during weeks 7–12. The LIM domain-binding factor 2 (LDB2) gene in this region had the strongest association with body weight for weeks 7–12 and with average daily gain for weeks 6–12. This GGA4 region was previously reported to contain body weight QTL. GGA1 and GGA18 had three SNP effects on body weight with genome-wide significance. Some of the SNP effects with the significance of “suggestive linkage” overlapped with previously reported results.

The Rose-comb Mutation in Chickens Constitutes a Structural Rearrangement Causing Both Altered Comb Morphology and Defective Sperm Motility

Rose-comb, a classical monogenic trait of chickens, is characterized by a drastically altered comb morphology compared to the single-combed wild-type. Here we show that Rose-comb is caused by a 7.4 Mb inversion on chromosome 7 and that a second Rose-comb allele arose by unequal crossing over between a Rose-comb and wild-type chromosome. The comb phenotype is caused by the relocalization of the MNR2 homeodomain protein gene leading to transient ectopic expression of MNR2 during comb development. We also provide a molecular explanation for the first example of epistatic interaction reported by Bateson and Punnett 104 years ago, namely that walnut-comb is caused by the combined effects of the Rose-comb and Pea-comb alleles. Transient ectopic expression of MNR2 and SOX5 (causing the Pea-comb phenotype) occurs in the same population of mesenchymal cells and with at least partially overlapping expression in individual cells in the comb primordium. Rose-comb has pleiotropic effects, as homozygosity in males has been associated with poor sperm motility. We postulate that this is caused by the disruption of the CCDC108 gene located at one of the inversion breakpoints. CCDC108 is a poorly characterized protein, but it contains a MSP (major sperm protein) domain and is expressed in testis. The study illustrates several characteristic features of the genetic diversity present in domestic animals, including the evolution of alleles by two or more consecutive mutations and the fact that structural changes have contributed to fast phenotypic evolution.

A genome-wide survey of copy number variation regions in various chicken breeds by array comparative genomic hybridization method.

The discovery of copy number variation (CNV) in the genome has provided new insight into genomic polymorphism. Studies with chickens have identified a number of large CNV segments using a 385k comparative genomic hybridization (CGH) chip (mean length >140 kb). We present a detailed CNV map for local Chinese chicken breeds and commercial chicken lines using an Agilent 400k array CGH platform with custom-designed probes. We identified a total of 130 copy number variation regions (CNVRs; mean length = 25.70 kb). Of these, 104 (80.0%) were novel segments reported for the first time in chickens. Among the 104 novel CNVRs, 56 (53.8%) of the segments were non-coding sequences, 65 (62.5%) showed the gain of DNA and 40 (38.5%) showed the loss of DNA (one locus showed both loss and gain). Overlapping with the formal selective sweep data and the quantitative trait loci data, we identified four loci that might be considered to be high-confidence selective segments that arose during the domestication of chickens. Compared with the CNVRs reported previously, genes for the positive regulation of phospholipase A2 activity were discovered to be significantly over-represented in the novel CNVRs reported here by gene ontology analysis. Availability of our results should facilitate further research in the study of the genetic variability in chicken breeds.

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.

A cis-regulatory mutation of PDSS2 causes silky-feather in chickens.

Silky-feather has been selected and fixed in some breeds due to its unique appearance. This phenotype is caused by a single recessive gene (hookless, h). Here we map the silky-feather locus to chromosome 3 by linkage analysis and subsequently fine-map it to an 18.9 kb interval using the identical by descent (IBD) method. Further analysis reveals that a C to G transversion located upstream of the prenyl (decaprenyl) diphosphate synthase, subunit 2 (PDSS2) gene is causing silky-feather. All silky-feather birds are homozygous for the G allele. The silky-feather mutation significantly decreases the expression of PDSS2 during feather development in vivo. Consistent with the regulatory effect, the C to G transversion is shown to remarkably reduce PDSS2 promoter activity in vitro. We report a new example of feather structure variation associated with a spontaneous mutation and provide new insight into the PDSS2 function.

Evaluation of SNPs in the chicken HMGA2 gene as markers for body weight gain

A QTL affecting body weight in chickens has been mapped to GGA1, between the markers GCT0006 and MCW0106. The gene HMGA2, which was previously identified as a candidate gene for determining body height in humans and mice, is also conspicuously close to the MCW0106 marker in chickens. Subsequently, 14 SNP markers of HMGA2 were genotyped in CAU chicken resource populations, and the associations between body weight and those SNP markers that displayed polymorphisms were analysed. Three SNPs (rs13849241, rs15231472 and rs13849381) were found to be significantly correlated with body weight in chickens (P < 0.05). Furthermore, haplotypes constructed based on these three SNPs were also discovered to be associated with body weight in chickens at the ages of 6, 7, 9 and 12 weeks. These results suggest that the chicken HMGA2 gene is indeed involved in body weight gain.

A Complex Structural Variation on Chromosome 27 Leads to the Ectopic Expression of HOXB8 and the Muffs and Beard Phenotype in Chickens.

Muffs and beard (Mb) is a phenotype in chickens where groups of elongated feathers gather from both sides of the face (muffs) and below the beak (beard). It is an autosomal, incomplete dominant phenotype encoded by the Muffs and beard (Mb) locus. Here we use genome-wide association (GWA) analysis, linkage analysis, Identity-by-Descent (IBD) mapping, array-CGH, genome re-sequencing and expression analysis to show that the Mb allele causing the Mb phenotype is a derived allele where a complex structural variation (SV) on GGA27 leads to an altered expression of the gene HOXB8. This Mb allele was shown to be completely associated with the Mb phenotype in nine other independent Mb chicken breeds. The Mb allele differs from the wild-type mb allele by three duplications, one in tandem and two that are translocated to that of the tandem repeat around 1.70 Mb on GGA27. The duplications contain total seven annotated genes and their expression was tested during distinct stages of Mb morphogenesis. A continuous high ectopic expression of HOXB8 was found in the facial skin of Mb chickens, strongly suggesting that HOXB8 directs this regional feather-development. In conclusion, our results provide an interesting example of how genomic structural rearrangements alter the regulation of genes leading to novel phenotypes. Further, it again illustrates the value of utilizing derived phenotypes in domestic animals to dissect the genetic basis of developmental traits, herein providing novel insights into the likely role of HOXB8 in feather development and differentiation.

Genomic regions associated with the sex-linked inhibitor of dermal melanin in Silkie chicken

A unique characteristic of the Silkie chicken is its fibromelanosis phenotype. The dermal layer of its skin, its connective tissue and shank dermis are hyperpigmen- ted. This dermal hyperpigmentation phenotype is con- trolled by the sex-linked inhibitor of dermal melanin gene (ID) and the dominant fibromelanosis allele. This study attempted to confirm the genomic region associated with ID. By genotyping, ID was found to be closely linked to the region between GGA_rs16127903 and GGA_rs14685542 (8406919 bp) on chromosome Z, which contains ten functional genes. The expression of these genes was characterized in the embryo and 4 days after hatching and it was concluded that MTAP, encoding methylthioadenosine- phosphorylase, would be the most likely candidate gene. Finally, target DNA capture and sequence analysis was performed, but no specific SNP(s) was found in the targeted region of the Silkie genome. Further work is necessary to identify the causal ID mutation located on chromosome Z.