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.

Generation of genome-edited chicken and duck lines by adenovirus-mediated in vivo genome editing.

Conventional avian genome editing is mediated by isolation, culture, and genome editing of primordial germ cells (PGCs); screening and propagating the genome-edited PGCs; and transplantation of the PGCs into recipient embryos. The PGC-mediated procedures, however, are technically difficult, and therefore, the conventional method has previously been utilized only in chickens. Here, we generated germline mosaic founder chicken and duck lines without the PGC-mediated procedures by injecting an adenovirus containing the CRISPR-Cas9 system into avian blastoderms. Genome-edited chicken and duck offspring produced from the founders carried different insertion or deletion mutations without mutations in the potential off-target sites. Our data demonstrate successful applications of the adenovirus-mediated method for production of genome-edited chicken and duck lines.

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.

A novel qPCR-based technique for identifying avian sex: An illustration within embryonic craniofacial bone.

Sex is a biological variable important to consider in all biomedical experiments. However, doing so in avian embryos can be challenging as sex can be morphologically indistinguishable. Unlike humans, female birds are the heterogametic sex with Z and W sex chromosomes. The female-specific W chromosome has previously been identified in chick using a species-specific polymerase chain reaction (PCR) technique. We developed a novel reverse transcription quantitative PCR (RT-qPCR) technique that amplifies the W chromosome gene histidine triad nucleotide-binding protein W (HINTW) in chick, quail, and duck. Accuracy of the HINTW RT-qPCR primer set was confirmed in all three species using species-specific PCR, including a novel quail-specific HINTW PCR primer set. Bone development-related gene expression was then analyzed by sex in embryonic lower jaws of duck and quail, as adult duck beak size is known to be sexually dimorphic while quail beak size is not. Trends toward sex differences were found in duck gene expression but not in quail, as expected. With these novel RT-qPCR and PCR embryo sexing methods, sex of chick, quail, and duck embryos can now be assessed by either/both RNA and DNA, which facilitates analysis of sex as a biological variable in studies using these model organisms.

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.

Genome-wide association study reveals serovar-associated genetic loci in Riemerella anatipestifer.

BACKGROUND: The disease caused by Riemerella anatipestifer (R. anatipestifer, RA) results in large economic losses to the global duck industry every year. Serovar-related genomic variation, such as the O-antigen and capsular polysaccharide (CPS) gene clusters, has been widely used for serotyping in many gram-negative bacteria. RA has been classified into at least 21 serovars based on slide agglutination, but the molecular basis of serotyping is unknown. In this study, we performed a pan-genome-wide association study (Pan-GWAS) to identify the genetic loci associated with RA serovars. RESULTS: The results revealed a significant association between the putative CPS synthesis gene locus and the serological phenotype. Further characterization of the CPS gene clusters in 11 representative serovar strains indicated that they were highly diverse and serovar-specific. The CPS gene cluster contained the key genes wzx and wzy, which are involved in the Wzx/Wzy-dependent pathway of CPS synthesis. Similar CPS loci have been found in some other species within the family Weeksellaceae. We have also shown that deletion of the wzy gene in RA results in capsular defects and cross-agglutination. CONCLUSIONS: This study indicates that the CPS synthesis gene cluster of R. anatipestifer is a serotype-specific genetic locus. Importantly, our finding provides a new perspective for the systematic analysis of the genetic basis of the R anatipestifer serovars and a potential target for establishing a complete molecular serotyping scheme.

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

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

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.

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.

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.

Integrated transcriptomics and metabolomics study of embryonic breast muscle of Jiaji ducks.

Because number of matured muscle fibers in poultry does not increase after birth, the meat yield is mainly determined during embryogenesis. We previously indicated breast muscle grew rapidly from 18th day after hatching (E18) to E27, and almost stopped from E27 to E34 of Jiaji ducks, while the mechanism is unclear. This study utilized RNA-seq to explore the related genes of muscle development and their relationship with small molecule metabolites at E18, E27 and E34 of Jiaji ducks. Several thousand differentially expressed genes (DEGs) were detected among E18, E27 and E34. DEGs expression profiles included 8 trend maps, among which trend 1 was opposite to and trend 6 was consistent with breast muscle development trend of Jiaji ducks. Through joint analysis between trend 1 of DEGs and trend 1 of differential metabolites (DEMs), protein digestion and absorption pathway stood out. The decrease of COL8A2 gene expression will lead to the decrease of arginine content, which will inhibit the development of breast muscle in embryonic Jiaji duck. Similarly, joint analysis between trend 6 of DEGs and trend 6 of DEMs indicated the increase of GAMT gene expression will cause the increase of proline content, and then promote the development of breast muscle of Jiaji duck in embryonic period. These results will be helpful for further understanding the mechanism of muscle yields of Jiaji ducks.

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.

Determining evolutionary origin and phylogenetic relationships of mallard-like ducks of Oceania, greater Indonesia, and the Philippines with ddRAD-seq data.

AIM: We aim to determine the evolutionary origins and population genetics of mallard-like ducks of Oceania, greater Indonesia, and the Philippines. LOCATION: Oceania, greater Indonesia, and the Philippines. TAXON: Mallard (Anas platyrhynchos), Pacific black duck (A. superciliosa spp.), and Philippine duck (A. luzonica) METHODS: Thousands of nuclear ddRAD-seq loci and the mitochondrial DNA control region were assayed across individuals representative of each species' range. We assessed population structure and phylogenetic relationships, as well as estimated demographic histories to reconstruct the biogeographical history of each species. RESULTS: Philippine and Pacific black ducks represent unique genetic lineages that diverged from the mallard 1-2 million years ago. We find no support for the Philippine duck representing a hybrid species as once posited; however, their low levels of genetic diversity requires further attention. We find a lack of substructure among Philippine ducks. However, we found pronounced differentiation between subspecies of Pacific black ducks, especially between A. s. superciliosa from New Zealand and A. s. rogersi from Australia, Papua New Guinea, and Timor-Leste, Indonesia. Anas superciliosa pelewensis gave mixed results; individuals from the Solomon Islands were differentiated from the other subspecies, but those from the island of Aunu'u, American Samoa, were genetically more similar to A. s. rogersi than A. s. pelewensis samples from the Solomon Islands. Finally, we find limited evidence of interspecific gene flow at evolutionary scales, and mallard introgression among contemporary samples. MAIN CONCLUSIONS: Mallard-like ducks radiated across Oceania, greater Indonesia, and the Philippines within the last 2 million years. Only the Pacific black duck showed unique sub-structuring that largely followed known sub-species ranges, except for A. s. pelewensis. We posit that the high interrelatedness among Solomon Island samples suggests that their genetic distinctiveness may simply be the result of high levels of genetic drift. In contrast, we conclude that mainland Australian Pacific black ducks were the most likely source for the recent colonization of American Samoa. As a result, our findings suggest that either the A. s. pelewensis subspecies designations and/or its geographical range may require re-evaluation. Continued re-evaluation of evolutionary and taxonomic relationships is necessary when attempting to reconstruct and understand biogeographical histories, with important implications towards any attempts to implement conservation strategies.

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.

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.

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-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.

Discover QTLs for the level of blood components in Shaoxing duck using GWAS and haplotype sharing analysis.

Blood composition is indicative of health-related traits such as immunity and metabolism. The use of molecular genetics to investigate alterations in these attributes in laying ducks is a novel approach. Our objective was to employ genome - wide association studies (GWAS) and haplotype - sharing analysis to identify genomic regions and potential genes associated with 11 blood components in Shaoxing ducks. Our findings revealed 35 SNPs and 1 SNP associated with low-density lipoprotein cholesterol (LDL) and globulin (GLB), respectively. We identified 36 putative candidate genes for the LDL trait in close proximity to major QTLs and key loci. Based on their biochemical and physiological properties, TRA2A, NPY, ARHGEF26, DHX36, and AADAC are the strongest putative candidate genes. Through linkage disequilibrium analysis and haplotype sharing analysis, we identified three haplotypes and one haplotype, respectively, that were significantly linked with LDL and GLB. These haplotypes could be selected as potential candidate haplotypes for molecular breeding of Shaoxing ducks. Additionally, we utilized a bootstrap test to verify the reliability of GWAS with small experimental samples. The test can be accessed at https://github.com/xuwenwu24/Bootstrap-test.

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 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.