Genomic scan revealed KIT gene underlying white/gray plumage color in Chinese domestic geese

RNA sequencing analysis reveals key genes and pathways associated with feather pigmentation in mule ducks

Feather color is an important morphological trait of poultry. At present, the reports on the inheritance of plumage color of mule ducks at the molecular level are few, and the regulatory mechanism in white plumage rates of different mule ducks remains unclear. This study aimed to broaden the understanding of the white plumage rates in mule ducks to improve their production value. We used RNA sequencing to analyze and compare the mRNA expression profiles in hair follicle tissues from 10-week-old mule ducks with black and white plumages, thereby revealing the temporal gene expression patterns and pathways associated with plumage color regulation. In total, 1672 annotated differentially expressed genes (DEGs) were identified in black and white plumages from different databases between mule ducks with the 2 plumage colors. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment in hair follicle tissues indicated that the aforementioned DEGs were mainly involved in the melanin signaling pathway. Concurrently, we use weighted gene co-expression network analysis to detect core modules and hub genes associated with melanin biosynthesis in feathers. The green module exhibited the strongest correlation with the phenotypic traits, encompassing a total of 1049 genes. Subsequent Kyoto Encyclopedia of Genes and Genomes enrichment analysis identified 11 genes as pivotal in the melanin biosynthetic pathway, including EDRNB2, TYR, KIT, EDNRB, and MC1R. The differential expression of eight selected DEGs was verified using quantitative reverse transcription-PCR, and the results were consistent with RNA-seq data. This study provides a basis for understanding the differences in plumage color development in mule ducks.

A novel codominant plumage color pattern of white breast patches in WugangTong geese was controlled by EDNRB2

Two basic plumage color patterns are observed in adult geese: solid grey (G) or colorless white (W). However, a Chinese indigenous breed, the Wugangtong goose (WGT), continues to be subject to selective breeding efforts as it displays segregation of plumage colors, including G, W, and a novel color pattern designated Wb (G with white breast circles). The underlying genetic mechanisms responsible for the Wb phenotype are yet to be determined. The current study employed the population differentiation index (FST) to analyze 90 geese exhibiting diverse plumage colors, identifying the fifth intron of EDNRB2 as a particularly noteworthy region with the highest FST values. Sanger sequencing of the region surrounding the EDNRB2 gene identified a 14-bp insertion within exon 3 as the causal mutation. The heterozygosity of this 14-bp insertion and wild-type alleles was completely associated with the Wb phenotype, thereby substantiating the codominant nature of the G and W phenotypes. An inter-species corroborated this finding cross between the graylag (no 14-bp insertion) and the swan goose (homozygous for the 14-bp insertion) breeds, as hybrids from this cross exhibited the Wb phenotype. Transcriptomes from white breast patches and gray dorsal skins of 4 Wb geese were compared. A significant downregulation of genes involved in melanin synthesis and melanocyte development was observed, including EDRNB2 and MLANA. The downregulation of MLANA indicated that the mutated EDNRB2 resulted in melanocyte loss in specific body regions, as MLANA is a marker gene for melanocytes. The findings were corroborated by melanin staining using the Mansson-Fontana method, which revealed no melanin particles deposited in the white breast patches. In summary, the gray plumage color was codominant to the white color in WGT geese, and plumage color variations were controlled by EDNRB2. The findings of our study offer valuable and practical guidance for the purification of plumage colors among WGT, whether through traditional phenotype selection or molecular breeding methods.

Genome-wide association study reveals 2 copy number variations associated with the variation of plumage color in the white duck hybrid population

Plumage color is an intuitive external poultry characteristic with rich manifestations and complex genetic mechanisms. In our previous study, we observed that there were more dark variations in plumage color in the F2 population derived from the hybridization of 2 white duck varieties. Therefore, based on the statistics of plumage color of 308 F2 populations, we further used the resequencing data of these individuals to detect copy number variations (CNVs) in the whole genome and conducted genome-wide association studies (GWAS) to determine the genetic basis related to plumage color traits. The CNV detection revealed 9,337 CNVs, with an average length of 15,950 bp and a total length of 142.02 MB, accounting for approximately 12.91% of the reference genome. The CNV distribution on the chromosomes was relatively uniform, and the number of CNVs on each chromosome positively correlated with the length of the chromosome. In the pure black plumage group, 2,101 CNVs were only identified, and 1,714 were specifically identified in the pure white plumage group. Ten CNVs were randomly selected for validation using quantitative real-time PCR, and 9 CNVs had the same CNV types as predicted, with an accuracy of 90%. Based on GWAS, we identified 2 CNVs potentially associated with plumage color variations, with the associated CNV regions covering 9 genes. Enrichment analysis of these 9 candidate genes showed significant enrichment of 3 pathways (ribosome biogenesis in eukaryotes, RNA transport, and protein export) and 17 gene ontology terms. Among these, VWA5A can downregulate MITF by binding to the regulatory factors SOX10. The occurrence of CNV may indirectly contribute to duck plumage color variation by affecting the regulatory factors of the switch gene MITF in the melanogenesis pathway. These findings have improved the understanding of the genetic basis of duck plumage color variation and have been beneficial for developing and using plumage color traits in subsequent poultry breeding.

Identification of Differentially Expressed Genes and microRNAs in the Gray and White Feather Follicles of Shitou Geese

The Shitou goose, a highly recognized indigenous breed with gray plumage originating from Chaozhou Raoping in Guangdong Province, China, is renowned for being the largest goose species in the country. Notably, during the pure breeding process of Shitou geese, approximately 2% of the offspring in each generation unexpectedly exhibited white plumage. To better understand the mechanisms underlying white plumage color formation in Shitou geese, we conducted a comparative transcriptome analysis between white and gray feather follicles, aiming to identify key genes and microRNAs that potentially regulate white plumage coloration in this unique goose breed. Our results revealed a number of pigmentation genes, encompassing TYR, TYRP1, EDNRB2, MLANA, SOX10, SLC45A2, GPR143, TRPM1, OCA2, ASIP, KIT, and SLC24A5, which were significantly down-regulated in the white feather follicles of Shitou geese. Among these genes, EDNRB2 and KIT emerged as the most promising candidate genes for white plumage coloration in Shitou geese. Additionally, our analysis also uncovered 46 differentially expressed miRNAs. Of these, miR-144-y may play crucial roles in the regulation of feather pigmentation. Furthermore, the expression of novel-m0086-5p, miR-489-y, miR-223-x, miR-7565-z, and miR-3535-z exhibits a significant negative correlation with the expression of pigmentation genes including TYRP1, EDNRB2, MLANA, SOX10, TRPM1, and KIT, suggesting these miRNAs may indirectly regulate the expression of these genes, thereby influencing feather color. Our findings provide valuable insights into the genetic mechanisms underlying white plumage coloration in Shitou geese and contribute to the broader understanding of avian genetics and coloration research.

Genome-wide association analysis of neck ring traits in NongHua ma male ducks

1. Male NongHua ma ducks have more colourful feathers than females, especially considering that the former have a distinctive neck ring that is different from that of females. This ring development might be influenced by sex selection, the environment, genetics and other elements.2. Genome-wide association analysis (GWAS) was used to locate candidate genes that affect the neck ring formation of male ducks to investigate the genetic basis of this phenomenon.3. In this study, the neck ring area and width of 180 male ducks were assessed at ages 80, 90, 100, 110 and 120 d. GWAS was used to identify associated genes. There were 0, 7, 14, 48 and 21 possible candidate genes annotated around the 0, 12, 25, 76 and 40 SNP loci n corresponding regions. A total of 13 candidate genes were identified around 21 SNP sites at the neck ring width of 120 d.4. These significant genes were annotated and GO and KEGG enrichment analyses were performed. All SNPs that exceeded the significance threshold were annotated and preliminarily screened as candidate genes affecting neck ring formation. From analysis of gene function and enriched KEGG pathways, genes such as THSD1, SLC6A4, DGAT2, PRKDC, B3GAT2, ROR1, GRK7, EXTL3, TXNDC12, COL4A2, PRKG1, ACTR3, were considered important candidate marker sites related to the neck ring. This provided a reference starting point for the genetic mechanism underlying duck feather colour.

Insights into genetic diversity and phenotypic variations in domestic geese through comprehensive population and pan-genome analysis

BACKGROUND

Domestic goose breeds are descended from either the Swan goose (Anser cygnoides) or the Greylag goose (Anser anser), exhibiting variations in body size, reproductive performance, egg production, feather color, and other phenotypic traits. Constructing a pan-genome facilitates a thorough identification of genetic variations, thereby deepening our comprehension of the molecular mechanisms underlying genetic diversity and phenotypic variability.

RESULTS

To comprehensively facilitate population genomic and pan-genomic analyses in geese, we embarked on the task of 659 geese whole genome resequencing data and compiling a database of 155 RNA-seq samples. By constructing the pan-genome for geese, we generated non-reference contigs totaling 612 Mb, unveiling a collection of 2,813 novel genes and pinpointing 15,567 core genes, 1,324 softcore genes, 2,734 shell genes, and 878 cloud genes in goose genomes. Furthermore, we detected an 81.97 Mb genomic region showing signs of genome selection, encompassing the TGFBR2 gene correlated with variations in body weight among geese. Genome-wide association studies utilizing single nucleotide polymorphisms (SNPs) and presence-absence variation revealed significant genomic associations with various goose meat quality, reproductive, and body composition traits. For instance, a gene encoding the SVEP1 protein was linked to carcass oblique length, and a distinct gene-CDS haplotype of the SVEP1 gene exhibited an association with carcass oblique length. Notably, the pan-genome analysis revealed enrichment of variable genes in the "hair follicle maturation" Gene Ontology term, potentially linked to the selection of feather-related traits in geese. A gene presence-absence variation analysis suggested a reduced frequency of genes associated with "regulation of heart contraction" in domesticated geese compared to their wild counterparts. Our study provided novel insights into gene expression features and functions by integrating gene expression patterns across multiple organs and tissues in geese and analyzing population variation.

CONCLUSION

This accomplishment originates from the discernment of a multitude of selection signals and candidate genes associated with a wide array of traits, thereby markedly enhancing our understanding of the processes underlying domestication and breeding in geese. Moreover, assembling the pan-genome for geese has yielded a comprehensive apprehension of the goose genome, establishing it as an indispensable asset poised to offer innovative viewpoints and make substantial contributions to future geese breeding initiatives.

Whole-genome sequencing identifies potential candidate genes for egg production traits in laying ducks (Anas platyrhynchos)

Egg production traits are economically important in laying ducks. Genetic molecular mechanisms and candidate genes underlying these traits remain unclear. In this study, whole genome variants were identified through whole-genome resequencing using three high-egg producing (HEN) and three low-egg producing (LEN) laying ducks. The gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genome (KEGG) pathways for the genes of common differential variants between HEN and LEN ducks were determined. Frizzled class receptor 6 (FZD6) was further genotyped using the Sequenom MassARRAY iPLEX platform. The association of FZD6 gene polymorphisms with 73 egg production and weight traits in 329 female ducks were estimated. A total of 65,535 single nucleotide polymorphisms (SNPs) and 4,702 indels were identified across the genome. Fourteen GO terms and 14 KEGG pathways were determined for the genes of common differential variants, including MAPK signaling, Wnt signaling, melanogenesis and calcium signaling pathways, which are key functional pathways for poultry egg production reported in previous reports. Further analysis showed that 27 SNPs of FZD6 were associated with three early egg production of duck and egg weight traits, including egg production at 17 weeks (EP17), 18 weeks (EP18) and 19 weeks (EP19) and egg weight at 59 weeks (EW59). The FZD6 should be considered a novel candidate gene for egg production traits in laying ducks.

Chromosome-level genome and population genomics reveal evolutionary characteristics and conservation status of Chinese indigenous geese

Geese are herbivorous birds that play an essential role in the agricultural economy. We construct the chromosome-level genome of a Chinese indigenous goose (the Xingguo gray goose, XGG; Anser cygnoides) and analyze the adaptation of fat storage capacity in the goose liver during the evolution of Anatidae. Genomic resequencing of 994 geese is used to investigate the genetic relationships of geese, which supports the dual origin of geese (Anser cygnoides and Anser anser). Chinese indigenous geese show higher genetic diversity than European geese, and a scientific conservation program can be established to preserve genetic variation for each breed. We also find that a 14-bp insertion in endothelin receptor B subtype 2 (EDNRB2) that determines the white plumage of Chinese domestic geese is a natural mutation, and the linkaged alleles rapidly increase in frequency as a result of genetic hitchhiking, leading to the formation of completely different haplotypes of white geese under strong artificial selection. These genomic resources and our findings will facilitate marker-assisted breeding of geese and provide a foundation for further research on geese genetics and evolution.

Using comparative genomics to detect mutations regulating plumage variations in graylag (A. anser) and swan geese (A. cygnoides)

Although graylag geese (A. anser) showed similar plumages of white, grey, and white with grey patches compared to those in swan geese (A. cygnoides), it was believed the substantial molecular mechanism for plumage variations were different. To date, studies on genes responsible for diverse plumages among graylag geese were limited and causal mutations remain unknown. In this study, genomes from 57 individuals belonging to six breeds showing different plumages were sequenced at ∼10X depth. Firstly, the allele frequency differences (AFD) of variants on the scaffold394 (NW_013185915.1) between grey and white goose breeds (A. anser) was calculated and a genomic region between 768,290-779,889 bp was detected to carry candidate variants associated with plumages, including one SNP (g. 775,151G > T, ∼18.6 kb upstream of EDNRB2) found to be fixed in white geese. This region was overlapped with the one detected by the haplotype-based sweep analysis, in which significant signals defined a candidate region of 736,610-820,622 bp on the same scaffold. Results from the transcriptomic data showed that expression levels of EDNRB2 and many other melanogenesis-related genes were significantly decreased among white geese compared to that in grey geese, especially at late embryonic stages (>E15). Modifications at transcriptional levels might result in abnormal melanocyte developments and thus the white plumages when they grow up. In addition, a frameshift mutation (C > -) in exon4 of MLANA gene on scaffold176 (NW_013185876.1) was suggested as the causal mutation for sex-linked dilution phenotype in graylag geese although this requires more demonstration experiments. Together with observed white plumages caused by EDNRB2 mutations in coding regions among swan geese and chicken, our study provided new examples to study the parallel evolution.

Whole Genome Sequencing Provides New Insights Into the Genetic Diversity and Coat Color of Asiatic Wild Ass and Its Hybrids

The diversity of livestock coat color results from human positive selection and is an indispensable part of breed registration. As an important biodiversity resource, Asiatic wild ass has many special characteristics, including the most visualized feature, its yellowish-brown coat color, and excellent adaptation. To explore the genetic mechanisms of phenotypic characteristics in Asiatic wild ass and its hybrids, we resequenced the whole genome of one Mongolian Kulan (a subspecies of Asiatic wild ass) and 29 Kulan hybrids (Mongolian Kulan ♂×Xinjiang♀), and the ancestor composition indicated the true lineage of the hybrids. XP-EHH (Cross Population Extended Haplotype Homozygosity), θπ-ratio (Nucleotide Diversity Ratio), CLR (Composite Likelihood Ratio) and θπ (Nucleotide Diversity) methods were used to detect the candidate regions of positive selection in Asiatic wild ass and its hybrids. Several immune genes (DEFA1, DEFA5, DEFA7, GIMAP4, GIMAP1, IGLC1, IGLL5, GZMB and HLA) were observed by the CLR and θπ methods. XP-EHH and θπ-ratio revealed that these genes are potentially responsible for coat color (KITLG) and meat quality traits (PDE1B and MYLK2). Furthermore, the heatmap was able to show the clear difference in the haplotype of the KITLG gene between the Kulan hybrids and Asiatic wild ass group and the Guanzhong black donkey group, which is a powerful demonstration of the key role of KITLG in donkey color. Therefore, our study may provide new insights into the genetic basis of coat color, meat quality traits and immunity of Asiatic wild ass and its hybrids.

Pooled Sequencing Analysis of Geese (Anser cygnoides) Reveals Genomic Variations Associated With Feather Color

During the domestication of the goose a change in its feather color took place, however, the molecular mechanisms responsible for this change are not completely understood. Here, we performed whole-genome resequencing on three pooled samples of geese (feral and domestic geese), with two distinct feather colors, to identify genes that might regulate feather color. We identified around 8 million SNPs within each of the three pools and validated allele frequencies for a subset of these SNPs using PCR and Sanger sequencing. Several genomic regions with signatures of differential selection were found when we compared the gray and white feather color populations using the F_ST and Hp approaches. When we combined previous functional studies with our genomic analyses we identified 26 genes (KITLG, MITF, TYRO3, KIT, AP3B1, SMARCA2, ROR2, CSNK1G3, CCDC112, VAMP7, SLC16A2, LOC106047519, RLIM, KIAA2022, ST8SIA4, LOC106044163, TRPM6, TICAM2, LOC106038556, LOC106038575, LOC106038574, LOC106038594, LOC106038573, LOC106038604, LOC106047489, and LOC106047492) that potentially regulate feather color in geese. These results substantially expand the catalog of potential feather color regulators in geese and provide a basis for further studies on domestication and avian feather coloration.