Genome-wide association analysis explores the genetic loci of amino acid content in duck's breast muscle.

BACKGROUND: Amino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat. RESULT: In the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55-56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2-21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways. CONCLUSION: In this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

Whole-genome sequences restore the original classification of dabbling ducks (genus Anas).

Anas, is a genus of dabbling ducks and encompasses a considerable number of species, among which some are the progenitors of domestic ducks. However, the taxonomic position of the Anas genus remains uncertain because several of its species, initially categorized as Anas based on morphological characteristics, were subsequently reclassified and grouped with the South American genus Tachyeres, primarily based on analysis of their mitochondrial gene sequences. Here, we constructed a phylogenetic tree using nine of our recently assembled Anas genomes, two Tachyeres genomes, and one Cairina genome that are publicly available. The results showed that the Northern shoveler (Anas clypeata) and Baikal teal (Anas formosa) clustered with the other Anas species at the whole-genome level rather than with the Steamer ducks (genus Tachyeres). Therefore, we propose to restore the original classification of the Anas genus, which includes the Northern shoveler and Baikal teal species, 47 species in total. Moreover, our study unveiled extensive incomplete lineage sorting and an ancient introgression event from Tachyeres to Anas, which has led to notable phylogenetic incongruence within the Anas genome. This ancient introgression event not only supports the theory that Anas originated in South America but also that it played a significant role in shaping the evolutionary trajectory of Anas, including the domestic duck.

Impacts of longitudinal water curtain cooling system on transcriptome-related immunity in ducks.

BACKGROUND: The closed poultry houses integrated with a longitudinal water curtain cooling system (LWCCS) are widely used in modern poultry production. This study showed the variations in environmental conditions in closed houses integrated with a longitudinal water curtain cooling system. We evaluated the influence of different environmental conditions on duck growth performance and the transcriptome changes of immune organs, including the bursa of Fabricius and the spleen. RESULT: This study investigated the slaughter indicators and immune organ transcriptomes of 52-day-old Cherry Valley ducks by analyzing the LWCC at different locations (water curtain end, middle position, and fan cooling end). The results showed that the cooling effect of the LWCCS was more evident from 10:00 a.m. -14:00. And from the water curtain end to the fan cooling end, the hourly average temperature differently decreased by 0.310℃, 0.450℃, 0.480℃, 0.520℃, and 0.410℃, respectively (P < 0.05). The daily and hourly average relative humidity decreased from the water curtain end to the fan cooling end, dropping by 7.500% and 8.200%, respectively (P < 0.01). We also observed differences in production performance, such as dressing weight, half-eviscerated weight, skin fat rate, and percentage of abdominal fat (P < 0.01), which may have been caused by environmental conditions. RNA-sequencing (RNA-seq) revealed 211 and 279 differentially expressed genes (DEGs) in the ducks' bursa of Fabricius and spleen compared between the water curtain end and fan cooling end, respectively. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the two organs showed the DEGs were mainly enriched in cytokine-cytokine receptor interaction, integral component of membrane, Retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) signaling pathway, etc. Our results implied that full-closed poultry houses integrated with LWCCS could potentially alter micro-environments (water curtain vs. fan cooling), resulting in ducks experiencing various stressful situations that eventually affect their immunity and production performance. CONCLUSION: In this study, our results indicated that uneven distributions of longitudinal environmental factors caused by LWCCS would affect the dressed weight, breast muscle weight, skin fat rate, and other product performance. Moreover, the expression of immune-related genes in the spleen and bursa of ducks could be affected by the LWCCS. This provides a new reference to optimize the use of LWCCS in conjunction with close duck houses in practical production.

Integrating miRNA and full-length transcriptome profiling to elucidate the mechanism of muscle growth in Muscovy ducks reveals key roles for miR-301a-3p/ANKRD1.

BACKGROUND: The popularity of Muscovy ducks is attributed not only to their conformation traits but also to their slightly higher content of breast and leg meat, as well as their stronger-tasting meat compared to that of typical domestic ducks. However, there is a lack of comprehensive systematic research on the development of breast muscle in Muscovy ducks. In addition, since the number of skeletal muscle myofibers is established during the embryonic period, this study conducted a full-length transcriptome sequencing and microRNA sequencing of the breast muscle. Muscovy ducks at four developmental stages, namely Embryonic Day 21 (E21), Embryonic Day 27 (E27), Hatching Day (D0), and Post-hatching Day 7 (D7), were used to isolate total RNA for analysis. RESULTS: A total of 68,161 genes and 472 mature microRNAs were identified. In order to uncover deeper insights into the regulation of mRNA by miRNAs, we conducted an integration of the differentially expressed miRNAs (known as DEMs) with the differentially expressed genes (referred to as DEGs) across various developmental stages. This integration allowed us to make predictions regarding the interactions between miRNAs and mRNA. Through this analysis, we identified a total of 274 DEGs that may serve as potential targets for the 68 DEMs. In the predicted miRNA‒mRNA interaction networks, let-7b, miR-133a-3p, miR-301a-3p, and miR-338-3p were the hub miRNAs. In addition, multiple DEMs also showed predicted target relationships with the DEGs associated with skeletal system development. These identified DEGs and DEMs as well as their predicted interaction networks involved in the regulation of energy homeostasis and muscle development were most likely to play critical roles in facilitating the embryo-to-hatchling transition. A candidate miRNA, miR-301a-3p, exhibited increased expression during the differentiation of satellite cells and was downregulated in the breast muscle tissues of Muscovy ducks at E21 compared to E27. A dual-luciferase reporter assay suggested that the ANKRD1 gene, which encodes a transcription factor, is a direct target of miR-301a-3p. CONCLUSIONS: miR-301a-3p suppressed the posttranscriptional activity of ANKRD1, which is an activator of satellite cell proliferation, as determined with gain- and loss-of-function experiments. miR-301a-3p functions as an inducer of myogenesis by targeting the ANKRD1 gene in Muscovy ducks. These results provide novel insights into the early developmental process of black Muscovy breast muscles and will improve understanding of the underlying molecular mechanisms.

Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits.

BACKGROUND: Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits. RESULTS: In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits. CONCLUSIONS: Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production traits.

Establishment of a hydrodynamic delivery system in ducks.

Chronic hepatitis B virus (HBV) poses a significant global health challenge as it can lead to acute or chronic liver disease and hepatocellular carcinoma (HCC). To establish a safety experimental model, a homolog of HBV-duck HBV (DHBV) is often used for HBV research. Hydrodynamic-based gene delivery (HGD) is an efficient method to introduce exogenous genes into the liver, making it suitable for basic research. In this study, a duck HGD system was first constructed by injecting the reporter plasmid pLIVE-SEAP via the ankle vein. The highest expression of SEAP occurred when ducks were injected with 5 µg/mL plasmid pLIVE-SEAP in 10% bodyweight volume of physiological saline for 6 s. To verify the distribution and expression of exogenous genes in multiple tissues, the relative level of foreign gene DNA and ß-galactosidase staining of LacZ were evaluated, which showed the plasmids and their products were located mainly in the liver. Additionally, ß-galactosidase staining and fluorescence imaging indicated the delivered exogenous genes could be expressed in a short time. Further, the application of the duck HGD model on DHBV treatment was investigated by transferring representative anti-HBV genes IFNα and IFNγ into DHBV-infected ducks. Delivery of plasmids expressing IFNα and IFNγ inhibited DHBV infection and we established a novel efficient HGD method in ducks, which could be useful for drug screening of new genes, mRNAs and proteins for anti-HBV treatment.

Complex N-glycans are important for interspecies transmission of H7 influenza A viruses.

Influenza A viruses (IAVs) can overcome species barriers by adaptation of the receptor-binding site of the hemagglutinin (HA). To initiate infection, HAs bind to glycan receptors with terminal sialic acids, which are either N-acetylneuraminic acid (NeuAc) or N-glycolylneuraminic acid (NeuGc); the latter is mainly found in horses and pigs but not in birds and humans. We investigated the influence of previously identified equine NeuGc-adapting mutations (S128T, I130V, A135E, T189A, and K193R) in avian H7 IAVs in vitro and in vivo. We observed that these mutations negatively affected viral replication in chicken cells but not in duck cells and positively affected replication in horse cells. In vivo, the mutations reduced virus virulence and mortality in chickens. Ducks excreted high viral loads longer than chickens, although they appeared clinically healthy. To elucidate why these viruses infected chickens and ducks despite the absence of NeuGc, we re-evaluated the receptor binding of H7 HAs using glycan microarray and flow cytometry studies. This re-evaluation demonstrated that mutated avian H7 HAs also bound to α2,3-linked NeuAc and sialyl-LewisX, which have an additional fucose moiety in their terminal epitope, explaining why infection of ducks and chickens was possible. Interestingly, the α2,3-linked NeuAc and sialyl-LewisX epitopes were only bound when presented on tri-antennary N-glycans, emphasizing the importance of investigating the fine receptor specificities of IAVs. In conclusion, the binding of NeuGc-adapted H7 IAV to tri-antennary N-glycans enables viral replication and shedding by chickens and ducks, potentially facilitating interspecies transmission of equine-adapted H7 IAVs.IMPORTANCEInfluenza A viruses (IAVs) cause millions of deaths and illnesses in birds and mammals each year. The viral surface protein hemagglutinin initiates infection by binding to host cell terminal sialic acids. Hemagglutinin adaptations affect the binding affinity to these sialic acids and the potential host species targeted. While avian and human IAVs tend to bind to N-acetylneuraminic acid (sialic acid), equine H7 viruses prefer binding to N-glycolylneuraminic acid (NeuGc). To better understand the function of NeuGc-specific adaptations in hemagglutinin and to elucidate interspecies transmission potential NeuGc-adapted viruses, we evaluated the effects of NeuGc-specific mutations in avian H7 viruses in chickens and ducks, important economic hosts and reservoir birds, respectively. We also examined the impact on viral replication and found a binding affinity to tri-antennary N-glycans containing different terminal epitopes. These findings are significant as they contribute to the understanding of the role of receptor binding in avian influenza infection.

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

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

Genome-wide association study for bone quality of ducks during the laying period.

The cage-rearing model of the modern poultry industry makes the bones of birds, especially egg-laying birds, more vulnerable to fracture, which poses serious damage to the health of birds. Research confirms that genetic material plays an important role in regulating bone growth, development, and remodeling. However, the genetic architecture underlying bone traits is not well understood. The objectives of this study are to identify valuable genes and genetic markers through a genome-wide association study (GWAS) for breeding to improve the duck bone quality. First, we quantified the tibia and femur quality traits of 260 laying ducks. Based on GWAS, a total of 75 SNP loci significantly associated with bone quality traits were identified, and 67 potential candidate genes were annotated. According to gene function analysis, genes P4HA2, WNT3A, and BST1 et al may influence bone quality by regulating bone cell activity, calcium and phosphate metabolism, or bone collagen maturation and cross-linking. Meanwhile, combined with the transcriptome results, we found that HOXB cluster genes are also important in bone growth and development. Therefore, our findings were helpful in further understanding the genetic architecture of the duck bone quality and provided a worthy theoretical basis and technological support to improve duck bone quality by breeding.

Interleukin-2 enhancer binding factor 2 negatively regulates the replication of duck hepatitis A virus type 1 by disrupting the RNA-dependent RNA polymerase activity of 3D polymerase.

The interaction between viral components and cellular proteins plays a crucial role in viral replication. In a previous study, we showed that the 3'-untranslated region (3'-UTR) is an essential element for the replication of duck hepatitis A virus type 1 (DHAV-1). However, the underlying mechanism is still unclear. To gain a deeper understanding of this mechanism, we used an RNA pull-down and a matrix-assisted laser desorption/ionization time-of-flight mass spectrometry assay to identify new host factors that interact with the 3'-UTR. We selected interleukin-2 enhancer binding factor 2 (ILF2) for further analysis. We showed that ILF2 interacts specifically with both the 3'-UTR and the 3D polymerase (3Dpol) of DHAV-1 through in vitro RNA pull-down and co-immunoprecipitation assays, respectively. We showed that ILF2 negatively regulates viral replication in duck embryo fibroblasts (DEFs), and that its overexpression in DEFs markedly suppresses DHAV-1 replication. Conversely, ILF2 silencing resulted in a significant increase in viral replication. In addition, the RNA-dependent RNA polymerase (RdRP) activity of 3Dpol facilitated viral replication by enhancing viral RNA translation efficiency, whereas ILF2 disrupted the role of RdRP in viral RNA translation efficiency to suppress DHAV-1 replication. At last, DHAV-1 replication markedly suppressed the expression of ILF2 in DEFs, duck embryo hepatocytes, and different tissues of 1 day-old ducklings. A negative correlation was observed between ILF2 expression and the viral load in primary cells and different organs of young ducklings, suggesting that ILF2 may affect the viral load both in vitro and in vivo.

Whole-transcriptome sequencing revealed the role of noncoding RNAs in susceptibility and resistance of Pekin ducks to DHAV-3.

As the most prevalent pathogen of duck viral hepatitis (DVH), duck hepatitis A virus genotype 3 (DHAV-3) has caused huge economic losses to the duck industry in China. Herein, we obtained whole-transcriptome sequencing data of susceptible (S) and resistant (R) Pekin duckling samples at 0 h, 12 h, and 24 h after DHAV-3 infection. We found that DHAV-3 infection induces 5,396 differentially expressed genes (DEGs), 85 differentially expressed miRNAs (DEMs), and 727 differentially expressed lncRNAs (DELs) at 24 hpi in S vs. R ducks, those upregulated genes were enriched in inflammation and cell communications pathways and downregulated genes were related to metabolic processes. Upregulated genes showed high connectivity with the miR-33, miR-193, and miR-11591, and downregulated genes were mainly regulated by miR-2954, miR-125, and miR-146b. With the construction of lncRNA-miRNA-mRNA axis, we further identified a few aberrantly expressed lncRNAs (e.g., MSTRG.36194.1, MSTRG.50601.1, MSTRG.34328.7, and MSTRG.29445.1) that regulate expression of hub genes (e.g., THBD, CLIC2, IL8, ACOX2, GPHN, SMLR1, and HAO1) by sponging those highly connected miRNAs. Altogether, our findings defined a dual role of ncRNAs in immune and metabolic regulation during DHAV-3 infection, suggesting potential new targets for treating DHAV-3 infected ducks.

Identification of candidate genes affecting the tibia quality in Nonghua duck.

The skeleton is a vital organ providing structural support in poultry. Weakness in bone structure can lead to deformities, osteoporosis, cage fatigue, and fractures, resulting in economic losses. Research has substantiated that genetic factors play a significant role in influencing bone quality. The discovery of genetic markers associated with bone quality holds paramount importance for enhancing genetic traits related to the skeletal system in poultry. This study analyzed nine phenotypic indicators of tibia quality in 120-day-old ducks. The phenotypic correlation revealed a high correlation among diameter, Perimeter, and weight (0.69-0.78), and a strong correlation was observed between toughness and breaking strength (0.62). Then, we conducted a genome-wide association analysis of the phenotypic indicators to elucidate the genetic basis of tibial quality in Nonghua ducks. Among the 11 candidate genes that were annotated, TAPT1, BST1, and STIM2 were related to the diameter indicator, ZNF652, IGF2BP1, CASK, and GREB1L were associated with the weight and toughness indicators. RFX8, GLP1R, and DNAAF5 were identified for ash, calcium, and phosphorus content, respectively. Finally, KEGG and GO analysis for annotated genes were performed. STIM2 and BST1 were enriched into the Calcium signalling pathway and Niacin and nicotinamide metabolic pathway, which may be key candidate genes affecting bone quality phenotypes. Gene expression analysis of the candidate genes, such as STIM2, BST1, TAPT1, and CASK showed higher expression levels in bones compared to other tissues. The obtained results can contribute to new insights into tibial quality and provide new genetic biomarkers that can be employed in duck breeding.

An in vivo duck hepatitis B virus model recapitulates key aspects of nucleic acid polymer treatment outcomes in chronic hepatitis B patients.

Nucleic acid polymers (NAPs) are an attractive treatment modality for chronic hepatitis B (CHB), with REP2139 and REP2165 having shown efficacy in CHB patients. A subset of patients achieve functional cure, whereas the others exhibit a moderate response or are non-responders. NAP efficacy has been difficult to recapitulate in animal models, with the duck hepatitis B virus (DHBV) model showing some promise but remaining underexplored for NAP efficacy testing. Here we report on an optimized in vivo DHBV duck model and explore several characteristics of NAP treatment. REP2139 was efficacious in reducing DHBV DNA and DHBsAg levels in approximately half of the treated ducks, whether administered intraperitoneally or subcutaneously. Intrahepatic or serum NAP concentrations did not correlate with efficacy, nor did the appearance of anti-DHBsAg antibodies. Furthermore, NAP efficacy was only observed in experimentally infected ducks, not in endogenously infected ducks (vertical transmission). REP2139 add-on to entecavir treatment induced a deeper and more sustained virological response compared to entecavir monotherapy. Destabilized REP2165 showed a different activity profile with a more homogenous antiviral response followed by a faster rebound. In conclusion, subcutaneous administration of NAPs in the DHBV duck model provides a useful tool for in vivo evaluation of NAPs. It recapitulates many aspects of this class of compound's efficacy in CHB patients, most notably the clear division between responders and non-responders.

Multiplex-PCR method application to identify duck blood and its adulterated varieties.

This study applied and validated the Multiplex-PCR method to identify the authenticity of duck blood and four common adulterated animal blood varieties. To this end, the genomic DNAs of duck blood and its counterfeit products were extracted using an efficient high-throughput extraction method. Specific primers were designed using the cytochrome b gene. The reaction system and conditions of a multiplex (namely, Five-plex) PCR were optimized, and the proposed methodology was verified, proving its good specificity, repeatability, and sensitivity. The Five-plex PCR system detected nine duck blood samples sold in the local market, revealing the adulteration of duck blood products. The Multiplex-PCR system can accurately and quickly detect adulterated animal blood in duck blood products, effectively finding counterfeits and identifying the authenticity of genuine duck blood products.

A new chromosome-scale duck genome shows a major histocompatibility complex with several expanded multigene families.

BACKGROUND: The duck (Anas platyrhynchos) is one of the principal natural hosts of influenza A virus (IAV), harbors almost all subtypes of IAVs and resists to many IAVs which cause extreme virulence in chicken and human. However, the response of duck's adaptive immune system to IAV infection is poorly characterized due to lack of a detailed gene map of the major histocompatibility complex (MHC). RESULTS: We herein reported a chromosome-scale Beijing duck assembly by integrating Nanopore, Bionano, and Hi-C data. This new reference genome SKLA1.0 covers 40 chromosomes, improves the contig N50 of the previous duck assembly with highest contiguity (ZJU1.0) of more than a 5.79-fold, surpasses the chicken and zebra finch references in sequence contiguity and contains a complete genomic map of the MHC. Our 3D MHC genomic map demonstrated that gene family arrangement in this region was primordial; however, families such as AnplMHCI, AnplMHCIIß, AnplDMB, NKRL (NK cell receptor-like genes) and BTN underwent gene expansion events making this area complex. These gene families are distributed in two TADs and genes sharing the same TAD may work in a co-regulated model. CONCLUSIONS: These observations supported the hypothesis that duck's adaptive immunity had been optimized with expanded and diversified key immune genes which might help duck to combat influenza virus. This work provided a high-quality Beijing duck genome for biological research and shed light on new strategies for AIV control.

Analysis of Transcriptomic Differences in the Ovaries of High- and Low-Laying Ducks.

The egg-laying performance of Shan Ma ducks (Anas Platyrhynchos) is a crucial economic trait. Nevertheless, limited research has been conducted on the egg-laying performance of this species. We examined routine blood indicators and observed higher levels of metabolic and immune-related factors in the high-egg-production group compared with the low-egg-production group. Furthermore, we explored the ovarian transcriptome of both high- and low-egg-production groups of Shan Ma ducks using Illumina NovaSeq 6000 sequencing. A total of 1357 differentially expressed genes (DEGs) were identified, with 686 down-regulated and 671 up-regulated in the high-egg-production (HEP) ducks and low-egg-production (LEP) ducks. Several genes involved in the regulation of ovarian development, including neuropeptide Y (NPY), cell cycle protein-dependent kinase 1 (CDK1), and transcription factor 1 (E2F1), exhibited significant differential expressions at varying stages of egg production. Pathway functional analysis revealed that the DEGs were primarily associated with the steroid biosynthesis pathway, and the neuroactive ligand-receptor interaction pathway exhibited higher activity in the HEP group compared to the LEP group. This study offers valuable information about and novel insights into high egg production.

Comparative analysis of the hypothalamus transcriptome of laying ducks with different residual feeding intake.

Feed costs account for approximately 60 to 70% of the cost of poultry farming, and feed utilization is closely related to the profitability of the poultry industry. To understand the causes of the differences in feeding in Shan Partridge ducks, we compared the hypothalamus transcriptome profiles of 2 groups of ducks using RNA-seq. The 2 groups were: 1) low-residual feed intake (LRFI) group with low feed intake but high feed efficiency, and 2) high-residual feed intake (HRFI) group with high feed intake but low feed efficiency. We found 78 DEGs were enriched in 9 differential Kyoto Encyclopedia of Genes and Genome (KEGG) pathways, including neuroactive ligand-receptor interaction, GABAergic synapse, nitrogen metabolism, cAMP signaling pathway, calcium signaling pathway, nitrogen metabolism, tyrosine metabolism, ovarian steroidogenesis, and gluconeogenesis. To further identify core genes among the 78 DEGs, we performed protein-protein interaction and coexpression network analyses. After comprehensive analysis and experimental validation, 4 core genes, namely, glucagon (GCG), cholecystokinin (CCK), gamma-aminobutyric acid type A receptor subunit delta (GABRD), and gamma-aminobutyric acid type A receptor subunit beta1 (GABRB1), were identified as potential core genes responsible for the difference in residual feeding intake between the 2 breeds. We also investigated the level of cholecystokinin (CCK), neuropeptide Y (NPY), peptide YY (PYY), ghrelin, and glucagon-like peptide1 (GLP-1) hormones in the sera of Shan Partridge ducks at different feeding levels and found that there was a difference between the 2 groups with respect to GLP-1 and NPY levels. The findings will serve as a reference for future research on the feeding efficiency of Shan Partridge ducks and assist in promoting their genetic breeding.

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

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

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.

Characteristics of the a sequence of the duck Plague virus genome and specific cleavage of the viral genome based on the a sequence.

During the replication process, the herpesvirus genome forms the head-to-tail linked concatemeric genome, which is then cleaved and packaged into the capsid. The cleavage and packing process is carried out by the terminase complex, which specifically recognizes and cleaves the concatemeric genome. This process is governed by a cis-acting sequence in the genome, named the a sequence. The a sequence and genome cleavage have been described in some herpesviruses, but it remains unclear in duck plague virus. In this study, we analysed the location, composition, and conservation of a sequence in the duck plague virus genome. The structure of the DPV genome has an a sequence of (DR4)m-(DR2)n-pac1-S termini (32 bp)-L termini (32 bp)-pac2, and the length is 841 bp. Direct repeat (DR) sequences are conserved in different DPV strains, but the number of DR copies is inconsistent. Additionally, the typical DR1 sequence was not found in the DPV a sequence. The Pac1 and pac2 motifs are relatively conserved between DPV and other herpesviruses. Cleavage of the DPV concatemeric genome was detected, and the results showed that the DPV genome can form a concatemer and is cleaved into a monomer at a specific site. We also established a sensitive method, TaqMan dual qRT‒PCR, to analyse genome cleavage. The ratio of concatemer to total viral genome was decreased during the replication process. These results will be critical for understanding the process of DPV genome cleavage, and the application of TaqMan dual qRT‒PCR will greatly facilitate more in-depth research.