Whole-genome resequencing provides insights into the population structure and domestication signatures of ducks in eastern China

Comparative genome-wide association study on body weight in Chinese native ducks using four models

The Jinling White duck represents a newly developed breed characterized by a rapid growth rate and a superior meat quality, offering significant economic value and research potential; however, the genetic basis underlying their body weight traits remains less understood. Here, we performed whole-genome resequencing for 201 diverse Jinling White male ducks and conducted population genomic analyses, suggesting a rich genetic diversity within the Jinling White duck population. Equipped with our genomic resources, we applied genome-wide association analysis for body weight on birth (BWB), body weight on 1 wk (BW1), body weight on 3 wk (BW3), body weight on 5 wk (BW5) and body weight on 7 wk (BW7) using 4 statistical models. Comparative studies indicated that factored spectrally transformed linear mixed models (FaST-LMM) demonstrated the most superior efficiency, yielding more results with the minimal false positives. We discovered that PUS7, FBXO11, FOXN2, MSH6, and SLC4A4 were associated with BWB. RAG2, and TMEFF2 were candidate genes for BW1, and STARD13, Klotho, ZAR1L are likely candidates for BW3 and BW5. PLXNC1, ATP1A1, CD58, FRYL, OCIAD1, and OCIAD2 were linked to BW7. These findings provide a genetic reference for the selection and breeding of Jinling White ducks, while also deepened our understanding of Growth and development phenotypic in ducks.

Unique feather color characteristics and transcriptome analysis of hair follicles in Liancheng White ducks

Avian feather color is a fascinating trait, and the genetic mechanism of duck plumage formation is still in the preliminary stage. In this study, feather color of Liancheng White ducks was analyzed by determination of melanin content and RNA-seq analysis. In this research, 9 ducks from Mallards (n = 3), Liancheng White (n = 3) and Pekin ducks (n = 3) were used by high performance liquid chromatography (HPLC) and Masson-Fontana staining to reveal the difference of feather melanin content. RNA-seq from 11 hair follicle tissues (1- and 8-wk-old) of Liancheng White ducks (n = 5) and Pekin ducks (n = 7) was used to analyze the candidate genes for the feather melanin synthesis, and Immunofluorescence experiment was used to show the protein expression in 6 black- and white-feathered ducks. Pectorale, skin, liver, fat, brain, heart, kidney, lung, spleen of an 8-wk-old black-feathered Mallard were collected for candidate gene expression. The results showed that the contents of feathers, beak, web melanin in Liancheng White ducks were higher than in Pekin ducks (p < 0.05). Melanin within hair follicles was located in the barb ridge and hair matrix of black feather duck, also we found that TYRP1, TYR, SOX10 genes were differentially expressed between Liancheng White and Pekin ducks (p < 0.05), and these genes were mainly expressed showed in duck skin tissues. This study revealed the unique feather color phenotype of Liancheng White duck shedding light on the transcriptome that underlies it.

Whole-genome resequencing revealed the population structure and selection signal of 4 indigenous Chinese laying ducks

The assessment of animal genetic structure had significant importance for the preservation and breeding of animal germplasm resources. Selection signals are genotype markers generated during the process of biological evolution, and the detection of selection signals could reveal the direction of species evolution. The aim of this study was to generate a whole-genome resequencing data from Jinding duck, Shanma duck, Youxian Partridge duck, and Taiwan Brown tsaiya duck to reveal their population structure and selection signals. The population structure analysis revealed significant genetic differences among the 4 indigenous laying ducks, indicating their independent lineage. Specifically, Shanma duck and Youxian partridge duck were closely and likely originated from a common ancestor. In addition, selection sweep analysis was performed using the population genetic differentiation coefficient (Fst) and nucleotide diversity ratio (π ratio). The top 5% was used as the threshold for the Fst and π ratio, and the 2 thresholds were combined to identify selected genomic regions. In the selected regions of the 3 comparison groups, 136, 143, and 268 candidate genes were detected. Further screening of all candidate genes revealed that 35 candidate genes appeared simultaneously in 3 comparative groups, with 16 genes annotated. The 16 genes were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The results revealed 5 functional genes (AQP3, PIK3C3, NOL6, RPP25, and DCTN3) that may be related to important economic traits in laying ducks and involved mainly invasopressin-regulated water reabsorption, ribosome biogenesis, and the PI3K signaling pathway. The results provide insights into the protection and exploitation of genetic resources of Chinese indigenous laying ducks.

Whole-genome selection signature differences between Chaohu and Ji'an red ducks

Assessing the genetic structure of local varieties and understanding their genetic data are crucial for effective management and preservation. However, the genetic differences among local breeds require further explanation. To enhance our understanding of their population structure and genetic diversity, we conducted a genome-wide comparative study of Chaohu and Ji'an Red ducks using genome sequence and restriction site-associated DNA sequencing technology. Our analysis revealed a distinct genetic distinction between the two breeds, leading to divided groups. The phylogenetic tree for Chaohu duck displayed two branches, potentially indicating minimal impact from artificial selection. Additionally, our ROH (runs of homozygosity) analysis revealed that Chaohu ducks had a lower average inbreeding coefficient than Ji'an Red ducks. We identified several genomic regions with high genetic similarity in these indigenous duck breeds. By conducting a selective sweep analysis, we identified 574 candidate genes associated with muscle growth (BMP2, ITGA8, MYLK, and PTCH1), fat deposits (ELOVL1 and HACD2), and pigmentation (ASIP and LOC101797494). These results offer valuable insights for the further enhancement and conservation of Chinese indigenous duck breeds.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Molecular genetic characterization and meat-use functional gene identification in Jianshui yellow-brown ducks through combined resequencing and transcriptome analysis

The Jianshui yellow-brown duck is a unique country-specific waterfowl species in Yunnan Province, well known for its tender meat. However, there is a lack of comprehensive systematic research on the molecular genetic characteristics, especially germplasm resources and economic traits, of the Jianshui yellow-brown ducks. This study investigated the molecular genetic characteristics of Jianshui yellow-brown ducks, compared their selection signals with those of ancestral mallard and meat-type Pekin ducks, and identified genes specific to their meat-use performance. Furthermore, this study also evaluated the breeding potential for its meat performance. In this study, phylogenetic trees, PCA and Admixture analysis were used to investigate the population genetic structure among local duck breeds in China; population genetic differentiation index (Fst), nucleotide diversity and Tajima's D were used to detect selected loci and genes in the population of Jianshui yellow-brown ducks; and transcriptome technology was used to screen for differentially expressed genes in the liver, sebum and breast muscle tissues, and finally, the results of the genome selection signals and transcriptome data were integrated to excavate functional genes affecting the meat performance of the Jianshui yellow-brown ducks. The results of the genetic structure of the population showed that Jianshui yellow-brown ducks were clustered into a separate group. Selection signal analysis indicated significant selection pressure on certain genes related to meat characteristics (ELOVL2, ELOVL3, GDF10, VSTM2A, PHOSPHO1, and IGF2BP1) in both Jianshui yellow-brown ducks and mallards. Transcriptomic data analysis suggested that ELOVL3, PHOSPHO1, and GDF10 are vital candidate genes influencing meat production and quality in Jianshui yellow-brown ducks. A comparison of selection signals between Jianshui yellow-brown ducks and Pekin ducks revealed only 21 selected genes in the Jianshui yellow-brown duck population, and no significant genes were related to meat traits. Moreover, whole-genome resequencing data suggested that the Jianshui yellow-brown duck represents a unique category with distinct genetic mechanisms. Through selection signaling and transcriptomic approaches, we successfully screened and identified important candidate genes affecting meat traits in Jianshui yellow-brown ducks. Furthermore, the Jianshui yellow-brown duck has good potential for improved meat performance, highlighting the need for further improvement.

The meaning of wild: Genetic and adaptive consequences from large-scale releases of domestic mallards

The translocation of individuals around the world is leading to rising incidences of anthropogenic hybridization, particularly between domestic and wild congeners. We apply a landscape genomics approach for thousands of mallard (Anas platyrhynchos) samples across continental and island populations to determine the result of over a century of supplementation practices. We establish that a single domestic game-farm mallard breed is the source for contemporary release programs in Eurasia and North America, as well as for established feral populations in New Zealand and Hawaii. In particular, we identify central Europe and eastern North America as epicenters of ongoing anthropogenic hybridization, and conclude that the release of game-farm mallards continues to affect the genetic integrity of wild mallards. Conversely, self-sustaining feral populations in New Zealand and Hawaii not only show strong differentiation from their original stock, but also signatures of local adaptation occurring in less than a half-century since game-farm mallard releases have ceased. We conclude that 'wild' is not singular, and that even feral populations are capable of responding to natural processes. Although considered paradoxical to biological conservation, understanding the capacity for wildness among feral and feral admixed populations in human landscapes is critical as such interactions increase in the Anthropocene.

Whole-genome resequencing to unveil genetic characteristics and selection signatures of specific pathogen-free ducks

Specific pathogen-free ducks are important high-grade laboratory animals, with a key role in research related to poultry biosecurity, production, and breeding. However, the genetic characteristics of experimental duck varieties remain poorly explored. Herein we performed whole-genome resequencing to construct a single nucleotide polymorphism genetic map of the genomes of 3 experimental duck varieties [Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)] to determine their genetic characteristics and identify selection signatures. Subsequent analyses of population structure and genetic diversity revealed that each duck variety formed a monophyletic group, with SM showing richer genetic diversity than JD and SX. Further, on exploring shared selection signatures, we found 2 overlapping genomic regions on chromosome Z of all experimental ducks, which comprised immune response-related genes (IL7R and IL6ST). Moreover, growth and skeletal development (IGF1R and GDF5), meat quality (FoxO1), and stress resistance (HSP90B1 and Gpx8-b) candidate gene loci were identified in strongly selected signatures specific to JD, SM, and SX, respectively. Our results identified the population genetic basis of experimental ducks at the whole-genome level, providing a framework for future molecular investigations of genetic variations and phenotypic changes. We believe that such studies will eventually contribute to the management of experimental animal resources.

Genomic characteristics and selection signals of Zhongshan ducks

The Zhongshan duck (ZSP) is a duck breed indigenous to China and is known for its moderate body size, strong disease resistance, tender meat, and little subcutaneous fat. However, the genomic basis of such excellent breeding characteristics remains poorly understood. Therefore, we generated whole-genomes of 58 ZSPs and 180 other indigenous Chinese ducks (60 Jinding ducks, 60 Shan Partridge ducks, and 60 Liancheng white ducks) and identified 10 560 032 single nucleotide polymorphisms and 1 334 893 structural variants. Based on genetic diversity and population structure indices, our results confirm that the ZSP is a unique germplasm resource. In addition, three reproduction-related genes (i.e., OAZ, AMH, and RLF) were located in highly differentiated regions between the ZSP and the other three duck breeds (Jinding duck; Liancheng White duck; Shan Partridge duck), suggesting that these genes may have a strong influence on egg production. Among these genes, AMH may have introgressed from an unknown species of the Anatidae family. We also identified other significant genes in the significantly differentiated window (i.e., 1% cut-off), some of which are responsible for growth and development (SEMA5B and MIB1), metabolism (EDEM3 and Xylb), skeletal system morphogenesis (bglap and MGP), and egg shape (ITPR2). These findings highlight the genetic characteristics of the ZSP that shape an array of its morphological traits. Overall, this study should facilitate a more fine-scale approach towards improving the ZSP and other indigenous ducks in China and even all over the world.

Chinese Domestic Ducks Evolved from Mallard Duck (Anas platyrhynchos) and Spot-Billed Duck (A. zonorhyncha)

China has a rich genetic resource in its 34 domestic duck breeds. In order to detect the genetic diversity and explore the origin of these indigenous duck populations, the mitochondrial DNA (mtDNA) control region was sequenced and analyzed for 208 individual ducks, including 22 domestic breeds, wild Mallards ducks, Eastern spot-billed ducks, White Muscovy ducks, and Black Muscovy ducks. The haplotype diversity (Hd) was 0.653 and the average nucleotide diversity (Pi) was 0.005, indicating moderate genetic diversity. Sixty haplotypes were detected, and the maximum-likelihood (ML) phylogenetic tree and median-joining (MJ) network were generated from the sequence analyses. In this study, haplotypes from the Mallard duck (Anas platyrhynchos) were detected in most of the Chinese domestic duck breeds. In addition, the Eastern spot-billed duck (A. zonorhyncha) H8 haplotype was detected in two duck breeds. Only two haplotypes were found in Muscovy ducks, suggesting low genetic diversity within this population. The sequence and haplotype analyses revealed that both A. platyrhynchos and A. zonorhyncha contributed to the evolution of domestic ducks in China.

Genomic diversity and signals of selection processes in wild and farm-reared red-legged partridges (Alectoris rufa)

The genetic dynamics of wild populations with releases of farm-reared reinforcements are very complex. These releases can endanger wild populations through genetic swamping or by displacing them. We assessed the genomic differences between wild and farm-reared red-legged partridges (Alectoris rufa) and described differential selection signals between both populations. We sequenced the whole genome of 30 wild and 30 farm-reared partridges. Both partridges had similar nucleotide diversity (π). Farm-reared partridges had a more negative Tajima's D and more and longer regions of extended haplotype homozygosity than wild partridges. We observed higher inbreeding coefficients (F_IS and F_ROH) in wild partridges. Selective sweeps (Rsb) were enriched with genes that contribute to the reproductive, skin and feather colouring, and behavioural differences between wild and farm-reared partridges. The analysis of genomic diversity should inform future decisions for the preservation of wild populations.

The study of selection signature and its applications on identification of candidate genes using whole genome sequencing data in chicken - a review

Chicken is a major source of protein for the increasing human population and is useful for research purposes. There are almost 1,600 distinct regional breeds of chicken across the globe, among which a large body of genetic and phenotypic variations has been accumulated due to extensive natural and artificial selection. Moreover, natural selection is a crucial force for animal domestication. Several approaches have been adopted to detect selection signatures in different breeds of chicken using whole genome sequencing (WGS) data including integrated haplotype score (iHS), cross-populated extend haplotype homozygosity test (XP-EHH), fixation index (F_ST), cross-population composite likelihood ratio (XP-CLR), nucleotide diversity (Pi), and others. In addition, gene enrichment analyses are utilized to determine KEGG pathways and gene ontology (GO) terms related to traits of interest in chicken. Herein, we review different studies that have adopted diverse approaches to detect selection signatures in different breeds of chicken. This review systematically summarizes different findings on selection signatures and related candidate genes in chickens. Future studies could combine different selection signatures approaches to strengthen the quality of the results thereby providing more affirmative inference. This would further aid in deciphering the importance of selection in chicken conservation for the increasing human population.

Effects of Compound Probiotics on Growth Performance, Serum Biochemical and Immune Indices, Antioxidant Capacity, and Intestinal Tissue Morphology of Shaoxing Duck

This experiment was conducted to investigate the effects of compound probiotics on growth performance, serum biochemical and immune indices, antioxidant capacity, and the intestinal tissue morphology of Shaoxing ducks. A total of 640 1-day-old healthy Shaoxing ducks of similar body weight were randomly divided into two treatment groups with eight replicates each and forty ducks per replicate. The ducks were fed a basal diet (Ctrl) or a basal diet supplemented with 0.15% compound probiotics (CP) for 125 d. The results revealed that the live body weight (BW; day 85 and 125) and the average daily gain (ADG; 28−85 and 85−125 d) of the CP group were significantly higher (p < 0.05) than those of the Ctrl group. In the CP group, total protein and total cholesterol contents were significantly increased (p < 0.05) on days 28 and 85, while triglyceride and low-density lipoprotein contents were significantly decreased (p < 0.05) on day 85. Furthermore, interferon-γ content was significantly increased (p < 0.05) in the CP group on days 28, 85, and 125. Interleukin-2 content was significantly increased (p < 0.05) in the CP group on days 28 and 85. Interleukin-4 content was significantly decreased (p < 0.05) in the CP group on day 85. Moreover, in the CP group, superoxide dismutase content was significantly increased (p < 0.05) on days 28 and 125, and glutathione peroxidase content was significantly increased (p < 0.05) on day 125. The crypt depth (CD) in the duodenum of the CP group was significantly decreased (p < 0.05) on days 28 and 125, whereas the villus height (VH) in the jejunum of the CP group was significantly increased (p < 0.05) on days 85 and 125. The VH/CD ratio in the ileum of the CP group was significantly increased (p < 0.05) on days 28 and 85. The VH in the ileum of the CP group was significantly increased (p < 0.05) on day 28. The CD in the ileum of the CP group was significantly decreased (p < 0.05) on day 28. In summary, the compound probiotics improved the growth performance, increased serum biochemical and immune indices, increased antioxidant capacity, and improved the intestinal tissue morphology of Shaoxing ducks.

Sex Identification of Feather Color in Geese and the Expression of Melanin in Embryonic Dorsal Skin Feather Follicles

In production practice, we have found that the gray and black down on the backs of the Holdobaggy goslings is usually darker in females than in males. Melanin is the key pigment affecting the color of poultry plumage. Therefore, to determine whether the darkness of the dorsal plumage of the Holdobaggy goslings is related to sex, we study the melanin in the feather follicles of the dorsal skin during the embryonic period. The feather follicle structure and melanin distribution on the dorsal surface of the goose embryo is observed by HE staining and melanin-specific staining. The melanin content in the feather follicles of the dorsal skin of goslings is determined by ELISA. The results showed that the melanin content is higher in female geese than in males (p < 0.05). In addition, we also analyze the mRNA and protein expression levels of melanin-related genes (TYRP1 and ASIP) by quantitative real-time PCR and Western blotting analysis. The results show that the mRNA expression level of TYRP1 is significantly higher in the females’ dorsal skin feather follicles (p < 0.05), while the mRNA expression level of ASIP is significantly higher in the dorsal skin feather follicles of male geese (p < 0.05). In conclusion, the difference between males and females in the color of the black feathers on the dorsal track of the Holdobaggy goslings is verified, and it is feasible to identify the sex by the initial plumage color.

Tissue Specific Transcriptome Changes Upon Influenza A Virus Replication in the Duck

Ducks are the natural host and reservoir of influenza A virus (IAV), and as such are permissive to viral replication while being unharmed by most strains. It is not known which mechanisms of viral control are globally regulated during infection, and which are specific to tissues during infection. Here we compare transcript expression from tissues from Pekin ducks infected with a recombinant H5N1 strain A/Vietnam 1203/04 (VN1203) or an H5N2 strain A/British Columbia 500/05 using RNA-sequencing analysis and aligning reads to the NCBI assembly ZJU1.0 of the domestic duck (Anas platyrhynchos) genome. Highly pathogenic VN1203 replicated in lungs and showed systemic dissemination, while BC500, like most low pathogenic strains, replicated in the intestines. VN1203 infection induced robust differential expression of genes all three days post infection, while BC500 induced the greatest number of differentially expressed genes on day 2 post infection. While there were many genes globally upregulated in response to either VN1203 or BC500, tissue specific gene expression differences were observed. Lungs of ducks infected with VN1203 and intestines of birds infected with BC500, tissues important in influenza replication, showed highest upregulation of pattern recognition receptors and interferon stimulated genes early in the response. These tissues also appear to have specific downregulation of inflammatory components, with downregulation of distinct sets of proinflammatory cytokines in lung, and downregulation of key components of leukocyte recruitment and complement pathways in intestine. Our results suggest that global and tissue specific regulation patterns help the duck control viral replication as well as limit some inflammatory responses in tissues involved in replication to avoid damage.

Functional Haplotypes and Evolutionary Insight into the Granule-Bound Starch Synthase II (GBSSII) Gene in Korean Rice Accessions (KRICE_CORE)

Granule-bound starch synthase 2 (GBSSII), a paralogous isoform of GBSSI, carries out amylose biosynthesis in rice. Unlike GBSSI, it mainly functions in transient organs, such as leaves. Despite many reports on the starch gene family, little is known about the genetics and genomics of GBSSII. Haplotype analysis was conducted to unveil genetic variations (SNPs and InDels) of GBSSII (OS07G0412100) and it was also performed to gain evolutionary insight through genetic diversity, population genetic structure, and phylogenetic analyses using the KRICE_CORE set (475 rice accessions). Thirty nonsynonymous SNPs (nsSNPs) were detected across the diverse GBSSII coding regions, representing 38 haplotypes, including 13 cultivated, 21 wild, and 4 mixed (a combination of cultivated and wild) varieties. The cultivated haplotypes (C_1-C_13) contained more nsSNPs across the GBSSII genomic region than the wild varieties. Nucleotide diversity analysis highlighted the higher diversity values of the cultivated varieties (weedy = 0.0102, landrace = 0.0093, and bred = 0.0066) than the wild group (0.0045). The cultivated varieties exhibited no reduction in diversity during domestication. Diversity reduction in the japonica and the wild groups was evidenced by the negative Tajima's D values under purifying selection, suggesting the domestication signatures of GBSSII; however, balancing selection was indicated by positive Tajima's D values in indica. Principal component analysis and population genetics analyses estimated the ambiguous evolutionary relationships among the cultivated and wild rice groups, indicating highly diverse structural features of the rice accessions within the GBSSII genomic region. FST analysis differentiated most of the classified populations in a range of greater FST values. Our findings provide evolutionary insights into GBSSII and, consequently, a molecular breeding program can be implemented for select desired traits using these diverse nonsynonymous (functional) alleles.