Pekin duck products have become popular in Asia over recent decades and account for an increasing market share. However, the genetic mechanisms affecting carcass growth in Pekin ducks remain unknown. This study aimed to identify quantitative trait loci affecting body size and carcass yields in Pekin ducks. We measured 18 carcass traits in 639 Pekin ducks and performed genotyping using genotyping-by-sequencing (GBS). Loci-based association analysis detected 37 significant loci for the 17 traits. Thirty-seven identified candidate genes were involved in many biological processes. One single nucleotide polymorphism (SNP) (Chr1_140105435 A > T) located in the intron of the ATPase phospholipid transporting 11A gene (ATP11A) attained genome-wide significance associated with five weight traits. Eight SNPs were significantly associated with three body size traits, including the candidate gene plexin domain containing 2 (PLXDC2) associated with breast width and tensin 3 (TNS3) associated with fossil bone length. Only two SNPs were significantly associated with foot weight and four SNPs were significantly associated with heart weight. In the gene-based analysis, three genes (LOC101791418, TUBGCP3 (encoding tubulin gamma complex-associated protein 3), and ATP11A) were associated with four traits (42-day body weight, eviscerated weight, half-eviscerated weight, and leg muscle weight percentage). However, no loci were significantly associated with leg muscle weight in this study. The novel results of this study improve our understanding of the genetic mechanisms regulating body growth in ducks and thus provide a genetic basis for breeding programs aimed at maximizing the economic potential of Pekin ducks.

Domestic ducks are raised for meat, eggs and feather down, and almost all varieties are descended from the Mallard (Anas platyrhynchos). Here, we report chromosome-level high-quality genome assemblies for meat and laying duck breeds, and the Mallard. Our new genomic databases contain annotations for thousands of new protein-coding genes and recover a major percentage of the presumed “missing genes” in birds. We obtain the entire genomic sequences for the C-type lectin (CTL) family members that regulate eggshell biomineralization. Our population and comparative genomics analyses provide more than 36 million sequence variants between duck populations. Furthermore, a mutant cell line allows confirmation of the predicted anti-adipogenic function of NR2F2 in the duck, and uncovered mutations specific to Pekin duck that potentially affect adipose deposition. Our study provides insights into avian evolution and the genetics of oviparity, and will be a rich resource for the future genetic improvement of commercial traits in the duck.

Backfat thickness is strongly associated with meat quality, fattening efficiency, reproductive performance, and immunity in pigs. Fat storage and fatty acid synthesis mainly occur in adipose tissue. Therefore, we used a high-throughput massively parallel sequencing approach to identify transcriptomes in adipose tissue, and whole-genome differences from three full-sibling pairs of pigs with opposite (high and low) backfat thickness phenotypes. We obtained an average of 38.69 million reads for six samples, 78.68% of which were annotated in the reference genome. Eighty-nine overlapping differentially expressed genes were identified among the three pair comparisons. Whole-genome resequencing also detected multiple genetic variations between the pools of DNA from the two groups. Compared with the animal quantitative trait loci (QTL) database, 20 differentially expressed genes were matched to the QTLs associated with fatness in pigs. Our technique of integrating transcriptome, whole-genome resequencing, and QTL database information provided a rich source of important differentially expressed genes and variations. Associate analysis between selected SNPs and backfat thickness revealed that two SNPs and one haplotype of ME1 significantly affected fat deposition in pigs. Moreover, genetic analysis confirmed that variations in the differentially expressed genes may affect fat deposition.

Mitochondrial remodeling during the peri-implantation stage is the hallmark event essential for normal embryogenesis. Among the changes, enhanced oxidative phosphorylation is critical for supporting high energy demands of postimplantation embryos, but increases mitochondrial oxidative stress, which in turn threatens mitochondrial DNA (mtDNA) stability. However, how mitochondria protect their own histone-lacking mtDNA, during this stage remains unclear. Concurrently, the mitochondrial genome gain DNA methylation by this stage. Its spatiotemporal coincidence with enhanced mitochondrial stress led us to ask if mtDNA methylation has a role in maintaining mitochondrial genome stability. Herein, we report that mitochondrial genome undergoes de novo mtDNA methylation that can protect mtDNA against enhanced oxidative damage during the peri-implantation window. Mitochondrial genome gains extensive mtDNA methylation during transition from blastocysts to postimplantation embryos, thus establishing relatively hypermethylated mtDNA from hypomethylated state in blastocysts. Mechanistic study revealed that DNA methyltransferase 3A (DNMT3A) and DNMT3B enter mitochondria during this process and bind to mtDNA, via their unique mitochondrial targeting sequences. Importantly, loss- and gain-of-function analyses indicated that DNMT3A and DNMT3B are responsible for catalyzing de novo mtDNA methylation, in a synergistic manner. Finally, we proved, in vivo and in vitro, that increased mtDNA methylation functions to protect mitochondrial genome against mtDNA damage induced by increased mitochondrial oxidative stress. Together, we reveal mtDNA methylation dynamics and its underlying mechanism during the critical developmental window. We also provide the functional link between mitochondrial epigenetic remodeling and metabolic changes, which reveals a role for nuclear-mitochondrial crosstalk in establishing mitoepigenetics and maintaining mitochondrial homeostasis.