Mitochondrial genome undergoes de novo DNA methylation that protects mtDNA against oxidative damage during the peri-implantation window.

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.

Genome-wide association analysis identifies candidate genes for carcass yields in Peking ducks.

Poultry meat, particularly Peking ducks, holds a significant global market share, prized for their high meat yield and fat content. However, understanding of the molecular genetic mechanisms influencing carcass yield in ducks is limited. This research aims to use genome-wide association analysis to uncover single-nucleotide polymorphisms influencing carcass yield in Peking ducks, followed by identifying candidate genes linked to carcass traits. In this study, we analyzed seven traits of 643 Peking ducks at age 42 days and identified novel loci associated with these traits. A total of 35 significant loci were detected, with eight SNPs reaching genome-wide significance. KIF20B, AGBL5, SGSM1, MRO, PLAG1, XKR4, and TGS1 were considered as important candidate genes influencing carcass yield in ducks. This study adds to the list of genes affecting Peking duck body traits, aiding marker-assisted breeding and enhancing economic yield.

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.

Genome-wide association study identified the candidate genes associated with angel wing trait in Pekin duck.

Angel wing is a developmental wing deformity that can influence breeding and reproduction in the commercial duck industry. The nutrition foundation of angel wing trait was initially explored, but the genetic basic remains poorly understood. In this study, we identified candidate genes and single-nucleotide polymorphisms (SNPs) associated with angel wing trait in Pekin ducks using a genome-wide association study (GWAS) and selective sweep analysis. The GWAS results showed that nine SNPs across five chromosomes were significantly correlated with the angel wing trait. In total, 468 selection signals were shown between the angel wing ducks and normal ducks, and these signals harbored 154 genes, which were enriched in the nervous system and metabolism. This study provides the new insights into the genetic factors that may influence duck angel wing.

Chromosome-level genome assembly of the Muscovy duck provides insight into fatty liver susceptibility.

The Muscovy duck (Cairina moschata) is an economically important poultry species, which is susceptible to fatty liver. Thus, the Muscovy duck may serve as an excellent candidate animal model of non-alcoholic fatty liver disease. However, the mechanisms underlying fatty liver development in this species are poorly understood. In this study, we report a chromosome-level genome assembly of the Muscovy duck, with a contig N50 of 11.8 Mb and scaffold N50 of 83.16 Mb. The susceptibility of Muscovy duck to fatty liver was mainly attributed to weak lipid catabolism capabilities (fatty acid ß-oxidation and lipolysis). Furthermore, conserved noncoding elements (CNEs) showing accelerated evolution contributed to fatty liver formation by down-regulating the expression of genes involved in hepatic lipid catabolism. We propose that the susceptibility of Muscovy duck to fatty liver is an evolutionary by-product. In conclusion, this study revealed the potential mechanisms underlying the susceptibility of Muscovy duck to fatty liver.

Genome-wide association study reveals novel loci associated with feeding behavior in Pekin ducks.

BACKGROUND: Feeding behavior traits are an essential part of livestock production. However, the genetic base of feeding behavior traits remains unclear in Pekin ducks. This study aimed to determine novel loci related to feeding behavior in Pekin ducks. RESULTS: In this study, the feeding information of 540 Pekin ducks was recorded, and individual genotype was evaluated using genotyping-by-sequencing methods. Genome-wide association analysis (GWAS) was conducted for feeding behavior traits. Overall, thirty significant (P-value < 4.74E-06) SNPs for feeding behavior traits were discovered, and four of them reached the genome-wide significance level (P-value < 2.37E-07). One genome-wide significance locus associated with daily meal times was located in a 122.25 Mb region on chromosome 2, which was within the intron of gene ubiquitin-conjugating enzyme E2 E2 (UBE2E2), and could explain 2.64% of the phenotypic variation. This locus was also significantly associated with meal feed intake, and explained 2.72% of this phenotypic variation. CONCLUSIONS: This study is the first GWAS for feeding behavior traits in ducks. Our results provide a list of candidate genes associated with feeding behavior, and also help to better understand the genetic mechanisms of feeding behavior patterns in ducks.

The transcriptome landscapes of ovary and three oviduct segments during chicken (Gallus gallus) egg formation.

The avian embryo develops within a specialized biological container (eggshell) that contains crucial nutritional compartments (albumen, yolk). We analyzed the transcriptome of ovary and three segments of oviduct, including magnum, isthmus and uterus in the chicken during egg formation. RNA-Seq libraries (42 in total) for ovary and three different parts of the oviduct were sequenced for two different phases of egg formation. We obtained 8365 novel transcripts with an mRNA length longer than 200 bp; of these, 6832 were long intergenic non-coding RNA transcripts. We identified 547 differentially expressed genes in magnum (actively secreting albumen versus inactive) and 585 in uterus (active eggshell calcification versus quiescent). By combining QTL, transcriptome and proteome data, we obtained high quality gene lists for chicken egg formation. This is the first study to describe the ovary and oviduct transcriptomes by mRNA sequencing, and to elucidate the global repertoire of functional genes involved in egg formation.

Genome-wide association study of the level of blood components in Pekin ducks.

Blood components are considered to reflect nutrient metabolism and immune activity in both humans and animals. In this study, we measured 12 blood components in Pekin ducks and performed genome-wide association analysis to identify the QTLs (quantitative trait locus) using a genotyping-by-sequencing strategy. A total of 54 QTLs were identified for blood components. One genome-wide significant QTL for alkaline phosphatase was identified within the intron-region of the OTOG gene (P = 1.31E-07). Moreover, 21 genome-wide significant SNPs for the level of serum cholinesterase were identified on six different scaffolds. In addition, for serum calcium, one genome-wide significant QTL was identified in the upstream region of gene RAB11B. These results provide new markers for functional studies in Pekin ducks, and several candidate genes were identified, which may provide additional insights into specific mechanisms for blood metabolism in ducks and their potential application for duck breeding programs.

Genome-wide association study of bone quality and feed efficiency-related traits in Pekin ducks.

Feeding and bone traits are vital for breeding and reproduction in the commercial duck industry. In this study, we performed a genome-wide association study for feeding and bone traits in a population of 540 lean-type Pekin ducks, followed by genotyping-by-sequencing procedures. The genetic parameters of feeding and bone traits were also estimated using genomic information. In total, seventy-eight significant SNPs were determined, and eleven of them reached the genome-wide significant level for 7 traits except for body weight at 42-day old. A peak of genome-wide significant SNPs was detected on chromosome 2 for feed conversion ratio (P-value = 7.46E-11), and the top SNP (P-value = 2.23E-08) for bone-breaking strength was also obtained in the upstream of gene RAPGEF5. This study provided a list of novel markers and candidate genes associated with feeding and bone traits in Pekin ducks, which could contribute to the genetic selection in duck breeding.

Integrating transcriptome, proteome and QTL data to discover functionally important genes for duck eggshell and albumen formation.

Duck egg quality improvement is an essential target for Asian poultry breeding. In total, 15 RNA-Seq libraries (magnum, isthmus, and uterus at two different physiological states) were sequenced from 48 weeks old Pekin ducks. De novo assembly and annotation methods were utilized to generate new reference transcripts. Our results revealed that 1264 and 2517 genes were differentially expressed in magnum and uterus in the presence versus absence of an egg, respectively. We identified 1089 genes that were differentially expressed in isthmus compared to uterus (in both presence and absence of a calcifying egg). We observed that 11 common DEGs were detected in the egg white proteomes of 6 different bird species including domestic Chicken, Duck, Goose, Turkey, Quail, and Pigeon. On the other hand, only one of the top five most highly expressed genes in duck isthmus was in this category for the chicken isthmus (SPINK7). Among the large number of DEGs during eggshell formation in ducks, only 41 genes showed a similar differential expression pattern in both duck and chicken. By combining chicken QTL database, chicken oviduct transcriptome and egg proteome data for five bird species, we have obtained high-quality gene lists for egg formation. This is the first study to elucidate the transcriptomic changes in different duck oviduct segments during egg formation, and to integrate QTL, proteome and transcriptome data to probe the functional genes associated with albumen secretion and eggshell mineralization.

iTRAQ-Based Quantitative Proteomic Analysis of Duck Eggshell During Biomineralization.

The avian egg is a valuable model for the calcitic biomineralization process as it is the fastest calcification process occurring in nature and is a clear example of biomineralization. In this study, iTRAQ MS/MS is used to detect and study for the first time: 1) the overall duck eggshell proteome; 2) regional differences in the proteome between the inner and outer portions of the duck eggshell. The new reference protein datasets allow us to identify 179 more eggshell proteins than solely using the current release of Ensembl duck annotations. In total, 484 proteins are identified in the entire duck eggshell proteome. Twenty-eight novel proteins of unknown function that are involved in eggshell formation are also identified. Among the identified eggshell proteins, 54 proteins show differential abundances between the inner, partially mineralized eggshell (obtained 16 h after ovulation) compared to the overall complete eggshell (normally expulsed eggshell). At least 64 of the abundant matrix proteins are common to eggshell of 4 different domesticated bird species (chicken, duck, quail, turkey) and zebra finch. This study provides a new resource for avian eggshell proteomics, and augments the inventory of eggshell matrix proteins that will lead to a deeper understanding of calcitic biomineralization.

Dynamics of transcriptome changes during subcutaneous preadipocyte differentiation in ducks.

BACKGROUND: Pekin duck is an important animal model for its ability for fat synthesis and deposition. However, transcriptional dynamic regulation of adipose differentiation driven by complex signal cascades remains largely unexplored in this model. This study aimed to explore adipogenic transcriptional dynamics before (proliferation) and after (differentiation) initial preadipocyte differentiation in ducks. RESULTS: Exogenous oleic acid alone successfully induced duck subcutaneous preadipocyte differentiation. We explored 36 mRNA-seq libraries in order to study transcriptome dynamics during proliferation and differentiation processes at 6 time points. Using robust statistical analysis, we identified 845, 652, 359, 2401 and 1933 genes differentially expressed between -48 h and 0 h, 0 h and 12 h, 12 h and 24 h, 24 h and 48 h, 48 h and 72 h, respectively (FDR < 0.05, FC > 1.5). At the proliferation stage, proliferation related pathways and basic cellular and metabolic processes were inhibited, while regulatory factors that initiate differentiation enter the ready-to-activate state, which provides a precondition for initiating adipose differentiation. According to weighted gene co-expression network analysis, pathways positively related to adipogenic differentiation are significantly activated at the differentiation stage, while WNT, FOXO and other pathways that inhibit preadipocyte differentiation are negatively regulated. Moreover, we identified and classified more than 100 transcription factors that showed significant changes during differentiation, and found novel transcription factors that were not reported to be related to preadipoctye differentiation. Finally, we manually assembled a proposed regulation network model of subcutaneous preadipocyte differentiation base on the expression data, and suggested that E2F1 may serve as an important link between the processes of duck subcutaneous preadipocyte proliferation and differentiation. CONCLUSIONS: For the first time we comprehensively analyzed the transcriptome dynamics of duck subcutaneous preadipocyte proliferation and differentiation. The current study provides a solid basis for understanding the synthesis and deposition of subcutaneous fat in ducks. Furthermore, the information generated will allow future investigations of specific genes involved in particular stages of duck adipogenesis.

Genome-wide association study reveals novel loci associated with body size and carcass yields in Pekin ducks.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Revisiting avian 'missing' genes from de novo assembled transcripts.

BACKGROUND: Argument remains as to whether birds have lost genes compared with mammals and non-avian vertebrates during speciation. High quality-reference gene sets are necessary for precisely evaluating gene gain and loss. It is essential to explore new reference transcripts from large-scale de novo assembled transcriptomes to recover the potential hidden genes in avian genomes. RESULTS: We explored 196 high quality transcriptomic datasets from five bird species to reconstruct transcripts for the purpose of discovering potential hidden genes in the avian genomes. We constructed a relatively complete and high-quality bird transcript database (1,623,045 transcripts after quality control in five birds) from a large amount of avian transcriptomic data, and found most of the presumed missing genes (83.2%) could be recovered in at least one bird species. Most of these genes have been identified for the first time in birds. Our results demonstrate that 67.94% genes have GC content over 50%, while 2.91% genes are AT-rich (AT% > 60%). In our results, 239 (53.59%) genes had a tissue-specific expression index of more than 0.9 in chicken. The missing genes also have lower Ka/Ks values than average (genome-wide: Ka/Ks = 0.99; missing gene: Ka/Ks = 0.90; t-test = 1.25E-14). Among all presumed missing genes, there were 135 for which we did not find any meaningful orthologues in any of the 5 species studied. CONCLUSION: Insufficient reference genome quality is the major reason for wrongly inferring missing genes in birds. Those presumably missing genes often have a very strong tissue-specific expression pattern. We show multi-tissue transcriptomic data from various species are necessary for inferring gene family evolution for species with only draft reference genomes.

Impaired imprinted X chromosome inactivation is responsible for the skewed sex ratio following in vitro fertilization.

Dynamic epigenetic reprogramming occurs during normal embryonic development at the preimplantation stage. Erroneous epigenetic modifications due to environmental perturbations such as manipulation and culture of embryos during in vitro fertilization (IVF) are linked to various short- or long-term consequences. Among these, the skewed sex ratio, an indicator of reproductive hazards, was reported in bovine and porcine embryos and even human IVF newborns. However, since the first case of sex skewing reported in 1991, the underlying mechanisms remain unclear. We reported herein that sex ratio is skewed in mouse IVF offspring, and this was a result of female-biased peri-implantation developmental defects that were originated from impaired imprinted X chromosome inactivation (iXCI) through reduced ring finger protein 12 (Rnf12)/X-inactive specific transcript (Xist) expression. Compensation of impaired iXCI by overexpression of Rnf12 to up-regulate Xist significantly rescued female-biased developmental defects and corrected sex ratio in IVF offspring. Moreover, supplementation of an epigenetic modulator retinoic acid in embryo culture medium up-regulated Rnf12/Xist expression, improved iXCI, and successfully redeemed the skewed sex ratio to nearly 50% in mouse IVF offspring. Thus, our data show that iXCI is one of the major epigenetic barriers for the developmental competence of female embryos during preimplantation stage, and targeting erroneous epigenetic modifications may provide a potential approach for preventing IVF-associated complications.

Ancestral resurrection of anthropoid estrogen receptor ß demonstrates functional consequences of positive selection.

Anthropoid primates arose during the Eocene approximately 55 million years ago (mya), and extant anthropoids share a most recent common ancestor ∼40mya. Paleontology has been very successful at describing the morphological phenotypes of extinct anthropoids. Less well understood is the molecular biology of these extinct species as well as the phenotypic consequences of evolutionary variation in their genomes. Here we resurrect the most recent common ancestral anthropoid estrogen receptor ß gene (ESR2) and demonstrate that the function of this ancestral estrogen receptor has been maintained during human descent but was altered during early New World monkey (NWM) evolution by becoming a more potent transcriptional activator. We tested hypotheses of adaptive evolution in the protein coding sequences of ESR2, and determined that ESR2 evolved via episodic positive selection on the NWM stem lineage. We separately co-transfected ESR2 constructs for human, NWM, and the anthropoid ancestor along with reporter gene vectors and performed hormone binding dose response experiments that measure transactivation activity. We found the transactivation potentials of the ancestral and human sequences to be significantly lower (p<0.0001 in each comparison) than that of the NWM when treated with estradiol, the most prevalent estrogen. We conclude the difference in fold activation is due to positive selection in the NWM ERß ligand binding domain. Our study validates inferential methods for detecting adaptive evolution that predict functional consequences of nucleotide substitutions and points a way toward examining the functional consequences of positive Darwinian selection.

De novo transcriptome assembly and identification of genes associated with feed conversion ratio and breast muscle yield in domestic ducks.

Breast muscle yield and feed conversion efficiency are the major breeding aims in duck breeding. Understanding the role of specific transcripts in the muscle and small intestine might lead to the elucidation of interrelated biological processes. In this study, we obtained jejunum and breast muscle samples from two strains of Peking ducks that were sorted by feed conversion ratio (FCR) and breast muscle percentage into two-tailed populations. Ten RNA-Seq libraries were developed from the pooled samples and sequenced using the Hiseq2000 platform. We created a reference duck transcript database using de novo assembly methods, which included 16 663 irredundant contigs with an N50 length of 1530 bp. This new duck reference cDNA dataset significantly improved the mapping rate for RNA-Seq data, from 50% to 70%. Mapping and annotation were followed by Gene Ontology analysis, which showed that numerous genes were differentially expressed between the low and high FCR groups. The differentially expressed genes in the jejunum were enriched in biological processes related to immune response and immune response activation, whereas those in the breast muscle were significantly enriched in biological processes related to muscle cell differentiation and organ development. We identified new candidate genes, that is, PCK1, for improving the FCR and breast muscle yield of ducks and obtained much better reference duck transcripts. This study suggested that de novo assembly is essential when applying transcriptome analysis to a species with an incomplete genome.

Dynamic proteomic profiles of in vivo- and in vitro-produced mouse postimplantation extraembryonic tissues and placentas.

As the interface between the mother and the developing fetus, the placenta is believed to play an important role in assisted reproductive technology (ART)-induced aberrant intrauterine and postnatal development. However, the mechanisms underlying aberrant placentation remain unclear, especially during extraembryonic tissue development and early stages of placental formation. Using a mouse model, this investigation provides the first comparative proteomic analysis of in vivo (IVO) and in vitro-produced (IVP) extraembryonic tissues and placentas after IVO fertilization and development, or in vitro fertilization and culture, respectively. We identified 165 and 178 differentially expressed proteins (DEPs) between IVO and IVP extraembryonic tissues and placentas on Embryonic Day 7.5 (E7.5) and E10.5, respectively. Many DEPs were functionally associated with genetic information processing, such as impaired de novo DNA methylation, as well as posttranscriptional, translational and posttranslational dysregulation. These novel findings were further confirmed by global hypomethylation, and a lower level of correlation was found between the transcriptome and proteome in the IVP groups. In addition, numerous DEPs were involved in energy and amino acid metabolism, cytoskeleton organization and transport, and vasculogenesis and angiogenesis. These disturbed processes and pathways are likely to be associated with embryonic intrauterine growth restriction, an enlarged placenta, and impaired labyrinth morphogenesis. This study provides a direct and comprehensive reference for the further exploration of the placental mechanisms that underlie ART-induced developmental aberrations.

A genome-wide association study identifies novel single nucleotide polymorphisms associated with dermal shank pigmentation in chickens.

Shank color of domestic chickens varies from black to blue, green, yellow, or white, which is controlled by the combination of melanin and xanthophylls in dermis and epidermis. Dermal shank pigmentation of chickens is determined by sex-linked inhibitor of dermal melanin (Id), which is located on the distal end of the long arm of Z chromosome, through controlling dermal melanin pigmentation. Although previous studies have focused on the identification of Id and the linear relationship with barring and recessive white skin, no causal mutations have yet been identified in relation to the mutant dermal pigment inhibiting allele at the Id locus. In this study, we first used the 600K Affymetrix Axiom HD genotyping array, which includes ~580,961 SNP of which 26,642 SNP were on the Z chromosome to perform a genome-wide association study on pure lines of 19 Tibetan hens with dermal pigmentation shank and 21 Tibetan hens with yellow shank to refine the Id location. Association analysis was conducted by the PLINK software using the standard chi-squared test, and then Bonferroni correction was used to adjust multiple testing. The genome-wide study revealed that 3 SNP located at 78.5 to 79.2 Mb on the Z chromosome in the current assembly of chicken genome (galGal4) were significantly associated with dermal shank pigmentation of chickens, but none of them were located in known genes. The interval we refined was partly converged with previous results, suggesting that the Id gene is in or near our refined genome region. However, the genomic context of this region was complex. There were only 15 SNP markers developed by the genotyping array within the interval region, in which only 1 SNP marker passed quality control. Additionally, there were about 5.8-Mb gaps on both sides of the refined interval. The follow-up replication studies may be needed to further confirm the functional significance for these newly identified SNP.