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.

Identification of alternative splicing events related to fatty liver formation in duck using full-length transcripts.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is one of most common diseases in the world. Recently, alternative splicing (AS) has been reported to play a key role in NAFLD processes in mammals. Ducks can quickly form fatty liver similar to human NAFLD after overfeeding and restore to normal liver in a short time, suggesting that ducks are an excellent model to unravel molecular mechanisms of lipid metabolism for NAFLD. However, how alternative splicing events (ASEs) affect the fatty liver process in ducks is still unclear. RESULTS: Here we identify 126,277 unique transcripts in liver tissue from an overfed duck (77,237 total transcripts) and its sibling control (69,618 total transcripts). We combined these full-length transcripts with Illumina RNA-seq data from five pairs of overfed ducks and control individuals. Full-length transcript sequencing provided us with structural information of transcripts and Illumina RNA-seq data reveals the expressional profile of each transcript. We found, among these unique transcripts, 30,618 were lncRNAs and 1,744 transcripts including 155 lncRNAs and 1,589 coding transcripts showed significantly differential expression in liver tissues between overfed ducks and control individuals. We also detected 27,317 ASEs and 142 of them showed significant relative abundance changes in ducks under different feeding conditions. Full-length transcript profiles together with Illumina RNA-seq data demonstrated that 10 genes involving in lipid metabolism had ASEs with significantly differential abundance in normally fed (control) and overfed ducks. Among these genes, protein products of five genes (CYP4F22, BTN, GSTA2, ADH5, and DHRS2 genes) were changed by ASEs. CONCLUSIONS: This study presents an example of how to identify ASEs related to important biological processes, such as fatty liver formation, using full-length transcripts alongside Illumina RNA-seq data. Based on these data, we screened out ASEs of lipid-metabolism related genes which might respond to overfeeding. Our future ability to explore the function of genes showing AS differences between overfed ducks and their sibling controls, using genetic manipulations and co-evolutionary studies, will certainly extend our knowledge of genes related to the non-pathogenic fatty liver process.

Multiscale 3D genome organization underlies duck fatty liver with no adipose inflammation or serious injury.

Non-alcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease. Little is known about how gene expression and chromatin structure are regulated in NAFLD due to lack of suitable model. Ducks naturally develop fatty liver similar to serious human non-alcoholic fatty liver (NAFL) without adipose inflammation and liver fibrosis, thus serves as a good model for investigating molecular mechanisms of adipose metabolism and anti-inflammation. Here, we constructed a NAFLD model without adipose inflammation and liver fibrosis in ducks. By performing dynamic pathological and transcriptomic analyses, we identified critical genes involving in regulation of the NF-κB and MHCII signaling, which usually lead to adipose inflammation and liver fibrosis. We further generated dynamic three-dimensional chromatin maps during liver fatty formation and recovery. This showed that ducks enlarged hepatocyte cell nuclei to reduce inter-chromosomal interaction, decompress chromatin structure, and alter strength of intra-TAD and loop interactions during fatty liver formation. These changes partially contributed to the tight control the NF-κB and the MHCII signaling. Our analysis uncovers duck chromatin reorganization might be advantageous to maintain liver regenerative capacity and reduce adipose inflammation. These findings shed light on new strategies for NAFLD control.

Functional divergence of oligoadenylate synthetase 1 (OAS1) proteins in Tetrapods.

OASs play critical roles in immune response against virus infection by polymerizing ATP into 2-5As, which initiate the classical OAS/RNase L pathway and induce degradation of viral RNA. OAS members are functionally diverged in four known innate immune pathways (OAS/RNase L, OASL/IRF7, OASL/RIG-I, and OASL/cGAS), but how they functionally diverged is unclear. Here, we focus on evolutionary patterns and explore the link between evolutionary processes and functional divergence of Tetrapod OAS1. We show that Palaeognathae and Primate OAS1 genes are conserved in genomic and protein structures but differ in function. The former (i.e., ostrich) efficiently synthesized long 2-5A and activated RNase L, while the latter (i.e., human) synthesized short 2-5A and did not activate RNase L. We predicted and verified that two in-frame indels and one positively selected site in the active site pocket contributed to the functional divergence of Palaeognathae and Primate OAS1. Moreover, we discovered and validated that an in-frame indel in the C-terminus of Palaeognathae OAS1 affected the binding affinity of dsRNA and enzymatic activity, and contributed to the functional divergence of Palaeognathae OAS1 proteins. Our findings unravel the molecular mechanism for functional divergence and give insights into the emergence of novel functions in Tetrapod OAS1.