Origin and genetic diversity of Chinese domestic ducks

Inter and intra population genetic variability in ducks under conservation programs

This study focuses on estimation of the inter and intra population genetic variability of 6 duck populations. Microsatellite loci were used to assess the genetic variation and population structure of 6 duck populations under a conservation program in Poland. DNA polymorphism was assessed using 24 microsatellite markers and 50 individuals from each population. Polymorphism information content (PIC), heterozygosity with 2 estimators of genetic differentiation (FST and GST), and Nei's standard genetic distance were calculated. The results showed that these 6 endangered duck populations showed high genetic polymorphism. Observed heterozygosity within populations ranged from 0.14 to 0.83, with the average value of 0.58. PIC within populations ranged from 0.038 (P-8 and P-9 lines) to 0.89 (LsA line). Average number of alleles in the studied populations ranged from 4.5 (KhO-1 line) to 7.3 (LsA line). Based on the results, the pairs of lines LsA: P-33 and P-8: P-9 were found to be the most related; and the most genetically distant group was KhO-1 line, which originated as a cross between Khaki Campbell with Orpington duck.

Duck pan-genome reveals two transposon insertions caused bodyweight enlarging and white plumage phenotype formation during evolution

Structural variations (SVs) are a major source of domestication and improvement traits. We present the first duck pan-genome constructed using five genome assemblies capturing ∼40.98 Mb new sequences. This pan-genome together with high-depth sequencing data (∼46.5×) identified 101,041 SVs, of which substantial proportions were derived from transposable element (TE) activity. Many TE-derived SVs anchoring in a gene body or regulatory region are linked to duck's domestication and improvement. By combining quantitative genetics with molecular experiments, we, for the first time, unraveled a 6945 bp Gypsy insertion as a functional mutation of the major gene IGF2BP1 associated with duck bodyweight. This Gypsy insertion, to our knowledge, explains the largest effect on bodyweight among avian species (27.61% of phenotypic variation). In addition, we also examined another 6634 bp Gypsy insertion in MITF intron, which triggers a novel transcript of MITF, thereby contributing to the development of white plumage. Our findings highlight the importance of using a pan-genome as a reference in genomics studies and illuminate the impact of transposons in trait formation and livestock breeding.

Evolution, Inheritance, and Strata Formation of the W Chromosome in Duck (Anas platyrhynchos)

The nonrecombining female-limited W chromosome is predicted to experience unique evolutionary processes. Difficulties in assembling W chromosome sequences have hindered the identification of duck W-linked sequences and their evolutionary footprint. To address this, we conducted three initial contig-level genome assemblies and developed a rigorous pipeline by which to successfully expand the W-linked data set, including 11 known genes and 24 newly identified genes. Our results indicate that the W chromosome expression may not be subject to female-specific selection; a significant convergent pattern of upregulation associated with increased female-specific selection was not detected. The genetic stability of the W chromosome is also reflected in the strong evolutionary correlation between it and the mitochondria; the complete consistency of the cladogram topology constructed from their gene sequences proves the shared maternal coevolution. By detecting the evolutionary trajectories of W-linked sequences, we have found that recombination suppression started in four distinct strata, of which three were conserved across Neognathae. Taken together, our results have revealed a unique evolutionary pattern and an independent stratum evolutionary pattern for sex chromosomes.

Comparison of the Gut Microbial Communities of Domestic and Wild Mallards (Anas platyrhynchos) Based on High-Throughput Sequencing Technology

Mallards (Anas platyrhynchos) are currently one of the most popular species in rare bird breeding in several southern provinces of China, but there have been no studies comparing the gut microbial communities of domestic and wild mallards. In this study, 16S rRNA gene high-throughput sequencing technology was used to compare the composition and diversity of gut microbial communities in domestic and wild mallards. Alpha diversity analysis showed significant differences in gut microbial communities between the two groups of mallards, and the diversity and richness of gut microbial communities were significantly higher in wild mallards than in domestic mallards. Beta diversity analysis showed that the two groups of stool samples were mostly separated on the principal coordinate analysis (PCoA) plot. In domestic mallards, Firmicutes (68.0% ± 26.5%) was the most abundant bacterial phylum, followed by Proteobacteria (24.5% ± 22.9%), Bacteroidetes (3.1% ± 3.2%), Fusobacteria (2.2% ± 5.9%), and Actinobacteria (1.1% ± 1.8%). The dominant bacterial phyla in wild mallards were Firmicutes (79.0% ± 10.2%), Proteobacteria (12.9% ± 9.5%), Fusobacteria (3.4% ± 2.5%), and Bacteroidetes (2.8% ± 2.4%). At the genus level, a total of 10 dominant genera (Streptococcus, Enterococcus, Clostridium, Lactobacillus, Soilbacillus, Bacillus, Acinetobacter, Comamonas, Shigella, and Cetobacterium) with an average relative abundance greater than 1% were detected in the fecal samples of both groups. The average relative abundance of five potential pathogenic genera (Streptococcus, Enterococcus, Acinetobacter, Comamonas, and Shigella) was higher in domestic mallards than in wild mallards. The enrichment of pathogenic bacteria in the intestinal tract of domestic mallards should be of sufficient concern.

Mutations, sex, and genetic diversity: New arguments for partially directed evolution

A review of the theories of the existence of sexes, genetic diversity, and the distribution of mutations among organisms shows that all these concepts are not a product of random evolution and cannot be explained within the framework of Darwinism. Most mutations are the result of the genome acting on itself. This is an organized process that is implemented very differently in different species, in different places in the genome. Because of the fact that it is not random, this process must be directed and regulated, albeit with complex and not fully understood laws. This means that an additional reason must be included in order to model such mutations during evolution. The assumption of directionality must not only be explicitly included in evolutionary theory but must also occupy a central place in it. In this study an updated model of partially directed evolution is constructed, which is capable of qualitatively explaining the indicated features of evolution. Experiments are described that can confirm or disprove the proposed model.

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.

Genome-Wide Population Structure Analysis and Genetic Diversity Detection of Four Chinese Indigenous Duck Breeds from Fujian Province

Simple Summary

The aim of this study was to conduct a genome-wide comparative analysis of four indigenous Chinese duck breeds (Jinding, Liancheng white, Putian black, and Shanma ducks) from Fujian Province, to understand their genetic diversity and population structure. Population parameters showed that the four indigenous breeds were separated groups. Five genomic regions are presented as hotspots of autozygosity among these indigenous duck breeds, with candidate genes involved in muscle growth, pigmentation, and neuroregulation. Genomic information may play a vital role in the improvement of conservation strategies.

Abstract

The assessment of population genetic structure is the basis for understanding the genetic information of indigenous breeds and is important for the protection and management of indigenous breeds. However, the population genetic differentiation of many local breeds still remains unclear. Here, we performed a genome-wide comparative analysis of Jinding, Liancheng white, Putian black, and Shanma ducks based on the genomic sequences using RAD sequencing to understand their population structure and genetic diversity. The population parameters showed that there were obvious genetic differences among the four indigenous breeds, which were separated groups. Among them, Liancheng white and Shanma ducks may come from the same ancestor because the phylogenetic tree forms three tree trunks. In addition, during the runs of homozygosity (ROH), we found that the average inbreeding coefficient of Liancheng white and Putian black ducks was the lowest and the highest, respectively. Five genomic regions were considered to be the hotspots of autozygosity among these indigenous duck breeds, and the candidate genes involved a variety of potential variations, such as muscle growth, pigmentation, and neuroregulation. These findings provide insights into the further improvement and conservation of Fujian duck breeds.

Integrated multi-omic data reveal the potential molecular mechanisms of the nutrition and flavor in Liancheng white duck meat

The Liancheng white (LW) duck is one of the most valued Chinese indigenous poultry breeds. Its meat is rich in nutrients and has distinct flavors, but the molecular mechanisms behind them are unknown. To address this issue, we measured and compared multi-omic data (genome, transcriptome, and metabolome) of breast meat from LW ducks and the Mianyang Shelduck (MS) ducks. We found that the LW duck has distinct breed-specific genetic features, including numerous mutant genes with differential expressions associated with amino acid metabolism and transport activities. The metabolome driven by genetic materials was also seen to differ between the two breeds. For example, several amino acids that are beneficial for human health, such as L-Arginine, L-Ornithine, and L-lysine, were found in considerably higher concentrations in LW muscle than in MS duck muscle (p < 0.05). SLC7A6, a mutant gene, was substantially upregulated in the LW group (p < 0.05), which may lead to excessive L-arginine and L-ornithine accumulation in LW duck meat through transport regulation. Further, guanosine monophosphate (GMP), an umami-tasting molecule, was considerably higher in LW muscle (p < 0.05), while L-Aspartic acid was significantly abundant in MS duck meat (p < 0.05), showing that the LW duck has a different umami formation. Overall, this study contributed to our understanding of the molecular mechanisms driving the enriched nutrients and distinct umami of LW duck meat, which will provide a useful reference for duck breeding.

The complete mitochondrial genome of Anas zonorhyncha (Swinhoe, 1866) (Anatidae: Anas)

The complete mitochondrial genome of Anas zonorhyncha was first reported. The length of the entire mitochondrial genome was 16,605 base pairs, including 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a D-loop region. A phylogenetic tree of A. zonorhyncha was constructed with a group of related species in the family of Anatidae, indicating a close genetic relationship between A. zonorhyncha and A. poecilorhyncha.

Positive selection of skeleton-related genes during duck domestication revealed by whole genome sequencing

BACKGROUND

Domestication alters several phenotypic, neurological, and physiological traits in domestic animals compared to those in their wild ancestors. Domestic ducks originated from mallards, and some studies have shown that spot-billed ducks may have also made minor genetic contributions to domestication. Compared with the two ancestral species, domestic ducks generally differ in body size and bone morphology. In this study, we performed both genomic and transcriptomic analyses to identify candidate genes for elucidating the genetic mechanisms underlying phenotypic variation.

METHODS

In this study, the duck genome data from eight domestic breeds and two wild species were collected to study the genetic changes during domestication. And the transcriptome data of different tissues from wild ducks and seven domestic ducks were used to reveal the expression difference between wild and domestic ducks.

RESULTS

Using fixation index (Fst) algorithm and transcriptome data, we found that the genes related to skeletal development had high Fst values in wild and domestic breeds, and the differentially expressed genes were mainly enriched in the ossification pathway. Our data strongly suggest that the skeletal systems of domestic ducks were changed to adapt to artificial selection for larger sizes. In addition, by combining the genome and transcriptome data, we found that some Fst candidate genes exhibited different expression patterns, and these genes were found to be involved in digestive, immune, and metabolic functions.

CONCLUSIONS

A wide range of phenotypic differences exists between domestic and wild ducks. Through both genome and transcriptome analyses, we found that genes related to the skeletal system in domestic ducks were strongly selected. Our findings provide new insight into duck domestication and selection effects during the domestication.

Whole-genome resequencing reveals genetic characteristics of different duck breeds from the Guangxi region in China

Distinctive indigenous duck (Anas platyrhynchos) populations of Guangxi, China, evolved due to the geographical, cultural, and environmental variability of this region. To investigate the genetic diversity and population structure of the indigenous ducks of Guangxi, 78 individuals from eight populations were collected and sequenced by whole-genome resequencing with an average depth of ∼9.40×. The eight indigenous duck populations included four breeds and four resource populations. Moreover, the genome data of 47 individuals from two typical meat-type breeds and two native egg-type breeds were obtained from a public database. Calculation of heterozygosity, nucleotide diversity (π), Tajima’s D, and F_ST indicated that the Guangxi populations were characterized by higher genetic diversity and lower differentiation than meat-type breeds. The highest diversity was observed in the Xilin-Ma ducks. Principal component, structure, and phylogenetic tree analyses revealed the relationship between the indigenous duck populations of Guangxi. A mild degree of differentiation was observed among the Guangxi populations, although three populations were closer to the meat or egg breeds. Indigenous populations are famous for their special flavor, small body size, and slow growth rates. Selective sweep analysis revealed the candidate genes and pathways associated with these growth traits. Our findings provide a valuable source of information regarding genetic diversity, population conservation, and genome-associated breeding of ducks.

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

Background

Duck is an ancient domesticated animal with high economic value, used for its meat, eggs, and feathers. However, the origin of indigenous Chinese ducks remains elusive. To address this question, we performed whole-genome resequencing to first explore the genetic relationship among variants of these domestic ducks with their potential wild ancestors in eastern China, as well as understand how the their genomes were shaped by different natural and artificial selective pressures.

Results

Here, we report the resequencing of 60 ducks from Chinese spot-billed ducks (Anas zonorhyncha), mallards (Anas platyrhnchos), Fenghua ducks, Shaoxing ducks, Shanma ducks and Cherry Valley Pekin ducks of eastern China (ten from each population) at an average effective sequencing depth of ~ 6× per individual. The results of population and demographic analysis revealed a deep phylogenetic split between wild (Chinese spot-billed ducks and mallards) and domestic ducks. By applying selective sweep analysis, we identified that several candidate genes, important pathways and GO categories associated with artificial selection were functionally related to cellular adhesion, type 2 diabetes, lipid metabolism, the cell cycle, liver cell proliferation, and muscle functioning in domestic ducks.

Conclusion

Genetic structure analysis showed a close genetic relationship of Chinese spot-billed ducks and mallards, which supported that Chinese spot-billed ducks contributed to the breeding of domestic ducks. During the long history of artificial selection, domestic ducks have developed a complex biological adaptation to captivity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07710-2.

Analysis of whole-genome re-sequencing data of ducks reveals a diverse demographic history and extensive gene flow between Southeast/South Asian and Chinese populations

BACKGROUND

The most prolific duck genetic resource in the world is located in Southeast/South Asia but little is known about the domestication and complex histories of these duck populations.

RESULTS

Based on whole-genome resequencing data of 78 ducks (Anas platyrhynchos) and 31 published whole-genome duck sequences, we detected three geographic distinct genetic groups, including local Chinese, wild, and local Southeast/South Asian populations. We inferred the demographic history of these duck populations with different geographical distributions and found that the Chinese and Southeast/South Asian ducks shared similar demographic features. The Chinese domestic ducks experienced the strongest population bottleneck caused by domestication and the last glacial maximum (LGM) period, whereas the Chinese wild ducks experienced a relatively weak bottleneck caused by domestication only. Furthermore, the bottleneck was more severe in the local Southeast/South Asian populations than in the local Chinese populations, which resulted in a smaller effective population size for the former (7100-11,900). We show that extensive gene flow has occurred between the Southeast/South Asian and Chinese populations, and between the Southeast Asian and South Asian populations. Prolonged gene flow was detected between the Guangxi population from China and its neighboring Southeast/South Asian populations. In addition, based on multiple statistical approaches, we identified a genomic region that included three genes (PNPLA8, THAP5, and DNAJB9) on duck chromosome 1 with a high probability of gene flow between the Guangxi and Southeast/South Asian populations. Finally, we detected strong signatures of selection in genes that are involved in signaling pathways of the nervous system development (e.g., ADCYAP1R1 and PDC) and in genes that are associated with morphological traits such as cell growth (e.g., IGF1R).

CONCLUSIONS

Our findings provide valuable information for a better understanding of the domestication and demographic history of the duck, and of the gene flow between local duck populations from Southeast/South Asia and China.

Genome-wide association studies demonstrate that TASP1 contributes to increased muscle fiber diameter

Muscle fiber diameter is an economically important trait because it affects meat yield and quality. However, the genetic basis underlying muscle fiber diameter has not been determined. In this study, we collected THREE muscular histological phenotypes in 479 ducks from an F₂ segregating population generated by mallard × Pekin duck crosses. We performed genome-wide association studies (GWAS) and identified a quantitative trait locus (QTL) significantly associated with muscle fiber diameter on chromosome 3. Then, we discovered the selection signatures using the fixation index among 40 mallards and 30 Pekin ducks in this QTL region. Furthermore, we characterized the recombination event in this QTL region and identified a 6-kb block located on TASP1 that was significantly associated with muscle fiber diameter. Finally, five SNPs were screened as potential causative mutations within the 6-kb block. In conclusion, we demonstrated that TASP1 contributes to an increase in muscle fiber diameter, which helps to characterize muscle development and contributes to the genetic improvement of meat yield and quality in livestock.

Revisiting the evolutionary history of domestic and wild ducks based on genomic analyses

Although domestic ducks have been important poultry species throughout human history, their origin remains enigmatic, with mallards and/or Chinese spot-billed ducks being proposed as the direct wild ancestor(s) of domestic ducks. Here, we analyzed 118 whole genomes from mallard, Chinese spot-billed, and domestic ducks to reconstruct their evolutionary history. We found pervasive introgression patterns among these duck populations. Furthermore, we showed that domestic ducks separated from mallard and Chinese spot-billed ducks nearly 38 thousand years ago (kya) and 54 kya, respectively, which is considerably outside the time period of presumed duck domestication. Thus, our results suggest that domestic ducks may have originated from another wild duck population that is currently undefined or unsampled, rather than from present-day mallard and/or Chinese spot-billed ducks, as previously thought. Overall, this study provides new insight into the complex evolution of ducks.

A joint analysis strategy reveals genetic changes associated with artificial selection between egg-type and meat-type ducks

Egg-type ducks and meat-type ducks are predominantly commercial or indigenous and have been subjected to artificial directional selection. These two duck types differ substantially in body shape, production performance and reproductivity. However, the genetic changes associated with phenotypic differences remain unclear. Here, we compared the two duck types at the genomic and transcriptomic levels. We identified a large number of SNPs and genes in genomic divergent regions in terms of F_ST and θπ values. The corresponding genes were mainly enriched in embryonic development function and metabolic pathway. RNA-seq analysis also revealed differential gene expression in the liver and gonads. The differentially expressed genes were functionally associated with signal transmission and substance metabolism respectively. Furthermore, we found that seven genes were related to differentiation between the two types by both g genome and transcriptome analysis and were plausible candidate genes. These genes were annotated to GO categories of cell development and disease immunity. These findings will enable a better understanding of the artificial selection history of meat and egg ducks and provide a valuable resource for future research on the breeding of these two lineages.

Genetic signature of hybridization between Chinese spot-billed ducks and domesticated ducks

In this study, we analyzed 93 whole genomes from Chinese spot-billed ducks (CSB), meat-type ducks (MET), and egg and dual purpose-type ducks (EDT) to characterize the genetic material flowing between the CSB and modern ducks. Using a frequency of shared identical-by-descent method, approximately 10.9 Mb introgression segments containing 140 genes were identified showing the signatures of introgression between CSB and EDT. Meanwhile, nearly 10.6 M introgression regions containing 149 genes were identified between CSB and MET. Based on the haplotypes tree of each segment, we found that the introgression between CSB and domesticated ducks was asymmetric with a high level of gene flow from domestic to CSB and a low level of migration in the opposite direction. Moreover, we identified several genes that were introgressions from CSB and showed the signature of positive selection, which may contribute to the breeding of modern ducks. Our results provide new insight into the evolution and breeding history of domestic ducks and may be useful for the future management of wild and domestic duck populations.

Population genomic transformations induced by isolation of wild bird avian influenza viruses (Orthomyxoviridae) in embryonated chicken eggs

Isolation and cultivation of wild-type viruses in model organism cells or tissues is standard practice in virology. Oftentimes, the virus host species is distantly related to the species from which the culture system was developed. Thus, virus culture in these tissues and cells basically constitutes a host jump, which can lead to genomic changes through genetic drift and/or adaptation to the culture system. We directly sequenced 70 avian influenza virus (Orthomyxoviridae) genomes from oropharyngeal/cloacal swabs collected from wild bird species and paired virus isolates propagated from the same samples following isolation in specific-pathogen-free embryonated chicken eggs. The data were analyzed using population genetic approaches including evaluation of single nucleotide polymorphism (SNP) frequencies and divergence with pooled-sequencing analyses, consensus sequence placement in neighbor-joining trees, and haplotype reconstruction and networks. We found that propagation of virus in eggs leads to skewed SNP mutation spectra with some SNPs going to fixation. Both synonymous and nonsynonmous SNP frequencies shifted. We found multiple consensus sequences that differed between the swabs and the isolates, with some sequences from the same sample falling into divergent genetic clusters. Twenty of 23 coinfections detected had different dominant subtypes following virus isolation, thus sequences from both the swab and isolate were needed to obtain full subtype data. Haplotype networks revealed haplotype frequency shifts and the appearance or loss of low-frequency haplotypes following isolation. The results from this study revealed that isolation of wild bird avian influenza viruses in chicken eggs leads to skewed populations that are different than the input populations. Consensus sequence changes from virus isolation can lead to flawed phylogenetic inferences, and subtype detection is biased. These results suggest that for genomic studies of wild bird influenza viruses the biological field should move away from chicken egg isolation towards directly sequencing the virus from host samples.

Understanding the Genetic Domestication History of the Jianchang Duck by Genotyping and Sequencing of Genomic Genes Under Selection

The Jianchang duck is mainly distributed in Southwest China, and has the characteristics of fast growth rate and strong abilities in lipid deposition in the liver. In order to investigate the effects of domestication process on formation of the unique characteristics of Jianchang duck, the whole genome of sixteen individuals and three pooling of Jianchang duck were re-sequenced, and genome data of 70 mallards and 83 domestic ducks from thirteen different places in China were obtained from NCBI. The population stratification and evolution analysis showed gene exchanges existed between the Jianchang and other domestic duck populations, as well as Jianchang ducks and mallards. Genomic comparison between mallards and Jianchang ducks showed genes, including CNTN1, CHRNA9, and SHANK2, which is involved in brain and nerve development, experienced strong positive selection in the process of Jianchang duck domestication. The genomic comparison between Jianchang and domestic duck populations showed that HSD17B12 and ESM1, which affect lipid metabolism, experienced strong positive selection during the domestication process. FST analysis among populations of Jianchang duck with different plumage colors indicated that MITF was related to the phenotype of a white feather, while MC1R was related to the phenotype of hemp feather. Our results provided a base for the domestication process of Jianchang duck and the genomic genes for unique traits.

Genetic characteristics of Jiaji Duck by whole genome re-sequencing

Jiaji Duck (JJ) is a Muscovy duck species that possesses many superior characteristics, and it has become an important genetic resource in China. However, to date, its genetic characteristics and genetic relationship with other duck breeds have not been explored yet, which greatly limits the utilization of JJ. In the present study, we investigated the genome sequences of 15 individual ducks representing five different duck populations, including JJ, French Muscovy duck (FF), mallard (YD), hong duck (HD) and Beijing duck (BD). Moreover, we investigated the characteristics of JJ-specific single nucleotide polymorphisms (SNPs) and compared the genome sequences of JJ vs. YD and JJ vs. BD using integrated strategies, including mutation detection, selective screening, and Gene Ontology (GO) analysis. More than 40 Gb of clean data were obtained for each population (mean coverage of 13.46 Gb per individual). A total number of 22,481,367 SNPs and 4,156,829 small insertion-deletions (Indels) were identified for the five duck populations, which could be used as molecular markers in breeding and utilization of JJ. Moreover, we identified 1,447,932 JJ-specific SNPs, and found that genes covering at least one JJ-specific SNP mainly involved in protein phosphorylation and dephosphorylation, as well as DNA modification. Phylogenetic tree and principal components analysis (PCA) revealed that the genetic relationship of JJ was closest to FF, while it was farthest to BD. A total of 120 and 111 genes were identified as positive selection genes for JJ vs. BD and JJ vs. YD, respectively. GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the positive selection genes for JJ vs. BD ducks mainly involved in pigmentation, muscle contraction and stretch, gland secretion, and immunology, while the positive selection genes obtained from JJ vs. YD ducks mainly involved in embryo development, muscle contraction and stretch, and gland secretion. Taken together, our findings enabled us to better understand the characteristics of JJ and provided a molecular basis for the breeding and hybrid utilization of JJ in the future.

Analysis of Anasplatyrhynchos genome resequencing data reveals genetic signatures of artificial selection

Three artificially selected duck populations (AS), higher lean meat ratios (LTPD), higher fat ratios (FTPD) and higher quality meat (CMD), have been developed in China, providing excellent populations for investigation of artificial selection effects. However, the genetic signatures of artificial selection are unclear. In this study, we sequenced the genome sequences of these three artificially selected populations and their ancestral population (mallard, M). We then compared the genome sequences between AS and M and between LTPD and FTPD using integrated strategies such as anchoring scaffolds to pseudo-chromosomes, mutation detection, selective screening, GO analysis, qRT-PCR, and protein multiple sequences alignment to uncover genetic signatures of selection. We anchored duck scaffolds to pseudo-chromosomes and obtained 28 pseudo-chromosomes, accounting for 84% of duck genome in length. Totally 78 and 99 genes were found to be under selection between AS and M and between LTPD and FTPD. Genes under selection between AS and M mainly involved in pigmentation and heart rates, while genes under selection between LTPD and FTPD involved in muscle development and fat deposition. A heart rate regulator (HCN1), the strongest selected gene between AS and M, harbored a GC deletion in AS and displayed higher mRNA expression level in M than in AS. IGF2R, a regulator of skeletal muscle mass, was found to be under selection between FTPD and LTPD. We also found two nonsynonymous substitutions in IGF2R, which might lead to higher IGF2R mRNA expression level in FTPD than LTPD, indicating the two nonsynonymous substitutions might play a key role for the regulation of duck skeletal muscle mass. Taken together, these results of this study provide valuable insight for the genetic basis of duck artificial selection.

An intercross population study reveals genes associated with body size and plumage color in ducks

Comparative population genomics offers an opportunity to discover the signatures of artificial selection during animal domestication, however, their function cannot be directly revealed. We discover the selection signatures using genome-wide comparisons among 40 mallards, 36 indigenous-breed ducks, and 30 Pekin ducks. Then, the phenotypes are fine-mapped based on resequencing of 1026 ducks from an F₂ segregating population generated by wild × domestic crosses. Interestingly, the two key economic traits of Pekin duck are associated with two selective sweeps with fixed mutations. A novel intronic insertion most possibly leads to a splicing change in MITF accounted for white duck down feathers. And a putative long-distance regulatory mutation causes continuous expression of the IGF2BP1 gene after birth which increases body size by 15% and feed efficiency by 6%. This study provides new insights into genotype-phenotype associations in animal research and constitutes a promising resource on economically important genes in fowl.

Evolutionary history of Mexican domesticated and wild Meleagris gallopavo

Background

The distribution of the wild turkey (Meleagris gallopavo) extends from Mexico to southeastern Canada and to the eastern and southern regions of the USA. Six subspecies have been described based on morphological characteristics and/or geographical variations in wild and domesticated populations. In this paper, based on DNA sequence data from the mitochondrial D-loop, we investigated the genetic diversity and structure, genealogical relationships, divergence time and demographic history of M. gallopavo populations including domesticated individuals.

Results

Analyses of 612 wild and domesticated turkey mitochondrial D-loop sequences, including 187 that were collected for this study and 425 from databases, revealed 64 haplotypes with few mutations, some of which are shared between domesticated and wild turkeys. We found a high level of haplotype and nucleotide diversity, which suggests that the total population of this species is large and stable with an old evolutionary history. The results of genetic differentiation, haplotype network, and genealogical relationships analyses revealed three main genetic groups within the species: mexicana as a population relict (C1), merriami (C2), and mexicana/intermedia/silvestris/osceola (C3). Haplotypes detected in domesticated turkeys belong to group C3. Estimates of divergence times agree with range expansion and diversification events of the relict population of M. gallopavo in northwestern Mexico during the Pliocene–Pleistocene and Pleistocene–Holocene boundaries. Demographic reconstruction showed that an expansion of the population occurred 110,000 to 130,000 years ago (Kya), followed by a stable period 100 Kya and finally a decline ~ 10 Kya (Pleistocene–Holocene boundary). In Mexico, the Trans-Mexican Volcanic Belt may be responsible for the range expansion of the C3 group. Two haplotypes with different divergence times, MGMDgoB/MICH1 and MICH2, are dominant in domesticated and commercial turkeys.

Conclusions

During the Pleistocene, a large and stable population of M. gallopavo covered a wide geographic distribution from the north to the center of America (USA and Mexico). The mexicana, merriami, and mexicana/intermedia/silvestris/osceola genetic groups originated after divergence and range expansion from northwestern Mexico during the Pliocene–Pleistocene and Pleistocene–Holocene boundaries. Old and new maternal lines of the mexicana/intermedia/silvestris/osceola genetic group were distributed within the Trans-Mexican Volcanic Belt where individuals were captured for domestication. Two haplotypes are the main founder maternal lines of domesticated turkeys.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0388-8) contains supplementary material, which is available to authorized users.

The genetic diversity of chicken breeds from Jiangxi, assessed with BCDO2 and the complete mitochondrial DNA D-loop region

The Jiangxi Province of China has numerous native domestic chicken breeds, including some black skin breeds. The genetic diversity of Jiangxi native chickens is largely unknown, and specifically, the genetic contribution of the grey junglefowl to black skin chickens is not well understood. To address these questions, the complete D-loop region of the mitochondrial DNA (mtDNA) and beta-carotene dioxygenase 2(BCDO2)gene was sequenced in a total of 209 chickens representing seven Jiangxi native breeds. Thirty-one polymorphic sites were identified across the complete mtDNA D-loop region sequence. Twenty-three haplotypes were observed in the seven breeds, which belonged to four distinct mitochondrial clades (A, B, C and E). Clade A and B were dominant in the chickens with a frequency of approximately 67.9%. There were five SNPs that defined two haplotypes, W and Y in BCDO2. Four breeds had one haplotype and three breeds had two. We conclude that Jiangxi native chicken breeds have relatively low genetic diversity and likely share four common maternal lineages from two different maternal ancestors of junglefowl. Furthermore, some Jiangxi chicken populations may have been mixed with chickens with exotic lineage. Further research should be established to protect these domestic chicken resources.

PCK1 expression is correlated with the plasma glucose level in the duck

Phosphoenolpyruvate carboxykinase 1 (soluble) (PCK1) is a key gene in gluconeogenesis and glyceroneogenesis. Although its functions have been extensively studied in mice, bats and humans, little is known in ducks. Here, PCK1 functions were studied using a duck domestication model and a 48-h fasting experiment. We found PCK1 expression significantly decreased in two breeds of domestic ducks (Jinyun Pockmark ducks and Cherry Valley ducks) as compared with wild ducks (Anas platyrhynchos). Simultaneously, plasma levels of glucose, triglycerides and free fatty acid in domestic ducks were lower than in wild ducks. When compared with fed ducks, the plasma triglyceride level was observed to be significantly decreased, while the glucose and free fatty acid levels remained constant in 48-h fasting ducks. The expression analysis of gluconeogenic genes revealed that fructose-1,6-bisphosphatase genes (FBP1 and FBP2) and the glucose-6-phosphatase gene (G6PC2) were not changed, whereas PCK1 was significantly upregulated. In addition, the reported regulators of PCK1, including forkhead box A2 (FOXA2) gene and orphan nuclear receptor NR4A family genes (NR4A1, NR4A2 and NR4A3), exhibited similar expression levels between 48-h fasting ducks and fed ducks, suggesting that PCK1 is not regulated by these genes in the duck under fasting conditions. In conclusion, PCK1 expression may affect plasma levels of glucose, triglycerides and free fatty acid during the duck domestication process. This work demonstrates for the first time in duck that PCK1 is a key gene in maintaining plasma glucose homeostasis during fasting and that the upregulated expression of PCK1 may be responsible for constant plasma free fatty acid level by the glyceroneogenesis process.

SNP discovery and genotyping using Genotyping-by-Sequencing in Pekin ducks

Genomic selection and genome-wide association studies need thousands to millions of SNPs. However, many non-model species do not have reference chips for detecting variation. Our goal was to develop and validate an inexpensive but effective method for detecting SNP variation. Genotyping by sequencing (GBS) can be a highly efficient strategy for genome-wide SNP detection, as an alternative to microarray chips. Here, we developed a GBS protocol for ducks and tested it to genotype 49 Pekin ducks. A total of 169,209 SNPs were identified from all animals, with a mean of 55,920 SNPs per individual. The average SNP density reached 1156 SNPs/MB. In this study, the first application of GBS to ducks, we demonstrate the power and simplicity of this method. GBS can be used for genetic studies in to provide an effective method for genome-wide SNP discovery.

Transcriptome analysis of adiposity in domestic ducks by transcriptomic comparison with their wild counterparts

Excessive adiposity is a major problem in the duck industry, but its molecular mechanisms remain unknown. Genetic comparisons between domestic and wild animals have contributed to the exploration of genetic mechanisms responsible for many phenotypic traits. Significant differences in body fat mass have been detected between domestic and wild ducks. In this study, we used the Peking duck and Anas platyrhynchos as the domestic breed and wild counterpart respectively and performed a transcriptomic comparison of abdominal fat between the two breeds to comprehensively analyze the transcriptome basis of adiposity in ducks. We obtained approximately 350 million clean reads; assembled 61 250 transcripts, including 23 699 novel ones; and identified alternative 5' splice sites, alternative 3' splice sites, skipped exons and retained intron as the main alternative splicing events. A differential expression analysis between the two breeds showed that 753 genes exhibited differential expression. In Peking ducks, some lipid metabolism-related genes (IGF2, FABP5, BMP7, etc.) and oncogenes (RRM2, AURKA, CYR61, etc.) were upregulated, whereas genes related to tumor suppression and immunity (TNFRSF19, TNFAIP6, IGSF21, NCF1, etc.) were downregulated, suggesting adiposity might closely associate with tumorigenesis in ducks. Furthermore, 280 576 single-nucleotide variations were found differentiated between the two breeds, including 8641 non-synonymous ones, and some of the non-synonymous ones were found enriched in genes involved in lipid-associated and immune-associated pathways, suggesting abdominal fat of the duck undertakes both a metabolic function and immune-related function. These datasets enlarge our genetic information of ducks and provide valuable resources for analyzing mechanisms underlying adiposity in ducks.

The investigation of genetic diversity and evolution of Daweishan Mini chicken based on the complete mitochondrial (mt)DNA D-loop region sequence

This study evaluated the genetic diversity and origin of Daweishan Mini chickens using mtDNA sequence polymorphism. Blood samples from 30 Daweishan Mini chickens were collected. The complete D-loop was PCR amplified, sequenced and compared with the DNA data of five Red Junglefowl (Gallus gallus) subspecies. Eighteen variable sites that defined six haplotypes were observed. The six haplotypes were clustered into four clades (A, B, D and E), of which clade A and B were dominant. Clades Aand B were clustered with G.g. spadiceus, indicating these two clades may have originated from this subspecies. These results show there is diversity in the middle of the mtDNA D-loop, and indicate there are multiple maternal origins for Daweishan Mini chickens. It appears that G.g. spadiceus contributed more to the evolution of the Daweishan Mini chickens breed than the other four subspecies tested here.

Phylogenetic analysis between domestic and wild duck species in Korea using mtDNA D-loop sequences

Recently, the consumption of duck meat has increased; therefore, we need to reveal the origin and gene flow of domestic ducks in Korea. In order to discriminate between duck species, D-loop variations in mitochondrial DNA (mtDNA) have been investigated. In this study, 45 individuals from seven species of wild and domestic ducks in Korea were considered for the D-loop region sequences. With the participation of all the sequences, a phylogenetic neighbor-joining tree was constructed to differentiate between the wild and domestic duck species. In consideration of these sequences, a total 66 haplotypes were obtained (indel included) with an average haplotype of 76.9%, and a haplotype and nucleotide diversity of 0.91 and 0.01, respectively. Also, an estimation of the sequence divergence within and between species was measured in 0.045 and 0.013-0.095, respectively. Meanwhile, the lowest distances of 0.024, 0.013 and 0.018 were observed in three species, including the Mallard, Spot-billed and domestic duck, respectively, which have relatively close genetic relationships. All haplotypes were used for the median-joining network analysis to differentiate all duck species, while three duck species were closely related. Moreover, 26 indel polymorphisms were identified which could be used for the discrimination among the duck species. Based on our results, duck species were effectively discriminated in a D-loop region, which could then be used for an appropriate genetic conservation program for the wild duck and domestic duck breeds in Korea.

Genetic structure of Eurasian and North American mallard ducks based on mtDNA data

To elucidate the origin and genetic structure of the domesticated duck in Eurasia and North America, we sequenced 114 duck D-loop sequences and retrieved 489 D-loop sequences from GenBank. In total, 603 ducks including 50 duck breeds/populations from eight countries (China, France, Russia, India, Kazakhstan, Mongolia, Thailand and USA) were used in this study. One hundred and thirty-four haplotypes and 81 variable sites were detected. H49 was the predominant haplotype, which was considered to be the same dominant haplotype found in the previous studies, and was found in 309 birds. The smallest values for both genetic differentiation index (F(ST), 0.04156) and the number of the net nucleotide substitutions between two populations (D(A), 0.00018) were observed between Eurasian domestic ducks and Eurasian mallards. No geography, breed or population clusters were observed in the Eurasian domestic ducks and mallards. Five haplotypes were shared by USA mallards and Eurasian domestic duck/Eurasian mallards. Only one haplotype (H49) was shared by Eurasian domestic ducks and China spot-billed ducks. By combining phylogenetic analyses, haplotype network profile, genetic distances and shared haplotypes, we can draw two major conclusions: (i) Eurasian and North American mallards show a clear geographic distribution pattern; (ii) Eurasian domestic ducks are derived from the Eurasian mallards, not from the spot-billed ducks.