A single nucleotide polymorphism variant located in the cis-regulatory region of the ABCG2 gene is associated with mallard egg colour

Genome-wide association study reveals the candidate genes of humerus quality in laying duck

Bone plays a crucial role in poultry's health and production. However, during the selection and cage farming, there has been a decline in bone quality. As the development of breeding theory, researchers find that it's possible to enhance bone quality through selective breeding.This study measure 8 humerus quality in 260 samples of the 350-day-old female duck. By descripting the basic characteristic traits, mechanical property traits we found that all the bone quality traits had a large variable coefficient, especially mechanical properties trait (20-70%), indicating that there was a large difference in bone health status among laying ducks. The phenotypic correlations showed a high correlation between weight and density, diameter and perimeter, breaking and toughness (r = 0.52-0.68). And then, we performed the Genome-wide association study (GWAS) to reveal the candidate genes of humerus quality in ducks. Seven candidate protein-coding genes were identified with perimeter trait, and 52 protein-coding genes were associated with toughness trait. We also analysed the candidate region and performed KEGG and GO analyse for 75 candidate genes. Furthermore, the expression analyse of the above candidate genes in different stage of humerus and different tissues were performed. Finally, AP2A2, SMAD3, SMNDC1, NFIA, EPHB2, PMEPA1, UNC5C, ESR1, VAV3, NFATC2 deserve further focus. The obtained results can contribute to new insight into bone quality and provide new genetic biomarkers for application in duck breeding programs.

Identification and Characterization of the miRNA Transcriptome Controlling Green Pigmentation of Chicken Eggshells

Green eggs are mainly caused by inserting an avian endogenous retrovirus (EVA-HP) fragment into the SLCO1B3 gene. Although the genotypes for this insertion allele are consistent, eggshell color (ESC) may vary after a peak laying period; light-colored eggs are undesired by consumers and farmers and result in financial loss, so it is necessary to resolve this problem. miRNAs are small non-coding RNAs that exert essential functions in animal development and diseases. However, the regulatory miRNAs and detailed molecular mechanisms regulating eggshell greenness remain unclear. In the present study, we determined the genotype of green-eggshell hens through the detection of a homozygous allele insertion in the SLCO1B3 gene. The shell gland epithelium was obtained from green-eggshell hens that produced white and green shell eggs to perform transcriptome sequencing and investigate the important regulatory mechanisms that influence the ESC. Approximately 921 miRNAs were expressed in these two groups, which included 587 known miRNAs and 334 novel miRNAs, among which 44 were differentially expressed. There were 22 miRNAs that were significantly upregulated in the green and white groups, respectively, which targeted hundreds of genes, including KIT, HMOX2, and several solute carrier family genes. A Gene Ontology enrichment analysis of the target genes showed that the differentially expressed miRNA-targeted genes mainly belonged to the functional categories of homophilic cell adhesion, gland development, the Wnt signaling pathway, and epithelial tube morphogenesis. A KEGG enrichment analysis showed that the Hedgehog signaling pathway was significantly transformed in this study. The current study provides an overview of the miRNA expression profiles and the interaction between the miRNAs and their target genes. It provides valuable insights into the molecular mechanisms underlying green eggshell pigmentation, screening more effective hens to produce stable green eggs and obtaining higher economic benefits.

Genome-wide characterization and comparison of endogenous retroviruses among 3 duck reference genomes

Endogenous retroviruses (ERV) are viral genomes integrated into the host genome and can be stably inherited. Although ERV sequences have been reported in some avian species' genome, the duck endogenous retroviruses (DERV) genome has yet to be quantified. This study aimed to identify ERV sequences and characterize genes near ERVs in the duck genome by utilizing LTRhavest and LTRdigest tools to forecast the duck genome and analyze the distribution of ERV copies. The results revealed 1,607, 2,031, and 1,908 full-length ERV copies in the Pekin duck (ZJU1.0), Mallard (CAU_wild_1.0), and Shaoxing duck (CAU_laying_1.0) genomes, respectively, with average lengths of 7,046, 7,027, and 6,945 bp. ERVs are mainly distributed on the 1, 2, and sex chromosomes. Phylogenetic analysis demonstrated the presence of Betaretrovirus in 3 duck genomes, whereas Alpharetrovirus was exclusively identified in the Shaoxing duck genome. Through screening, 596, 315, and 343 genes adjacent to ERV were identified in 3 duck genomes, respectively, and their functions of ERV neighboring genes were predicted. Functional enrichment analysis of ERV-adjacent genes revealed enrichment for Focal adhesion, Calcium signaling pathway, and Adherens junction in 3 duck genomes. The overlapped genes were highly expressed in 8 tissues (brain, fat, heart, kidney, liver, lung, skin, and spleen) of 8-wk-old Mallard, revealing their important expression in different tissues. Our study provides a new perspective for understanding the quantity and function of DERVs, and may also provide important clues for regulating nearby genes and affecting the traits of organisms.

Genome-wide association study identifies multiple loci influencing duck serum biochemical indicators in the laying period

1. Laying performance is an important economic trait in poultry. The blood is essential in transporting nutrients to the yolk and albumen and is necessary for egg formation.2. This study calculated the phenotypic relationships of duck egg quality, egg production efficiency and 22 serum parameters in the egg-laying stage. Using a variety of methodologies, a genome-wide association study (GWAS) was carried out to uncover the genetic foundations of the 22 serum biochemical markers of laying ducks.3. Spearman correlation coefficients between the egg production (226-329 per day) and the serum parameters were all weak, being less than 0.3. This analysis was done on 22 serum parameters, with total protein (TP), total triglycerides (TG), calcium (Ca) and phosphorous (P) having the highest correlation coefficients (r = 0.56-0.88). The coefficients for blood markers, such as total cholesterol (CHOL), total bilirubin (TBIL), low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) varied from 0.70-0.94.4. Based on single-marker single-trait genome-wide analyses by a mixed linear model program of EMMAX, nine candidate genes were associated with enzyme traits (AST/ALT aspartate transaminase/glutamic-pyruvic transaminase, creatine kinase) and 19 candidate genes were associated with metabolism and protein-related serum parameters (glucose, total bile acid, uric acid (UA), albumin (ALB).5. The mvLMM (multivariate linear mixed model) of GEMMA software was used to carry out multiple trait integrated GWAS. Two candidate genes were found in the TP-TG-CA-P analysis and seven candidate genes in the CHOL_LDL-C_HDL-C_TBIL study. There was a high genetic correlation between the two groups.

Comparative analysis of hepatic transcriptomes and metabolomes of Changshun green-shell laying hens based on different green eggshell color intensities

The eggshell color of avian species is an important trait that is predominantly determined by the pigments biliverdin and protoporphyrin. Various factors affect eggshell pigment deposition and coloration; however, the underlying mechanisms remain unclear. We analyzed the hepatic transcriptomes and metabolomes of Changshun green-shell hens laying dark green and light green eggs to investigate the potential role of the liver in regulating the intensity of the green eggshell color. In total, 350 differentially expressed genes and 211 differentially altered metabolites were identified. Gene set enrichment analysis revealed that the enriched pathways and Gene Ontology (GO) terms were mainly associated with energy, immunity, and nutrient metabolism. Metabolite set enrichment analysis revealed that the enriched pathways were mainly associated with amino acid, vitamin, bile acid, and lipid metabolism. Moreover, gene-metabolite interaction network analysis revealed 1 subnetwork. Most genes and metabolites in this subnetwork were determined to be related to melanin metabolism and transport. In conclusion, our results suggest that hepatic melanin metabolism and transport are critical for eggshell coloration. Six candidate genes (CDKN2B, DDC, PYCR1, ABCG5, SLC3A1, and P2RX2) and 7 candidate metabolites (serotonin, 5-hydroxyindoleacetic acid, ornithine, acetylcholine, L-tryptophan, D-ornithine, and ADP) were suggested to play important roles in this process. Meanwhile, this study suggests that changes in hepatic energy metabolism, immune status, antioxidation activity, nutrient availability, and bile acid synthesis can impair eggshell coloration.

Selection on the promoter regions plays an important role in complex traits during duck domestication

BACKGROUND

Identifying the key factors that underlie complex traits during domestication is a great challenge for evolutionary and biological studies. In addition to the protein-coding region differences caused by variants, a large number of variants are located in the noncoding regions containing multiple types of regulatory elements. However, the roles of accumulated variants in gene regulatory elements during duck domestication and economic trait improvement are poorly understood.

RESULTS

We constructed a genomics, transcriptomics, and epigenomics map of the duck genome and assessed the evolutionary forces that have been in play across the whole genome during domestication. In total, 304 (42.94%) gene promoters have been specifically selected in Pekin duck among all selected genes. Joint multi-omics analysis reveals that 218 genes (72.01%) with selected promoters are located in open and active chromatin, and 267 genes (87.83%) with selected promoters were highly and differentially expressed in domestic trait-related tissues. One important candidate gene ELOVL3, with a strong signature of differentiation on the core promoter region, is known to regulate fatty acid elongation. Functional experiments showed that the nearly fixed variants in the top selected ELOVL3 promoter in Pekin duck decreased binding ability with HLF and increased gene expression, with the overexpression of ELOVL3 able to increase lipid deposition and unsaturated fatty acid enrichment.

CONCLUSIONS

This study presents genome resequencing, RNA-Seq, Hi-C, and ATAC-Seq data of mallard and Pekin duck, showing that selection of the gene promoter region plays an important role in gene expression and phenotypic changes during domestication and highlights that the variants of the ELOVL3 promoter may have multiple effects on fat and long-chain fatty acid content in ducks.

The Pigments in Eggshell with Different Color and the Pigment Regulatory Genes Expression in Corresponding Chicken’s Shell Gland

Eggshell colour is the unique appearance and economically valuable trait of eggs, whereas the colour is often short of uniformity, especially in the blue-shelled breeds, hence, their pigment differences and molecular mechanism need clarity. To investigate the relationship between the pigment content of eggshells and related gene expression in the eggshell glands of chickens, four subtypes of blue-shelled eggs ('Olive', 'Green', 'Blue', and 'Light') from the same blue-eggshell chicken line were selected; Hy-Line 'White' and 'Brown'-shelled eggs were used as control groups. The L*, a*, b* values, and protoporphyrin-IX and biliverdin contents in each group of eggshells were measured. In addition, the shell glands of the corresponding hens were collected to detect SLCO1B3 genotype and mRNA expression, and ABCG2 and HMOX1 transcription and protein expression. Eggshell colour L* values were negatively correlated with protoporphyrin-IX, b* values were positively correlated with total pigment content (P < 0.001), and a* values were positively correlated with protoporphyrin-IX (P < 0.001) but negatively with biliverdin. Moreover, all four blue-eggshell subtypes were SLCO1B3 homozygous, with SLCO1B3 mRNA expression in shell glands being significantly higher than in the White and Brown groups. ABCG2 and HMOX1 mRNA expression were highest in the Brown and Green groups, respectively (P < 0.05), and were positively correlated with protoporphyrin-IX (P < 0.001) and biliverdin contents in eggshells, respectively. Western blot and immunohistochemical results demonstrated that the Brown group had the highest ABCG2 expression (P < 0.05), followed by the Green and Olive groups. HMOX1 protein expression was higher in the Olive and Green groups (P < 0.05), and lowest in the White group. This study suggests that ABCG2 and HMOX1 have important regulatory roles in the production and transport of protoporphyrin-IX and biliverdin in blue-shelled chicken eggs, respectively.

Genome-wide association studies demonstrate the genes associated with perimysial thickness in ducks

The thickness of the perimysium has an essential effect on the tenderness of the meat. However, the genetic basis underlying perimysial thickness has not been determined. The objective of this study was to explore the quantitative trait loci (QTL) that influence perimysial thickness in an F₂ segregating population generated by Mallard × Pekin duck using the genome-wide association study (GWAS) method. Two QTL identified in chromosomes 27 and 13 displayed significant associations with perimysial thickness traits at the genome-wide level. The strongest association was the QTL located in chromosome 27, and this region had an effect on perimysial thickness and contained a promising candidate gene MAGI3 (Membrane-associated guanylate kinase, WW and PDZ domain containing 3). Meanwhile, association analysis showed that the top SNP within the MAGI3 gene was also associated with intramuscular fat content traits, which showed that perimysial thickness was positively correlated with intramuscular fat content. The second strongest association was the QTL region of chromosome 13. SUCLG2 (Succinate-CoA ligase GDP-forming subunit beta) is proximal to the top SNP and stood out as another candidate gene. Furthermore, the Transposase-Accessible Chromatin using Sequencing result showed that some key transcription factors (MYF5, MYOD1, KLF11) related to muscle development or energy metabolism might bind to the open regions of MAGI3 and SUCLG2. By analyzing the expression of different genotypes of the candidate gene, we speculate that different genotypes of MAGI3 may have an effect on breast muscle development, and then affect the thickness of the perimysium. This study maps two major genes of the duck breast muscle perimysial thickness trait, which helps to characterize muscle development and contributes to the genetic improvement of meat yield and quality in livestock.

Genomic divergence during artificial selection by feed conversion ratio in Pekin ducks

Pekin ducks are world-famous for its fast growth and have become the majority of breeds rearing in duck industry. Feed conversion ratio (FCR) is an important trait in Pekin ducks breeding and production, and the underlying biological processes are complex. To gain an insight to the possible biological mechanism underlying the FCR in Pekin ducks, an artificial selection population (S) and a natural population (Z7) were used in this study. The FCR of S line decreased from 2.184 ± 0.057 in the first generation to 1.886 ± 0.063 in the eighth generation, which displays significantly low FCR (p = 0.0032) than that of the Z7 line (2.23 ± 0.046). Then, 9 samples from eighth generation of S line and 10 samples from Z7 were used for whole-genome resequencing. Analyses of F_ST, θπ and XP-EHH revealed 450, 479 and 356 candidate genes, which involved in 1,955, 1,933 and 1,964 candidate divergent regions (CDRs), respectively. And the integration of three approaches resulted in 30 overlapping genes. Functional analysis of 30 candidate genes revealed that variants of KCNQ1 and ADCY7, which were involved in the pancreatic secretion signal pathway, could be important molecular markers for high feed conversion efficiency in S line breeding.

Genome-wide association analysis reveals 6 copy number variations associated with the number of cervical vertebrae in Pekin ducks

As a critical developmental stage in vertebrates, the vertebral column formation process is under strict control; however, we observed variations in the number of cervical vertebrae in duck populations in our previous study. Here, we further explored the variations in the number of vertebrae in two duck populations: 421 Pekin duck × mallard F2 ducks and 850 Pekin ducks. Using resequencing data of 125 Pekin ducks with different numbers of cervical vertebrae and 352 Pekin duck × mallard F2 ducks with different numbers of thoracic vertebrae, we detected whole-genome copy number variations (CNVs) and implemented a genome-wide association study (GWAS) to identify the genetic variants related to the traits. The findings verified the existence of variations in the number of cervical vertebrae in duck populations. The number of cervical vertebrae in most ducks was 15, while that in a small number of the ducks was 14 or 16. The number of cervical vertebrae had a positive influence on the neck production, and one cervical vertebra addition could increase 11 g or 2 cm of duck neck. Genome-wide CNV association analysis identified six CNVs associated with the number of cervical vertebrae, and the associated CNV regions covered 15 genes which included WNT10A and WNT6. These findings improve our understanding of the variations in the number of vertebrae in ducks and lay a foundation for future duck breeding.

Genome-wide association analysis of the primary feather growth traits of duck: identification of potential Loci for growth regulation

The feather is an important epidermal appendage, plays an important role in the life activities of avian specie, and has important economic value. Revealing the molecular regulation mechanism of feather growth has a significant meaning in studying adaptive evolution, physiology, and mating of avian species and also provides a theoretical reference for poultry breeding. In this study, the genome-wide association analysis (GWAS) of 358 ducks was based on primary feather length phenotypic data (28-60 d), length growth rates (LGRs), and maturity scores (60 d) to explore the genetic basis affecting feather growth and maturation. The results showed that, among the primary feather 1 to 5 in ducks, the mean LGR of primary feather 2 was the fastest, with the longest length. The primary feathers in males grew and matured slightly faster than in females. The mean maturity scores of primary feather 10∼7 were higher than primary feather 1 to 3 in ducks. GWAS further showed 116 SNPs associated with feather length traits. In addition, 2 candidate regions (Chr1: 127,407,230-127,524,879 bp and Chr21: 182,061,707-183,616,298 bp) were associated with LGR, which contain total 13 candidate genes (The extremely significant SNPs were mainly located in 2 genes: Chr1: REPS2 and Chr21: PTPRT). Four candidate regions (Chr1: 29,113,036-28,675,018 bp, Chr2: 18,253,612-149,111,290 bp, Chr15: 6,489,774 to 12,138,221 bp and Chr21: 6,578,021-8,472,904 bp) were associated with feather maturity, which contain total 24 candidate genes (The extremely significant SNPs were mainly located in 4 genes: Chr1: IMMP2L, DOCK4 and DDX10, Chr2: LDLRAD4). In conclusion, sex factors influence feather growth and maturity, and the genetic basis of the growth /maturity trait between different feathers is similar. REPS2, PTPRT genes, and IMMP2L, DOCK4, DDX10, and LDLRAD4 are important candidate genes that influence feather growth and maturity, respectively.

The study of candidate genes in the improvement of egg production in ducks – a review

Duck is the second-largest poultry species aside from chicken. The rate of egg production is a major determinant of the economic income of poultry farmers. Among the reproductive organs, the ovary is a major part of the female reproductive system which is highly important for egg production. Based on the importance of this organ, several studies have been carried out to identify candidate genes at the transcriptome level, and also the expression level of these genes at different tissues or egg-laying conditions, and single nucleotide polymorphism (SNPs) of genes associated with egg production in duck. In this review, expression profile and association study analyses at SNPs level of different candidate genes with egg production traits of duck were highlighted. Furthermore, different studies on transcriptome analysis, Quantitative Trait Loci (QTL) mapping, and Genome Wide Association Study (GWAS) approach used to identify potential candidate genes for egg production in ducks were reported. This review would widen our knowledge on molecular markers that are associated or have a positive correlation to improving egg production in ducks, for the increasing world populace.

Genome-wide association and selective sweep analyses reveal genetic loci for FCR of egg production traits in ducks

Background

As a major economic trait in poultry, egg production efficiency attracts widespread interest in breeding and production. However, limited information is available about the underlying genetic architecture of egg production traits in ducks. In this paper, we analyzed six egg production-related traits in 352 F₂ ducks derived from reciprocal crosses between mallard and Pekin ducks.

Results

Feed conversation ratio (FCR) was positively correlated with feed intake but negatively correlated with egg-related traits, including egg weight and egg production, both phenotypically and genetically. Estimates of pedigree-based heritability were higher than 0.2 for all traits investigated, except hip-width. Based on whole-genome sequencing data, we conducted genome-wide association studies to identify genomic regions associated with these traits. In total, 11 genomic regions were associated with FCR. No genomic regions were identified as significantly associated with hip-width, total feed intake, average daily feed intake, and total egg production. Analysis of selective sweeps between mallard and Pekin ducks confirmed three of these genomic regions on chromosomes 13, 3 and 6. Within these three regions, variants in candidate genes that were in linkage disequilibrium with the GWAS leader single nucleotide polymorphisms (SNPs) (Chr13:2,196,728, P = 7.05 × 10^–14; Chr3:76,991,524, P = 1.06 × 10^–12; Chr6:20,356,803, P = 1.14 × 10^–10) were detected. Thus, we identified 31 potential candidate genes associated with FCR, among which the strongest candidates are those that are highly expressed in tissues involved in reproduction and nervous system functions of ducks: CNTN4, CRBR, GPR63, KLHL32, FHL5, TRNT1, MANEA, NDUFAF4, and SCD.

Conclusions

For the first time, we report the identification of genomic regions that are associated with FCR in ducks and our results illustrate the genomic changes that occurred during their domestication and are involved in egg production efficiency.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00684-5.

Avian Coloration Genetics: Recent Advances and Emerging Questions

The colorful phenotypes of birds have long provided rich source material for evolutionary biologists. Avian plumage, beaks, skin, and eggs-which exhibit a stunning range of cryptic and conspicuous forms-inspired early work on adaptive coloration. More recently, avian color has fueled discoveries on the physiological, developmental, and-increasingly-genetic mechanisms responsible for phenotypic variation. The relative ease with which avian color traits can be quantified has made birds an attractive system for uncovering links between phenotype and genotype. Accordingly, the field of avian coloration genetics is burgeoning. In this review, we highlight recent advances and emerging questions associated with the genetic underpinnings of bird color. We start by describing breakthroughs related to 2 pigment classes: carotenoids that produce red, yellow, and orange in most birds and psittacofulvins that produce similar colors in parrots. We then discuss structural colors, which are produced by the interaction of light with nanoscale materials and greatly extend the plumage palette. Structural color genetics remain understudied-but this paradigm is changing. We next explore how colors that arise from interactions among pigmentary and structural mechanisms may be controlled by genes that are co-expressed or co-regulated. We also identify opportunities to investigate genes mediating within-feather micropatterning and the coloration of bare parts and eggs. We conclude by spotlighting 2 research areas-mechanistic links between color vision and color production, and speciation-that have been invigorated by genetic insights, a trend likely to continue as new genomic approaches are applied to non-model species.

Genome-wide association analysis reveals that EDNRB2 causes a dose-dependent loss of pigmentation in ducks

Background

Birds have various plumage color patterns, and spot is a common phenotype. Herein, we conducted genome-wide association studies (GWAS) in a population of 225 ducks with different sized black spots to reveal the genetic basis of this phenomenon.

Results

First, we quantified the black spot phenotype within the duck population. The results showed that the uncolored area of the body surface first appeared on the ventral side. With increasing duck age, the area of the black spots was highly conserved across the whole body surface. The GWAS results identified a 198 kb (Chr4: 10,149,651 bp to 10,348,068 bp) genetic region that was significantly associated with the black spot phenotype. The conditional GWAS and linkage disequilibrium (LD) analysis further narrowed the ultimate candidate region to 167 kb (Chr4: 10,180,939 bp to 10,348,068 bp). A key gene regulating melanoblast migration and differentiation, EDNRB2 (Endothelin B receptor-like), was found in the candidate region and having significant mRNA expression level changes in embryonic duck skin tissue with different spot sizes. The significant SNPs (single nucleotide polymorphisms) associated with the EDNRB2 gene were annotated, and two mutations (Chr4: 10,180,939 T > C and Chr4: 10,190,671 A > T) were found to result in the loss of binding sites for two trans-factors, XBP1 and cMYB. The phenotypic effect of these two mutations suggested that they can regulate the size of black spots in a dose-dependent manner, and Chr4: 10,180,939 T > C was the major allele locus.

Conclusions

Our results revealed that EDNRB2 was the gene responsible for the variation in duck body surface spot size. Chr4: 10,180,939 T > C was the major allele that explained 49.5 % (dorsal side) and 32.9 % (ventral side) of the variation in duck body surface spot size, while 32.1 % (dorsal side) and 19.1 % (ventral side) of the variation could be explained by Chr4: 10,190,671 A > T. The trans-factor prediction also suggested that XBP1 and cMYB have the potential to interact with EDNRB2, providing new insights into the mechanism of action of these genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07719-7.

Genetic variations for egg internal quality of ducks revealed by genome-wide association study

Egg internal quality traits are important traits related to egg production in poultry industry. To better understand the genetic architecture of egg internal quality traits in ducks, we performed genetic parameters estimates and a genome-wide association study (GWAS). The phenotypic values of egg weight, yolk color, albumin height (AH), yolk weight, and Haugh unit (HU) were collected individually from 352 F₂ laying ducks produced by reciprocal crosses between mallards and Pekin ducks, and their genotypes were assayed by whole genome re-sequencing. The results showed that the AH and HU traits have a clear coefficient of variance, around 15% for both mallards and Pekin ducks. The pedigree-based genetic parameters estimates rane from 0.26 to 0.71 for all eight egg quality traits, while the highest heritability was 0.71 for egg weight. The GWAS showed that a clear signal was associated with AH and HU traits. The locus zoom analysis and conditional GWAS helped to narrow the candidate region to ~5.8-Mb spanning from 14.7 to 20.5 Mb on Chromosome 5, which harbored 111 candidate genes. MUC6 and LDLRAD3 were finally promised as the major candidate genes affecting albumen composition. Our data revealed the egg internal quality traits for the first time in ducks, which provides a theoretical basis and technological support for improving duck egg internal quality.